R

將非嵌套模型與 AIC 進行比較

  • September 26, 2014

假設我們必須使用 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 的其他版本在以下論文中進行了討論和比較:

引用自:https://stats.stackexchange.com/questions/116935

comments powered by Disqus