找回密码
 立即注册
查看: 3046|回复: 0

使用limma包时遇到的问题

[复制链接]
发表于 2016-8-17 11:53:26 | 显示全部楼层 |阅读模式
按照一个例子编写的时候遇到了这样一个问题,求大神解答,具体代码如下:
> library(gcrma)
> library(CLL)
> data(CLLbatch)
> data(disease)
> CLLbatch<-CLLbatch[,-match(c("CLL10.CEL","CLL1.CEL","CLL13.CEL"),sampleNames(CLLbatch))]
> CLLgcrma<-gcrma(CLLbatch)
> sampleNames(CLLgcrma)<-gsub(".CEL$","",sampleNames(CLLgcrma))
> disease<-disease[match(sampleNames(CLLgcrma),disease[,"SampleID"]),]
> eset<-exprs(CLLgcrma)
> disease<-factor(disease[,"Disease"])
> design<-model.matrix(~-1+disease)
> contrast.matrix<-makeContrasts(contrasts="diseaseprogres.-diseasestable",levels=design)
> fit<-lmFit(eset,design)
Error in lmFit(eset, design) :
  row dimension of design doesn't match column dimension of data object

请问我该如何解决这个error,谢谢啦!
回复

使用道具 举报

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

本版积分规则

Archiver|手机版|小黑屋|R语言中文网

GMT+8, 2024-11-24 19:28 , Processed in 0.024202 second(s), 16 queries .

Powered by Discuz! X3.5

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表