合并GEO数据后分析差异基因报错

[iridescentii]

> #write table
> diffSig <- allDiff[with(allDiff, (abs(logFC)>logFoldChange & adj.P.Val < adjustP )), ]
> write.table(diffSig,file="diff.xls",sep="\t",quote=F,row.names=F)
> diffUp <- allDiff[with(allDiff, (logFC>logFoldChange & adj.P.Val < adjustP )), ]
> write.table(diffUp,file="up.xls",sep="\t",quote=F,row.names=F)
> diffDown <- allDiff[with(allDiff, (logFC<(-logFoldChange) & adj.P.Val < adjustP )), ]
> write.table(diffDown,file="down.xls",sep="\t",quote=F,row.names=F)
> 
> #write expression level of diff gene
> hmExp=rt[as.vector(diffSig[,1]),]
Error in rt[as.vector(diffSig[, 1]), ] : 只有负下标里才能有零
> diffExp=rbind(id=colnames(hmExp),hmExp)
Error in is.data.frame(x) : 找不到对象'hmExp'
> write.table(diffExp,file="diffExp.txt",sep="\t",quote=F,col.names=F)
Error in is.data.frame(x) : 找不到对象'diffExp'
> 
> #volcano
> pdf(file="vol.pdf")
> xMax=max(-log10(allDiff$adj.P.Val))
> yMax=max(abs(allDiff$logFC))
> plot(-log10(allDiff$adj.P.Val), allDiff$logFC, xlab="-log10(adj.P.Val)",ylab="logFC",
+      main="Volcano", xlim=c(0,xMax),ylim=c(-yMax,yMax),yaxs="i",pch=20, cex=0.8)
Error in plot.window(...) : 'xlim'值不能是无限的
> diffSub=subset(allDiff, adj.P.Val<adjustP & logFC>logFoldChange)
> points(-log10(diffSub$adj.P.Val), diffSub$logFC, pch=20, col="red",cex=0.8)
> diffSub=subset(allDiff, adj.P.Val<adjustP & logFC<(-logFoldChange))
> points(-log10(diffSub$adj.P.Val), diffSub$logFC, pch=20, col="green",cex=0.8)
> abline(h=0,lty=2,lwd=3)
> dev.off()
null device 
          1 
请问有大神知道这个报错是什么意思吗 要怎么解决呢

[CSDN-Ada助手]

您可以前往 CSDN问答-其他 发布问题, 以便更快地解决您的疑问