Artigo Revisado por pares

Statistical models for autocorrelated count data

2005; Wiley; Volume: 25; Issue: 8 Linguagem: Inglês

10.1002/sim.2274

ISSN

1097-0258

Autores

Kerrie P. Nelson, Brian G. Leroux,

Tópico(s)

Statistical Distribution Estimation and Applications

Resumo

Statistics in MedicineVolume 25, Issue 8 p. 1413-1430 Research Article Statistical models for autocorrelated count data Kerrie P. Nelson, Corresponding Author Kerrie P. Nelson [email protected] Department of Statistics, University of South Carolina, U.S.A.Department of Statistics, University of South Carolina, 1523 Greene Street, Columbia, SC 29208, U.S.A.Search for more papers by this authorBrian G. Leroux, Brian G. Leroux Department of Biostatistics, University of Washington, Seattle, U.S.A.Search for more papers by this author Kerrie P. Nelson, Corresponding Author Kerrie P. Nelson [email protected] Department of Statistics, University of South Carolina, U.S.A.Department of Statistics, University of South Carolina, 1523 Greene Street, Columbia, SC 29208, U.S.A.Search for more papers by this authorBrian G. Leroux, Brian G. Leroux Department of Biostatistics, University of Washington, Seattle, U.S.A.Search for more papers by this author First published: 30 September 2005 https://doi.org/10.1002/sim.2274Citations: 17AboutPDF ToolsRequest permissionExport citationAdd to favoritesTrack citation ShareShare Give accessShare full text accessShare full-text accessPlease review our Terms and Conditions of Use and check box below to share full-text version of article.I have read and accept the Wiley Online Library Terms and Conditions of UseShareable LinkUse the link below to share a full-text version of this article with your friends and colleagues. Learn more.Copy URL Share a linkShare onEmailFacebookTwitterLinkedInRedditWechat Abstract 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. REFERENCES 1 Zeger SL. A regression model for time series of counts. Biometrika 1988; 75: 621–629. 2 Centers for Disease Control. Morbidity and Mortality Weekly Report Annual Summary: Reported Morbidity and Mortality in the U.S. 1970–1983; 18–31(54). 3 Chan KS, Ledolter J. Monte-Carlo EM estimation for time-series models involving counts. Journal of the American Statistical Association 1995; 92(429): 242–252. 4 McCulloch CE. Maximum likelihood algorithms for generalized linear mixed models. Journal of the American Statistical Association 1997; 92(437): 162–170. 5 Kuk AYC, Cheng YW. The Monte Carlo Newton–Raphson algorithm. Journal of Statistical Computation and Simulation 1997; 59: 233–350. 6 Kuk AYC, Cheng YW. Pointwise and functional approximations in Monte Carlo maximum likelihood estimation. Statistics and Computing 1999; 9: 91–99. 7 Breslow NE, Clayton DG. Approximate inference in generalized linear mixed models. Journal of the American Statistical Association 1993; 88(421): 9–25. 8 Oh MS, Lim YB. Bayesian analysis of time series Poisson data. Journal of Applied Statistics 2001; 28(2): 259–271. 9 Pinheiro JC, Bates DM. Approximations to the loglikelihood function in the nonlinear mixed effects model. Journal of Computational and Graphical Statistics 1995; 4: 12–35. 10 Kuk AYC. Asymptotically unbiased estimation in generalized linear models with random effects. Journal of the Royal Statistical Society, Series B 1995; 57: 395–407. 11 Nelson KP. Generalized linear mixed models: development and comparison of different estimation methods. Ph.D. Thesis, Department of Statistics, University of Washington, Seattle, Washington, 2002. 12 Liang K-Y, Zeger SL. Longitudinal data analysis using generalized linear mixed models. Biometrika 1986; 73: 13–22. 13 McCullagh P. Quasi-likelihood functions. Annals of Statistics 1983; 11: 59–67. 14 McGilchrist CA. Estimation in generalized mixed models. Journal of the Royal Statistical Society, Series B 1994; 56: 61–69. 15 Engel B, Keen A. A simple approach for the analysis of generalized linear mixed models. Statistica Neerlandica 1994; 48: 1–22. 16 Schall R. Estimation in generalized linear models with random effects. Biometrika 1991; 40: 719–727. 17 Wolfinger R, O'Connell M. Generalized linear mixed models: a pseudo-likelihood approach. Journal of Statistical Computation and Simulation 1993; 48: 233–243. 18 Centre for Multilevel Modelling. MLwiN. Institute of Education: London, UK. 19 Hoeschele I, Tier B. Estimation of variance components of threshold characters by marginal posterio modes and means via Gibbs sampling. Genetics Selection Evolution 1995; 27: 519–540. 20 Lin X, Breslow NE. Bias correction in generalized linear mixed models with multiple components of dispersion. Journal of the American Statistical Association 1996; 91: 1007–1016. Citing Literature Volume25, Issue830 April 2006Pages 1413-1430 ReferencesRelatedInformation

Referência(s)
Altmetric
PlumX