(y|u)} --></p>
<p><!-- %% which depends on the fixed effect and covariance parameters, -->
  <!-- %% \eqn{\beta} and \eqn{\theta}.  We define the likelihood itself as: --></p>
<p><!-- %% \deqn{f_{\theta,\beta}(y) = \int f_{\theta,\beta}(y|u) f(u)}{f(y) = integral(f(y|u)f(u))} --></p>
<p><!-- %% where \eqn{f(u)} is an independent multivariate normal.  It is not -->
  <!-- %% typically possible to evaluate this integral for GLMMs.  In -->
  <!-- %% \code{lme4} we approximate it using the Laplace approximation