R
將非嵌套模型與 AIC 進行比較
假設我們必須使用 GLMM
mod1 <- glmer(y ~ x + A + (1|g), data = dat) mod2 <- glmer(y ~ x + B + (1|g), data = dat)
這些模型不是通常意義上的嵌套:
a <- glmer(y ~ x + A + (1|g), data = dat) b <- glmer(y ~ x + A + B + (1|g), data = dat)
所以我們不能
anova(mod1, mod2)
像我們想要的那樣做anova(a ,b)
。我們可以用 AIC 來說明哪個是最好的模型嗎?
AIC 可以與非嵌套模型一起應用。事實上,這是關於 AIC 的最廣泛的神話(誤解?)之一。看:
您必須注意的一件事是包含所有規範化常量,因為這些對於不同的(非嵌套)模型是不同的:
也可以看看:
在 GLMM 的上下文中,一個更微妙的問題是用於比較此類模型的 AIC 的可靠性(另見@BenBolker’s)。AIC 的其他版本在以下論文中進行了討論和比較: