 doubly iterative) a weighted normal mixed model with a working variate, is implemented by various commercial and open source statistical programs.Fitting GLMMs via maximum likelihood (as via AIC) involves integrating over the random effects.In general, those integrals cannot be expressed in analytical form.