Inferring COVID-19 spreading rates and potential change points for case number forecasts

Introduction When faced with the outbreak of a novel epidemic like COVID-19, rapid response measures are required by individuals as well as by society as a whole to mitigate the spread of the virus. During this initial, time-critical period, neither the central epidemiological parameters, nor the effectiveness of measures like cancellation of public events, school closings, and social distancing are known.

Rationale As one of the key epidemiological parameters, we infer the spreading rate λ from confirmed COVID-19 case numbers at the example in Germany by combining Bayesian inference with an SIR (Susceptible-Infected-Recovered) model from compartmental epidemiology. Our analysis characterizes the temporal change of the spreading rate and, importantly, allows us to identify potential change points and to provide short-term forecast scenarios based on various degrees of social distancing. A detailed, educational description is provided in the accompanying paper, and the model, inference, and prediction are available on github. While we apply it to Germany, our approach can be readily adapted to any other country or region.

Results In Germany, political interventions to contain the outbreak were implemented in three steps over three weeks: Around March 9, large public events like soccer matches were cancelled. On March 16, schools and other educational institutions as well as many non-essential stores were closed. One week later, on March 23, a farreaching contact ban (“Kontaktsperre”), which includes the prohibition of even small public gatherings as well as the further closing of restaurants and non-essential stores, was imposed by the government authorities.

From the observed case numbers of COVID-19, we can quantify the impact of these measures on the spread (Fig. 1). As of April 10, we have evidence of the first change point in the spreading rate from λ₀ = 0.40 (95 % Confidence interval (CI: [0.33,0.49]) to λ₁ = 0.24 (CI: [0.20,0.28]), which occurred around March 8 (CI: March 5 to March 10). Moreover, we have evidence for a second change point to λ₂ = 0.15 (CI: [0.12,0.19]), which occurred around March 16 (CI: March 15 to March 18). Both changes in λ slowed the spread of the virus, but still imply exponential growth (Fig. 1, red and orange traces). To contain the disease spread, and turn from exponential growth to a decline of novel cases, a further decrease in λ is necessary. We have first indications that this transition has been reached by the third change-point around March 23 (CI: March 21 to March 25).

With the start of this third change point, λ takes approximately the critical value where the spreading rate λ balances the recovery rate µ, i.e. the effective growth rate λ* = λ− µ≈ 0 (Fig. 1, green traces). The case numbers in the coming week will provide more information on its precise value. Importantly, λ* = 0 presents the watershed between exponential growth or decay. Together with the delay of approximately two weeks between infection and first inference of λ*, any future intervention such as lifting restrictions therefore warrants careful consideration.

Our detailed analysis shows that, in the current phase, reliable short- and long-term forecasts are very difficult, if not impossible: In Fig. 1C,D already the three example scenarios quickly diverge from each other, and consequently span a huge range of future case numbers. Thus, any uncertainty on the magnitude of our social distancing in the past two weeks can have a major impact on the case numbers in the next two weeks. Beyond two weeks, the case numbers depend on our future behavior, for which we have to make explicit assumptions. We illustrate how the precise magnitude and timing of potential change points impact the forecast of case numbers (see Fig. 2, main paper).

Conclusions We developed a Bayesian framework to infer the spreading rate λ and the timing and magnitude of change points. Thereby, the efficiency of political and individual measures for social distancing and containment can be assessed in a timely manner. We find first evidence for a successive decrease of the spreading rate in Germany around March 9 and around March 16, which significantly reduced the magnitude of exponential growth, but was not sufficient to turn growth into decay. The development in the coming week will reveal the efficiency of the contact ban initiated on March 23. In general, our analysis code may help to infer the efficiency of measures taken in other countries and inform policy makers about tightening, loosening and selecting appropriate rules for containment.