Disons que nous devons les GLMM
mod1 <- glmer(y ~ x + A + (1|g), data = dat)
mod2 <- glmer(y ~ x + B + (1|g), data = dat)
Ces modèles ne sont pas imbriqués dans le sens habituel de:
a <- glmer(y ~ x + A + (1|g), data = dat)
b <- glmer(y ~ x + A + B + (1|g), data = dat)
donc nous ne pouvons pas faire anova(mod1, mod2)
comme nous le ferions avec anova(a ,b)
.
Pouvons-nous utiliser AIC pour dire quel est le meilleur modèle à la place?