roduct matrix for -->
  <!-- %% the random effects (\code{\link{getME}(., "L")}).  Note that in the -->
  <!-- %% non-canonical link case, the above equation is technically an -->
  <!-- %% approximation to the Laplace approximation. --></p>
<p><!-- %% For models with a single scalar random effect one may also obtain an -->
  <!-- %% approximate log-likelihood using adaptive Gauss-Hermite quadrature as: --></p>
<p><!-- %% \deqn{-\log|\bm L| - \frac{q}{2}\log(2\pi) + \sum_j^q \log \left[\sum_k^{N_{AGQ}} w_k 