我正在使用affy软件包处理从http://carrier.gnf.org/publications/CellCycle/获得的恶性疟原虫的affymetrix数据作为R中的原始CEL数据 并通过rma方法进行标准化。芯片描述文件也可从上述链接获得。 cdf包是使用" makecdfenv"创建的。封装
处理时我遇到了错误
> mydata
AffyBatch object
size of arrays=712x712 features (12 kb)
cdf=tsrimalariacdf (5337 affyids)
number of samples=7
number of genes=5337
annotation=tsrimalariacdf
notes=
> eset <-rma(mydata)
Error in exprs(object)[index, , drop = FALSE] : subscript out of bounds
感谢任何帮助
谢谢