A Logistic Curve in the SEIR Model and the Basic Reproduction Number of COVID-19 in Japan ========================================================================================= * Takesi Saito * Kazuyasu Shigemoto ## Abstract The SEIR model is one of modified models of SIR, especially taken into account of exposed people. SEIR equations can be solved numerically, but it is hard to obtain analytically. Here, we propose some approximate solutions of SEIR equations, one of which is related with the logistic formula in Biology. As the second aim, the SEIR model is applied to the 7th-wave of COVID-19 in Japan. The basic reproduction number (*α*) in the SEIR model is estimated for the Omicron wave. We make use of data of the removed number *R*(*t*) rather than that of the infective number, because the latter seems to be ambiguous. This analysis gives *α* = 10 with *γ* = 1 and σ = 0.5. Keywords * Variants of SARS-COV-2 * SIR model * Epidemiology * Molecular biology ## 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 [3-16]. In this article we consider the SEIR model [2] which is one of modified models of SIR, especially taken into account of exposed people. SEIR equations can be solved numerically, but it is hard to obtain analytically. In order to analyze infections, analytic solutions are important. Here, therefore, we propose some approximate solutions of SEIR equations, one of which is related with the logistic formula in Biology. As our second aim, the SEIR model is applied to the 7th-wave of Omicron-COVID-19 in Japan. The basic reproduction number (*α*) in the SEIR model is estimated for the 7th - wave. We make use of data of the removed number *R*(*t*) rather than that of the infective number, because the latter seems to be ambiguous. According to this analysis we have *α* = 10 with *γ* = 1 and *σ* = 0.5, where *γ* and *σ* are the removed and the exposed ratios, respectively. ## 2. A logistic formula from the SEIR model The sequential SEIR equations are given by ![Formula][1] ![Formula][2] ![Formula][3] ![Formula][4] where *S*(*t*), *E*(*t*), *I*(*t*) and *R*(*t*) are numbers for susceptible, exposed, infectious and recovered, respectively, *α* = *β*/*γ* the basic reproduction number, *σ* the exposed ratio and *γ* the removed ratio. Four numbers are normalized as ![Formula][5] The SEIR equations can be solved numerically, when parameters *β, γ and σ* are given. In Fig.1 we give four curves for *S*(*t*), *E*(*t*), *I*(*t*) and *R*(*t*) with *β* = 10, *γ* = 1 *and σ* = 0.5. Here we propose some approximate solutions of SEIR equations, one of which is related with a logistic formula in Biology. From Eqs. (1) and (4) we get *dS*/*dR* = −*αS*, which is integrated to be ![Fig. 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/09/26/2022.09.18.22279896/F1.medium.gif) [Fig. 1.](http://medrxiv.org/content/early/2022/09/26/2022.09.18.22279896/F1) Fig. 1. Four curves for *S*(*t*), *E*(*t*), *I*(*t*) and *R*(*t*) with *β* = 10, *γ* = 1 *and σ* = 0.5. ![Formula][6] Hence we have from Eq. (5) ![Formula][7] Now, from Eq. (3) *E*(*t*) can be expressed as ![Formula][8] Substituting *I*(*t*) = *γ*−1*dR*(*t*)/*dt* of Eq. (4) into Eq. (8) above, we get ![Formula][9] so that Eq. (7) reduces to the second order differential equation of *R*(*t*). ![Formula][10] In order to find approximate solution of the exact equation (10), we consider three cases: Case 1: |*d*2*R*(*t*)/*dt*2 |≪ |*dR*(*t*)/*dt* |, *R*(*t*) In this case the first term in Eq. (10) can be neglected so that we have ![Formula][11] When *σ* ≫ *γ*, we find that this equation reduces to the SIR equation: ![Formula][12] Case 2: When *αR*(*t*) ≪ 1, we can use the approximate formula exp(−*αR*) ≅ 1− *αR* + *α*2*R*2/2, then Eq. (10) reduces to ![Formula][13] Case 3: A combination of Case 1 and Case 2 gives ![Formula][14] In Fig. 2, we draw curves of approximate solutions for Cases 1, 2 and 3, against the Exact equation (10). Here we find that the contribution of the second order differential part, *d*2*R*(*t*)/*dt*2, into any case is very small. The last equation (14), Case 3, is nothing but the logistic equation. ![Fig. 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/09/26/2022.09.18.22279896/F2.medium.gif) [Fig. 2.](http://medrxiv.org/content/early/2022/09/26/2022.09.18.22279896/F2) Fig. 2. Curves of approximate solutions for Cases 1, 2 and 3, against the Exact equation. solution. The Case 3 is for the logistic equation. ## 3. The basic reproduction number of COVID-19 in Japan Let us estimate the basic reproduction number in the SEIR model for the 7th-wave of COVI-19 in Japan. We make use of data [17] of the removed number *R*(*t*) rather than that of the infective number, because the latter seems to be ambiguous. In Fig. 3, the red curve of ![Fig. 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/09/26/2022.09.18.22279896/F3.medium.gif) [Fig. 3.](http://medrxiv.org/content/early/2022/09/26/2022.09.18.22279896/F3) Fig. 3. The red line shows a curve of *R*(*t*) with *α* = 10, *γ* = 1 *and σ* = 0.5. Blue bar charts stand for the 7th-wave data of *R*(*t*). *R*(*t*) for the 7th-wave is drawn with *α* = 10, *γ* = 1 *and σ* = 0.5 in the SEIR model. Blue bar charts stand for the 7th-wave data of *R*(*t*). The number on the horizon axis corresponds with the date for each in Table 1. View this table: [Table 1.](http://medrxiv.org/content/early/2022/09/26/2022.09.18.22279896/T1) Table 1. The number above corresponds with the date below for each. In Fig. 3, we see that a coincidence between both curves is fairly well. So, we conclude that the basic reproduction number in the SEIR model is *α* = 10 *withγ* = 1 *and σ* = 0.5 for the 7th-wave of COVID-19 in Japan. ## 4. Conclusion We have found some approximate solutions of SEIR equations. First, when *σ* ≫ *γ* in Case 1, the SEIR equation reduces to the SIR equation (12). Approximate solutions are given by Case1, Case 2 and Case 3, and their curves are drawn in Fig. 2, where Case 3 is known to be the logistic formula in Biology. Here we find that the contribution of the second order differential part, *d*2*R*(*t*)/*dt*2, into any case is very small. The logistic equation (14) can be rewritten as ![Formula][15] Then its explicit solution is easily obtained by ![Formula][16] Hence, we have useful formulas ![Formula][17] from Eq. (4) in Sec. 2, and ![Formula][18] from Eq. (9) in Sec. 2, by neglecting the second order differential part, *d*2*R*(*t*)/*dt*2. The second aim is the estimation of *α*, the basic reproduction number in the SEIR model for the 7th-wave of Omicron COVID-19 in Japan. We have made use of data of the removed number *R*(*t*) rather than that of the infective number, because the latter seems to be ambiguous. We have obtained *α* = 10 *with γ* = 1 *and σ* = 0.5. ## Data Availability All data produced are available online at Yahoo Japan News page "Summary of the New Coronavirus Infection"(in Japanese), https://news.yahoo.co.jp/pages/article/20200207 and Toyokeizai net page "Status of the Domestic New Coronavirus Infection"(in Japanese), https://toyokeizai.net/sp/visual/tko/covid19/ [https://news.yahoo.co.jp/pages/article/20200207](https://news.yahoo.co.jp/pages/article/20200207) [https://toyokeizai.net/sp/visual/tko/covid19/](https://toyokeizai.net/sp/visual/tko/covid19/) ## Acknowledgment We would like to express our deep gratitude to Kokado A. for valuable discussions. ## Footnotes * Title is changed; The reference [12] is supplemented. * Received September 18, 2022. * Revision received September 24, 2022. * Accepted September 26, 2022. * © 2022, 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 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]. M. Li, H. L. Smith and L. Wang, “Global Dynamics of an SEIR Epidemic Model with Vertical Transmission”, SIAM J. Appl. Math. 62(1), 58-69(2001) 3. [3]. W. D. Murray, “Epidemic models and the dynamics of infectious diseases”, Mathematical Biology 42, 610–650(1993). 4. [4]. 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%2F2022%2F09%2F26%2F2022.09.18.22279896.atom) 5. [5]. 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%2F2022%2F09%2F26%2F2022.09.18.22279896.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000272684600002&link_type=ISI) 6. [6]. S. Pathak, A. Maiti and G. P. Samanta, “Rich dynamics of an SIR epidemic model”, Nonlinear Analysis: Modelling and Control 15, 71–81(2010). 7. [7]. T. Harko, F. S. N. Lobo, M. K. Mak, “Exact analytical solutions of the SusceptibleInfected-Recovered (SIR) epidemic model and of the SIR model with equal death and birth rates”, Applied Mathematics and Computation 236, 184194(2014). 8. [8]. J. C. Miller, “Mathematical models of SIR disease spread with combined non-sexual and sexual transmission routes”, Infectious Disease Modelling 2, 3555(2017). 9. [9]. R. Sameni, “Mathematical Modeling of Epidemic Diseases; A Case Study of the COVID19 Coronavirus”, [arXiv:2003.11371 [q-bio.PE]](2017). 10. [10]. P-F. Verhulst, “Notice sur la loi que la population poursuit dans son accroissement”, Correspondance Mathmatique et Physique 10, 113–121(1838). 11. [11]. T. Saito and K. Shigemoto, “A Logistic Curve in the SIR Model and Its Application to Deaths by COVID-19 in Japan”, medRxiv doi:10.1101/2020.06.25.20139865. European Journal of Applied Sciences 10(5), 157-160(2022). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1101/2020.06.25.20139865&link_type=DOI) 12. [12]. T. Saito, “An Application of Logistic Formula to Deaths by COVID-19 in Japan”, medRxiv doi:10.1101/2020. 09.15.20195081. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1101/2020&link_type=DOI) 13. [13]. A. Kokado and T. Saito, “A Logistic Formula and Its Application to Deaths by the Third Wave of COVID-19 in Japan”, medRxiv doi:10.1101/2021.01.30.21250827, [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMS4wMS4zMC4yMTI1MDgyN3YxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDkvMjYvMjAyMi4wOS4xOC4yMjI3OTg5Ni5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 14. [14]. T. Saito, “Variants of SARS-COV-2 and the Death Toll”, Integ. Mol. Bio. Biotechnol. 1, 1–3(2021). 15. [15]. T. Saito, “Area Theorem of Infection Curves for Large Basic Reproduction Number”, European Journal of Applied Sciences 10(1), 482–486(2022). 16. [16]. T. Saito, “Analyses of COVID-19 in Japan by the Mathematical Model of SIR”, EJAS,10(3), 784–796(2022). 17. [17].Data published from the Ministry of Health, Labor and Welfare of Japan. See for example, “Summary of the New Coronavirus Infection”(in Japanese), [https://news.yahoo.co.jp/pages/article/20200207](https://news.yahoo.co.jp/pages/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-7.gif [7]: /embed/graphic-8.gif [8]: /embed/graphic-9.gif [9]: /embed/graphic-10.gif [10]: /embed/graphic-11.gif [11]: /embed/graphic-12.gif [12]: /embed/graphic-13.gif [13]: /embed/graphic-14.gif [14]: /embed/graphic-15.gif [15]: /embed/graphic-19.gif [16]: /embed/graphic-20.gif [17]: /embed/graphic-21.gif [18]: /embed/graphic-22.gif