Pǿlya–Aeppli regression model for overdispersed count data. (August 2019)

• Record Type: Journal Article
• Authors: Borges, Patrick; Godoi, Luciana G.
• Is Part Of: Statistical modelling. Volume 19:Number 4(2019)
• Page Start: 362
• Page End: 385
• Publication Date: 2019-08
• Subjects: bootstrap -- EM algorithm -- Generalized linear models (GLM) -- overdispersion -- zero-inflated models; Linear models (Statistics) -- Periodicals; Mathematical models -- Periodicals; Modèles linéaires (Statistique) -- Périodiques; Modèles mathématiques -- Périodiques; Modèle statistique; Modèle linéaire; Modélisation statistique; Périodique électronique (Descripteur de forme); Ressource Internet (Descripteur de forme); 519.5011
• DOI: 10.1177/1471082X18766797
• Languages: English
• ISSNs: 1471-082X

Abstract

The log-linear Poisson model, characterized by linear variance function and a logarithmic relation between means and covariates, embedded in the exponential family regression framework provided by generalized linear models (GLM) is still the standard approach for analyzing count data responses with regression models. In practice, however, count data are often overdispersed and, thus, not conducive to Poisson regression. Therefore, the main goal of this article is to introduce a log-linear model based on the P∅ ′ lya–Aeppli (PA) distribution, which is an extension of the Poisson distribution by including a dispersion parameter ρ, to address the problem of overdispersion. Maximum likelihood (ML) estimation procedure is discussed as well as a test for determining the need for a PA regression over a standard Poisson regression. In addition, a simple EM-type algorithm for iteratively computing ML estimates is presented. In order to study departures from the error assumption as well as the presence of outliers, we perform residual analysis based on the standardized Pearson residuals. Furthermore, for different parameter settings and sample sizes, various simulations are performed. Finally, we also illustrated the new method on three real datasets, two of them are from biological researches and the other is from a violence study.