Negative Binomial Modeling of District-Level Diarrhea Incidence in West Java Province

Nurul Nabilah, A'yunin Sofro

Abstract


This study investigated district-level diarrhea incidence in West Java Province during 2022–2024 and identified an appropriate count-data regression model for longitudinal epidemiological data. The response variable was the annual number of diarrhea cases, while explanatory variables included sanitation access, inpatient health centers, community health workers, nutritional indicators, drinking water facilities, and households implementing clean and healthy living behavior (PHBS). A population offset term was incorporated to account for differences in population exposure across districts/cities. Poisson Regression was initially fitted as a baseline model; however, severe overdispersion was detected. Therefore, Negative Binomial Regression and Generalized Linear Mixed Models (GLMM) were estimated and compared using Akaike Information Criterion (AIC), Bayesian Information Criterion (BIC), log-likelihood values, and simulation-based residual diagnostics. The results showed that Negative Binomial Regression provided the best overall fit, yielding substantially lower AIC and BIC values than competing models. Residual diagnostics indicated no evidence of remaining overdispersion, zero inflation, or serious outliers. Sanitation access was significantly associated with lower diarrhea incidence rates, with a one-standard-deviation increase corresponding to an estimated 15.7% reduction in incidence. These findings highlight the importance of accounting for overdispersion and suggest that Negative Binomial Regression is an appropriate framework for modeling regional diarrhea incidence data.

Keywords


diarrhea cases; Poisson regression; Poisson Mixed Regression; Regional random effects; Akaike Information Criterion

Full Text:

PDF

References


[1] Karen L. Kotloff et al. “The Burden and Etiology of Diarrheal Illness in Developing Countries”. The Lancet 382.9888 (2013), pp. 209–222. DOI: https://doi.org/10.1016/S0140-6736(13)60844-2.

[2] Frank Osei. “Spatial variation and hot-spots of district-level diarrhea incidences in Ghana: 2010–2014”. BMC Public Health 17 (2017). DOI: https://doi.org/10.1186/s12889-017-4541-z.

[3] A. C. Cameron and P. K. Trivedi. Regression Analysis of Count Data. 2nd ed. Cambridge University Press, 2013. DOI: https://doi.org/10.1017/CBO9780511814365.

[4] E. Oral. “Modified maximum likelihood estimation in Poisson regression”. Biometrics & Biostatistics International Journal 6.3 (2017). DOI: https://doi.org/10.15406/bbij.2017.06.00154.

[5] A. M. Alomair, F. Ayyaz, S. Tariq, and M. Ahsan-ul-Haq. “Discrete Poisson Haq distribution with mathematical properties and count data modeling”. Scientific Reports 15 (2025), p. 23281. DOI: https://doi.org/10.1038/s41598-025-07223-y.

[6] J. Hinde and C. G. B. Demétrio. “Overdispersion: Models and estimation”. Computational Statistics & Data Analysis 27.2 (1998), pp. 151–170. DOI: https://doi.org/10.1016/S0167-9473(98)00007-3.

[7] E. H. Payne, J. W. Hardin, L. E. Egede, V. Ramakrishnan, A. Selassie, and M. Gebregziabher. “Approaches for dealing with various sources of overdispersion in modeling count data: Scale adjustment versus modeling”. Statistical Methods in Medical Research 26.4 (2017), pp. 1802–1823. DOI: https://doi.org/10.1177/0962280215588569.

[8] A. M. Azevedo, Í. J. Silva, M. C. Nery, H. P. Rocha, and R. A. Santana. “Counting models for overdispersed data: A review with application to tuberculosis data”. Brazilian Journal of Biometrics 41.3 (2023), pp. 274–286. DOI: https://doi.org/10.28951/bjb.v41i3.635.

[9] Bijesh Yadav, Lakshmanan Jeyaseelan, Visalakshi Jeyaseelan, Jothilakshmi Durairaj, Sebastian George, K. G. Selvaraj, and Shrikant I. Bangdiwala. “Can Generalized Poisson model replace any other count data models? An evaluation”. Clinical Epidemiology and Global Health 11 (2021), p. 100774. DOI: https://doi.org/10.1016/j.cegh.2021.100774.

