 and -->
  <!-- %% adaptive Gauss-Hermite quadrature. --></p>
<p><!-- %% For canonical links, the Laplace approximation to the log-likelihood -->
  <!-- %% is: --></p>
<p><!-- %% \deqn{log f_{\theta,\beta}(y | u)-0.5|u|^2 - log|L_{\theta}|} -->
  <!-- %% {log f(y|u) - 0.5|u|^2 - log|L|} --></p>
<p><!-- %% where \eqn{|u|^2} is the squared length of the vector of spherical -->
  <!-- %% random effects, \eqn{u}, and \eqn{|L_{\theta}|}{|L|} is the -->
  <!-- %% determinant of the Cholesky factor for the cross-p