The reproductive index from SEIR model of Covid-19 epidemic in Asean ==================================================================== * Pongkeaw Udomsamuthirun * Grittichon Chanilkul * Pongkarn Tongkhonburi * Chatcharawan Meesubthong ## Abstract As we calculate analytic to link the coefficient of third-order polynomial equations from raw data of an Asean to the SEIR model. The Reproductive index depending on the average incubation period and the average infection period and the coefficient polynomial equations fitted from raw are derived. We also consider the difference of the average incubation period as 5 days and 3 days with the average infection period as 10 day of an Asean. We find that the value of *R* are Indonesia (7.97), Singapore (6.22), Malaysia (3.86), Thailand (2.48), respectively. And we also find that Singapore has 2 values of *R* as 1.54 (16 Feb to 37 March) and 6.22 (31 March-4 April).The peak of infection rate are not found for Singapore and Indonesia at the time of consideration. The model of external stimulus is added into raw data of Singapore and Indonesia to find the maximum rate of infection. We find that Singapore need more magnitude of external stimulus than Indonesia. And the external stimulus for 14 days can stimulate to occur the peak of infected daily case of both country. Keywords * SEIR model * The Reproductive index * Covid-19 ## 1. Introduction An Asean (The Association of Southeast Asian Nations) is consist of 10 ten member states as Indonesia, Malaysia, Philippines, Singapore, Thailand, Brunei Darussalam, Viet Nam, Lao PDR, Myanmar, and Cambodia. And we try to predict the number of coronavirus (Covid-19) victims as number of persons who caught the infection and got sick only in this area. The complicated mathematical models are necessary for long-time predictions. The SIR (Susceptible-Infectious-Removed) model are used to obtain the prediction values of the model parameters using the statistical approach for predication the number of infected, susceptible and removed persons [1,2]. The Susceptible-Infectious-Recovered/Death (SIRD) Model was used to formulate an optimal control problem with an expanded epidemic model to compute (Non-pharmaceutical) implementation strategy.[3] A modify SIR called SEIRUS model (Susceptible – Exposed – Infectious – Removed – Undetectable – Susceptible) is generated for evaluate the new deterministic pandemic Covid-19 endemic that originally developed for the control of the prevalence of HIV/AIDS in Africa.[4] The Susceptible-Exposed-Infectious-Removed (SEIR) model was adopt to be SEIRNDC [5] that the total population size N with two extra classes “D” mimicking the public perception of risk regarding the number of severe and critical cases and deaths; and “C” representing the number of cumulative cases. This model proposed the compartmental model that sustained human-to-human transmission of Covid-19 after December 2019 of Wuhan, China. In this research, we use the daily cases and total cases of Covid-19 infection population of an Asean from website: Worldometer [6]: [https://www.worldometers.info/coronavirus/](https://www.worldometers.info/coronavirus/) between 15/02/2020 to 15/04/2020. The raw data are fitted with the regression method of third-order polynomial formula. We calculate analytic to link the coefficient of polynomial to the Reproductive index (R) of the SEIR model (Susceptible-Exposed-Infected-Removed). ## 2. Model We apply the well-known SEIR compartmental model [7] (Susceptible-Exposed-Infected-Removed) for the prediction properties of how a disease spread. The variable description is *S*(*t*) is number of susceptible populations, E(t) is the number of exposed populations, *I*(*t*) is number of infected populations, *R*(*t*) is number of infected populations quarantined and expecting recovery at time. There is no emigration from the total population and there is no immigration into the population. A negligible proportion of individuals move in and out of at a given time that ![Formula][1] The people susceptible are able to get infected when they contact infectious people. Once infected, they move into the infectious compartment. People recovered are assumed to be immune within our study horizon. Then ![Formula][2] The dynamics non-linear differential equations are below ![Formula][3] ![Formula][4] ![Formula][5] ![Formula][6] Here, the *β, σ* and *γ* denoted the infection rate, the onset rate, and the removal rate. The ![Graphic][7] and ![Graphic][8] imply the average incubation period and the average infection period. Base on Ref.[8,9], the *σ* = 0.2 and *γ* = 0.1. The normalize parameters of model are defined by ![Formula][9] In our model, we assume that the rate of infection is in the form of third-order polynomial formula as ![Formula][10] Integrate with respect to time, we can get the total accumulation of infective population that ![Formula][11] Here *a* is the total cases found in the day before. *a*1, *a*2, *a*3 *and a*4 are constant that found from regression method on the raw data from website: Worldometer at [https://www.worldometers.info/coronavirus/](https://www.worldometers.info/coronavirus/). After take some calculate on solving the first order differential equation on Eqs.(1-6), we get ![Formula][12] ![Formula][13] ![Formula][14] Substitution Eq.(6), Eq.(7) and Eq.(8) in to Eq.(2), we get ![Formula][15] By setting t=0, and ![Graphic][16], the reproductive index in the parameters of raw data is derived as ![Formula][17] The reproductive index is defined as the dominant eigenvalue of the matric *FV* −1; the linearized infection subsystem F and V. In our model, we assume that there are two disease state but one way to create new infection, *X* = (*e,i*). The linearized infection subsystem F and V are in form ![Graphic][18] that R is the dominant eigenvalue of the matrix *FV* −1. From SEIR model, the linearized infection subsystem are ![Graphic][19] and ![Graphic][20]. Then, we can get ![Graphic][21]. Our R is agreed with Ref. [10,11]. The reproductive rate or the reproductive index, R is the course of an epidemic shape. It represents the number of further cases each new case will give rise to. For high value of R, the number of newly infected people climbs more quickly to a maximum than low value of *R*. The higher R, the higher infectious population are. The World Health Organization (WHO) reported an incubation period for COVID-19 between 1 and 14 days but most commonly around five days. The most infection period is thought to be 1 to 3 days before symptoms start, and in the first 7 days after symptoms begin. Mean incubation period observed are 3.0 days (0 - 24 days range case-study) and 5.2 days (4.1 - 7.0 days range case-study) [6]. Base on Ref.[8] and [9], the *σ* = 0.2 and *γ* = 0.1, or ![Graphic][22]. The Reproductive index becomes ![Formula][23] And base on WHO, we set ![Graphic][24], that ![Formula][25] Here *a* is the accumulation infected population before the date of calculation. *a*1 *and a*2 are constant values and the coefficient of term “t” of ![Graphic][26], respectively. ## 3. Results and Discussion When the spreading of viruses in a population occurred, the number of new cases rises rapidly, peaks, and then declines called the epidemiological curve. The spreading curve should be the flatten as the spreading the infections out of time. In this paper, we model the shape of spreading rate to be the third-order polynomial that we called “bell-shape “. According to the total accumulate infected cases, the shape of curve should be gotten the saturation value so that the rate of infectious should be the bell-shape at the end of disease spreading. After analysis the raw data, we find that the daily new cases of Asean are divided into 2 cases; the bell-shape and no bell-shape. The bell-shape cases mean that the infection is in the state of the beginning of the saturated status as we can find the maximum value and we cannot find the maximum value for no bell-shape case. However, both cases can be solve to find the *R*. An Asean is consist of 10 ten member states as Indonesia, Malaysia, Philippines, Singapore, Thailand, Brunei Darussalam, Viet Nam, Lao PDR and Myanmar, and Cambodia. Because of not enough data to calculate efficiency, we ignore Lao PDR and, Myanmar, Cambodia for our consideration. Here, the raw data fitted well with *R**fit* 2 > 0.80 are used to find the reproductive index so there are only 4 country suitable as Malaysia, Singapore, Thailand and Indonesia. From Eq.(10), the Reproductive index (*R*) depend on the average incubation period ![Graphic][27] and the average infection period ![Graphic][28]. Upon to the difference on the average incubation period as 5 days and 3 days with the average infection period as 10 day, Eq.(10 become to Eq.(11) and Eq.(12) as ![Graphic][29] and ![Graphic][30]. And *a* is the accumulation infected population before the date of calculation [6]. The others constant, *a*1 *and a*2, come from our fitting equation. The results are shown in Table 1. View this table: [Table 1.](http://medrxiv.org/content/early/2020/04/29/2020.04.24.20078287/T1) Table 1. The Reproductive index (*R*) of Malaysia, Thailand Singapore and Indonesia. As there are only 4 countries that have enough data to fit well with third-order polynomial; Malaysia, Thailand, Singapore, and Indonesia. There for, Philippines also has enough data but raw data show the swing-type behavior so the third-order polynomial equation cannot fit well, while the other country have not enough data to conform our condition, ![Graphic][31]. We find that Singapore having 2 values of as 1.54 (16 Feb to 37 March) and 6.22 (31 March-4 April). The higher are Indonesia (7.97), Singapore (6.22), Malaysia (3.86), Thailand (2.48), respectively. The higher *R* are corresponded to the higher total case as the total cases are found; Singapore (6588 cases), Indonesia (6575 cases), Malaysia (5384 cases), Thailand (2765 cases) (Worldmeter on 19 April 2020.). The more incubation period, the more value of are found. After taking the regression method to raw data, we find that Malaysia, Philippines, Thailand, Brunei Darussalam, Viet Nam and Cambodia show the peak of the spreading. We take the extrapolate on equations to reach the end of time so “Bell-shape” occurred. To compare with the equal peak, we normalized the results with the maximum of each country and shown in Figure 1. ![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/29/2020.04.24.20078287/F1.medium.gif) [Figure 1.](http://medrxiv.org/content/early/2020/04/29/2020.04.24.20078287/F1) Figure 1. The normalized peak of spreading curve of an Asean. Most of these countries found first-case on January; Thailand, Singapore, Viet Nam, Malaysia, Cambodia, Philippines and the other found on March; Indonesia, Brunei Darussalam, Lao PDR, Myanmar [Wikipedia]. It takes a time to be outbreak to global pandemic depending on condition of each country. In Figure. 1, the maximum of the spreading peak is found in chronological order as Brunei Darussalam, Cambodia, Philippines, Viet Nam, Thailand, and Malaysia, respectively. In this case Thailand and Malaysia have, ![Graphic][32], so the equation are well enough for consideration. The graph of Thailand has more slope than the graph of Malaysia and the total case of Malaysia (5384 cases) are higher than Thailand (2765 cases) (Worldometer on 19 April 2020.) [6]. This result is agree with more R of Malaysia than R of Thailand as Table 1. The raw data of Indonesia and Singapore have shown that these countries take more time than the other country to get a peak. So, their curves do not show the Bell-shape at the time of consideration. However, we can calculate the R of them that are 6.22 and 7.97 for Singapore and Indonesia, respectively. The high value of the R mean that they will take more time to get a peak and the infective population will be abundant if government do nothing. In this situation, the external stimulus is need for stimulate the system to exhibit the new character such as the FitzHugh - Nagumo equations [12,13]. The Eq.(5) is similar to the FitzHugh - Nagumo equations then we add the external stimulus, *I* *ext* into Eq.(5), ![Formula][33] The external stimulus will reset the variable of system to become new value. If the external stimulus exceeds a certain threshold value, the system will generate a new behavior. In our model, we need the external stimulus to reduce the exponential increase to be the Bell-shape form. So, we set ![Graphic][34], that is the condition for maximum value. The external stimulus is the policy of public-health organization of each country. The aim of public-health policy is to flatten the curve and decrease the maximum infected population of spreading over time. The basic protective measures against the Covid-19 such as wash your hands frequently and avoid touching eyes, nose and mouth (WHO). In Asean, aside from the basic suggestions from WHO, many policies are done for decrease the spreading of Covid-19. A common policy used are travel restriction, social distancing and work at home. The difference is the time of these policy beginning that depend on first case confirmed and the attention of government and people to solve this problem. The medical supply and living conditions are still the problem of this area [Wikipedia]. In this paper, we take the constant external stimulus as ![Graphic][35]. And this constant is equal to the latest value of daily infection per day. From Figure 2 and Figure 3, Singapore need more the magnitude of external stimulus than Indonesia. However, the external stimulus for 14 days can stimulate to occur the peak of infective daily case of both country. So, the effective policy for decrease the spreading of Covid-19 should be at least 14 days. In Singapore the preventive measures were collectively called a ‘circuit breaker’. But there is still some unsuccessfulness of this policy. The small territory of this country makes some problems such as living conditions at foreign worker dormitories and social distancing on public transport is not long enough. In Indonesia, there are many policies were done by central government such as travel restrictions, making medical facility specially equipped to handle the coronavirus, social distancing, work from home and having stimulus policy on economic sector. [Wikipedia]. ![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/29/2020.04.24.20078287/F2.medium.gif) [Figure 2.](http://medrxiv.org/content/early/2020/04/29/2020.04.24.20078287/F2) Figure 2. The dot “Normal” are raw data of daily case of Indonesia [6]. The “Cal” is from the regression. “Cal Act 3 day” is regression with the external stimulus for 3 days. “Cal Act 14 day” is regression with the external stimulus for 14 days. ![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/29/2020.04.24.20078287/F3.medium.gif) [Figure 3.](http://medrxiv.org/content/early/2020/04/29/2020.04.24.20078287/F3) Figure 3. The dot “Normal” are raw data of daily case of Singapore [6]. The “Cal” is from the regression. “Cal Act 7 day” is regression with the external stimulus for 7 days. “Cal Act 14 day” is regression with the external stimulus for 14 days. ## 4. Conclusions The Covid-19 raw data of infected rate of Asean are fitted with the regression method of third-order polynomial formula under the scope of SEIR model. We derived the reproductive index and set into 2 equations as ![Graphic][36] and ![Graphic][37] for the average incubation period as 5 days and 3 days with the average infection period as 10 day, respectively. The raw data show that there are only 4 country having enough data to fit well with third-order polynomial as Malaysia, Thailand, Singapore, and Indonesia. The higher *R* are Indonesia (7.97), Singapore (6.22), Malaysia (3.86), Thailand (2.48), respectively. And we also find that Singapore having 2 values of *R* as 1.54 (16 Feb to 37 March) and 6.22 (31 March-4 April). Our estimated R are in range that consistent with other preliminary estimates posted on public domains [14–17]. However, the peak of infected rate is not found for Singapore and Indonesia at the time of consideration. The model of external stimulus on raw data are added into polynomial equation of Singapore and Indonesia to control the Covid-19. We find that Singapore need more the magnitude of external stimulus than Indonesia. However, the external stimulus for 14 days can stimulate to occur the peak of infective daily case of both country. ## Data Availability Worldometer from [https://www.worldometers.info/coronavirus/](https://www.worldometers.info/coronavirus/) ## Acknowledgments We would like to thank to Srinakharinwirot University and Huachiew Chalermprakiet University for support at all the time of doing research. * Received April 24, 2020. * Revision received April 24, 2020. * Accepted April 29, 2020. * © 2020, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution 4.0 International), CC BY 4.0, as described at [http://creativecommons.org/licenses/by/4.0/](http://creativecommons.org/licenses/by/4.0/) ## References 1. [1].Nesteruk, I. (2020). Statistics-based predictions of coronavirus 2019-nCoV spreading in mainland China. Innov Biosyst Bioeng 4(1):13–18 doi: 10.20535/ibb.2020.4.1.195074. Taylor & Francis. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.20535/ibb.2020.4.1.195074&link_type=DOI) 2. [2].Nesteruk I. (2017). Statistics based models for the dynamics of Chernivtsi children disease. Naukovi Visti NTUU KPI. 5:26–34. doi: 10.20535/1810-0546.2017.5.108577 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.20535/1810-0546.2017.5.108577&link_type=DOI) 3. [3].Lin Feng.; Muthuraman K. and Lawley M. (2010). An optimal control theory approach to non-pharmaceutical interventions. BMC Infectious Diseases 10:32. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1471-2334-10-32&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20170501&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F04%2F29%2F2020.04.24.20078287.atom) 4. [4].Okhuese V.A. (2020). Mathematical predictions for Covid-19 as a global pandemic. [https://doi.org/10.1101/2020.03.19.20038794doi](https://doi.org/10.1101/2020.03.19.20038794doi) 5. [5].Lin Q., et al. (2020). A conceptual model for the coronavirus disease 2019 (COVID-19) outbreak in Wuhan, China with individual reaction and governmental action. International journal of infectious diseases 93: 211–216. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ijid.2020.02.058&link_type=DOI) 6. [6].Worldometer from [https://www.worldometers.info/coronavirus/](https://www.worldometers.info/coronavirus/) 7. [7].Inaba, H. (2017). Age-Structured Populatin Dynamics in Demography and Epidemiology; Springer: Berlin/Heidelberg, Germany. 8. [8].Linton, N.M.; Kobayashi, T.; Yang, Y.; Hayashi, K.; Akhmetzhanov, A.R.; Jung, S.; Yuan, B.; Kinoshita, R.; Nishiura, H. (2020). Incubation period and other epidemiological characteristics of 2019 novel coronavirus infections with right truncation: A statistical analysis of publicly available case data. J. Clin. Med. 9: 538. 9. [9].Sun, H.; Qiu, Y.; Yan, H.; Huang, Y.; Zhu, Y.; Chen, S.X. (2020). Tracking and predicting COVID-19 epidemic in China mainland. medRxive. 10. [10].Diekmann, O.; Heesterbeek, J.A.P.; Metz, J.A.J. (1990). On the definition and the computation of the basic reproduction ratio R0 in models for infectious diseases in heterogeneous populations. J. Math. Biol. 28: 365–382. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/BF00178324&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=2117040&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F04%2F29%2F2020.04.24.20078287.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1990DG35000001&link_type=ISI) 11. [11].van den Driessche, P.; Watmough, J. (2002). Reproduction numbers and sub- threshold endemic equilibria for compartmental models of disease transmission. Math.Biosci.180:29–48. 12. [12].FitzHugh R. (1961). Impulses and physiological states in theoretical models of nerve membrane. Biophysical journal 1: 445–466. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0006-3495(61)86902-6&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19431309&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F04%2F29%2F2020.04.24.20078287.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A19617285A00001&link_type=ISI) 13. [13].Nagumo, J.; Animoto, S.; Yoshizawa, S. (1962). An active pulse transmission line simulating nerve axon. Proc. Inst. Radio Engineers 50: 2061–2070. 14. [14].Zhao, S.; Ran, J.; Musa, S.S.; Yang, G.; Lou, Y.; Gao, D.; Yang, L.; He, D. (2020). Preliminary Estimation of the Basic Reproduction Number of Novel Coronavirus (COVID-19) in China, from 2019 to 2020: A Data-Driven Analysis in the Early Phase of the Outbreak. Available online: [https://www.biorxiv.org/content/10.1101/2020.01.23.916395v1](https://www.biorxiv.org/content/10.1101/2020.01.23.916395v1) (accessed on 29 January 2020). 15. [15].Imperial College London-MRC Centre for Global Infectious Disease Analysis. News/Wuhan Coronavirus. Available online: [https://www.imperial.ac.uk/mrc-global-infectious-disease-analysis/news-wuhan-coronavirus/](https://www.imperial.ac.uk/mrc-global-infectious-disease-analysis/news-wuhan-coronavirus/) (accessed on 28 January 2020). 16. [16].The University of Hong Kong-HKUMed WHO Collaborating Centre for Infectious Disease Epidemiology and Control. Nowcasting and Forecasting the Wuhan COVID-19 Outbreak. Available online: [https://www/med.hku.hk/news/press/HKUMed-WHO-Collaborating-Centre-for-Infectious-Disease-Epidemiologyand-Control-releases%20real-time-nowcast-on-the-likely-extent-of-the-Wuhan-coronavirus](https://www/med.hku.hk/news/press/HKUMed-WHO-Collaborating-Centre-for-Infectious-Disease-Epidemiologyand-Control-releases%20real-time-nowcast-on-the-likely-extent-of-the-Wuhan-coronavirus) (accessed on 28 January 2020). 17. [17].Read, J.M.; Bridgen, J.R.; Cummings, D.A.; Ho, A.; Jewell, C.P. (2020).Novel Coronavirus COVID-19: Early Estimation of Epidemiological Parameters and Epidemic Predictions. MedRxiv. [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/inline-graphic-1.gif [8]: /embed/inline-graphic-2.gif [9]: /embed/graphic-7.gif [10]: /embed/graphic-8.gif [11]: /embed/graphic-9.gif [12]: /embed/graphic-10.gif [13]: /embed/graphic-11.gif [14]: /embed/graphic-12.gif [15]: /embed/graphic-13.gif [16]: /embed/inline-graphic-3.gif [17]: /embed/graphic-14.gif [18]: /embed/inline-graphic-4.gif [19]: /embed/inline-graphic-5.gif [20]: /embed/inline-graphic-6.gif [21]: /embed/inline-graphic-7.gif [22]: /embed/inline-graphic-8.gif [23]: /embed/graphic-15.gif [24]: /embed/inline-graphic-9.gif [25]: /embed/graphic-16.gif [26]: /embed/inline-graphic-10.gif [27]: /embed/inline-graphic-11.gif [28]: /embed/inline-graphic-12.gif [29]: /embed/inline-graphic-13.gif [30]: /embed/inline-graphic-14.gif [31]: /embed/inline-graphic-15.gif [32]: /embed/inline-graphic-16.gif [33]: /embed/graphic-19.gif [34]: /embed/inline-graphic-17.gif [35]: /embed/inline-graphic-18.gif [36]: /embed/inline-graphic-19.gif [37]: /embed/inline-graphic-20.gif