You are here: MIMS > EPrints
MIMS EPrints

2007.169: Quasi-Monte Carlo estimation in generalized linear mixed models

2007.169: Jianxin Pan and Robin Thompson (2006) Quasi-Monte Carlo estimation in generalized linear mixed models. Computational Statistics and Data Analysis, 51 (12). pp. 5765-5775. ISSN 0167-9473

Full text available as:

PDF - Requires a PDF viewer such as GSview, Xpdf or Adobe Acrobat Reader
272 Kb

DOI: 10.1016/j.csda.2006.10.003


Generalized linear mixed models (GLMMs) are useful for modelling longitudinal and clustered data, but parameter estimation is very challenging because the likelihood may involve high-dimensional integrals that are analytically intractable. Gauss–Hermite quadrature (GHQ) approximation can be applied but is only suitable for low-dimensional random effects. Based on the Quasi-Monte Carlo (QMC) approximation, a heuristic approach is proposed to calculate the maximum likelihood estimates of parameters in the GLMM. The QMC points scattered uniformly on the high-dimensional integration domain are generated to replace the GHQ nodes. Compared to the GHQ approximation, the proposed method has many advantages such as its affordable computation, good approximation and fast convergence. Comparisons to the penalized quasi-likelihood estimation and Gibbs sampling are made using a real dataset and a simulation study. The real dataset is the salamander mating dataset whose modelling involves six 20-dimensional intractable integrals in the likelihood.

Item Type:Article
Uncontrolled Keywords:Generalized linear mixed models; High-dimensional integrals; Maximum likelihood estimates; Quasi-Monte Carlo approximation; Salamander data analysis
Subjects:MSC 2000 > 62 Statistics
MIMS number:2007.169
Deposited By:Mrs Louise Healey
Deposited On:19 November 2007

Download Statistics: last 4 weeks
Repository Staff Only: edit this item