# Poisson GLMM
count_mod <- galamm(
formula = y ~ lbas * treat + lage + v4 + (1 | subj),
data = epilep, family = poisson
)
# Model without interaction
count_mod0 <- galamm(
formula = y ~ lbas + treat + lage + v4 + (1 | subj),
data = epilep, family = poisson
)
# Model comparison
anova(count_mod, count_mod0)
Run the code above in your browser using DataLab