使用DESeq2进行差异表达分析时遇到如下报错
> res<-results(dds,contrast = c("condition",rev(levels(group_list))))
> resOrdered<-res[order(res$pavlue),]#按照P排序
Error in h(simpleError(msg, call)) :
在为'['函数选择方法时评估'i'参数出了错: argument 1 is not a vector
前面的代码如下
library(DESeq2)
colData<- data.frame(row.names = colnames(a),condition=group_list)
dds<-DESeqDataSetFromMatrix(countData = a,colData = colData,design=~condition)
dds<-DESeq(dds)
#两两比较
res<-results(dds,contrast = c("condition",rev(levels(group_list))))
resOrdered<-res[order(res$pavlue),]#按照P排序
DEG<-as.data.frame(resOrdered)
head(DEG)
```