Abstract:Compartmental models that describe infectious disease transmission across subpopulations are central for assessing the impact of non-pharmaceutical interventions, behavioral changes and seasonal effects on the spread of respiratory infections. We present a Bayesian workflow for such models, including four features: (1) an adjustment for incomplete case ascertainment, (2) an adequate sampling distribution of laboratory-confirmed cases, (3) a flexible, time-varying transmission rate, and (4) a stratification by age group. Within the workflow, we benchmarked the performance of various implementations of two of these features (2 and 3). For the second feature, we used SARS-CoV-2 data from the canton of Geneva (Switzerland) and found that a quasi-Poisson distribution is the most suitable sampling distribution for describing the overdispersion in the observed laboratory-confirmed cases. For the third feature, we implemented three methods: Brownian motion, B-splines, and approximate Gaussian processes (aGP). We compared their performance in terms of the number of effective samples per second, and the error and sharpness in estimating the time-varying transmission rate over a selection of ordinary differential equation solvers and tuning parameters, using simulated seroprevalence and laboratory-confirmed case data. Even though all methods could recover the time-varying dynamics in the transmission rate accurately, we found that B-splines perform up to four and ten times faster than Brownian motion and aGPs, respectively. We validated the B-spline model with simulated age-stratified data. We applied this model to 2020 laboratory-confirmed SARS-CoV-2 cases and two seroprevalence studies from the canton of Geneva. This resulted in detailed estimates of the transmission rate over time and the case ascertainment. Our results illustrate the potential of the presented workflow including stratified transmission to estimate age-specific epidemiological parameters. The workflow is freely available in the R package HETTMO, and can be easily adapted and applied to other infectious diseases. Mathematical models are a central tool for understanding the spread of infectious diseases. These models can frequently be fitted to surveillance data such as the number of laboratory-confirmed cases and seroprevalence over time. We identified that in these situations, four crucial features are required for a model to provide insightful information for managing an epidemic. These features relate to the adjustment for incomplete case ascertainment, to the choice of sampling distribution, to the variation of transmission over time and to the stratification by age. For each feature, we identify and compare several implementation options on simulated data. This structural comparison of methods results in a Bayesian workflow that is optimized for modeling the transmission of SARS-CoV-2 over a short period. We highlight the advantages and limitations of our approach in a real situation, using real-world SARS-CoV-2 data from the canton of Geneva. In addition to providing validated solutions to important technical points, such a comprehensive workflow helps to improve the reliability and the transparency of epidemic models.

Estimating epidemic dynamics with genomic and time series data

A simulation-based approach for estimating the time-dependent reproduction number from temporally aggregated disease incidence time series data

Incorporating epidemiological data into the genomic analysis of partially sampled infectious disease outbreaks

Integrating Transmission Dynamics and Pathogen Evolution Through a Bayesian Approach

Accounting for reporting delays in real-time phylodynamic analyses with preferential sampling

Uncovering heterogeneous inter-community disease transmission from neutral allele frequency time series

Forecasting Epidemics Through Nonparametric Estimation of Time-Dependent Transmission Rates Using the SEIR Model

Bayesian workflow for time-varying transmission in stratified compartmental infectious disease transmission models

Real time analysis of epidemic data

Refining Estimation of the Instantaneous Reproduction Number During Early Pandemic Stages: Addressing Case-Reporting Variability and Serial Interval Uncertainty

Tracking Disease Outbreaks from Sparse Data with Bayesian Inference

Bayesian inference for the onset time and epidemiological characteristics of emerging infectious diseases

A New Framework and Software to Estimate Time-Varying Reproduction Numbers During Epidemics

Bayesian reconstruction of transmission trees from genetic sequences and uncertain infection times

Detecting changes in the transmission rate of a stochastic epidemic model

Fitting stochastic epidemic models to gene genealogies using linear noise approximation

Bayesian Inference of Reproduction Number from Epidemiological and Genetic Data Using Particle MCMC

Online Bayesian phylodynamic inference in BEAST with application to epidemic reconstruction

Bayesian estimation of transmission networks for infectious diseases

Model design for non-parametric phylodynamic inference and applications to pathogen surveillance

Estimation of the Basic Reproductive Number and Mean Serial Interval of a Novel Pathogen in a Small, Well-Observed Discrete Population