 Where formula_7 and formula_8 are the fixed effects design matrix, and fixed effects; formula_9 and formula_1 are the random effects design matrix and random effects.The complete likelihood,
has no general closed form, and integrating over the random effects is usually extremely computationally intensive.In addition to numerically approximating this integral(e.g.