Sagen wir, wir müssen GLMMs
mod1 <- glmer(y ~ x + A + (1|g), data = dat)
mod2 <- glmer(y ~ x + B + (1|g), data = dat)
Diese Modelle sind nicht im üblichen Sinne verschachtelt:
a <- glmer(y ~ x + A + (1|g), data = dat)
b <- glmer(y ~ x + A + B + (1|g), data = dat)
Also können wir nicht so machen, anova(mod1, mod2)
wie wir es wollten anova(a ,b)
.
Können wir mit AIC sagen, welches das beste Modell ist?