全部版块 我的主页
论坛 站务区 十一区 新手入门区
1287 0
2021-04-13
求大神解答,差异分析遇到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)}
}
该怎么修改?谢谢谢谢!!!!

二维码

扫码加我 拉你入群

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

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

栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群