 Fitting GLMMs via maximum likelihood (as via AIC) involves integrating over the random effects.In general, those integrals cannot be expressed in analytical form.Various approximate methods have been developed, but none has good properties for all possible models and data sets (e.g.