A Bayesian modelling framework with model comparison for epidemics with super-spreading

The transmission dynamics of an epidemic are rarely homogeneous. Super-spreading events and super-spreading individuals are two types of heterogeneous transmissibility. Inference of super-spreading is commonly carried out on secondary case data, the expected distribution of which is known as the offspring distribution. However, this data is seldom available. Here we introduce a multi-model framework fit to incidence time-series, data that is much more readily available. The framework consists of five discrete-time, stochastic, branching-process models of epidemics spread through a susceptible population. The framework includes a baseline model of homogeneous transmission, a unimodal and a bimodal model for super-spreading events, as well as a unimodal and a bimodal model for super-spreading individuals. Bayesian statistics is used to infer model parameters using Markov Chain Monte-Carlo methods. Model comparison is conducted by computing Bayes factors, with importance sampling used to estimate the marginal likelihood of each model. This estimator is selected for its consistency and lower variance compared to alternatives. Application to simulated data from each model identifies the correct model for the majority of simulations and accurately infers the true parameters, such as the basic reproduction number. We also apply our methods to incidence data from the 2003 SARS outbreak and the Covid-19 pandemic caused by SARS-CoV-2. Model selection consistently identifies the same model and mechanism for a given disease, even when using different time series. Our estimates are consistent with previous studies based on secondary case data. Quantifying the contribution of super-spreading to disease transmission has important implications for infectious disease management and control. Our modelling framework is disease-agnostic and implemented as an R package, with potential to be a valuable tool for public health.