Powiedzmy, że mamy GLMM
mod1 <- glmer(y ~ x + A + (1|g), data = dat)
mod2 <- glmer(y ~ x + B + (1|g), data = dat)
Te modele nie są zagnieżdżone w zwykłym sensie:
a <- glmer(y ~ x + A + (1|g), data = dat)
b <- glmer(y ~ x + A + B + (1|g), data = dat)
więc nie możemy zrobić tego, anova(mod1, mod2)
co byśmy zrobili anova(a ,b)
.
Czy zamiast AIC możemy powiedzieć, który model jest najlepszy?