Mathematical Model with Social Distancing Parameter for Early Estimation of COVID-19 Spread =========================================================================================== * Saroj Kumar Chandra * Avaneesh Singh * Manish Kumar Bajpai ## Abstract COVID-19 is well known to everyone in the world. It has spread around the world. No vaccine or antiviral treatment is available till now. COVID-19 patients are increasing day by day. All countries have adopted social distancing as a preventive measure to reduce spread. It becomes necessary to estimate the number of peoples going to be affected with COVID-19 in advance so that necessary arrangements can be done. Mathematical models are used to provide early disease estimation based on limited parameters. In the present manuscript, a novel mathematical model with a social distancing parameter has been proposed to provide early COVID-19 spread estimation. The model has been validated with real data set. It has been observed that the proposed model is more accurate in spread estimation. Keywords * Corona Virus * COVID-19 * Mathematical Modelling * Epidemic ## 1. Introduction COVID-19 has emerged as a life-threatening outbreak disease. World Health Organization has declared as a pandemic in January 2020 [1]. The first case has been reported in Wuhan city of Hubei Province in south China on 31, December 2019 as unidentified pneumonia [2, 3]. COVID-19 disease has been identified as a member of the Severe Acute Respiratory Syndrome (SARS) that outbroke also in South China in 2002-2003 [4]. This disease is developed in the human body in the presence of a coronavirus. The coronavirus has been described in 1966 by Tyrell and Bynoe [5]. The virus is divided into four categories namely alpha, beta, gamma, and delta. Bats are a major source of alpha and beta virus category of the coronavirus. While gamma and delta originate from pigs and birds. Among these viruses, beta can infect human beings. The most common symptoms of COVID-19 are fever, tiredness, and dry cough. Some patients may have aches and pains, nasal congestion, runny nose, sore throat or diarrhea. Around 1 out of every 6 people who get COVID-19 becomes seriously ill and develops difficulty breathing [6]. COVID-19 is a transmissible disease. Hence, peoples get infected with coronavirus from others (COVID patients) by coughs or sneezing [7]. No vaccine or antiviral treatment is available till now. The mortality rate is increasing day by day. The rapid spread of the COVID-19 may be due to multiple causes. One cause is the lacking of information transparency at the early stage of the epidemic outbreak. Releasing the epidemic information in a timely and accurate way is extremely important for the anti-epidemic response of the public. The authentic and transparent information could have prohibited the spread of the COVID-19 at the early stage. The other cause is the lacking of the scientific diagnostic criterion for the COVID-19. As preventive measures, all the major countries have canceled events, classes in schools and colleges, businesses have pushed work from home policies. All of these measures are adopted to slow the spread of the disease. These measures are broadly referred to as social distancing. In the present manuscript, a mathematical model is being developed for estimating COVID-19 disease. Also, the proposed model incorporates the social distancing parameter for better spread estimation. Mathematical modelling has gained more attention and awareness in epidemiology and the medical sciences [8, 9, 10, 11]. It has been used for cancer detection, segmentation and classification [12, 13]. These models are useful in cases where disease dynamics are not unclear. It estimates the number of cases in worst and best-case scenarios. Susceptible-Infected-Removed (SIR) and Susceptible-Exposed-Infected-Removed (SEIR) models are the most studied and used models to estimate disease spread. These models are helpful in estimating the effect of preventive measures on disease spread. The SIR model originated from the study of the plague almost one hundred years ago. This model estimates spread by using contact rate and infection period only. It has been observed that this model is not suitable for estimating the spread in which the incubation period is involved such as COVID-19. Hence, the SEIR model is a suitable model for estimating COVID-19 spread since it incorporates the incubation period. However, the SEIR model failed to estimate the spread where preventive measures are adopted such as social distancing. In the present manuscript, a modified SEIR model is presented estimating COVID-19 spread. The proposed model has incorporated the social distancing parameter more accurate estimation of COVID-19 disease. The present work is organized as follows. The mathematical model for COVID-19 spread estimation is presented in section 2. Section 3 discusses the results obtained by the proposed mathematical model. ## 2. Proposed Methodology with Social Distancing In this section, a modified SEIR mathematical model is being presented for estimating COVID-19 spread. Proposed model uses basic SEIR model to design novel mathematical model. The SEIR model is a compartmental model for estimating disease spread in the population. Its an acronym for Susceptible, Exposed, Infected, Recovered. When a disease is introduced to a population, the people move from one compartment to another compartment. When a people has reached to R state, then the people either survived the disease or out of the population. The classical SEIR is defined be equations ![Formula][1] ![Formula][2] ![Formula][3] ![Formula][4] Here, equations (1), (2), (3) and (4) represent four Ordinary Differential Equations (ODE’s). The three parameters are *α*, *β* and *γ*. These are defined as * *α*: is the inverse of the incubation period * *β* is the average contact rate in the population * *γ* is the inverse of the mean infectious period Equation (1) is the change in people susceptible to the disease and is moderated by the number of infected people and their contact with the infected. Equation (2) gives the people who have been exposed to the disease. It grows based on the contact rate and decreases based on the incubation period whereby people then become infected. Equation (3) gives us the change in infected people based on the exposed population and the incubation period. It decreases based on the infectious period, so the higher *γ* is, the more quickly people die/recover and move on to the final stage in Equation (4). The proposed mathematical model is shown in the Figure 1 ![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/05/2020.04.30.20086611/F1.medium.gif) [Figure 1:](http://medrxiv.org/content/early/2020/05/05/2020.04.30.20086611/F1) Figure 1: Proposed Mathematical Model With Social Distancing In the proposed model a novel factor *ρ* has been introduced as social distancing parameter. Social distancing is maintained in transmissible diseases like COVID-19. It includes avoiding large gatherings, physical contact, and other efforts to mitigate the spread of infectious disease. This parameter controls contact rate, *β*. Social distancing parameter *ρ* lies in the range 01, where 0 indicates everyone is locked down and quarantined while 1 is equivalent to our base case above. The equations (1) and (2) are being multiplied with parameter to have social distancing effect. The modified equations are defined as ![Formula][5] ![Formula][6] ## 3. Results and Discussion All the experimental studies have been performed on PYTHON. The machine used for performing simulation work has CPU clock speed 1.60 GHz, 8GB RAM, 256 KB L1 cache, 1.0 MB L2 cache, and 6.0 MB L3 cache hardware configuration. The proposed mathematical model has been validated on cases of COVID-19 in India [14]. The day-wise COVID-19 cases are shown in the Figure 2 State-wise active COVID-19 cases in India is shown in Figure 2 [15]. It can be analyzed from Figure 2 that COVID-19 cases are increasing day by day. ![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/05/2020.04.30.20086611/F2.medium.gif) [Figure 2:](http://medrxiv.org/content/early/2020/05/05/2020.04.30.20086611/F2) Figure 2: Day-wise active COVID-19 cases in India The state-wise COVID-19 cases are Tabulated in the Table 1. It can be easily analyzed from the Table 1 that Maharashtra, Delhi, Tamil Nadu, Uttar Pradesh, Madhya Pradesh, Telengana and Rajasthan are major affected states with COVID-19 diseases. The COVID-19 spread behavior with different measurements such as before and after applying social distancing and spread due to breakage in maintaining social distancing have been shown in Figure 3. View this table: [Table 1:](http://medrxiv.org/content/early/2020/05/05/2020.04.30.20086611/T1) Table 1: Total number of confirmed state wise cases in India till 17-March-2020 ![Figure 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/05/2020.04.30.20086611/F3.medium.gif) [Figure 3:](http://medrxiv.org/content/early/2020/05/05/2020.04.30.20086611/F3) Figure 3: Active COVID-19 cases spread scenarios in India Figure 3 (a) shows polynomial spread behavior with degree 2 before applying social distancing. The effect of social distancing can be seen in Figure 3 (b) which shows a reduction of spread behavior from polynomial to linear. Exponential behavior has been analyzed by breaking the social distancing as can be seen in Figure 3 (c). The overall spread behavior up to 10 – *April –* 2020 is shown in Figure 3 (d). Figure 3 (d) shows how spread behavior has been changed from polynomial of degree 2 to polynomial of degree 3. The overall population has been considered as susceptible cases. The experimental studies have been performed on 10000 population. The cases which do not belong to susceptible cases either belong to exposed, infected or recovered. The peoples who have contacted recently to infected people move to exposed cases. The exposed cases move to infected after completion of the incubation period. The infected people’s moves recover cases after successful completion of coronavirus cycle, clinical process or death. The value of incubation period a has been fixed with 0.2 and an infection period of 0.5 days has been considered in all modes for validating proposed methodology [16, 17]. The performance of the proposed mathematical model works in two different modes. The first mode works on the classical SEIR model in which the value of the social distancing rate *ρ* is fixed to 1. In this mode, the spread of disease depends only on the contact rate *β*. The result obtained by proposed mathematical model with *ρ* =1 and *β* = 1.5, *β* = 1.75 and *β* = 2.0 has been shown in the Figure 4. Greenline shows COVID-19 spread with *β* = 1.5, similarly red and blue line shows spread with *β* = .1.75 and *β* = 2.0 respectively. It can be easily analyzed from the green dotted line that spread is about 3% in 21 day which is approximately the same as COVID-19 disease spread before applying social distancing in India. It can be also analyzed that if this spread has been followed then the maximum 35% population could have been suffered from COVID-19 disease in 30 days. ![Figure 4:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/05/2020.04.30.20086611/F4.medium.gif) [Figure 4:](http://medrxiv.org/content/early/2020/05/05/2020.04.30.20086611/F4) Figure 4: COVID-19 Spread before applying social distancing with *ρ* =1 The second mode is with different social distancing rates. The results obtained have been shown in Figure 5. It has been analyzed with about 300 active cases. It can be easily analyzed from the green dotted line that spread is about 10% in 10 days which is approximately the same as COVID-19 disease spread after applying social distancing in India. Hence, it can be also analyzed that spread behavior is linear in nature. Also, it can be said that if this spread has been followed then the maximum 30% population could have been suffered from COVID-19 disease 30 days. ![Figure 5:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/05/2020.04.30.20086611/F5.medium.gif) [Figure 5:](http://medrxiv.org/content/early/2020/05/05/2020.04.30.20086611/F5) Figure 5: COVID-19 Spread after applying different social distancing rates The sudden increase has been analyzed after 30 – *March –* 2020. This is due to the breakage of social distancing as a gathering of peoples. It can be easily observed from Figure 6 that spread is about to 60% peoples in 10 days which closely resembles COVID-19 spread after gathering. Also, It can be said that due to this gathering the spreading behavior has been changed from linear to a polynomial of degree 3. Also, it is found that 60 days lockdown is required in the current scenario to complete recovery of COVID-19 disease. ![Figure 6:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/05/2020.04.30.20086611/F6.medium.gif) [Figure 6:](http://medrxiv.org/content/early/2020/05/05/2020.04.30.20086611/F6) Figure 6: COVID-19 Spread due to breakage in social distancing ## 4. Conclusion The present work has investigated the problem of COVID-19 spread in India in current scenarios. A mathematical model has been established, which follows the actual data trend of COVID-19 spread in India. It has been proved from analytical (based on mathematical modeling) and simulation results that social distancing plays an important role in spread estimation. The effect of social distancing has been discussed with different social distancing rates. It has been found that social distancing can reduce spread from polynomial to linear. It has been also observed breakage in social distancing can rise spread from linear to exponentially. The maximum cases and recovery periods are also analyzed. It has been found that in the current spread scenario 60 days lockdown is required for complete recovery. ## Data Availability Code will be available on request ## Footnotes * *Email addresses:* saroj.chandra08{at}gmail.com (Saroj Kumar Chandra), avaneesh.singh{at}iiitdmj.ac.in (Avaneesh Singh), mkbajpai{at}iiitdmj.ac.in (Manish Kumar Bajpai) * Received April 30, 2020. * Revision received April 30, 2020. * Accepted May 5, 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].W. H. Organization, Novel Coronavirus (2019-nCoV) Advice for the Public, (2020), [https://www.who.int/emergencies/diseases/novel-coronavirus-2019/advice-for-public](https://www.who.int/emergencies/diseases/novel-coronavirus-2019/advice-for-public). 2. [2]. C. W. S. Hongzhou Lu1, Y. Tang, Outbreak of pneumonia of unknown etiology in Wuhan, China: The mystery and the miracle, Journal of Medical Virology (2020) 401–402. 3. [3]. L. Zhong, L. Mu, J. Li, J. Wang, Z. Yin, D. Liu, Early prediction of the 2019 novel coronavirus outbreak in the mainland china based on simple mathematical model, IEEE Access 8 (2020) 51761–51769. doi:10.1109/ACCESS.2020.2979599. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1109/ACCESS.2020.2979599&link_type=DOI) 4. [4].H. M. Zou L, Ruan F, Sars-cov-2 viral load in upper respiratory specimens of infected patients, New England Journal of Medicine 382 (2020) 11771179. 5. [5]. G. Williamson, Covid-19 epidemic editorial, The Open Nursing Journal 14 (2020) 37–38. doi:10.2174/1874434602014010037. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2174/1874434602014010037&link_type=DOI) 6. [6]. Z. A. Mizumoto K, Kagaya K, Estimating the asymptomatic proportion of coronavirus disease 2019 (covid-19) cases on board the diamond princess cruise ship, yokohama, japan, 2020. eurosurveillance, Epub ahead of print 25. 7. [7]. K. Kwok, F. Lai, W. Wei, S. Wong, J. Tang, Herd immunity estimating the level required to halt the covid-19 epidemics in affected countries, Journal of Infection doi:10.1016/j.jinf.2020.03.027. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jinf.2020.03.027&link_type=DOI) 8. [8]. C. L. Blackwood JC, An introduction to compartmental modeling for the budding infectious disease modeler, Letters in Biomathematics 5 (2018) 195–221. 9. [9]. A. Huppert, G. Katriel, Mathematical modelling and prediction in infectious disease epidemiology, Clinical Microbiology and Infection 19 (11) (2013) 999—1005. doi:[https://doi.org/10.1111/1469-0691.12308](https://doi.org/10.1111/1469-0691.12308). URL [http://www.sciencedirect.com/science/article/pii/S1198743X14630019](http://www.sciencedirect.com/science/article/pii/S1198743X14630019) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/1469-0691.12308&link_type=DOI) 10. [10]. D. P. Lizarralde-Bejarano, S. Arboleda-Snchez, M. E. Puerta-Yepes, Understanding epidemics from mathematical models: Details of the 2010 dengue epidemic in bello (antioquia, colombia), Applied Mathematical Modelling 43 (2017) 566—578. doi:[https://doi.org/10.1016/j.apm.2016.11.022](https://doi.org/10.1016/j.apm.2016.11.022). URL [http://www.sciencedirect.com/science/article/pii/S0307904X16306278](http://www.sciencedirect.com/science/article/pii/S0307904X16306278) 11. [11]. I. M. Longini, A mathematical model for predicting the geographic spread of new infectious agents, Mathematical Biosciences 90 (1) (1988) 367—383. doi:[https://doi.org/10.1016/0025-5564(88)90075-2](https://doi.org/10.1016/0025-5564(88)90075-2). URL [http://www.sciencedirect.com/science/article/pii/0025556488900752](http://www.sciencedirect.com/science/article/pii/0025556488900752) [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/0025-5564(88)90075-2&link_type=DOI) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1988P802900022&link_type=ISI) 12. [12]. S. K. Chandra, M. K. Bajpai, Mesh free alternate directional implicit method based three dimensional super-diffusive model for benign brain tumor segmentation, Computers Mathematics with Applications doi:[https://doi.org/10.1016/j.camwa.2019.02.009](https://doi.org/10.1016/j.camwa.2019.02.009). URL [http://www.sciencedirect.com/science/article/pii/S089812211930077X](http://www.sciencedirect.com/science/article/pii/S089812211930077X) 13. [13]. S. K. Chandra, M. K. Bajpai, Fractional mesh-free linear diffusion method for image enhancement and segmentation for automatic tumor classification, Biomedical Signal Processing and Control 58 (2020) 101841. doi:[https://doi.org/10.1016/j.bspc.2019.101841](https://doi.org/10.1016/j.bspc.2019.101841). URL [http://www.sciencedirect.com/science/article/pii/S1746809419304227](http://www.sciencedirect.com/science/article/pii/S1746809419304227) 14. [14].Worldometer, Coronavirus Update (Live): 1,521,090 Cases, 88,565 Deaths and 331,354 recovered cases from COVID-19 Virus Outbreak, [https://www.worldometers.info/coronavirus/](https://www.worldometers.info/coronavirus/). 15. [15].Ministry of Health and Family Welfare: Goverment of India, [https://www.mohfw.gov.in/](https://www.mohfw.gov.in/). 16. [16].J. H. et al., Feasibility of controlling covid-19 outbreaks by isolation of cases and contacts, The Lancet Global Health 8 (4) (2020) 488–496. 17. [17]. P. Liangrong, W. Yang, D. Zhang, C. Zhuge, L. Hong, Epidemic analysis of covid-19 in china by dynamical modeling doi:10.1101/2020.02.16.20023465. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMC4wMi4xNi4yMDAyMzQ2NXYxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDUvMDUvMjAyMC4wNC4zMC4yMDA4NjYxMS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) [1]: /embed/graphic-1.gif [2]: /embed/graphic-2.gif [3]: /embed/graphic-3.gif [4]: /embed/graphic-4.gif [5]: /embed/graphic-6.gif [6]: /embed/graphic-7.gif