Laplace's approximation for nonlinear mixed models

An approximation to Laplace's method for integrals is applied to marginal distributions of data arising from models in which both fixed and random effects enter nonlinearly. The approach provides alternative derivations of some recent algorithms for fitting such models, and it has direct ties with Gaussian restricted maximum likelihood and the accompanying mixed model equations.