特定于使用R 2.15.2的limma的错误消息

Question

我最近更新了R，之前使用limma包编写的代码现在返回一个错误。

代码是

classes = c(rep('Phenotype1',), rep('Phenotype 2',))
classes= as.factor(classes)
design=model.matrix(~classes-1)
colnames(design)=levels(classes) 
fit1= lmFit(exset, design)
cm=makeContrasts(Phenotype 1 - Phenotype 2, levels=design)
fit2= contrasts.fit(fit1,cm) 
fit3= eBayes(fit2)

对于新版本的R，可以使用以下命令

fit2= contrasts.fit(fit1,cm)

返回此错误消息...

Error in .Call("La_chol", as.matrix(x), PACKAGE = "base") : 
Incorrect number of arguments (1), expecting 2 for 'La_chol'

Answer 1

BioConductor是一个来自BioConductor项目的包，R的每个版本都有一个关联的limma版本。如果你更新了R而没有更新你的BioConductor包(包括limma)，可能会有一些不兼容的地方。您确定您使用的是最新版本的limma和BioConductor吗？