搜索
查看: 109|回复: 0

[expression-profile] R语言中在设计matrix时怎么进行多组分组,并分析差异基因...

[复制链接]

2

主题

4

帖子

28

积分

新手上路

Rank: 1

积分
28
发表于 2021-1-7 21:12:45 | 显示全部楼层 |阅读模式


  #f <-factor(pdata$GROUP)
  #design<-model.matrix(~0+f)
  #design<-model.matrix(~0+as.factor(pdata$GROUP))
  #colnames(design)<-c("CO","D15","D30","D5")
  #rownames(design)<-rownames(pdata)


fit<-lmFit(eset,design)
  cont.matrix<-makeContrasts(X,levels = design)
  fit2<-contrasts.fit(fit,cont.matrix)
  fit2e<-eBayes(fit2)


tested<-topTable(fit2e,adjust="fdr",number = Inf,sort.by = "B",coef = 1)

X那里应该怎么改呢?


请大神们指点~~

回复

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

QQ|手机版|小黑屋|生信技能树 ( 粤ICP备15016384号  

GMT+8, 2021-1-17 23:13 , Processed in 0.023471 second(s), 29 queries .

Powered by Discuz! X3.2

© 2001-2013 Comsenz Inc.