悬赏 2 个论坛币 未解决
> for(i in row.names(TF1)){
+ if(sd(TF1[i,])>1){
+ for(j in row.names(immuneGene1)){
+ x=as.numeric(TF1[i,])
+ y=as.numeric(immuneGene1[j,])
+ corT=cor.test(x,y)
+ cor=corT$estimate
+ pvalue=corT$p.value
+ if((cor>corFilter) & (pvalue<pvalueFilter)){
+ outTab=rbind(outTab,cbind(TF=i,immuneGene=j,cor,pvalue,Regulation="postive"))
+ }
+ if((cor< -corFilter) & (pvalue<pvalueFilter)){
+ outTab=rbind(outTab,cbind(TF=i,immuneGene=j,cor,pvalue,Regulation="negative"))
+ }
+ }
+ }
+ }
Error in if (sd(TF1[i, ]) > 1) { : 需要TRUE/FALSE值的地方不可以用缺少值