A simple model to fit the time evolution of the daily death rate of Covid-19 in European Union countries ======================================================================================================== * Tristan Beau * Julien Browaeys * Olivier Dadoun ## Abstract We suggest a minimal model to describe the evolution of daily deaths due to Covid-19 in the European Union (EU) and UK, without any epidemiological hypothesis. Assuming current lockdown conditions were to remain in place, as of May 3rd 2020, this *ad hoc* fitting model is forecasting a total of 204 586 deaths in the EU. We could currently be at 2/3 of the total casualty count, and reach an overall death rate of one over 2 500 people. ## I. — Introduction The current Covid-19 pandemia is a major worldwide concern since its outbreak in December 2019[1]. Although all countries are publishing their statistics in term of deaths, confirmed or recovered cases, the epidemiological modeling remains complex. In this study, we propose an modelization of data without any epidemiological assumptions. The realization of a major epidemic crisis happening in the European Union (EU) came in March 2020, and most EU countries consequently instituted some form of lockdown in a matter of days. Since then, at least in Europe, the epidemics has receded enough so that these countries have planned to lift the lockdown, at least partially, in the second week of May. We wish to contribute to the evaluation of the incidence of such a decision based on current data. We fit an effective physical model to public Covid-19 deaths data in all 28 EU countries (including UK) in order to forecast the evolution of the epidemic if lockdown was maintained. In that case, the total number of deaths is estimated. ## II. — Data and method The John Hopkins University database[5] provides a daily report of deaths, confirmed and recovered Covid-19 patients all over the world, since January 22nd, 2020; this date is considered to be our time reference *t* = 0. Because of large variations in testing policies in different countries, in this analysis we limit ourselves to the study of daily deaths as a function of time. We are aware that this remains an imperfect proxy of the epidemic as counting policies can also vary from one country to the other. In order to be able to compare countries, we normalize data by country population. This information is taken from the latest *United Nations Population Division* estimates [4]. Databases extraction of data and their initial analysis is ensured by the open source software called CoCoA[2, 3]. ## III. — Analysis and extrapolation The effective modeling of data is done by a standard least-square fit of a Crystal Ball function[6]. ![Formula][1] The Crystal Ball function consists of a smoothly joined Gaussian part and a power-law tail at a certain cut-off. It has 5 parameters: *N* is the full integral of the function over the time, *n* the power law parameter, *α* a normalized cutoff where the power-law replaces the Gaussian part, *t* the position of the peak, *σ* the time scale. We are only reporting the useful parameters *N* and *t*, which respectively describe the total number of death (assuming the daily death rate to follow this function) and the time at which the peak of the epidemic occurred. Commonly used in high-energy physics, this function is named after a neutral particle detector initially used at the Stanford Linear Accelerator Center. The reason for choosing this function is two-fold. First we needed to describe properly the asymmetrical bell shape curve of the data (see figures below). Second we wanted a minimum number of adjustable parameters to avoid over-fitting. For every EU countries (and UK), we represent daily deaths as a function of time since *t*=0, together with the fitted Crystal Ball function. It is represented over 200 days to better visualize our forecast. The best fit parameters *N* and *t* are given on each figure. ![Figure1](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/10/2020.05.06.20093062/F1.medium.gif) [Figure1](http://medrxiv.org/content/early/2020/05/10/2020.05.06.20093062/F1) ![Figure2](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/10/2020.05.06.20093062/F2.medium.gif) [Figure2](http://medrxiv.org/content/early/2020/05/10/2020.05.06.20093062/F2) ![Figure3](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/10/2020.05.06.20093062/F3.medium.gif) [Figure3](http://medrxiv.org/content/early/2020/05/10/2020.05.06.20093062/F3) The same procedure has also been applied to European Union (including UK) as a whole. ![Figure4](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/10/2020.05.06.20093062/F4.medium.gif) [Figure4](http://medrxiv.org/content/early/2020/05/10/2020.05.06.20093062/F4) ## IV. — Discussion Our modeling yields an estimation of the expected overall number of deaths for each country. We provide these numbers in the table below, and compare them with the total number deaths reported until May 3rd, 2020, to infer the progress of the epidemic at this date. We normalize the figures by population to facilitate comparison between countries. View this table: [Table1](http://medrxiv.org/content/early/2020/05/10/2020.05.06.20093062/T1) With merely five free parameters, our model fits reasonably well all data for the 28 EU countries (and UK). Fitting seems all the more accurate that the peak can be easily identified, and that the total number of daily death is high. Some countries plan lifting the lockdown in early May 2020, such as Germany (May 4th), or France (May 11th). One can assume that without lifting the lockdown, the evolution of the daily rate of death could have followed the trend we have exhibited. On May 3rd, EU countries (and UK) have recorded a total of 136 391 deaths due to Covid-19. Our effective model is projecting 204 586 deaths if all conditions remain as current. In this hypothesis, 2/3 of the total number of Covid-19 victims have already perished. It will be interesting to compare these figures with the total number of death at the end of the epidemic. We should stress that our proposal is prospective and not predictive, in the sense that it does *not* rely on epidemiological modeling. We merely describe data through *ad hoc* fitting of mathematical functions, a methodology similar to what is used in economic forecasts. ## V. — Conclusion We provide an minimal model, within the CoCoA framework, that reasonably fits the daily rate of death as a function of time with only five free parameters. Assuming lockdown is not lifted, we estimate that at the beginning of May 2020 about two third of deaths have already occurred in Europe before the end of the epidemic. This implies a final total death rate of 1 over 2 500 people. In the future, lockdown lifting policies could be evaluated with regard to our projection. ## Data Availability The data that support the findings of this study are available in CSSEGISandData/COVID-19 by the Center for Systems Science and Engineering (CSSE) at Johns Hopkins University, at [[https://github.com/CSSEGISandData/COVID-19](https://github.com/CSSEGISandData/COVID-19)]. And also in Population by Country (2020) - Worldometer, at [[https://www.worldometers.info/world-population/population-by-country/](https://www.worldometers.info/world-population/population-by-country/)]. These data were derived from the following resources available in the public domain: [[https://github.com/tjbtjbtjb/CoCoA](https://github.com/tjbtjbtjb/CoCoA)]. * Received May 6, 2020. * Revision received May 6, 2020. * Accepted May 10, 2020. * © 2020, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution 4.0 International), CC BY 4.0, as described at [http://creativecommons.org/licenses/by/4.0/](http://creativecommons.org/licenses/by/4.0/) ## References 1. [1].“About the Virus.” World Health Organization, May 4, 2020. [http://www.euro.who.int/en/health-topics/health-emergencies/coronavirus-covid-19/novel-coronavirus-2019-ncov](http://www.euro.who.int/en/health-topics/health-emergencies/coronavirus-covid-19/novel-coronavirus-2019-ncov). 2. [2]. Tristan Beau, Julien Browaeys, Olivier Dadoun. CoCoA, a free Collaborative Covid Analysis tool. Preprint under preparation. 3. [3]. Tristan Beau, Olivier Dadoun, Julien Browaeys. GitHub CoCoA Official Repository (version 0.1). Python, 2020. [https://github.com/tjbtjbtjb/CoCoA](https://github.com/tjbtjbtjb/CoCoA). 4. [4].“Population by Country (2020) - Worldometer.” Accessed May 4, 2020. [https://www.worldometers.info/world-population/population-by-country/](https://www.worldometers.info/world-population/population-by-country/). 5. [5].CSSEGISandData/COVID-19 by the Center for Systems Science and Engineering (CSSE) at Johns Hopkins University, 2020. [https://github.com/CSSEGISandData/COVID-19](https://github.com/CSSEGISandData/COVID-19). 6. [6].Skwarnicki, T. A study of the radiative CASCADE transitions between the Upsilon-Prime and Upsilon resonances. (PhD Thesis, Cracow, INP, 1986). [https://inspirehep.net/literature/230779](https://inspirehep.net/literature/230779) [1]: /embed/graphic-1.gif