models and multilevel generalized linear mixed models for a binary response. Likelihood based inferences are hindered by the … SR method. Compared to Pinheiro and Chao's multilevel Adaptive Gaussian quadrature, our proposed method has an improved …

parameters in the multilevel model are estimated by a combination of EM and Newton-Raphson algorithms. Even for samples of 50 …