Forecasting COVID-19 cases and deaths in epidemic-mitigating European countries by Richards function-based regression analyses ============================================================================================================================== * Cheng Long * Xinmiao Fu ## Abstract The COVID-19 pandemic has hit many countries, and in some European countries it has been mitigated since April. Here we applied Richards function to simulate and forecast the course of COVID-19 epidemics in Italy, Spain, France, Germany, Turkey, Belgium, Ireland, Netherlands, Portugal and Switzerland. Potential total COVID-19 confirmed cases in these countries were estimated to be 240400±1300, 294100±4000, 178500±800, 176900±700, 155400±1000, 57900±400, 24000±200, 46200±300, 30000±300 and 30700±100 respectively. Most of these countries are predicted to approach ending stage between late May and early June such that daily new cases will become minimal, which may guide societal and economic restorations. In addition, total COVID-19 deaths were estimated to be 33500 ±300, 28200±200, 27800±200, 8740±80, 4500 ±30, 9250±70, 1530±20, 6240±50, 1380±10 and 1960±8, respectively. To our best knowledge, this is the first study forecasting the COVID-19 epidemic by applying the Richard function-based regression analysis. Keywords * COVID-19 * SARS-CoV-2 * Richards function * coronavirus * pandemic * total confirmed cases * total COVID-19 deaths ## Introduction The COVID-19 pandemic caused by the novel coronavirus (SARS-CoV-2) has hit many countries [1], resulting in more than 4 million cases and over 300,000 deaths in the world as of May 14, 2020. In European countries such as Italy and Spain, daily new cases are declining since April, yet they are still as many as one thousand per day. In these epidemic-mitigating countries, governments are attempting to restore normal societal and economic activity, however, concerns are raised about the risk of COVID-19 resurgence. It is of significance to precisely forecast when the epidemic in each country will become minimal such that societal and economic restorations will not lead in the resurgence of the outbreak. Further, such forecasting may also inform governments to get better preparedness for making control measures and supplying medical resources. Estimations by traditional epidemiological models apparently require much detailed data for analysis [2, 3]. Here we applied a simple data-driven, Richards function-based approach for epidemic forecasting only based on the daily cumulative confirmed COVID-19 cases. This approach has been explored for simulating 2003 SARS outbreaks in several regions/countries by others [4, 5] and also for estimating COVID-19 deaths in China by us recently [6]. ## Methods The Richards function is expressed as follows [4, 5]: ![Formula][1] where *y* is the daily cumulative number of confirmed COVID-19 cases and *x* is day; *a*, *d, xc* and *k* are constants. Specifically, total number of COVID-19 cases are given by parameter *a* in the equation. Simulation was performed using Microcal Origin software with the Richard function. ## Results ### Estimation of potential total COVID-19 confirmed cases We decided to collect the data of the cumulative confirmed cases from the website of Worldometer [7] for 10 most affected and epidemic-mitigating European countries, i.e., Italy, Spain, France, Germany, Turkey, Belgium, Ireland, Netherlands, Portugal and Switzerland (for detail, refer to the supporting information file **original data.xlsx**). Regression analyses indicate that the data of each country were well fitted with the Richards function (all *R2* values being close to 0.999; **Table 1** and **Figure 1**). Potential total numbers of confirmed cases in the above countries were estimated to be 240400±1300, 294100±4000, 178500±800, 176900±700, 155400±1000, 57900±400, 24000±200, 46200±300, 30000±300 and 30700±100 respectively. All parameters from regression analyses are presented in **Table S1** and can be utilized to predict daily new cases in each country since May 15 according to equation 1. View this table: [Table 1](http://medrxiv.org/content/early/2020/05/26/2020.05.18.20106146/T1) Table 1 Estimation of potential total COVID-19 confirmed cases and deaths in European countries. ![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/26/2020.05.18.20106146/F1.medium.gif) [Figure 1.](http://medrxiv.org/content/early/2020/05/26/2020.05.18.20106146/F1) Figure 1. Fitting the cumulative confirmed COVID-19 cases of each European country to the Richards function (**A, B**) Plots of the cumulative number of COVID-19 confirmed cases of the indicated European countries as of May 14, 2020, with the simulation results being plotted as Magenta lines. Simulation was performed using Microcal Origin software with the Richard function. All parameters from regression analyses are shown in Table S1. In addition, we estimated the key date, on which daily new cases are lower than 0.1% of the potential total confirmed cases as defined by us subjectively. It appears that most countries will approach such late stage between late May and early June (refer to **Table 1**) such that daily new cases will become minimal in comparison with potential total confirmed cases. In our opinion, nationwide societal and economic restorations should be initiated at least not earlier than the estimated key dates, although local restorations may be operated earlier according to the local epidemic situation. ### Estimation of potential total COVID-19 deaths In several countries, crude fatality ratios are as high as over 10% such that people have raised grave concerns about how many patients will die eventually. Assuming that the number of deaths is proportional to the number of confirmed cases for the COVID-19 outbreak under specific circumstances, we speculated that the cumulative COVID-19 deaths would also obey the Richards function. Regression analyses indicate that the data of cumulative COVID-19 deaths in each country were well fitted with the Richards function (all *R* values being close to 0.999; **Table 1** and **Figure 2**). All parameters of regression analyses are displayed in **Table S2** and can be utilized to predict daily new deaths in each country since May 15 according to equation 1. Potential total COVID-19 deaths in Italy, Spain, France, Germany, Turkey, Belgium, Ireland, Netherlands, Portugal and Switzerland were estimated to be 33500 ±300, 28200±200, 27800±200, 8740±80, 4500 ±30, 9250±70, 1530±20, 6240±50, 1380±10 and 1960±8, respectively. Crude fatality ratios (total deaths among total confirmed cases) of these countries could be estimated to be 13.9±0.9%, 9.6±0.9%, 15.6±0.7%, 4.9±0.9%. 2.9±0.7%, 16.0±0.9%, 6.4±1.3%, 13.5±0.9%, 4.6±0.8% and 6.4±0.4%, respectively. ![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/26/2020.05.18.20106146/F2.medium.gif) [Figure 2.](http://medrxiv.org/content/early/2020/05/26/2020.05.18.20106146/F2) Figure 2. Fitting the cumulative COVID-19 deaths of each European country to the Richards function (**A, B**) Plots of the cumulative number of COVID-19 confirmed cases of the indicated European countries as of May 14, 2020, with the simulation results being plotted as Magenta lines. Simulation was performed using Microcal Origin software with the Richard function. All parameters from regression analyses are shown in Table S2. ## Discussion To our best knowledge, this is the first study forecasting the COVID-19 epidemic by applying the Richard function-based regression analysis. It should be pointed out that all the above estimates are based on the assumption that there are no significant changes regarding the current interventions, including non-pharmacological interventions such as quarantine and other prevention measures and also diagnosis and treatment approaches. In addition, UK and Sweden, two most affected European countries with more than 240,000 and 29,000 COVID-19 confirmed cases as of May 15, respectively, were not subjected to analysis because their daily new cases have not yet started to significantly decline (**Figure. S1**). Collectively, our forecasting on the course of COVID-19 outbreaks in these European countries may provide valuable guidance for governments and also the public to get better preparedness and optimize their efforts to contain this unprecedented crisis at current critical stage. The estimated key date for each country may serve as a reference for operating societal and economic restorations such that the resurgence of COVID-19 can be avoided, minimized or under control according to COVID-19 containment in China [8, 9]. In addition, because thousands of COVID-19 patients have died in these countries, our prediction on the course of COVID-19 deaths (refer to **Figure 2** and **Table S2**) may benefit the mental health service that needs to be timely provided to the families of passed patients [10]. As a matter of fact, we precisely estimated total COVID-19 deaths in different areas of China (including Hubei Province and Wuhan city, the epicenter of the COVID-19 outbreaks in China [1]) by applying both Boltzmann and Richards functions-based regression analyses [6]. Such estimation help governments to evaluate the severity of the COVID-19 epidemics. ## Data Availability Data are included in the manuscript ## Author contribution X.F. designed the study; C.L. organized data and performed analysis; X.F. wrote the manuscript. ## Conflict of interest none. ## Acknowledgments This work is support by the National Natural Science Foundation of China (No. 31972918 and 31770830 to XF). We thank graduate students (Zuqin Zhang, Youhui Zeng, Shuting Shi and Zhijie Huang, all from Professor Fu’s research group at Fujian Normal University) for data collection. * Received May 18, 2020. * Revision received May 18, 2020. * Accepted May 26, 2020. * © 2020, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/) ## References 1. (1).WHO: Coronavirus disease (COVID-19) outbreak, [https://www.who.int/emergencies/diseases/novel-coronavirus-2019](https://www.who.int/emergencies/diseases/novel-coronavirus-2019). 2. (2).Wu JT, Leung K, Leung GM. Nowcasting and forecasting the potential domestic and international spread of the 2019-nCoV outbreak originating in Wuhan, China: a modelling study. Lancet 2020: DOI: 10.1016/S0140-6736(1020)30260-30269. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(1020)30260-30269&link_type=DOI) 3. (3).Yang Y, et al. Epidemiological and clinical features of the 2019 novel coronavirus outbreak in China. [https://doiorg/101101/2020021020021675](https://doiorg/101101/2020021020021675) 2020. 4. (4).Wang XS, Wu J, Yang Y. Richards model revisited: Validation by and application to infection dynamics. J Theor Biol 2012; 313: 12-19. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jtbi.2012.07.024&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22889641&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F26%2F2020.05.18.20106146.atom) 5. (5).Hsieh YH. Richards Model: A Simple Procedure for Real-time Prediction of Outbreak Severity. Modeling and Dynamics of Infectious Diseases 2009: p. 216–236 [https://doi.org/210.1142/9789814261265_9789814260009](https://doi.org/210.1142/9789814261265_9789814260009). 6. (6).Gao Y, et al. Forecasting the cumulative number of COVID-19 deaths in China: a Boltzmann function-based modeling study. Infect Control Hosp Epidemiol 2020: 1-3. 7. (7).Worldometer: COVID-19 CORONAVIRUS PANDEMIC, [https://www.worldometers.info/coronavirus/](https://www.worldometers.info/coronavirus/). 8. (8).Report of the WHO-China Joint Mission on Coronavirus Disease 2019 (COVID-19), [https://www.who.int/docs/default-source/coronaviruse/who-china-ioint-mission-on-covid-19-final-report,pdf](https://www.who.int/docs/default-source/coronaviruse/who-china-ioint-mission-on-covid-19-final-report.pdf): 2020. 9. (9).National Health Commission of China: Update of the COVID-19 Outbreak Data (in Chinese), [http://www.nhc.gov.cn/xcs/yqfkdt/gzbdindex.shtml](http://www.nhc.gov.cn/xcs/yqfkdt/gzbdindex.shtml). 10. (10).Xiang YT, et al. Timely mental health care for the 2019 novel coronavirus outbreak is urgently needed. Lancet Psychiatry 2020; 7(3): 228–229.158 [1]: /embed/graphic-1.gif