Abstract
The instantaneous reproduction number (Rt) is a widely used measure of the rate of spread of an infectious disease. Correct quantification of the uncertainty of Rt estimates is crucial for making well-informed decisions. Popular methods for estimating Rt leverage smoothing techniques to distinguish signal from noise. Examples include EpiEstim and EpiFilter, each are controlled by a single "smoothing parameter", which is traditionally chosen by the user. We demonstrate that the values of these smoothing parameters are unknown and vary markedly with epidemic dynamics. We argue that data-driven smoothing choices are crucial for accurately representing uncertainty about Rt estimates. We derive model likelihoods for the smoothing parameters in both EpiEstim and EpiFilter. Adopting a flexible Bayesian framework, we use these likelihoods to automatically marginalise out the relevant smoothing parameters from these models when fitting to incidence time-series. Applying our methods, we find that the default parameterisations of these models can negatively impact inferences of Rt, delaying detection of epidemic growth, and misrepresenting uncertainty (typically by producing overconfident estimates), with substantial implications for public health decision-making. Our extensions mitigate these issues, provide a principled approach to uncertainty quantification, and improve the robustness of inference of Rt in real-time.
Competing Interest Statement
The authors have declared no competing interest.
Funding Statement
N.S. acknowledges support from the Oxford-Radcliffe Scholarship from University College, Oxford, the EPSRC CDT in Modern Statistics and Statistical Machine Learning (Imperial College London and University of Oxford), and A. Maslov for studentship support.
Author Declarations
I confirm all relevant ethical guidelines have been followed, and any necessary IRB and/or ethics committee approvals have been obtained.
Yes
I confirm that all necessary patient/participant consent has been obtained and the appropriate institutional forms have been archived, and that any patient/participant/sample identifiers included were not known to anyone (e.g., hospital staff, patients or participants themselves) outside the research group so cannot be used to identify individuals.
Yes
I understand that all clinical trials and any other prospective interventional studies must be registered with an ICMJE-approved registry, such as ClinicalTrials.gov. I confirm that any such study reported in the manuscript has been registered and the trial registration ID is provided (note: if posting a prospective study registered retrospectively, please provide a statement in the trial ID field explaining why the study was not registered in advance).
Yes
I have followed all appropriate research reporting guidelines, such as any relevant EQUATOR Network research reporting checklist(s) and other pertinent material, if applicable.
Yes
Data Availability
All data used in this study is publicly available from the original sources. Additionally, all data and code used in our analysis is available at https://github.com/nicsteyn2/RobustRtEstimators.
