Modeling the COVID-19 pandemic - parameter identification and reliability of predictions ======================================================================================== * Klaus Hackl ## Abstract In this paper, we try to identify the parameters in an elementary epidemic model, the so-called SI-model, via non-linear regression using data of the COVID-19 pandemic. This is done based on the data for the number of total infections and daily infections, respectively. Studying the convergence behavior of the two parameter sets obtained this way, we attempt to estimate the reliability of predictions made concerning the future course of the epidemic. We validate this procedure using data for the case numbers in China and South Korea. Then we apply it in order to find predictions for Germany, Italy and the United States. The results are encouraging, but no final judgment on the validity of the procedure can yet be made. Keywords * COVID-19 * modeling * parameter identification * predictions ## 1 Introduction In most countries, social distancing measures are in effect now in order to fight the covid-19 pandemic. Considering the serious effects of these measures on the affected societies and the ensuing political discussions on their intensity and duration, it would be highly desirable to be able to make modeling based predictions on the future timeline of the epidemic, as long as the measures are upheld. Of course, many attempts are made in this direction. However, most of them require very detailed data that are laborious and time-consuming to generate. In this work, we try to study the possibility to base predictions on data sets readily available, namely the number of reported infections. We are aware, that these numbers depend strongly on the intensity of testing done in the various countries and the reliability of the reported numbers. In this work we presume that there is a factor, country-specific, but constant in time, between the reported and the actual number of cases. If this assumption were valid, the total number of infected individuals would be off by this very factor. However, other parameters, like the point in time when the peak in the numbers of daily infections would occur, or the following rate of decay of these numbers, would not be affected. Finally, we would like to stress, that we intend this work to be the starting point of a discussion and maybe further research. By no means, having a background in engineering and not in virology or epidemiology, we are claiming any medical expertise. The paper should be rather seen as a general exercise in modeling and interpretation of data. ## 2 An elementary model Our aim is to model a situation where social distancing measures are in effect, as currently is the case in most countries. This means, that only a small portion of the population is affected, which is well but not completely isolated from the rest. As starting point, we refer to the compartmental model by Kermack, McKendrick and Walker, [3]. It is defined by the differential equations ![Formula][1] Here *I*(*t*) is the number of individuals in the infectious population and *S*(*t*) denotes the number of individuals in the susceptible population, in our case those who can get infected because they are not protected by social distancing. This formulation is also called the SIR-model, where the dependent variable *R*(*t*) stands for the removed (by recovery of death) population, and we have ![Graphic][2]. The parameter *α* is related to the basic reproduction number by ![Formula][3] where *N* is the initially susceptible population and *T*inf = 1/*γ* is the time period during which an individual is infectious. For Sars-Cov-2, no definite value for *T*inf has yet been reported. The parameter *γ* denotes the rate at which individuals are removed from the infected population because of an outcome (recovery or death). In our study, we are going to neglect the term *γI*, which will have an effect only in late stages of the epidemic when *S* ≈ *γ*/*α*. As a result, we obtain the so-called SI-model, [4]. Later, we will see, that it is basically impossible to identify the parameter *γ* from the data available unless the epidemic is in a very late stage. So, for our purposes, this simplification amounts to a necessity. Employing the assumption above, Eqs. (1) become equivatent to the so-called logistic differential equation having the closed form solution ![Formula][4] where ![Formula][5] and *t*max, marking the peak of the epidemic, is defined by ![Formula][6] ## 3 Parameter identification In order to achieve a more robust parameter identification, we precondition our solution by introducing new parameters *a, b* given by ![Formula][7] Note, that Eq. (6) implies ![Formula][8] After substitution of Eqs. (6) into Eqs. (3), the number of total infections is then given as ![Formula][9] and the rate of daily infections becomes ![Formula][10] We determine the three parameters {*a, b, t*max} of our model via non-linear regression. The data taken from the worldometer web page, [1], which essentially uses the data from the Johns Hopkins University Center for Systems Science and Engineering (JHU CCSE). For the parameter identification done in this paper, we have used the available data up to including Apr. 4, 2020. The data are provided in form of lists ![Graphic][11] for the total number of infections up to day *t**i*, and ![Graphic][12] for the number of daily infections. Time is measured in days, starting on Jan. 1, 2020. Hence, *t* = 1 d corresponds to Jan. 1, *t* = 32 d to Feb. 1, *t* = 61 d to Mar. 1, 2020, and so on. Obviously, we have ![Formula][13] Let us define errors *e*(*a, b, t*max) with respect to the total cases and *e*1(*a, b, t*max) with respect to the daily cases by ![Formula][14] In order to judge the accuracy of the modeling, let us define the data norms ![Formula][15] and the relative errors ![Formula][16] Finally, we find the parameters *a, b, t*max by minimizing the errors: ![Formula][17] ![Formula][18] Minimization is done using the computer algebra system *Mathematica*, [2]. For our purposes, the *simulated annealing* global minimization algorithm works best. Attention has to be given, though, to choosing appropriate initial intervals for the parameters in order to achieve convergence. ## 4 Results In Figs. 1 to 5, the numbers of daily cases (left) and total cases (right) are plotted versus time in days. In order to get an estimation of the variability of the predictions, we use both parameter identification schemes defined in Eqs, (14) and (15). The results obtained by fitting the number of total cases according to Eq. (14) are shown in red color and those obtained obtained by fitting the number of daily cases according to Eq. (15) are shown in magenta. The corresponding data are shown in blue color. ![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/11/2020.04.07.20056937/F1.medium.gif) [Figure 1:](http://medrxiv.org/content/early/2020/04/11/2020.04.07.20056937/F1) Figure 1: China, left: daily cases (![Graphic][19] red, ![Graphic][20] magenta), right: total cases (![Graphic][21] in red, ![Graphic][22] in magenta), data in blue ![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/11/2020.04.07.20056937/F2.medium.gif) [Figure 2:](http://medrxiv.org/content/early/2020/04/11/2020.04.07.20056937/F2) Figure 2: South Korea, left: daily cases (![Graphic][23] red, ![Graphic][24] magenta), right: total cases (![Graphic][25] in red, ![Graphic][26] in magenta), data in blue ![Figure 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/11/2020.04.07.20056937/F3.medium.gif) [Figure 3:](http://medrxiv.org/content/early/2020/04/11/2020.04.07.20056937/F3) Figure 3: Germany, left: left: daily cases (![Graphic][27] red, ![Graphic][28] magenta), right: total cases (![Graphic][29] in red, ![Graphic][30] in magenta), data in blue ![Figure 4:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/11/2020.04.07.20056937/F4.medium.gif) [Figure 4:](http://medrxiv.org/content/early/2020/04/11/2020.04.07.20056937/F4) Figure 4: Italy, left: left: daily cases (![Graphic][31] red, ![Graphic][32] magenta), right: total cases (![Graphic][33] in red, ![Graphic][34] in magenta), data in blue ![Figure 5:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/11/2020.04.07.20056937/F5.medium.gif) [Figure 5:](http://medrxiv.org/content/early/2020/04/11/2020.04.07.20056937/F5) Figure 5: United States, left: left: daily cases (![Graphic][35] red, ![Graphic][36] magenta), right: total cases (![Graphic][37] in red, ![Graphic][38] in magenta), data in blue In Fig. 1 and Fig. 2 the data for China and South Korea are displayed. Both countries can be considered to be in a late stage of the epidemic and the data are matched well by the model. Especially for China, the predictions obtained by both parameter identification schemes are close together. The pronounced spike in the number of daily cases is due to a change of the procedure how infections are counted, and is averaged out by the model. It is apparent that the model cannot fit the remaining almost constant level of daily infections around 100 and the corresponding ongoing rise in the total cases in the South Korea data. This causes the predictions generated by both procedures to lie a little further apart. In Figs. 3, 4 and 5, the numbers of daily new cases are plotted for Germany, Italy and the United States. These countries can be considered to be in earlier stages of the epidemic. For Germany and Italy, the predictions generated by both procedures agree closely in the time range where data are already available and are divergent for later points in time. This divergence is especially strong in case of the United States data, indicating a very dynamic epidemic process of exponential growth taking place there. Some key data provided by the model are given in Table 1. They are defined as follows: View this table: [Table 1:](http://medrxiv.org/content/early/2020/04/11/2020.04.07.20056937/T1) Table 1: values of key parameters on Apr. 4, 2020 ![Graphic][39]: This is the total number of population getting infected, if the social distancing measures taken are upheld until the epidemic has completely subsided. Note, that for China and South Korea,these numbers correspond closely to the total number of cases reported. ![Graphic][40]: This is the time predicted when the number of daily new cases will drop below 100 using both regression procedures, indicating a point in time when social distancing measures might be loosened. For Germany, Italy and the United States the values are in the range from mid of April to mid of May. *T*data: The dates, when the number of daily new cases dropped below 100 in China and South Korea, agreeing very well with the model data. (For South Korea, the number of daily new cases is fluctuating strongly. So we took the date, when the number dropped below 100 the second time.) *e*0,rel, *e*1,rel: The is relative errors produced by both parameter identification procedures, as defined in Eq. (13). Note, that due to the fluctuating nature of the number of daily cases, *e*1,rel is much larger than *e*0,rel. ## 5 Reliability of predictions In the initial stages of an epidemic, the number of cases grows exponentially. This is easy to see by considering the limit *t* → −∞ in Eqs. (8) and (9), giving ![Formula][41] From Eq. (16), we see that during an early stage, it is impossible to identify the parameters *a* and *t*max independently, because the occur via the common factor ![Graphic][42]. Hence, there are infinitely many pairs of *a* and *t*max giving the same behavior and thus being minimizers in Eqs. (14) and (15). By virtue, it is even harder to fit all parameters in the SIR-model or one of the many existing extensions of it. Only past this phase of exponential growth, it is possible to identify all three parameters and thus arrive at viable predictions. But how to identify this point in time from the available data? Our suggestion for a solution of this problem is to monitor the two different parameter sets ![Graphic][43] and ![Graphic][44] defined in Eqs. (14) and (15). Theoretically, the values should be close to each other. However, during a phase of exponential growth, the mentioned ill-posedness of the minimization problems given by Eqs. (14) and (15) will give results lying substantially apart. Lets test this hypothesis: In Figs. 6 to 10, we display the parameters ![Graphic][45] (in red) and ![Graphic][46] (in magenta) to the left and *b* (in red) and *b*1 (in magenta) to the right, respectively, versus time in days. ![Figure 6:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/11/2020.04.07.20056937/F6.medium.gif) [Figure 6:](http://medrxiv.org/content/early/2020/04/11/2020.04.07.20056937/F6) Figure 6: China, left: ![Graphic][47] (in red) and ![Graphic][48] (in magenta), right: *b* (in red) and *b*1 (in magenta) Looking at the timeline for China, we can state stable behavior starting between day 50 and 60, agreeing with the converged behavior in Fig. 1. For South Korea, we have convergence around day 66. Afterwards, the graphs diverge slightly again, which is likely due to the constant number of daily infections occurring in the later stage of the epidemic, already mentioned above. From this observation we can deduce with some caution, that the predictions based on the present model will be reliable to a certain extent as soon as the parameters identified by the two procedures stated in Eqs. (14) and (15) approach each other. Applying this reasoning to the data for Italy, Fig. 9, we can assume reliable predictions starting on day 83. This is supported by the close values for *T*100 in Table 1. Less confidence can be put into the predictions for Germany. In Fig. 8, the graphs for ![Graphic][49] and ![Graphic][50] seem to have converged, but for *b* and *b*1, this cannot be stated with certainty. We will have to wait for the development during the upcoming few days. No convergence can be observed up to now for the United States data, Fig. 10, where we are likely still in a phase of rapid growths of the case numbers. ![Figure 7:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/11/2020.04.07.20056937/F7.medium.gif) [Figure 7:](http://medrxiv.org/content/early/2020/04/11/2020.04.07.20056937/F7) Figure 7: South Korea, left: ![Graphic][51] (in red) and ![Graphic][52] (in magenta), right: *b* (in red) and *b*1 (in magenta) ![Figure 8:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/11/2020.04.07.20056937/F8.medium.gif) [Figure 8:](http://medrxiv.org/content/early/2020/04/11/2020.04.07.20056937/F8) Figure 8: Germany, left: ![Graphic][53] (in red) and ![Graphic][54] (in magenta), right: *b* (in red) and *b*1 (in magenta) ![Figure 9:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/11/2020.04.07.20056937/F9.medium.gif) [Figure 9:](http://medrxiv.org/content/early/2020/04/11/2020.04.07.20056937/F9) Figure 9: Italy, left: ![Graphic][55] (in red) and ![Graphic][56] (in magenta), right: *b* (in red) and *b*1 (in magenta) ![Figure 10:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/11/2020.04.07.20056937/F10.medium.gif) [Figure 10:](http://medrxiv.org/content/early/2020/04/11/2020.04.07.20056937/F10) Figure 10: United States, left: ![Graphic][57] (in red) and ![Graphic][58] (in magenta), right: *b* (in red) and *b*1 (in magenta) ## 6 Conclusion We have identified the parameters in an elementary epidemic model via non-linear regression using data of the covid-19 pandemic. Furthermore, we have attempted to get an insight into the reliability of predictions based on this procedure by observing the timeline of the parameters calculated by two different schemes. Our results indicate, that this approach might work. However, more detailed studies will be necessary in order to establish this method as valid. So caution is required when interpreting the results stated here. In the future, it would be desirable, too, to identify more complex models. It is uncertain, though, if this will be possible at all without more detailed data available. ## Data Availability Data for the COVID-19 pandemic available at the following web sites have been used. [https://www.worldometers.info/coronavirus/](https://www.worldometers.info/coronavirus/) [https://data.humdata.org/dataset/novel-coronavirus-2019-ncov-cases](https://data.humdata.org/dataset/novel-coronavirus-2019-ncov-cases) * Received April 7, 2020. * Revision received April 7, 2020. * Accepted April 11, 2020. * © 2020, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial 4.0 International), CC BY-NC 4.0, as described at [http://creativecommons.org/licenses/by-nc/4.0/](http://creativecommons.org/licenses/by-nc/4.0/) ## References 1. [1].worldometer web page. [https://www.worldometers.info/coronavirus/](https://www.worldometers.info/coronavirus/). 2. [2].Wolfram Research, Inc. Mathematica, Version 12.1. Champaign, IL, 2020. 3. [3]. William Ogilvy Kermack, A. G. McKendrick, and Gilbert Thomas Walker. A con-tribution to the mathematical theory of epidemics. Proceedings of the Royal Society of London. Series A, Containing Papers of a Mathematical and Physical Character, 115(772):700–721, 1927. 4. [4].James D. Murray. Mathematical Biology. Springer, 2002. [1]: /embed/graphic-1.gif [2]: /embed/inline-graphic-1.gif [3]: /embed/graphic-2.gif [4]: /embed/graphic-3.gif [5]: /embed/graphic-4.gif [6]: /embed/graphic-5.gif [7]: /embed/graphic-6.gif [8]: /embed/graphic-7.gif [9]: /embed/graphic-8.gif [10]: /embed/graphic-9.gif [11]: /embed/inline-graphic-2.gif [12]: /embed/inline-graphic-3.gif [13]: /embed/graphic-10.gif [14]: /embed/graphic-11.gif [15]: /embed/graphic-12.gif [16]: /embed/graphic-13.gif [17]: /embed/graphic-14.gif [18]: /embed/graphic-15.gif [19]: F1/embed/inline-graphic-4.gif [20]: F1/embed/inline-graphic-5.gif [21]: F1/embed/inline-graphic-6.gif [22]: F1/embed/inline-graphic-7.gif [23]: F2/embed/inline-graphic-8.gif [24]: F2/embed/inline-graphic-9.gif [25]: F2/embed/inline-graphic-10.gif [26]: F2/embed/inline-graphic-11.gif [27]: F3/embed/inline-graphic-12.gif [28]: F3/embed/inline-graphic-13.gif [29]: F3/embed/inline-graphic-14.gif [30]: F3/embed/inline-graphic-15.gif [31]: F4/embed/inline-graphic-16.gif [32]: F4/embed/inline-graphic-17.gif [33]: F4/embed/inline-graphic-18.gif [34]: F4/embed/inline-graphic-19.gif [35]: F5/embed/inline-graphic-20.gif [36]: F5/embed/inline-graphic-21.gif [37]: F5/embed/inline-graphic-22.gif [38]: F5/embed/inline-graphic-23.gif [39]: /embed/inline-graphic-24.gif [40]: /embed/inline-graphic-25.gif [41]: /embed/graphic-22.gif [42]: /embed/inline-graphic-26.gif [43]: /embed/inline-graphic-27.gif [44]: /embed/inline-graphic-28.gif [45]: /embed/inline-graphic-29.gif [46]: /embed/inline-graphic-30.gif [47]: F6/embed/inline-graphic-31.gif [48]: F6/embed/inline-graphic-32.gif [49]: /embed/inline-graphic-33.gif [50]: /embed/inline-graphic-34.gif [51]: F7/embed/inline-graphic-35.gif [52]: F7/embed/inline-graphic-36.gif [53]: F8/embed/inline-graphic-37.gif [54]: F8/embed/inline-graphic-38.gif [55]: F9/embed/inline-graphic-39.gif [56]: F9/embed/inline-graphic-40.gif [57]: F10/embed/inline-graphic-41.gif [58]: F10/embed/inline-graphic-42.gif