A Contribution to the Mathematical Modeling of the Corona/COVID-19 Pandemic =========================================================================== * Günter Bärwolff ## Abstract The responsible estimation of parameters is a main issue of mathematical pandemic models. Especially a good choice of *β* as the number of others that one infected person encounters per unit time (per day) influences the adequateness of the results of the model. For the example of the actual COVID-19 pandemic some aspects of the parameter choice will be discussed. Because of the incompatibility of the data of the Johns-Hopkins-University [3] to the data of the German Robert-Koch-Institut we use the COVID-19 data of the European Centre for Disease Prevention and Control [2] (ECDC) as a base for the parameter estimation. Two different mathematical methods for the data analysis will be discussed in this paper and possible sources of trouble will be shown. Parameters for several countries like UK, USA, Italy, Spain, Germany and China will be estimated and used in W. O. Kermack and A. G. McKendrick’s SIR model[1]. Strategies for the commencing and ending of social and economic shutdown measures are discussed. The numerical solution of the ordinary differential equation system of the modified SIR model is being done with a Runge-Kutta integration method of fourth order [4]. At the end the applicability of the SIR model could be shown. Suggestions about appropriate points in time at which to commence with lockdown measures based on the acceleration rate of infections conclude the paper. This paper is an improved sequel of [5]. ## 1 The mathematical SIR model At first we will describe the model. *I* denotes the infected people, *S* stands for the susceptible and *R* denotes the recovered people. The dynamics of infections and recoveries can be approximated by the ODE system ![Formula][1] ![Formula][2] ![Formula][3] We understand *β* as the number of others that one infected person encounters per unit time (per day). *γ* is the reciprocal value of the typical time from infection to recovery. *N* is the total number of people involved in the epidemic disease and there is *N* = *S* + *I* + *R*. The empirical data currently available suggests that the corona infection typically lasts for some 14 days. This means *γ* = 1*/*14 ≈ 0,07. The choice of *β* is more complicated and will be considered in the next section. ## 2 The estimation of *β* based on real data We use the European Centre for Disease Prevention and Control [2] as a data for the COVID-19 infected people for the period from December 31st 2019 to April 8th 2020. At the beginning of the pandemic the quotient *S/N* is nearly equal to 1. Also, at the early stage no-one has yet recovered. Thus we can describe the early regime by the equation ![Formula][4] with the solution ![Formula][5] We are looking for periods in the spreadsheets of infected people per day where the course can be described by a function of type (4). Starting with a spreadsheet like View this table: [Table1](http://medrxiv.org/content/early/2020/04/19/2020.04.01.20050229/T1) for a certain country and a chosen period [*t*1, *t**k*] with my favored method We search for the minimum of the functional ![Formula][6] i.e. ![Formula][7] We solved this non-linear minimum problem with the damped Gauss-Newton method (see [4]). After some numerical tests we found the subsequent results for the considered countries. Thereby we chose different periods for the countries with the aim to approximate the infection course in a good quality. The following figures show the graphs and the evaluated parameter. It must be said that evaluated *β*-values are related to the stated period. For the iterative Gauss-Newton method we guessed the respective periods for every country by a visual inspection of the graphs of the infected people over days. Especially in medicine, psychology and other life sciences the logarithm behavior of data was readily considered. Instead of the above table of values the following logarithmic one was used. View this table: [Table2](http://medrxiv.org/content/early/2020/04/19/2020.04.01.20050229/T2) The logarithm of (4) leads to ![Formula][8] and based on the logarithmic table the functional ![Formula][9] is to minimize. The solution of this linear optimization problem is trivial and it is available in most of computer algebra systems as a “block box” of the logarithmic-linear regression. The following figures show the results for the same periods as above for Spain, the UK, the USA and Italy. Figures 7-14 show that the logarithmic-linear regression implies poor results. Thus, the non-linear optimization problem (5) is to choose as the favored method for the estimation of *I* and *β*. We found some notes on the parameters of Italy in the literature, for example *β* = 0 25, and we are afraid that this is a result of the logarithmic-linear regression. ![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/19/2020.04.01.20050229/F1.medium.gif) [Figure 1:](http://medrxiv.org/content/early/2020/04/19/2020.04.01.20050229/F1) Figure 1: German course from January 31st 2020 to March 20th 2020 ![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/19/2020.04.01.20050229/F2.medium.gif) [Figure 2:](http://medrxiv.org/content/early/2020/04/19/2020.04.01.20050229/F2) Figure 2: Italian course from January 31st 2020 to March 20th 2020 ![Figure 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/19/2020.04.01.20050229/F3.medium.gif) [Figure 3:](http://medrxiv.org/content/early/2020/04/19/2020.04.01.20050229/F3) Figure 3: Spanish course from January 31st 2020 to March 20th 2020 ![Figure 4:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/19/2020.04.01.20050229/F4.medium.gif) [Figure 4:](http://medrxiv.org/content/early/2020/04/19/2020.04.01.20050229/F4) Figure 4: UK course from February 20th 2020 to March 20th 2020 ![Figure 5:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/19/2020.04.01.20050229/F5.medium.gif) [Figure 5:](http://medrxiv.org/content/early/2020/04/19/2020.04.01.20050229/F5) Figure 5: USA course from February 10th 2020 to April 4th 2020 ![Figure 6:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/19/2020.04.01.20050229/F6.medium.gif) [Figure 6:](http://medrxiv.org/content/early/2020/04/19/2020.04.01.20050229/F6) Figure 6: Chines course from December 31st 2019 to January 28th 2020 ![Figure 7:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/19/2020.04.01.20050229/F7.medium.gif) [Figure 7:](http://medrxiv.org/content/early/2020/04/19/2020.04.01.20050229/F7) Figure 7: log-lin-result for Spain (January 31st 2020 to March 20th 2020) ![Figure 8:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/19/2020.04.01.20050229/F8.medium.gif) [Figure 8:](http://medrxiv.org/content/early/2020/04/19/2020.04.01.20050229/F8) Figure 8: Logarithm of the Spanish result (January 31st 2020 to March 20th 2020) ![Figure 9:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/19/2020.04.01.20050229/F9.medium.gif) [Figure 9:](http://medrxiv.org/content/early/2020/04/19/2020.04.01.20050229/F9) Figure 9: log-lin-result of the UK (February 20th 2020 to March 20th 2020) ![Figure 10:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/19/2020.04.01.20050229/F10.medium.gif) [Figure 10:](http://medrxiv.org/content/early/2020/04/19/2020.04.01.20050229/F10) Figure 10: Logarithm of the UK result (February 20th 2020 to March 20th 2020) ![Figure 11:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/19/2020.04.01.20050229/F11.medium.gif) [Figure 11:](http://medrxiv.org/content/early/2020/04/19/2020.04.01.20050229/F11) Figure 11: log-lin-result of the USA (February 10th 2020 to April 4th 2020) ![Figure 12:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/19/2020.04.01.20050229/F12.medium.gif) [Figure 12:](http://medrxiv.org/content/early/2020/04/19/2020.04.01.20050229/F12) Figure 12: Logarithm of the USA result (February 10th 2020 to April 4th 2020) ![Figure 13:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/19/2020.04.01.20050229/F13.medium.gif) [Figure 13:](http://medrxiv.org/content/early/2020/04/19/2020.04.01.20050229/F13) Figure 13: log-lin-result of Italy (January 31st 2020 to March 20th 2020) ![Figure 14:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/19/2020.04.01.20050229/F14.medium.gif) [Figure 14:](http://medrxiv.org/content/early/2020/04/19/2020.04.01.20050229/F14) Figure 14: Logarithm of the Italian result (January 31st 2020 to March 20th 2020) ## 3 Some numerical computations for Germany and Spain With the choice of *β*-value 0,215 (see fig. 1) which was evaluated on the basis of the real data of ECDC and *γ* = 0,07 one gets the course of the pandemic dynamics pictured in fig. 15.1. *R* is the basis reproduction number of persons, infected by the transmission of a pathogen from one infected person during the infectious time (*R* = *β/γ*) in the following figures. ![Figure 15:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/19/2020.04.01.20050229/F15.medium.gif) [Figure 15:](http://medrxiv.org/content/early/2020/04/19/2020.04.01.20050229/F15) Figure 15: German course of one year, starting end of January 2020, *S*-green, *I*-red, *R*-blue Neither data from ECDC nor the data from the German Robert-Koch-Institut and the data from the Johns Hopkins University are correct, for we have to reasonably assume that there are a number of unknown cases. It is guessed that the data covers only 15% of the real cases. Considering this we get a slightly changed results and in the subsequent computations we will include estimated number of unknown cases to the initial values of *I*. For Spain we use the *β*-value 0,249 (see fig. 2) and *γ* = 0,07 we get the course pictured in fig. 16. *N* was set to 40 millions. ![Figure 16:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/19/2020.04.01.20050229/F16.medium.gif) [Figure 16:](http://medrxiv.org/content/early/2020/04/19/2020.04.01.20050229/F16) Figure 16: Spanish course of one year, starting end of January 2020, *S*-green, *I*-red, *R*-blue ## 4 Influence of a temporary lockdown and extensive social distancing In all countries concerned by the Corona pandemic a lockdown of the social life is discussed. In Germany the lockdown started at March 16th 2020. The effects of social distancing to decrease the infection rate can be modeled by a modification of the SIR model. The original ODE system (1)-(3) was modified to ![Formula][10] ![Formula][11] ![Formula][12] *κ* is a function with values in [0,1]. For example ![Formula][13] means for example a reduction of the infection rate of 50% in the period [*t*, *t*1] (Δ*t* = *t*1 *t* is the duration of the temporary lockdown in days). A good choice of *t* and *t**k* is going to be complicated. If we respect the chosen starting day of the German lockdown, March 16th 2020 (this conforms the 46th day of the concerned year), and we work with ![Formula][14] we got the result pictured in fig. 17. ![Figure 17:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/19/2020.04.01.20050229/F17.medium.gif) [Figure 17:](http://medrxiv.org/content/early/2020/04/19/2020.04.01.20050229/F17) Figure 17: German course of one year, starting end of January 2020, *S*-green, *I*-red, *R*-blue, 30 days lockdown, starting time March 16th 2020 The numerical tests showed that a very early start of the lockdown resulting in a reduction of the infection rate *β* results in the typical Gaussian curve to be delayed by *I*; however, the amplitude (maximum value of *I*) doesn’t really change. One knows that the development of the infected people looks like a Gaussian curve. The interesting points in time are those where the acceleration of the numbers of infected people increases or decreases, respectively. These are the points in time where the curve of *I* was changing from a convex to a concave behavior or vice versa. The convexity or concavity can be controlled by the second derivative of *I*(*t*). Let us consider equation (2). By differentiation of (2) and the use of (1) we get ![Formula][15] With that the *I*-curve will change from convex to concave if the relation ![Formula][16] is valid. For the switching time follows ![Formula][17] A lockdown starting at *t* (assigning *β**∗* = *κβ, κ*∈ [0,1[) up to a point in time *t*1 = *t* + Δ*t*, with Δ*t* as the duration of the lockdown in days, will be denoted as a **dynamical lockdown** (for *t* > *t*1 *β**∗* was reset to the original value *β*). *t* means the point in time up to which the growth rate increases and from which on it decreases. Fig. 18 shows the result of such a computation of a dynamical lockdown. We got *t* = 108 (*κ* = 0,2)-The result is significant. In fig. 20 a typical behavior of ![Graphic][18] is plotted. ![Figure 18:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/19/2020.04.01.20050229/F18.medium.gif) [Figure 18:](http://medrxiv.org/content/early/2020/04/19/2020.04.01.20050229/F18) Figure 18: German course of one year, starting end of January 2020, dynamical lockdown, *S*-green, *I*-red, *R*-blue The result of a dynamical 30 days lockdown for Spain is shown in fig. 19, where we found *t* = 106 (*κ* = 0,2). ![Figure 19:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/19/2020.04.01.20050229/F19.medium.gif) [Figure 19:](http://medrxiv.org/content/early/2020/04/19/2020.04.01.20050229/F19) Figure 19: Spanish course of one year, starting end of March 2020, dynamical lockdown, *S*-green, *I*-red, *R*-blue Data from China and South Korea suggests that the group of infected people with an age of 70 or more is of magnitude 10%. This group has a significant higher mortality rate than the rest of the infected people. Thus we can presume that *α*=10% of *I* must be especially sheltered and possibly medicated very intensively as a highrisk group. Fig. 21 shows the German time history of the above defined high-risk group with a dynamical lockdown with *κ* = 0,2 compared to regime without social distancing. The maximum number of infected people decreases from approximately 1,7 millions of people to 0,8 millions in the case of the lockdown (30 days lockdown). ![Figure 20:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/19/2020.04.01.20050229/F20.medium.gif) [Figure 20:](http://medrxiv.org/content/early/2020/04/19/2020.04.01.20050229/F20) Figure 20: History of the second derivative of *I* (de) ![Figure 21:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/19/2020.04.01.20050229/F21.medium.gif) [Figure 21:](http://medrxiv.org/content/early/2020/04/19/2020.04.01.20050229/F21) Figure 21: German history of the infected people of high-risk groups depending on a dynamical lockdown This result proves the usefulness of a lockdown or a strict social distancing during an epidemic disease. We observe a flattening of the infection curve as requested by politicians and health professionals. With a strict social distancing for a limited time one can save time to find vaccines and time to improve the possibilities to help high-risk people in hospitals. To see the influence of a social distancing we look at the spanish situation without a lockdown and a dynamical lockdown of 30 days with fig. 22 (*κ* = 0,2) for the 10% high-risk people. ![Figure 22:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/19/2020.04.01.20050229/F22.medium.gif) [Figure 22:](http://medrxiv.org/content/early/2020/04/19/2020.04.01.20050229/F22) Figure 22: Spanish history of the infected people of high-risk groups depending on a dynamical lockdown The computations with the SIR model show, that the social distancing with a lockdown will only be successful with a start behind the time greater or equal to *t*, found by the evaluation of the second derivative of *I* (formula (10)). If the lockdown is started at a time less then *t* the effect of such a social distancing is not significant. ## 5 Closing remarks If we write (2) or (7) resp. in the form ![Formula][19] we realize that the number of infected people decreases if ![Formula][20] is complied. The relation (11) shows that there are two possibilities for the rise of infected people to be inverted and the medical burden to be reduced. 1. The reduction of the stock of the species *S*. This can be obtained by immunization or vaccination. Another possibility is the isolation of high-risk people (70 years and older). Positive tests for antibodies reduce the stock of susceptible persons. 2. A second possibility is the reduction of the infection rate *κβ*. This can be achieved by strict lockdowns, social distancing at appropriate times, or rigid sanitarian moves. The results are pessimistic in total with respect to a successful fight against the COVID-19-virus. Hopefully the reality is a bit more merciful than the mathematical model. But we rather err on the pessimistic side and be surprised by more benign developments. Note again that the parameters *β* and *κ* are guessed very roughly. Also, the percentage *α* of the group of high-risk people is possibly overestimated. Depending on the capabilities and performance of the health system of the respective countries, those parameters may look different. The interpretation of *κ* as a random variable is thinkable, too. ## Data Availability I use actual data if coronavirus infected people given by the Johns-Hopkins-University [https://www.jhu.edu](https://www.jhu.edu) ## Footnotes * 1 *I*0 denotes the initial value of the *I* species, that is January 31th 2020. *Imax* stands for the maximum of *I*. The total number *N* for Germany is guessed to be 75 millions. * Received April 1, 2020. * Revision received April 15, 2020. * Accepted April 19, 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]. W.O. Kermack and A.G. McKendrick, A contribution to the mathematical theory of epidemics. Proc. R. Soc. London A 115(1927)700. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1098/rspa.1927.0118&link_type=DOI) 2. [2].Bulletins of the European Centre for Disease Prevention and Control ([https://www.ecdc.europa.eu/en/geographical-distribution-2019-ncov-cases](https://www.ecdc.europa.eu/en/geographical-distribution-2019-ncov-cases)) 2020. 3. [3].Bulletins of the John Hopkins University of world-wide Corona data ([https://www.jhu.edu](https://www.jhu.edu)) 2020. 4. [4]. G. Bärwolff, Numerics for engineers, physicists and computer scientists (3rd ed., in German). Springer-Spektrum 2020. 5. [5]. G. Bärwolff, A Contribution to the Mathematical Modeling of the Corona/COVID-19 Pandemic. medRxiv.preprint 2020, doi: [https://doi.org/10.1101/2020.04.01.20050229](https://doi.org/10.1101/2020.04.01.20050229). 6. [6]. Toshihisa Tomie, Understandig the present status and forcasting of COVID-19 in Wuhan. medRxiv.preprint 2020. [1]: /embed/graphic-1.gif [2]: /embed/graphic-2.gif [3]: /embed/graphic-3.gif [4]: /embed/graphic-4.gif [5]: /embed/graphic-5.gif [6]: /embed/graphic-7.gif [7]: /embed/graphic-8.gif [8]: /embed/graphic-10.gif [9]: /embed/graphic-11.gif [10]: /embed/graphic-28.gif [11]: /embed/graphic-29.gif [12]: /embed/graphic-30.gif [13]: /embed/graphic-31.gif [14]: /embed/graphic-32.gif [15]: /embed/graphic-34.gif [16]: /embed/graphic-35.gif [17]: /embed/graphic-36.gif [18]: /embed/inline-graphic-1.gif [19]: /embed/graphic-42.gif [20]: /embed/graphic-43.gif