GEO多芯片差异基因分析时候出现以下错误:Error in rt[as.vector(diffSig[, 1]), ] :
only 0’s may be mixed with negative subscripts。合并了两个数据,数据数值都小于10,应该是已经取了log的。但是为什么会出现这个情况,求大神告知!!!!!!!!
logFoldChange=1
> adjustP=0.05
>
> library(limma)
> setwd("C:\\Users\\acer\\Desktop\\DEGs\\07.diff")
> rt=read.table("normalize.txt",sep="\t",header=T,check.names=F)
> rt=as.matrix(rt)
> rownames(rt)=rt[,1]
> exp=rt[,2:ncol(rt)]
> dimnames=list(rownames(exp),colnames(exp))
> rt=matrix(as.numeric(as.matrix(exp)),nrow=nrow(exp),dimnames=dimnames)
>
> #differential
> modType=c(rep("normal",6),rep("Endometrium"