Generalized linear mixed models using adaptive Gaussian quadrature
