按照一个例子编写的时候遇到了这样一个问题,求大神解答,具体代码如下:
> 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,谢谢啦!
|