使用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

limma是BioConductor项目的一个软件包，每个版本的R都有一个BioConductor的相关版本。如果您在不更新BioConductor包（包括limma）的情况下更新了R，则可能存在一些不兼容性。你确定你使用的是最新版本的limma和BioConductor吗？