Premium
Statistical models for autocorrelated count data
Author(s) -
Nelson Kerrie P.,
Leroux Brian G.
Publication year - 2006
Publication title -
statistics in medicine
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.996
H-Index - 183
eISSN - 1097-0258
pISSN - 0277-6715
DOI - 10.1002/sim.2274
Subject(s) - generalized linear mixed model , statistics , markov chain monte carlo , random effects model , autoregressive model , autocorrelation , count data , mathematics , quasi likelihood , generalized linear model , linear model , mixed model , computer science , monte carlo method , poisson distribution , medicine , meta analysis
A generalized linear mixed model is an increasingly popular choice for the modelling of correlated, non‐normal responses in a regression setting. A number of methods are currently available for fitting a generalized linear mixed model including Monte‐Carlo Markov‐Chain maximum likelihood algorithms, approximate maximum likelihood (PQL), iterative bias correction, and others. Of interest in this paper is to compare the parameter estimation of the various methods in the modelling of a count data set, the incidence of polio in the USA over the period 1970–1983, using a longlinear generalized linear mixed model with an autoregressive correlation structure. Despite the fact that all of these methods are considered valid modelling techniques, we find that parameter estimates and standard errors differ substantially between analyses, particularly in the estimation of the parameters describing the random effects distribution. A small simulation study is helpful in understanding some of these differences. The methods lead to reasonably similar predictions for future observations, with small differences observed in some monthly counts. Copyright © 2005 John Wiley & Sons, Ltd.