[10] E. Espinoza, U. Saputri, F. H. Fadilah, and D. Devianto. “Modeling the count data of public health service visits with overdispersion problem by using negative binomial regression”. Journal of Physics: Conference Series 1940.1 (2021), p. 012021. DOI: https://doi.org/10.1088/1742-6596/1940/1/012021.

[11] S. Q. Aina, Darnah, M. Fauziyah, and W. P. Nurmayanti. “Comparison of generalized Poisson regression and negative binomial regression models based on Akaike information criterion values”. STATISTIKA: Journal of Theoretical Statistics and Its Applications 25.1 (2025). DOI: https://doi.org/10.29313/statistika.v25i1.5402.

[12] Antonio J. Sáez-Castillo and Antonio Conde-Sánchez. “A hyper-Poisson regression model for overdispersed and underdispersed count data”. Computational Statistics & Data Analysis 61 (2013), pp. 148–157. DOI: https://doi.org/10.1016/j.csda.2012.12.009.

[13] N. E. Breslow and D. G. Clayton. “Approximate inference in generalized linear mixed models”. Journal of the American Statistical Association 88.421 (1993), pp. 9–25. DOI: https://doi.org/10.2307/2290687.

[14] Youngjo Lee and John A. Nelder. “Hierarchical Generalized Linear Models”. Journal of the Royal Statistical Society: Series B (Methodological) 58.4 (Nov. 1996), pp. 619–656. DOI: https://doi.org/10.1111/j.2517-6161.1996.tb02105.x.

[15] B. M. Bolker, M. E. Brooks, C. J. Clark, S. W. Geange, J. R. Poulsen, M. H. H. Stevens, and J. S. White. “Generalized linear mixed models: A practical guide for ecology and evolution”. Trends in Ecology & Evolution 24.3 (2009), pp. 127–135. DOI: https://doi.org/10.1016/j.tree.2008.10.008.

[16] Bertho Tantular, I Gede Nyoman Mindra Jaya, and Zulhanif. “Multilevel Models Approach on Longitudinal Studies for Disease Data”. Journal of Physics: Conference Series 1025.1 (2018), p. 012108. DOI: https://doi.org/10.1088/1742-6596/1025/1/012108.

[17] B. E. Yirdaw, L. K. Debusho, and A. Samuel. “Application of longitudinal multilevel zero inflated Poisson regression in modeling of infectious diseases among infants in Ethiopia”. BMC Infectious Diseases 24 (2024), p. 927. DOI: https://doi.org/10.1186/s12879-024-09820-0.

[18] Andrew B. Lawson. Bayesian Disease Mapping: Hierarchical Modeling in Spatial Epidemiology. 3rd ed. Boca Raton: Chapman and Hall/CRC, 2018. DOI: https://doi.org/10.1201/9781351271769.

[19] S. Xu, M. A. R. Ferreira, E. M. Porter, and C. T. Franck. “Bayesian model selection for generalized linear mixed models”. Biometrics 79.4 (2023), pp. 3266–3278. DOI: https://doi.org/10.1111/biom.13896.

[20] Peter J. Diggle, Patrick Heagerty, Kung-Yee Liang, and Scott L. Zeger. Analysis of Longitudinal Data. 2nd ed. Oxford: Oxford University Press, 2002. DOI: https://doi.org/10.1093/acprof:oso/9780198524847.001.0001.

[21] Geert Molenberghs and Geert Verbeke. Models for Discrete Longitudinal Data. 1st ed. Springer Series in Statistics. New York: Springer, 2005, p. 687. DOI: https://doi.org/10.1007/0-387-28980-1.

[22] H. Goldstein. Multilevel Statistical Models. John Wiley & Sons, 2011. DOI: https://doi.org/10.1002/9780470973394.

