Authors
Haoyu Wu, David A Stephens, Erica E M Moodie
Published in
Biometrics. Volume 82. Issue 4. Oct 07, 2026.
Abstract
Statistical modeling of infectious disease outbreaks with incidence data is crucial to understanding epidemic dynamics and informing public health interventions. Stochastic epidemic models (SEMs) account for the inherent stochasticity in epidemic evolution and are particularly useful in small outbreaks or low-prevalence settings. In SEMs, the epidemic is naturally represented as a continuous-time Markov jump process (MJP), while surveillance data are collected at discrete times and subject to underdetection. Fitting SEMs to partially observed epidemic trajectories requires integrating over all possible transition paths, which is often intractable. Individual-level data augmentation (DA) methods have been developed to address these challenges. However, existing approaches either rely on per-subject updates with limited scalability or impose typically unrealistic assumptions to simplify the complex dimensional change problem for state-space SEMs. We present an efficient DA Markov chain Monte Carlo (MCMC) framework for fitting continuous-time state-space SEMs to incidence data. Our method uses a block sampler to jointly update individual trajectories and proposes transition times using a discrete-time approximation to the latent MJP. Through simulation studies, we demonstrate that our method significantly improves MCMC mixing and scalability compared to the previous DA approach and remains computationally robust in scenarios where particle MCMC algorithms struggle. We apply the method to the 2014 Ebola outbreak in a rural village in Sierra Leone and the 2009 H1N1 outbreak in a remote community in Canada.
PMID:
42979708
Bibliographic data and abstract were imported from PubMed on 04 Nov 2026.
Read full publication at:
Please sign in
to see all details.
Advertisement
Stats
- Recommendations n/a n/a positive of 0 vote(s)
- Views 23
- Comments 0