A spatiotemporal negative binomial model with dynamic dispersion: An application to Tuberculosis infections
Abstract
Tuberculosis (TB) remains a critical public health concern in Brazil, characterized by pronounced spatial heterogeneity and fluctuating temporal volatility. In this paper, we study monthly TB notifications across 61 microregions of Sao Paulo state from 2001 to 2024. To do this, we introduce a negative binomial spatial integer-valued generalized autoregressive conditional heteroskedastic (INGARCH) model featuring jointly dynamic conditional means and time-varying dispersion. To capture inter-regional spillovers, we incorporate both discrete adjacency structures and a novel continuous distance-based formulation leveraging the Matern correlation function. Parameter estimation via conditional maximum likelihood employs a two-step profile-likelihood iterative scheme, demonstrating solid finite-sample performance in simulation studies. Applied to the Sao Paulo TB surveillance data, the framework substantially outperforms standard Poisson and fixed-dispersion spatiotemporal baselines in empirical fit and uncertainty quantification, maintaining nominal 95% predictive coverage across both dense metropolitan centers and rural microregions. Our results reveal marked spatial heterogeneity in baseline incidence, dynamic overdispersion driven by localized outbreaks, and short-range spatial interaction decay. By accurately modeling spatiotemporal volatility, the proposed methodology provides a robust statistical tool to support public health surveillance, policy-making, and resource allocation.
Related
- Scaling-up tuberculosis preventive therapy to HIV-negative contacts of bacteriologically confirmed TB cases aged 5 years and above – Lessons from a tuberculosis preventive therapy surge activity in Zambia
- Singularities of Non-negative Matrix Factorization and their application to Bayesian inference
- Nonparametric Estimation of Negative Cumulative Extropy for Linear Consecutive Systems: Reliability Properties and Image-Processing Application
- DISPERSION OF MAGNETIC CORIOLIS WAVES
- THE ZERO-DISPERSION LIMIT OF THE PERIODIC BENJAMIN-ONO EQUATION FOR SQUARE-INTEGRABLE DATA
BibTeX
@misc{silva2026spatiotemporal,
title = {{A spatiotemporal negative binomial model with dynamic dispersion: An application to Tuberculosis infections}},
author = {Rodrigo B. Silva and Luiza S. C. Piancastelli and Wagner Barreto-Souza},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2609.10174},
eprint = {2609.10174},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2609.10174},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - A spatiotemporal negative binomial model with dynamic dispersion: An application to Tuberculosis infections
AU - Rodrigo B. Silva
AU - Luiza S. C. Piancastelli
AU - Wagner Barreto-Souza
PY - 2026
DA - 2026/09/09
PB - arXiv
DO - 10.48550/arxiv.2609.10174
UR - https://arxiv.org/abs/2609.10174
AB - Tuberculosis (TB) remains a critical public health concern in Brazil, characterized by pronounced spatial heterogeneity and fluctuating temporal volatility. In this paper, we study monthly TB notifications across 61 microregions of Sao Paulo state from 2001 to 2024. To do this, we introduce a negative binomial spatial integer-valued generalized autoregressive conditional heteroskedastic (INGARCH) model featuring jointly dynamic conditional means and time-varying dispersion. To capture inter-regional spillovers, we incorporate both discrete adjacency structures and a novel continuous distance-based formulation leveraging the Matern correlation function. Parameter estimation via conditional maximum likelihood employs a two-step profile-likelihood iterative scheme, demonstrating solid finite-sample performance in simulation studies. Applied to the Sao Paulo TB surveillance data, the framework substantially outperforms standard Poisson and fixed-dispersion spatiotemporal baselines in empirical fit and uncertainty quantification, maintaining nominal 95% predictive coverage across both dense metropolitan centers and rural microregions. Our results reveal marked spatial heterogeneity in baseline incidence, dynamic overdispersion driven by localized outbreaks, and short-range spatial interaction decay. By accurately modeling spatiotemporal volatility, the proposed methodology provides a robust statistical tool to support public health surveillance, policy-making, and resource allocation.
ER -