Abstract
To understand the transmissibility and spread of infectious diseases, epidemiologists turn to estimates of the instantaneous reproduction number. While many estimation approaches exist, their utility may be limited. Challenges of surveillance data collection, model assumptions that are unverifiable with data alone, and computationally inefficient frameworks are critical limitations for many existing approaches. We propose a discrete spline-based approach that solves a convex optimization problem—Poisson trend filtering—using the proximal Newton method. It produces a locally adaptive estimator for instantaneous reproduction number estimation with heterogeneous smoothness. Our methodology remains accurate even under some process misspecifications and is computationally efficient, even for large-scale data. The implementation is easily accessible in a lightweight R package rtestim.
Competing Interest Statement
The authors have declared no competing interest.
Funding Statement
Daniel J. McDonald is partially supported by Centers for Disease Control and Prevention under U011P001121 and 75D30123C15907, and National Sciences and Engineering Research Council under ALLRP 581756-23 and RGPIN 2021-02618. Paul Gustafson is partially supported by National Sciences and Engineering Research Council Discovery Grant under RGPIN-2019-03957.
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
Footnotes
Here we highlight some of the major changes we have made to the manuscript. The first is to add a number of experiments to compare the performance of our method with alternatives under several scenarios: 1. investigating misspecified serial interval distributions (under mild and more major misspecification); 2. including all degrees ("k") of the discrete splines in our method across all settings, (thus investigating potential mismatches); 3. using relatively lower incident cases (<1000) in addition to the previous settings; 4. examining epidemics that are shorter (with 50 timepoints); and 5. directly investigating interval coverage. All experiments are now based on serial interval distributions for real epidemics (rather than synthetic distributions). We have included experimental results from EpiFilter and EpiNow2 wherever possible. Another key change is that we clarify that the "instantaneous reproduction number" is our focus. With this in mind, we have modified the title to misstating our intentions and also to be aligned more closely with existing related work. We also add more discussion on how the hyperparameters (degree "k" and tuning parameter "lambda") should be chosen in practice, the coverage of the approximate confidence interval, and a few future directions. We now use "rtestim" to refer to both of the method and software throughout the manuscript for consistency. Other changes include using time-varying serial interval distributions for Covid-19 incidence in Canada, discussing the level of overdispersion in the experiments, applying alternative methods to real epidemics for comparison to ours, wording clarification, and additional references.
Data Availability
All data produced are available online at the "jiapivialiu/rt-est-manuscript" Github repository with the following url: https://github.com/jiapivialiu/rt-est-manuscript/tree/main