[23] Clifford M. Hurvich and Chih-Ling Tsai. “Regression and time series model selection in small samples”. Biometrika 76.2 (1989), pp. 297–307. DOI: https://doi.org/10.1093/biomet/76.2.297.

[24] N. M. Mangan, J. N. Kutz, S. L. Brunton, and J. L. Proctor. “Model selection for dynamical systems via sparse regression and information criteria”. Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences 473.2204 (2017). DOI: https://doi.org/10.1098/rspa.2017.0009.

[25] A. F. Zuur, E. N. Ieno, N. Walker, A. A. Saveliev, and G. M. Smith. Mixed Effects Models and Extensions in Ecology with R. Springer, 2009. DOI: https://doi.org/10.1007/978-0-387-87458-6.

[26] Jae H. Kim. “Multicollinearity and misleading statistical results”. Korean Journal of Anesthesiology 72.6 (2019), pp. 558–569. DOI: https://doi.org/10.4097/kja.19087.

[27] Christopher G. Thompson, R. S. Kim, Ariel M. Aloe, and Betsy J. Becker. “Extracting the variance inflation factor and other multicollinearity diagnostics from typical regression results”. Basic and Applied Social Psychology 39.2 (2017), pp. 81–90. DOI: https://doi.org/10.1080/01973533.2016.1277529.

[28] Román Salmerón, Catalina Garcia, and Jose Pérez. “Overcoming the inconsistences of the variance inflation factor: A redefined VIF and a test to detect statistical troubling multicollinearity”. arXiv (2020). DOI: https://doi.org/10.48550/arXiv.2005.02245.

[29] M. Edmondson, C. Luo, M. N. Islam, N. Sheils, J. Buresh, Z. Chen, J. Bian, and Y. Chen. “Distributed quasi-Poisson regression algorithm for modeling multi-site count outcomes in distributed data networks”. Journal of Biomedical Informatics 132 (2022), p. 104097. DOI: https://doi.org/10.1016/j.jbi.2022.104097.

[30] Abobaker Mohamed Jaber, Mohamed Amraja Mohamed, Aissa Omar Assrhani, Hanadi Abdullah Amhimmid, and Kasem Abdinibi Farag. “Robust Modeling of Over Dispersed Count Data Using an Outlier-Weighted Poisson Regression Approach”. Communications on Applied Nonlinear Analysis 32.10s (2025), pp. 3191–3200. DOI: https://doi.org/10.52783/cana.v32.5959.

[31] A. R. Jamaludin, F. Yusof, and Suhartono. “Modelling asthma cases by count analysis approach: Poisson INGARCH and negative binomial INGARCH”. Matematika 36.2 (2020), pp. 173–184. DOI: https://doi.org/10.11113/matematika.v36.n2.1217.

[32] D. Handayani, K. Notodiputro, K. Sadik, and A. Kurnia. “A comparative study of approximation methods for maximum likelihood estimation in generalized linear mixed models (GLMM)”. AIP Conference Proceedings. Vol. 1827. 2017, p. 020033. DOI: https://doi.org/10.1063/1.4979449.

[33] S. N. Wood. “Fast stable restricted maximum likelihood and marginal likelihood estimation of semiparametric generalized linear models”. Journal of the Royal Statistical Society: Series B (Statistical Methodology) 73.1 (2011), pp. 3–36. DOI: https://doi.org/10.1111/j.1467-9868.2010.00749.x.




DOI: https://doi.org/10.18860/cauchy.v11i2.42057

Refbacks

  • There are currently no refbacks.


Copyright (c) 2026 Nurul Nabilah, A'yunin Sofro

Creative Commons License
This work is licensed under a Creative Commons Attribution-ShareAlike 4.0 International License.

Editorial Office
Mathematics Department,
Maulana Malik Ibrahim State Islamic University of Malang
Gajayana Street 50 Malang, East Java, Indonesia 65144
e-mail: cauchy@uin-malang.ac.id

Creative Commons License
CAUCHY: Jurnal Matematika Murni dan Aplikasi is licensed under a Creative Commons Attribution-ShareAlike 4.0 International License.