Quasi-Monte Carlo estimation in generalized linear mixed models

Authors:
Jianxin Pan;Robin Thompson
Affiliations:
School of Mathematics, University of Manchester, P.O. Box 88, Manchester M60 1QD, UK;Biomathematics and Bioinformatics Division, Rothamsted Research, Harpenden, Hertfordshire AL5 2JQ, UK
Venue:
Computational Statistics & Data Analysis
Year:
2007

Citing 1
Cited 1

Latin supercube sampling for very high-dimensional simulations

ACM Transactions on Modeling and Computer Simulation (TOMACS) - Special issue on uniform random number generation

Restricted likelihood inference for generalized linear mixed models

Statistics and Computing

Quantified Score

Hi-index	0.03

Visualization

Abstract

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.