Сообщение об ошибке, относящееся к лимме, используемой с использованием R 2.15.2

Я недавно обновил 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'

1 ответ

Решение

Limma - это пакет из проекта BioConductor, и каждая версия R имеет связанную версию BioConductor. Если вы обновили R без обновления пакетов BioConductor (включая limma), возможны некоторые несовместимости. Вы уверены, что используете последнюю версию Limma и BioConductor?

Другие вопросы по тегам