Сообщение об ошибке dredge (пакет MuMIn): нет подходящего метода для logLik
Я столкнулся с ошибкой, о которой я не смог найти никакой информации, и я пока не смог ее выяснить.
Вот пример:
x1 <- rnorm(100,0,1); x2 <- rnorm(100,5,1); x3 <- rnorm(100,8,4)
x4 <- rnorm(100,3,3); x5 <- rnorm(100,1,2); x6 <- rnorm(100,98,4)
y <- x1 + 2*x2 + x3/4 + 4.5*x4 + 2*x5 + 4*x6 + rnorm(100,1,50)
library(glmmADMB)
gm<-glmmadmb(y~x1+x2+x3+x4+x5+x6, family="gaussian")
library(MuMIn)
test.admb.dredge <- dredge(gm, nmax=2)
Error in AICc(global.model, nmax = 2) :
no applicable method for 'logLik' applied to an object of class "c('double', 'numeric')"
Я попытался предоставить не ранжированный аргумент по умолчанию (например, "BIC" или "QAIC"), но получил то же сообщение об ошибке.
Вот мой sessionInfo
в случае, если это полезно:
R version 3.0.3 (2014-03-06)
Platform: x86_64-w64-mingw32/x64 (64-bit)
locale:
[1] LC_COLLATE=English_United States.1252
[2] LC_CTYPE=English_United States.1252
[3] LC_MONETARY=English_United States.1252
[4] LC_NUMERIC=C
[5] LC_TIME=English_United States.1252
attached base packages:
[1] stats graphics grDevices utils datasets methods base
other attached packages:
[1] glmmADMB_0.7.7 R2admb_0.7.10 MASS_7.3-30 MuMIn_1.9.26
loaded via a namespace (and not attached):
[1] grid_3.0.3 lattice_0.20-27 Matrix_1.1-2-2 nlme_3.1-115
[5] tools_3.0.3
1 ответ
Решение
Я думаю, что это "Thinko", который на самом деле не имеет ничего общего с glmmADMB
; скорее (1) вы должны использовать m.max=2
скорее, чем nmax=2
и (2) nmax=2
не дает очень четкое сообщение об ошибке.
Я получаю точно такую же ошибку, если я заменяю
gm2 <- lm(y~x1+x2+x3+x4+x5+x6)
dredge(gm2,nmax=2)
И все, кажется, работает нормально, если я использую dredge(gm,m.max=2)