be approximated. The function works under the
assumption of normally distributed random effects with mean zero and variance-covariance
matrix \(D\). These integrals are approximated numerically using an adaptive
Gauss-Hermite quadrature rule. Using the control argument <code>nAGQ</code>, the user can
specify the number of quadrature points used in the approximation.</p>
<p><b>User-defined family:</b> The user can define its own family object; for an example,
see the help page of <code><a href="negative_bino