A Logistic Curve in the SIR Model and Its Application to Deaths by COVID-19 in Japan ==================================================================================== * Takesi Saito * Kazuyasu Shigemoto ## Abstract Approximate solutions of SIR equations are given, based on a logistic growth curve in the Biology. These solutions are applied to fix the basic reproduction number *α* and the removed ratio *c*, especially from data of accumulated number of deaths in Japan COVID-19. We then discuss the end of the epidemic. These logistic curve results are compared with the exact results of the SIR model. ## 1 Introduction The SIR model [1] in the theory of infection is powerful to analyze an epidemic about how it spreads and how it ends [2–8]. The SIR model is composed of three equations for *S, I* and *R*, where they are numbers for susceptibles, infectives and removed, respectively. Three equations can be solved completely by means of MATHEMATICA, if two parameters *α* and *c* are given, where *α* is the basic reproduction number and *c* the removed ratio. In Sec. 2 we would like to summarize some exact solutions of the SIR equations. These exact solutions are applied to COVID-19 in Japan. Here, our policy is a little use of data of cases. In Sec. 3 we propose approximate solutions of SIR equations, based on the logistic growth curve in the Biology [9]. These approximate solutions have simple forms, so that they are very useful to discuss an epidemic. The final section is devoted to concluding remarks. The logistic approach is compared with exact solutions of the SIR model to our epidemic. ## 2 The SIR model in the theory of infection Equations of the SIR model are given by ![Formula][1] ![Formula][2] ![Formula][3] where *S, I* and *R* are numbers for susceptibles, infectives and removed, respectively, *b* the infection ratio and *c* the removed ratio. Our aim is to propose an approximate solutions of these equations, based on a logistic growth curve in the Biology. We first summarize some of exact solutions of the SIR model. From Eq(2.1) and Eq.(2.3) we get *dS*/*dR* = −*αS*, (*α* = *b*/*c*), which is integrated to be ![Formula][4] where *α* stands for the basic reproduction number. In the same way, from Eq.(2.2) and Eq.(2.4), we have *dI*/*dR* = *αS* − 1 = *α* exp(−*αR*) − 1, which is integrated to be ![Formula][5] The solutions Eq.(2.4) and Eq.(2.5) satisfy boundary conditions *S* = 1 and *I* = 0 at *R* = 0. We normalize the total number to be unity, i.e., *S* + *I* + *R* = 1. Since *I*(*t*) is 0 when *t ⟶* ∞, we have 1 − *R*(∞) − exp [−*αR*(∞)] = 0 from Eq.(2.5). Hence, it follows a useful formula ![Formula][6] Some exact formulas at the peak *t* = *T* are summarized as follows: ![Formula][7] ![Formula][8] The first one is derived as follows: Since the peak point is given by *dI*/*dt* = *c*(*αS* − 1)*I* = 0 at *t* = *T*, we get *S*(*T*) = exp[−*αR*(*T*)] = 1/*α*. This proves Eq.(2.7). The second equation (2.8) directly follows by *I*(*T*) = 1 − *R*(*T*) − *S*(*T*) = 1 − ln *α*/*α* − 1/*α*. Now let us consider an application of the above exact results in the SIR model to COVID-19 in Japan. In order to fix *T, α* and *c*, we use the data of deaths [10]. At May 2nd, the accumulated number of deaths *D* takes 492 and the new increased number of deaths Δ*.D*(*t*)/Δ*.t* takes the maximum value 34. Then we find *T* to be May 2nd. Here, our policy is a little use of data of cases. Then we connect *D*(*t*) with *R*(*t*) at *t* = *T* by the formula *D*(*t*) = *rR*(*t*), where *r* is the death rate, so that ![Formula][9] Since Δ*.D*(*t*)/Δ*.t* is fluctuating, we take 5-day average from April 30th to May 4th. The 5-day values for Δ*.D*(*t*)/Δ*.t* are (April 30th: 17), (May 1st: 26), (May 2nd: 34), (May 3rd :18), (May 4th: 11), and the 5-day average around May 2nd gives Δ*.D*/Δ*.t*|5-day average = 21.2. By using *D*(*T*) = 492, we have (Δ*.D*/Δ*.t*)/*D*|*T* = 21.2/492 = 0.043. Next, let us fix *c* from the formula *c* = (Δ*.R*(*t*)/Δ*.t*)/*I*(*t*). We take the total average from Feb. 16 to June 9, in order to compensate the fluctuation. Thus we get *c* = (Δ*.R*(*t*)/Δ*.t*)/*I*(*t*)|total average = 0.041, which means that almost all cases are rejected approximately after 25 days. Substituting *c* = 0.041 into Eq.(2.9), we have an equation for *α* ![Formula][10] which gives a solution *α*exact = 3.66. To sum up, we have exact values of parameters: *T* =May 2nd, *α*exact = 3.66 and *c* = 0.041. We then draw curves of *S, I* and *R* by means of MATHEMATICA in Fig. 1. ![Figure 1:](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/06/27/2020.06.25.20139865/F1.medium.gif) [Figure 1:](http://medrxiv.org/content/early/2020/06/27/2020.06.25.20139865/F1) Figure 1: Graph of *S, I* and *R* for α = 3.66 and *c* = 0.041 ## 3 A logistic curve from the SIR model Let us consider approximate formulas for SIR’s functions. The third equation Eq.(2.3) can be written as ![Formula][11] with *t* the true time. Let us expand the exponential factor in the second order of *x* = *αR*, ![Formula][12] Then we have ![Formula][13] with *A* = *α*2/2, *B* = 2(*α* − 1)/*α*2. This equation is a type of the logistic growth curve in the Biology [9], easily solved as ![Formula][14] where *AB* = *α* − 1 = *a* and *B* = *R*(∞) = 2(*α* − 1)/*α*2 = 2*R*(*T*). Inserting *R*(*t*) into *I*(*t*) = (1/*c*)(*dR*(*t*)/*dt*), we have ![Formula][15] where the peak value of *I*(*t*) is given by *I*(*T*) = *aR*(∞)/4 at *t* = *T* . Here we have useful formulas ![Formula][16] In the following, we consider an application of the logistic curve to COVID-19 in Japan. By using Eq.(2.9), Eq.(3.4) and Eq.(3.5), we have ![Formula][17] According to the data of accumulated number of deaths, we have (Δ*.D*/Δ*.t*)/*D*|*T* = 21.2/492 = 0.043, which gives *α*appr = *a* + 1 = 3.10. In this way, we have fixed parameters in the logistic curve: *T* =May 2nd, *α*appr = 3.10 and *c* = 0.041. The value *α*appr = 3.10 should be compared with *α*exact = 3.66. The error of our approximation for *α* is about 15%. ## 4 Concluding remarks We have proposed logistic formulas in the theory of infection, which are approximate formulas driven from the SIR model. These formulas appear to be simple forms, therefore, very useful to analyze an epidemic. In logistic formulas we have determined the basic reproduction number *α*appr = 3.10 with the removed ratio *c* = 0.041, from data of the accumulated number of deaths in Japan COVID-19. One can see *D*(∞) = 984 from the approximate formula *D*(∞) = 2*D*(*T*) = 2 × 492, that is, the final accumulated number of deaths is 984. We have also fixed exactly *α* and *c*, in the SIR model for Japan COVID-19. Once having the basic reproduction number *α*exact = 3.66 and the removed ratio *c* = 0.041, we can draw curves of *S, I* and *R* by means of MATHEMATICA in Fig. 1. The peak day of Δ*.D*(*t*)/Δ*.t* is *T* =May 2nd from the data, which is the 77th day from Feb.15. Since Δ*.D*(*t*)/Δ*.t* is almost left-right symmetric with respect to the peak day, the epidemic ends on 77 days after May 2nd, a middle of July. Finally if we compare *α*appr = 3.10 with *α*exact = 3.66, the error of our approximation for *α* is about 15%. ## Data Availability “Summary of the New Coronavirus Infection ”(in Japanese), [https://hazard.yahoo.co.jp/article/20200207](https://hazard.yahoo.co.jp/article/20200207); “Status of the Domestic New Coronavirus Infection”(in Japanese), [https://toyokeizai.net/sp/visual/tko/covid19/](https://toyokeizai.net/sp/visual/tko/covid19/) ## Footnotes * * tsaito{at}k7.dion.ne.jp * Received June 25, 2020. * Revision received June 26, 2020. * Accepted June 27, 2020. * © 2020, Posted by Cold Spring Harbor Laboratory The copyright holder has placed this preprint in the Public Domain. It is no longer restricted by copyright. Anyone can legally share, reuse, remix, or adapt this material for any purpose without crediting the original authors. ## References 1. [1]. W. O. Kermack and A.G. McKendrick, A.G. (1927), “A Contribution to the Mathematical Theory of Epidemics”, Proceedings of the Royal Society A115 (1927), 700–721 (1927). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1098/rspa.1927.0118&link_type=DOI) 2. [2]. W. D. Murray, “Epidemic models and the dynamics of infectious diseases”, Mathematical Biology 42, 610–650 (1993). 3. [3]. H. Hethcote, “The Mathematics of Infectious Diseases”, SIAM Review 42, 599–653 (2000). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1137/s0036144500371907&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=WOS:00016567&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F06%2F27%2F2020.06.25.20139865.atom) 4. [4]. M.J. Keeling and L. Danon, “Mathematical modelling of infectious diseases “, British Medical Bulletin 92, 33–42 (2009). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bmb/ldp038&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19855103&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F06%2F27%2F2020.06.25.20139865.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000272684600002&link_type=ISI) 5. [5]. S. Pathak, A. Maiti and G.P. Samanta, “Rich dynamics of an SIR epidemic model”, Nonlinear Analysis: Modelling and Control 15, 71–81 (2010). 6. [6]. T. Harko, F.S.N. Lobo, M.K. Mak, “Exact analytical solutions of the Susceptible-Infected-Recovered (SIR) epidemic model and of the SIR model with equal death and birth rates”, Applied Mathematics and Computation 236, 184194 (2014). 7. [7]. J.C. Miller, “Mathematical models of SIR disease spread with combined non-sexual and sexual transmission routes”, Infectious Disease Modelling 2, 3555 (2017). 8. [8]. R. Sameni, “Mathematical Modeling of Epidemic Diseases; A Case Study of the COVID-19 Coronavirus”, [arxiv:2003.11371 [q-bio.PE]](2017). 9. [9]. P-F. Verhulst, “Notice sur la loi que la population poursuit dans son accroissement”, Correspondance Mathmatique et Physique 10, 113–121 (1838). 10. [10].Data published from the Ministry of Health, Labor and Welfare of Japan. See for example, “Summary of the New Coronavirus Infection “(in Japanese), [https://hazard.yahoo.co.jp/article/20200207](https://hazard.yahoo.co.jp/article/20200207); “Status of the Domestic New Coronavirus Infection”(in Japanese), [https://toyokeizai.net/sp/visual/tko/covid19/](https://toyokeizai.net/sp/visual/tko/covid19/). [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-6.gif [7]: /embed/graphic-7.gif [8]: /embed/graphic-8.gif [9]: /embed/graphic-9.gif [10]: /embed/graphic-10.gif [11]: /embed/graphic-12.gif [12]: /embed/graphic-13.gif [13]: /embed/graphic-14.gif [14]: /embed/graphic-15.gif [15]: /embed/graphic-16.gif [16]: /embed/graphic-17.gif [17]: /embed/graphic-18.gif