Dynamic parameterization of a modified SEIRD model to analyze and forecast the dynamics of COVID-19 outbreaks in the United States

The rapid spread of the numerous outbreaks of the coronavirus disease 2019 (COVID-19) pandemic has fueled interest in mathematical models designed to understand and predict infectious disease spread, with the ultimate goal of contributing to the decision making of public health authorities. Here, we propose a computational pipeline that dynamically parameterizes a modified SEIRD (susceptible-exposed-infected-recovered-deceased) model using standard daily series of COVID-19 cases and deaths, along with isolated estimates of population-level seroprevalence. We test our pipeline in five heavily impacted states of the US (New York, California, Florida, Illinois, and Texas) between March and August 2020, considering two scenarios with different calibration time horizons to assess the update in model performance as new epidemiologic data become available. Our results show a median normalized root mean squared error (NRMSE) of 2.38% and 4.28% in calibrating cumulative cases and deaths in the first scenario, and 2.41% and 2.30% when new data are assimilated in the second scenario, respectively. Then, 2-week (4-week) forecasts of the calibrated model resulted in median NRMSE of cumulative cases and deaths of 5.85% and 4.68% (8.60% and 17.94%) in the first scenario, and 1.86% and 1.93% (2.21% and 1.45%) in the second. Additionally, we show that our method provides significantly more accurate predictions of cases and deaths than a constant parameterization in the second scenario (p < 0.05). Thus, we posit that our methodology is a promising approach to analyze the dynamics of infectious disease outbreaks, and that our forecasts could contribute to designing effective pandemic-arresting public health policies.