Flexible Bayesian modeling of non-equidispersed counts with penalized complexity priors in disease incidence studies

Abstract

Counts in epidemiology often deviate from equidispersion and exhibit spatial, temporal, and nonlinear structure that the Poisson model cannot accommodate. We introduce a gamma-count structured additive regression model that strategically integrates penalized complexity priors in two critical aspects: (i) a principled penalized complexity prior on the dispersion parameter of the gamma-count distribution, which naturally shrinks toward the base Poisson model when the data support equidispersion, and (ii) scale-dependent penalized complexity hyperpriors on the smoothing variances for nonlinear, spatial, and temporal effects. By formulating the model within a latent Gaussian framework, we enable efficient approximate Bayesian inference through integrated nested Laplace approximations. Simulation studies across under-, equi-, and over-dispersed regimes show that the penalized complexity prior for dispersion parameter combined with scale-dependent hyperpriors yields accurate estimation of dispersion and smooth effects, favorable predictive scores, and robust inference. In empirical applications to larynx cancer mortality in Germany, COVID-19 incidence in Georgia (USA), and lung and bronchus cancer in Iowa (USA), the gamma-count structured additive regression model exhibits competitive or enhanced fit relative to Poisson and negative binomial counterparts, while elucidating interpretable nonlinear and spatial structures. This framework delivers robust, spatially resolved estimates of disease burden in the presence of non-equidispersion, thereby facilitating evidence-based resource allocation, epidemiological surveillance, and monitoring of health disparities, contributing to Sustainable Development Goals (SDGs) 3 (Good Health and Well-Being) and 10 (Reduced Inequalities). For geographically targeted analyses, it further supports informed decision-making in urban and community planning, aligning with SDG 11 (Sustainable Cities and Communities).

Keywords

Gamma-count structured additive (GCSA) regression penalized complexity (PC) priors dispersion prior scale-dependent hyperpriors INLA non-equidispersed counts spatial epidemiology disease mapping health inequalities,Sustainable Development Goals (SDGs)

Get full access to this article

View all access options for this article.

References

Fahrmeir

Kneib

Lang

. Penalized structured additive regression for space-time data: a Bayesian perspective. Stat Sin 2004; 14: 731–761.

Kneib

Hothorn

Tutz

. Variable selection and model choice in geoadditive regression models. Biometrics 2009; 65: 626–634.

Breslow

Clayton

. Approximate inference in generalized linear mixed models. J Am Stat Assoc 1993; 88: 9–25.

Baetschmann

Winkelmann

. A dynamic hurdle model for zero-inflated count data: With an application to health care utilization. SOEPpapers on Multidisciplinary Panel Data Research 648, DIW Berlin, The German Socio-Economic Panel (SOEP). https://ideas.repec.org/p/diw/diwsop/diw_sp648.html, 2014.

Cameron

Trivedi

. Regression Analysis of Count Data. Cambridge: Cambridge University Press, 2013.

Ridout

Besbeas

. An empirical model for underdispersed count data. Stat Modell 2004; 4: 77–89.

Lord

Geedipally

Guikema

. Extension of the application of Conway-Maxwell-Poisson models: analyzing traffic crash data exhibiting underdispersion. Risk Anal: An Off Publication Soc Risk Anal 2010; 30: 1268–1276.

Lord

Guikema

Geedipally

. Application of the Conway–Maxwell–Poisson generalized linear model for analyzing motor vehicle crashes. Accident Anal Preven 2008; 40: 1123–1134.

Zhu

Joe

. Modelling heavy-tailed count data using a generalised Poisson-inverse Gaussian family. Stat Probab Lett 2009; 79: 1695–1703.

10.

Zeviani

Ribeiro Jr

Bonat

, et al. The gamma-count distribution in the analysis of experimental underdispersed data. J Appl Stat 2014; 41: 2616–2626.

11.

Cox

. Renewal Theory. London: Methuen, 1962.

12.

Winkelmann

. Duration dependence and dispersion in count-data models. J Bus Econ Stat 1995; 13: 467–474.

13.

Toft

Innocent

Mellor

, et al. The gamma-Poisson model as a statistical method to determine if micro-organisms are randomly distributed in a food matrix. Food Microbiol 2006; 23: 90–94.

14.

McShane

Adrian

Bradlow

, et al. Count models based on Weibull interarrival times. J Bus Econ Stat 2008; 26: 369–378.

15.

Gonzales-Barron

Butler

. Characterisation of within-batch and between-batch variability in microbial counts in foods using Poisson-gamma and Poisson-lognormal regression models. Food Control 2011; 22: 1268–1278.

16.

Nadifar

Baghishani

Fallah

. A flexible Bayesian non-confounding spatial model for analysis of dispersed count data. Biomet J 2022; 64: 758–770.

17.

Nadifar

Baghishani

Fallah

. A flexible generalized Poisson likelihood for spatial counts constructed by renewal theory, motivated by groundwater quality assessment. J Agric Biol Environ Stat 2023; 28: 726–748.

18.

Simpson

Rue

Riebler

, et al. Penalising model component complexity: a principled, practical approach to constructing priors. Stat Sci 2017; 32: 1–28.

19.

Klein

Kneib

