Analysis of deviance for GLMMs — anova.mm

With two or more fitted models, anova() performs a sequential likelihood-ratio comparison (like anova(glmer1, glmer2)). For a single model, fixed-effect tests are routed to drop1() (term LRTs), summary() (Wald z), or contrast() (custom Wald contrasts), which the GLMM contract supports directly.

Usage

# S3 method for class 'mm_glmm'
anova(object, ...)

Arguments

object: A fitted mm_glmm.
...: Additional fitted models to compare.

Value

A data frame of class c("mm_anova_comparison", "mm_glmm_comparison", "anova", "data.frame") shaped like anova(glmer1, glmer2): columns npar, AIC, BIC, logLik, -2*log(L) (lme4 before 2.0 called it deviance), Chisq, Df, Pr(>Chisq), rows named after the arguments. The underlying comparison (with the model formulas) is reachable as x$table.