Back

A mechanistic statistical model of dengue dynamics in an endemic region

Luna-Martinez, N.; Cruz-Rodriguez, E. X.; Bernal-Castro, E. A.

2026-09-03 epidemiology
10.64898/2026.09.01.26361961 medRxiv
Show abstract

Background Dengue is a major public health challenge, and predictive models are crucial for early warning systems. However, many current modeling practices rely exclusively on climatic factors or employ complex algorithms that lack the interpretability needed for informed public health decision-making. To address these shortcomings, we developed and validated a multidimensional, interpretable statistical model to predict monthly dengue incidence. Methodology/Principal Findings We used a Generalized Linear Mixed Model (GLMM) with a Negative Binomial distribution to analyze 14 years (2010-2023) of spatiotemporal data from 37 municipalities in Huila, Colombia, an endemic region. The model integrates non-linear and lagged effects of climatic, demographic, and socioeconomic factors. The final model underwent rigorous external validation on an independent test set (2021-2023). Our model demonstrated high predictive discrimination (R2 = 0.743, Spearman's {rho} = 0.657), accurately capturing the timing of epidemic outbreaks. Key findings include the identification of an optimal thermal window for transmission at 27-28{degrees}C, a threshold effect for precipitation above 800 mm, and a saturation dynamic in outbreak autocorrelation. Conclusions/Significance This mechanistically-informed statistical approach provides a robust and transparent tool for epidemiological surveillance, successfully balancing high predictive performance with the explanatory power needed for effective, data-driven public health interventions.

Matching journals

The top 6 journals account for 50% of the predicted probability mass.

50% of probability mass above

"Similar papers" are the closest papers from that journal in the model's embedding space. They show what the match is built on, but the ranking comes mostly from a classifier over the whole training set, not from these examples alone.