f(y, u_{jk}) e^{-\frac{1}{2}u_{jk}^2 - l_k}\right]} --></p>
<p><!-- %% where \eqn{N_{AGQ}} is the number of quadrature nodes, \eqn{q} is the -->
  <!-- %% number of levels of the random effect, the \eqn{w_k} and \eqn{l_k} are -->
  <!-- %% weights and log of the normal density associated with the \eqn{k}th -->
  <!-- %% quadrature node (\code{\link{GHrule}}), and \eqn{u_{jk}} is the value -->
  <!-- %% of the spherical random effect coefficient associated with the -->
  <!-- %% \eqn{j}th level of the groupi