TY - JOUR T1 - Research Article Summary: Inferring change points in the COVID-19 spreading reveals the effectiveness of interventions JF - medRxiv DO - 10.1101/2020.04.02.20050922 SP - 2020.04.02.20050922 AU - Jonas Dehning AU - Johannes Zierenberg AU - F. Paul Spitzner AU - Michael Wibral AU - Joao Pinheiro Neto AU - Michael Wilczek AU - Viola Priesemann Y1 - 2020/01/01 UR - http://medrxiv.org/content/early/2020/05/06/2020.04.02.20050922.abstract N2 - As COVID-19 is rapidly spreading across the globe, short-term modeling forecasts provide time-critical information for decisions on containment and mitigation strategies. A main challenge for short-term forecasts is the assessment of key epidemiological parameters and how they change when first interventions show an effect. By combining an established epidemiological model with Bayesian inference, we analyze the time dependence of the effective growth rate of new infections. Focusing on the COVID-19 spread in Germany, we detect change points in the effective growth rate that correlate well with the times of publicly announced interventions. Thereby, we can quantify the effect of interventions, and we can incorporate the corresponding change points into forecasts of future scenarios and case numbers. Our code is freely available and can be readily adapted to any country or region.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 interventions 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 of Germany by combining Bayesian inference based on Markov-Chain Monte-Carlo sampling with a class of SIR (Susceptible-Infected-Recovered) compartmental models from 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 description is provided in the accompanying paper, and the models, inference, and predictions are available on github. While we apply it to Germany, our approach can be readily adapted to other countries or regions.Results In Germany, 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 childcare facilities as well as many non-essential stores were closed. One week later, on March 23, a far-reaching contact ban (“Kontaktsperre”), which included 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 disease spread (Fig. 0). Based on our analysis, which includes data until April 21, we have evidence of three change points: the first changed the spreading rate from λ0 = 0.43 (95 % credible interval (CI: [0.35, 0.51])) to λ1 = 0.25 (CI: [0.20, 0.30]), and occurred around March 6 (CI: March 2 to March 9); the second change point resulted in λ2 = 0.15 (CI: [0.12, 0.20]), and occurred around March 15 (CI: March 13 to March 17). Both changes in λ slowed the spread of the virus, but still implied exponential growth (Fig. 0, red and orange traces). To contain the disease spread, and turn from exponential growth to a decline of new cases, a further decrease in λ was necessary. Our analysis shows that this transition has been reached by the third change point that resulted in λ3 = 0.09 (CI: [0.06, 0.12]) around March 23 (CI: March 20 to March 25).With this third change point, λ transitioned below the critical value where the spreading rate λ balances the recovery rate μ, i.e. the effective growth rate λ* = λ − μ ≈ 0 (Fig. 0, gray traces). Importantly, λ* = 0 presents the watershed between exponential growth or decay. Given the delay of approximately two weeks between an intervention and first inference of the induced changes in λ*, future interventions such as lifting restrictions warrant careful consideration.Our detailed analysis shows that, in the current phase, reliable short- and long-term forecasts are very difficult as they critically hinge on how the epidemiological parameters change in response to interventions: In Fig. 0 already the three example scenarios quickly diverge from each other, and consequently span a considerable 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. In the main paper we illustrate how the precise magnitude and timing of potential change points impact the forecast of case numbers (Fig. 2).Conclusions We developed a Bayesian framework to infer central epidemiological parameters and the timing and magnitude of intervention effects. Thereby, the efficiency of political and individual intervention measures for social distancing and containment can be assessed in a timely manner. We find evidence for a successive decrease of the spreading rate in Germany around March 6 and around March 15, which significantly reduced the magnitude of exponential growth, but was not sufficient to turn growth into decay. Our analysis also shows that a further decrease of the spreading rate occurred around March 23, turning exponential growth into decay. Future interventions and lifting of restrictions can be modeled as additional change points, enabling short-term forecasts for case numbers. 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.Competing Interest StatementThe authors have declared no competing interest.Funding StatementAll authors received support from the Max-Planck-Society. JD and PS acknowledge funding by SMARTSTART, the joint training program in computational neuroscience by the VolkswagenStiftung and the Bernstein Network. JZ received financial support from the Joachim Herz Stiftung. MW is employed at the Campus Institute for Dynamics of Biological Networks funded by the VolkswagenStiftung.Author DeclarationsAll relevant ethical guidelines have been followed; any necessary IRB and/or ethics committee approvals have been obtained and details of the IRB/oversight body are included in the manuscript.YesAll necessary patient/participant consent has been obtained and the appropriate institutional forms have been archived.YesI 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).YesI have followed all appropriate research reporting guidelines and uploaded the relevant EQUATOR Network research reporting checklist(s) and other pertinent material as supplementary files, if applicable.YesAll code is available on Github. The data is publicly available from the John Hopkins University. https://github.com/Priesemann-Group/covid19_inference_forecast https://github.com/CSSEGISandData/COVID-19 ER -