Reproducibility of COVID-era infectious disease models.

Infectious disease modelling has been prominent throughout the COVID-19 pandemic, helping to understand the virus' transmission dynamics and inform response policies. Given their potential importance and translational impact, we evaluated the computational reproducibility of infectious disease modelling articles from the COVID era. We found that four out of 100 randomly sampled studies released between January 2020 and August 2022 could be completely computationally reproduced using the resources provided (e.g., code, data, instructions) whilst a further eight were partially reproducible. For the 100 most highly cited articles from the same period we found that 11 were completely reproducible with a further 22 partially reproducible. Reflecting on our experience, we discuss common issues affecting computational reproducibility and how these might be addressed.

Replicating superspreader dynamics with compartmental models.

Infectious disease outbreaks often exhibit superspreader dynamics, where most infected people generate no, or few secondary cases, and only a small fraction of individuals are responsible for a large proportion of transmission. Although capturing this heterogeneity is critical for estimating outbreak risk and the effectiveness of group-specific interventions, it is typically neglected in compartmental models of infectious disease transmission-which constitute the most common transmission dynamic modeling framework. In this study we propose different classes of compartmental epidemic models that incorporate transmission heterogeneity, fit them to a number of real outbreak datasets, and benchmark their performance against the canonical superspreader model (i.e., the negative binomial branching process model). We find that properly constructed compartmental models can capably reproduce observed superspreader dynamics and we provide the pathogen-specific parameter settings required to do so. As a consequence, we also show that compartmental models parameterized according to a binary clinical classification have limited support.