. Scale-dependent priors for variance parameters in structured additive distributional regression. Bayesian Anal 2016; 11: 1071–1106.

20.

Rue

Martino

Chopin

. Approximate Bayesian inference for latent Gaussian models by using integrated nested Laplace approximations. J R Stat Soc Ser B 2009; 71: 319–392.

21.

Bakka

Rue

Fuglstad

, et al. Spatial modeling with R-INLA: a review. Wiley Interdiscip Rev: Comput Stat 2018; 10: e1443.

22.

Barmoudeh

Baghishani

Martino

. Bayesian spatial analysis of crash severity data with the INLA approach: assessment of different identification constraints. Accident Anal Preven 2022; 167: 106570.

23.

Muff

Riebler

Held

, et al. Bayesian analysis of measurement error models using integrated nested Laplace approximations. J R Stat Soc: Ser C (Appl Stat) 2015; 64: 231–252.

24.

Rue

Riebler

Sørbye

, et al. Bayesian computing with INLA: a review. Annu Rev Stat Appl 2017; 4: 395–421.

25.

Van Niekerk

Bakka

Rue

. Competing risks joint models using R-INLA. Stat Modell 2021; 21: 56–71.

26.

Winkelmann

. Econometric Analysis of Count Data. Berlin, Heidelberg: Springer-Verlag, 2013.

27.

Ong

Biswas

Peiris

, et al. Count distribution for generalized Weibull duration with applications. Commun Stat - Theory Methods 2015; 44: 4203–4216.

28.

Fahrmeir

Kneib

Lang

, et al. Regression: Models, Methods and Applications. Berlin, Heidelberg: Springer, 2013.

29.

Fahrmeir

Osuna Echavarría

. Structured additive regression for overdispersed and zero-inflated count data. Appl Stoch Models Bus Ind 2006; 22: 351–369.

30.

Gelman

. Analysis of variance—why it is more important than ever. Ann Stat 2005; 33: 1–53.

31.

Gelman

. Prior distributions for variance parameters in hierarchical models (comment on article by Browne and Draper). Bayesian Anal 2006; 1: 515–534.

32.

Hodges

. Richly parameterized linear models: additive, time series, and spatial models using random effects. USA: CRC Press, 2013.

33.

Abramowitz

Stegun

. Handbook of mathematical functions with formulas, graphs, and mathematical tables, Vol. 55. New York: US Government Printing Office, 1964.

34.

Youssef

Delpha

Diallo

. Analytical model of the KL divergence for gamma distributed data: application to fault estimation. In: 2015 23rd European signal processing conference (EUSIPCO), 2015, pp.2266–2270. IEEE.

35.

Klein

. sdprior: Scale-dependent hyperpriors in structured additive distributional regression. R package version 0.3, 2015.

36.

Sørbye

Rue

. Scaling intrinsic Gaussian Markov random field priors in spatial modelling. Spat Stat 2013; 8: 39–51.

37.

Gelman

Jakulin

Pittau

, et al. A weakly informative default prior distribution for logistic and other regression models. Ann Appl Stat 2008; 2: 1360–1383.

38.

Polson

Scott

. On the half-Cauchy prior for a global scale parameter. Bayesian Anal 2012; 7: 887–902.

39.

Spiegelhalter

Best

Carlin

, et al. Bayesian measures of model complexity and fit (with discussion). J R Stat Soc Ser B 2002; 64: 1–34.

40.

Gelman

Hwang

Vehtari

. Understanding predictive information criteria for Bayesian models. Stat Comput 2013; 24: 997–1016.

41.

Watanabe

. A widely applicable Bayesian information criterion. J Mach Learn Res 2012; 14: 867–897.

42.

Pettit

. The conditional predictive ordinate for the normal distribution. J R Stat Soc: Ser B (Methodol) 1990; 52: 175–184.

43.

Adrion

Mansmann

. Bayesian model selection techniques as decision support for shaping a statistical analysis plan of a clinical trial: an example from a vertigo phase III study with longitudinal count data as primary endpoint. BMC Med Res Methodol 2012; 12: 137.

44.

Klein

Kneib

Lang

. Bayesian generalized additive models for location, scale, and shape for zero-inflated and overdispersed count data. J Am Stat Assoc 2015; 110: 405–419.

45.

Sun

Tsutakawa

. Propriety of posteriors with improper priors in hierarchical linear mixed models. Stat Sin 2001; 77–95.

46.

Ghosal

van der Vaart

. Fundamentals of nonparametric Bayesian inference, Vol. 44. Cambridge: Cambridge University Press, 2017.

47.

Natário

Knorr-Held

. Non-parametric ecological regression and spatial variation. Biomet J: J Math Methods Biosci 2003; 45: 670–688.

48.

Besag

York

Mollié

. Bayesian image restoration, with two applications in spatial statistics. Ann Inst Stat Math 1991; 43: 1–20.

49.

Johnson

Ravi

Braneon

. Spatiotemporal associations between social vulnerability, environmental measurements, and COVID-19 in the conterminous united states. GeoHealth 2021; 5: e2021GH000423.

Supplementary Material

Please find the following supplemental material available below.

For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.

For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.

5.42 MB

0.00 MB