求大神解答，TCGA数据使用R语言差异分析遇到50+警告信息：

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

求大神解答，差异分析遇到50+警告信息：
1: In if (dim(GenesExpDiff) == 0) { ... : 条件的长度大于一，因此只能用其第一元素
2: In if (dim(GenesExpDiff) == 0) { ... : 条件的长度大于一，因此只能用其第一元素
3: In if (dim(GenesExpDiff) == 0) { ... : 条件的长度大于一，因此只能用其第一元素
4: In if (dim(GenesExpDiff) == 0) { ... : 条件的长度大于一，因此只能用其第一元素
5: In if (dim(GenesExpDiff) == 0) { ... : 条件的长度大于一，因此只能用其第一元素
6: In if (dim(GenesExpDiff) == 0) { ... : 条件的长度大于一，因此只能用其第一
代码如下：
library("limma")
for(i in row.names(GenesExp)){
  rt <- rbind(expression=GenesExp[i,],grade=grade)
  rt <- as.matrix(t(rt))
  wilcoxTest <- wilcox.test(expression ~ grade, data=rt)
  conGeneMeans <- mean(GenesExp[i,1:conNum])
  treatGeneMeans <- mean(GenesExp[i,(conNum+1):ncol(GenesExp)])
  logFC <- log2(treatGeneMeans)-log2(conGeneMeans)  
  pvalue <- wilcoxTest$p.value
  conMed <- median(GenesExp[i,1:conNum])
  treatMed <- median(GenesExp[i,(conNum+1):ncol(GenesExp)])
  diffMed <- treatMed-conMed
  exp <- data.frame(gene=i,conMean=conGeneMeans,treatMean=treatGeneMeans,logFC=logFC,pValue=pvalue)
  if (dim(GenesExpDiff) == 0){
    GenesExpDiff <- exp
  }else {GenesExpDiff <- rbind(GenesExpDiff,exp)}
}
该怎么修改？谢谢谢谢!!!!

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：警告信息 警告信 TCG CGA R语言