Abstract
We describe a novel approach to recovering the underlying parameters of the SIR dynamic epidemic model from observed data on case incidence or deaths. We formulate a discrete-time approximation to the original continuous-time model and search for the parameter vector that minimizes the standard least squares criterion function. We show that the gradient vector and matrix of second-order derivatives of the criterion function with respect to the parameters adhere to their own systems of difference equations and thus can be exactly calculated iteratively. Applying our new approach, we estimate four-parameter SIR models on two types of datasets: (1) daily reported cases of COVID-19 during the SARS-CoV-2 Omicron/BA.1 surges of December 2021 - March 2022 in New York City and Los Angeles County; and (2) weekly deaths from a plague outbreak on the Isle of Bombay during December 1905 - July 1906, originally studied by Kermack and McKendrick in their now-classic 1927 paper. The estimated parameters from the COVID-19 data suggest a duration of persistent infectivity beyond that reported in small-scale clinical studies of mostly symptomatic subjects. The estimated parameters from the plague data suggest that the Bombay outbreak was in fact driven by pneumonic rather than bubonic plague.
Competing Interest Statement
The authors have declared no competing interest.
Funding Statement
This study did not receive any funding.
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
Data Availability
The data and programs used in this work are publicly available on the Open Science Framework at https://osf.io/3b4hv/.