Integration of Kalman filter in the epidemiological model: a robust approach to predict COVID-19 outbreak in Bangladesh ======================================================================================================================= * Md. Shariful Islam * Md. Enamul Hoque * Mohammad Ruhul Amin ## Abstract As one of the most densely populated countries in the world, Bangladesh have been trying to contain the impact of a pandemic like COVID-19 since March, 2020. Although government announced an array of restricted measures to slow down the diffusion in the beginning of the pandemic, the lockdown has been lifted gradually by reopening all the industries, markets and offices with a notable exception of educational institutes. As the physical geography of Bangladesh is highly variable across the largest delta, the population of different regions and their lifestyle also differ in the country. Thus, to get the real scenario of the current pandemic across Bangladesh, it is essential to analyze the transmission dynamics over the individual districts. In this article, we propose to integrate the Unscented Kalman Filter (UKF) with classic SIRD model to explain the epidemic evolution of individual districts in the country. We show that UKF-SIRD model results in a robust prediction of the transmission dynamics for 1-4 months. Then we apply the robust UKF-SIRD model over different regions in Bangladesh to estimates the course of the epidemic. Our analysis demonstrate that in addition to the densely populated areas, industrial areas and popular tourist spots are in the risk of higher COVID-19 transmission. In the light of these outcomes, we provide a set of suggestions to contain the pandemic in Bangladesh. All the data and relevant codebase is available at [https://mjonyh.github.io](https://mjonyh.github.io). **Highlights** * We integrate the UKF with classic SIRD model for the better estimation of the COVID-19 diffusion of 64 districts in Bangladesh. * Nationwide analysis show the strong correlation between population density and the number of COVID-19 positive cases in the country. * Industrial zones and popular tourists spots are at greater risk of spreading the Coronavirus. * With the better assessment of the COVID-19 cases dynamics, the Government will find effective policies to contain the current pandemic. ## 1 Introduction The world has been braving the COVID-19 crisis by which more than 52 million people across 113 countries over the globe got infected as of *November* 11, 2020. The Corona virus disease was first identified in Wuhan, the capital of Central China’s Hubei province, on December 2019. The World Health Organization (WHO) declared it as the public health emergency on January 20, 2020 and later as Pandemic on March 11, 2020 [1, 2]. Although China, Italy, Spain, USA and UK were severely affected by COVID-19 at first, beginning from April 2020 onward, countries from Global South are also impacted by this devastating Corona virus. As of *August* 30, 2020, the lower middle income country Bangladesh has the second lowest number of tests conducted per confirmed cased of COVID-19. The novel Corona virus can be transmitted via three commonly accepted modes, namely direct contact with an infected person, contaminated objects, and larger respiratory droplets or small airborne droplets [3]. Usually the virus infected person shows symptoms like fever, shortness in breath, cough, diarrhea, pneumonia and loss of smell, whereas the corona positive patient with premedical conditions such as diabetics, pneumonia and high blood pressure might experience the fatal death if proper medical treatment cannot be sought for in time [4, 5]. Since there are no recommended medicines or remedies available for treating COVID-19 disease as of now, many countries are setting up an array of measures including wearing face mask, early detection of infected people, contact tracing, isolation of infected people, increasing the number of testing per day, and imposing lockdown to curb the COVID-19 outbreak [6, 7]. Recent studies also support the effectiveness of wearing face mask, social distancing, and imposing lockdown on the declining of both infection and death rate [8]. Bangladesh is one of the most densely populated (ranked 8*th* worldwide) and under developing countries in the world. The recent economic development of this country makes it very competent in the South Asian region, especially in terms of foreign exchange reserves, reaching as high as second position in South Asia [9]. The first COVID-19 positive case was identified in Bangladesh on March 8, 2020. The Government of Bangladesh imposed a set of strict measures and closed all educational institutions, government, non-government offices including garments and factories [10] on March 26th, 2020. As a result, the number of infection and fatality rates are less compared to India and Pakistan. In Bangladesh, the contamination numbers spilled over 425,353 confirmed cases with death toll of more than 6,127 and it is still growing each day [11]. To understand the dynamics of a pandemic, mathematical model can play a major role. The classic mean-field Susceptible-Infected-Recovery-Death (SIRD) model by Kermack and McKendrick [10, 12], is frequently used to illustrate a quantitative picture of COVID-19 outbreak for many countries [6]. However, spread of an epidemic is a sophisticated process that depends largely on the mutated strains of the virus and its principle vectors [13, 14]. Moreover, the recorded data of infected, recovered and death cases might miss significant amount of information due to many kinds of biases [15]. For instance, time series of recovered numbers are heavily unreliable as it is solely based on how a country traces its asymptomatic or mildly symptomatic patients. Therefore, usual mean-field, compartment-based models and stochastic spatial epidemic models might estimate the dynamics of an epidemic with higher error margin [16–19]. The SIRD model is associated with a set of differential equations which compute values for an instantaneous event, whereas the time series of COVID-19 cases that we observe everyday is a discrete event. Therefore, all the daily events of COVID-19 cases lack the instantaneous effect in the differential equation based SIRD model. Here, we use the prediction based Unscented Kalman Filter (UKF) model to derive the dynamics. UKF is a classic non-linear estimation algorithm that accurately and timely forecasts the dynamic state in a nonlinear system. UKF has been used in various areas such as, navigation, target tracking, structural dynamics and vehicle positioning due to its high accuracy and rapid convergence merits [20, 21]. In this study, we integrate the UKF with classic SIRD model to examine the underlying process of COVID-19 transmission in a more precise way. As Bangladesh has 64 districts which are uneven in terms of population density, economical significance and geographical locations, we estimate the transmission trend of 64 districts to evaluate the complete picture of the outbreak in Bangladesh. ## 2 Method We used the standard SIRD model to estimate the number of COVID-19 active cases, recovered cases and death cases in Bangladesh. A homogeneous immunity has been considered for the whole country with no transmission from animals and no significant difference between natural death and birth. We segmented the total population size, *N*, into four stages of disease: *S*, susceptible; *I*, infected; *R*, recovered; *D*, death and can be written by *N* = *S* + *I* + *R* + *D* [10, 12, 22, 23]. The interaction between these four stages is illustrated in the Figure 1. ![Fig 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/12/23/2020.10.14.20212878/F1.medium.gif) [Fig 1.](http://medrxiv.org/content/early/2020/12/23/2020.10.14.20212878/F1) Fig 1. Schematic diagram of the compartmental SIRD model. It represents the total populations (*N*) into four different health status, such as, susceptible (*S*), infected (*I*), recovered (*R*) and death (*D*). The *β, γ* and *µ* denote the rate of infection, rate of recovery and rate of death, respectively. The differential equations of SIRD model can be written as: ![Formula][1] ![Formula][2] ![Formula][3] ![Formula][4] Where *t* denotes the time duration, *S*(*t*) represents the number of susceptible people at time *t, I*(*t*) shows the number of people infected at time *t, R*(*t*) stands for the number of people who have recovered from the infection and *D*(*t*) indicates the number of deaths. The population, *N* = *S*(*t*) + *I*(*t*) + *R*(*t*) + *D*(*t*), is a conserved quantity for every time step [24, 25]. The constant *β, γ* and *µ* illustrate the rate of infection, recovery and deaths, respectively. The Eq. 1 - 4 have been discretized by forward Euler method as, ![Formula][5] ![Formula][6] ![Formula][7] ![Formula][8] An augmented state vector *X* is presented for simplicity as, ![Formula][9] and the discrete-time augmented SIRD model, Eq. 5 - 8, can be written as, ![Formula][10] where *F* is the nonlinear term and *W* is the zero mean Gaussian uncertainty with covariance *Q**f*. All types of reported cases, the cumulative confirmed, recovered, deaths and active cases, are incorporated with the model using the vector as, ![Formula][11] where *ν* is the uncertainties due to the SERS-CoV-2 test results. The considered zero mean Gaussian uncertainty has the known covariance *R**f*. Let us consider an estimated vector state ![Graphic][12] for the Unscented Kalman Filter (UKF), then, ![Formula][13] where ![Graphic][14] is the Jacobian matrix of *F*(*X**t*), given by, ![Formula][15] Hence the algorithm of UKF is given by, ### Predict ![Formula][16] ![Formula][17] ### Update ![Formula][18] ![Formula][19] ![Formula][20] ![Formula][21] Here, *P* (*t* |*t*) and *C* represent the posterior estimate covariance matrix and the data augmented matrix respectively. We integrated this updated formula of UKF in SIRD model to estimate the dynamics of the SERS-CoV-2 in Bangladesh. The dynamical behaviour of the class of infected people is described by the basic (*R*), effective (*R**e*) and time-varying (*R**t*) reproduction number. In the case of SIRD model, this are defined as, ![Formula][22] ![Formula][23] Using Eq. 2, one can find that *S*(*t*) = *N/R* and *R**e* = 1, at *I* = *I**max*. Besides, one can estimate the onset confirmed (*C**′*) cases as, ![Graphic][24], where *θ* = *γ**′*(*R**t* *−* 1) observes random walk and *γ**′* vary independently [26]. Thus, ![Formula][25] If this reproduction number becomes larger than the ratio between the total population and initial susceptible people, there will be a proper epidemic outbreak [27]. Otherwise, the disease will never lead to a proper outbreak. ## 3 Results As a state of art, we first evaluate whether UKF-SIRD model provides more robust estimation of COVID-19 active cases in Bangladesh compared to classic SIRD model prediction. To illustrate the active cases of the outbreak, the UKF-SIRD and SIRD model simulations were performed in the GNU octave (Version 5.2.0) and Python programming language environment. It is clearly noticeable in figure 2 that UKF-SIRD model delivers more accurate data estimation compared to normal SIRD model prediction. For instance, the UKF-SIRD model predicts around 85,000 corona active cases at 100 days (third week of June, 2020) relative to 90,000 real cases, whereas SIRD model estimates 120,000 active cases at the same time in Bangladesh. Therefore, these results illustrate the fact that integrated UKF in SIRD model is a robust method which is capable of exclusive prediction of COVID-19 transmission in Bangladesh. ![Fig 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/12/23/2020.10.14.20212878/F2.medium.gif) [Fig 2.](http://medrxiv.org/content/early/2020/12/23/2020.10.14.20212878/F2) Fig 2. The estimation of corona active cases in Bangladesh by UKF-SIRD and SIRD model. Dynamics of the active cases population has been shown as a function of time. The real data are fitted from March 8 to July 28, 2020. Bangladesh has 64 districts that are heterogeneous in the context of different population densities, socioeconomic activities and geographical parameters as shown in figure 3. For example, Dhaka, the capital of Bangladesh, has a population density of 23,234 people/sq. km compared to Bandarban (87 people/sq. km) district in Bangladesh. On the other hand, Narayanganj has more industries than Bogra district. Recently, Nadjat and colleagues showed the positive impact of population density on COVID-19 diffusion in Algeria [32]. Thus, considering the heterogeneous density in Bangladesh, general estimation of the COVID-19 dynamics will never match with the real outbreak scenario unless it is done on the district level. To evaluate the local transmission in Bangladesh, next we cluster whole Bangladesh into three groups according to its population density, industrial area and geographical variance. ![Fig 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/12/23/2020.10.14.20212878/F3.medium.gif) [Fig 3.](http://medrxiv.org/content/early/2020/12/23/2020.10.14.20212878/F3) Fig 3. Classification of Bangladesh map based on population density, industrial regions, agricultural land, high land, Beel and Haor ![Fig 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/12/23/2020.10.14.20212878/F4.medium.gif) [Fig 4.](http://medrxiv.org/content/early/2020/12/23/2020.10.14.20212878/F4) Fig 4. Epidemic evolution of COVID-19 based on confirmed cases, daily growth and *R**t* values in Bangladesh. The UKF-SIRD model illustrates the positive impact of population densities on COVID-19 confirmed cases and daily growth in major districts such as, Dhaka, Gazipur and Narayanganj. ![Fig 5.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/12/23/2020.10.14.20212878/F5.medium.gif) [Fig 5.](http://medrxiv.org/content/early/2020/12/23/2020.10.14.20212878/F5) Fig 5. Epidemic evolution of COVID-19 based on confirmed cases, daily growth and *R**t* values in Bangladesh. The UKF-SIRD model illustrates the positive impact of industrial zones on COVID-19 confirmed cases and daily growth.High Industrial districts such as Dhaka, Chattogram are at high risk for COVID-19 diffusion compared to low industrial area like, Bogra. ![Fig 6.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/12/23/2020.10.14.20212878/F6.medium.gif) [Fig 6.](http://medrxiv.org/content/early/2020/12/23/2020.10.14.20212878/F6) Fig 6. Epidemic evolution of COVID-19 based on confirmed cases, daily growth and *R**t* values in Bangladesh. The UKF-SIRD model illustrates that Beel, Haor areas (Sunamganj, Satkhira), Highland districts (Rangamati, Khagrachhari) and argicultural zones (Noakhali, Dinajpur) are at low risk for COVID-19 transmission in Bangladesh. ### 3.1 Estimation of COVID-19 diffusion in highly dense cities of Bangladesh In this article, we consider 12 districts (high to low dense areas) of Bangladesh to examine the spread of the Corona virus infection. Fig. 3.1 clearly shows that densely populated area, Dhaka is severely affected by COVID-19 outbreak as the confirmed cases reach ∼ 100, 000, whereas low dense area, Bandarban has only ∼ 1, 000 positive cases. In addition, the daily growth rate also follows a high trajectory trend for Dhaka compared to Narsingdi, Shariatpur and Bagerhat, as shown in the Fig. 3.1. The *R**t* values of Dhaka were over 1.0 starting from April, 2020 to end of July,2020, whereas Bagerhat had the *R**t* values *>* 1.0 for only 30 days (June,2020). Interestingly, we notice that the district of Cox’s Bazar has similar number of corona positive cases ∼ 10, 000 like highly dense district of Gazipur despite the fact that Cox’s Bazar is a moderately dense area in Bangladesh. Since Cox’s Bazar is the most popular tourists spot in Bangladesh, we hypothesize that high mobility in this area might cause the high number of corona positive cases. Strikingly, at the end of July, 2020 *R**t* value is approaching to *<* 1.0 for all twelve districts in Fig. 3.1, which could lead to the end of this pandemic if the current trend continues for the next 3-4 months. ### 3.2 COVID-19 transmission trend in the industrial areas of Bangladesh The year 2020 has seen the unprecedented downfall of the economy for both the national and global market, trade, and investment to reduce the infection rate of the novel Corona virus by ceasing the means of transportation and forcing the lockdown measures that shutdown many business organizations. Moreover, many garment factories were temporarily closed to avoid the risk of Corona virus outbreak. Despite the government-announced shutdown to prevent the COVID-19 diffusion in Bangladesh, industrialists reopened the factories for the economic salvation. However, health experts have shown their grave concerns about the proper guidelines for social distancing inside the local factories, as several ready made garment (RMG) workers have been tested Corona virus positive in the recent days. Therefore, it is important to determine the COVID-19 transmission trend in these industrial districts such as, Narayanganj, Chattogram, Khulna and Rajshahi. Our systematic UKF-SIRD model analysis illustrates a sharp daily growth of corona positive cases for the industrial area, Chattogram, compared to another industrial zone, Narayanganj. Moreover, small industrial region like Rajshahi has less number of cases (*∼* 1, 000) relative to Narayanganj (*∼* 10, 000) and Dhaka (*∼* 100, 000). Furthermore, we observe the uncertain values of *R**t* for Gazipur district, whereas Khulna show a growing trend of *R**t* from August, 2020. These analysis clearly forecast the fact that high industrial areas are at greater risk of Corona virus infection in Bangladesh. ### 3.3 COVID-19 outbreak in the crop land, hilly areas and low land of Bangladesh Bangladesh is known as a predominant agricultural country where more than 50% of the population rely directly or indirectly on employment in the agricultural and fisheries sectors, as indicated by the Food and Agricultural Organization of United Nation, 2020. These sectors are playing significant role in national economy by contributing 18% to country’s GDP as mentioned in the Bangladesh Economic Review. Thus, next we estimate the COVID-19 cases in several districts of Bangladesh that are highly associated with agricultural and fisheries sectors. The UKF-SIRD model estimates *∼* 1, 500 confirmed cases at the end of October for Noakhali district compared to ∼ 1, 000 cases for Patuakhali and Dinajpur region. In addition, haors (large shallow water body or backswamp) and wetlands area, such as Sunamganj has approximately ∼ 1, 200 cases, whereas Satkhira has ∼ 700 Corona virus positive cases in October, 2020. Notably, the model predicts lower number of confirmed cases (*<* 1, 000) and a smaller growth rate for high lands areas namely, Rangamati and Khagrachhari districts. Since these districts also have low *R**t* values (*<* 0.9) at the moment, there is a high chance that Rangamati and Khagrachhari will wipe out the COVID-19 outbreak if the current transmission trend continues for next 3-4 months. Taken together, this study represents that the COVID-19 pandemic is less serious in agricultural and high land areas compared to industrial areas in Bangladesh. ## Discussion Epidemiological models are the powerful tools that are designed for the mathematical representations of the disease transmission and their associated dynamics [28, 29]. In addition, models can be used to estimate the trend of an epidemic outbreak under certain conditions so that necessary actions can be taken at various scales to restrain the disease diffusion in a community [30]. However, to make a series of policy decisions, it is a critical prerequisite to have a robust epidemiological model that is suitable for the intended purpose and shows accurate and precise prediction of a pandemic [31]. In this study, we integrated UKF with classic SIRD model and estimate COVID-19 diffusion in different parts of Bangladesh. Dhaka is the 6*t**h* largest densely populated city in the World, and the key financial hub which holds the major portions of the Bangladesh’s economy [32, 33]. Unfortunately, we observe that Dhaka is the epicenter of the COVID-19 pandemic with the highest number of infected cases due to its high population density in Bangladesh. Although the government has implemented zone-coded lockdown in small areas within the Dhaka city [34], we observe that the mobility between zones could not be maintained strictly due to economic activities. If the government had the capacity to test and isolate each positive cases, only then the zone-coded lockdown system in Bangladesh could have been successful. Several reports show the association of population density with higher rates of transmission, infection, and mortality from COVID-19 [35, 36]. Strikingly, our nationwide analysis of COVID-19 dynamics also demonstrate the relationship between population density and the COVID-19 diffusion rates for Bangladesh. Industrial zones are at greater risk of transmitting the Corona virus as large groups of workers enter and leave the industrial area without maintaining social distance and proper hygiene in most of the cases. As a result, Narayanganj, Dhaka and Chattogram have been braving with COVID-19 for last couple of months. Notably, our model speculates the dire consequences for Chattogram if the COVID-19 diffusion continues with the current trend. Therefore, industrialists from Chattogram must ensure proper guideline for social distancing and the hygiene rules at the industrial areas to save their workers. Moreover, to reduce the public mobility, we suggest that public transportation between high risk areas should also be suspended for a time until the situation gets under control. To support this claim we would like to mention that the high rate of COVID-19 positive cases in the New York city of USA could have been controlled by implementing the proper lockdown and public mobility suspension [37]. According to our analysis, districts with the hilly areas and the wet land areas in Bangladesh show lower rate of infection than elsewhere. Besides the low population density in areas of Bandarban, Khagracchari and Rangamati, the government also control the transportation facilities in the area with some restrictions. On the other hand, the Cox’s Bazar, a popular tourist destination in Bangladesh has seen rise in COVID-19 positive cases despite having less population density. So it is evident that lower population density and restricted movements are the keys to slow down the diffusion rate in any areas. Therefore, we suggest to suspend or restrict the mass transportation between densely populated area, and industrial zone in a bit to contain the pandemic in Bangladesh. Since Bangladesh is in emergency situation, we deal with incomplete and raw data of the current outbreak. For instance, if a person from a family becomes COVID-19 positive, the whole family might turn out to be positive. But, in many families only one of the members were tested for COVID-19 and if positive then counted as one despite the necessity to test for the entire family. Moreover, there were many cases who died with COVID-19 symptoms, but could not be tested, and conversely all the death cases counted in the analysis due to COVID-19 might have died by other diseases. In addition, the recovered cases count only considered those who went through two negative COVID-19 test results. Therefore, our UKF-SIRD model might not exemplify the pandemic situation accurately. ## Conclusion Taken together, this study proposes a robust epidemiological method to assess the dynamics of COVID-19 pandemic in Bangladesh, as well as in any other countries which are still in the middle/late phase of the outbreak. Moreover, our estimation provides policymakers with a tool to evaluate the consequences of possible arrangements, for instance district wise lockdown, COVID-19 testing, social distancing and contact tracing. This data-driven mathematical modeling suggests that enforcing restricted social-distancing measures are obligatory and it needs to be maintained to contain the diffusion until the situation gets back to normal. We are optimistic that our UKF-SIRD model can serve as a base model for the policy suggestions and controlling the disease transmission in a heterogeneous region by adding more sophisticated factors under different environmental and political conditions in Bangladesh. ## Data Availability All the data and relevant codebase is available at [https://mjonyh.github.io](https://mjonyh.github.io). [https://mjonyh.github.io](https://mjonyh.github.io). ## Supporting information ## Declarations ### Ethics Approval and Consent to Participate Not applicable. ### Consent for Publication Not applicable. ### Competing interests The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest. ### Funding The authors declare that the research was conducted in the absence of any commercial or financial relationships. ### Author’s contributions MH, MI, MA made substantial contributions to the conception or design of the work. MH and MI completed the relevant studies and created plots related to the SIRD estimation using Kalman filter, while MA critically revised the manuscript. All the authors met regularly to discuss the outcome of each experiment. The data were acquired from the IEDCR website and each of the results were cross validated by at least two of the three authors. All authors are responsible for acquiring, analysing, and interpreting the data for this article. All the authors approved the version to be published. ## Acknowledgements We acknowledge the contribution of [www.Pipilika.com](http://www.Pipilika.com) software development team for collecting the daily data of total COVID-19 tests, total positive cases, and total deaths in Bangladesh. ## Footnotes * 1. We have added highlights of the research 2. Updated discussion section 3. Added conclusion section 4. Added supporting information, such as a. abbreviations b. author contributions c. acknowledgements 5. Updated references ## Abbreviations CFR : Case Fatality Rate COVID-19 : Coronavirus Disease 2019 IEDCR : Institute of Epidemiology, Disease Control and Research in Bangladesh SIRD : Susceptible, Infected, Recovered and Death UKF : Unscented Kalman Filter WHO : World Heath Organization * Received October 14, 2020. * Revision received December 23, 2020. * Accepted December 23, 2020. * © 2020, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/) ## References 1. 1.Organization WH, et al. Laboratory testing for coronavirus disease 2019 (COVID-19) in suspected human cases: interim guidance, 2 March 2020. World Health Organization; 2020. 2. 2.Organization WH, et al. Statement on the second meeting of the International Health Regulations (2005) Emergency Committee regarding the outbreak of novel coronavirus (2019-nCoV); 2020. 3. 3.Organization WH, et al. Transmission of SARS-CoV-2: implications for infection prevention precautions: scientific brief, 09 July 2020. World Health Organization; 2020. 4. 4.Huang C, Wang Y, Li X, Ren L, Zhao J, Hu Y, et al. Clinical features of patients infected with 2019 novel coronavirus in Wuhan, China. The lancet. 2020;395(10223):497–506. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0146736(20)30183-5&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F23%2F2020.10.14.20212878.atom) 5. 5.Figueiredo A, Daponte Codina A, Figueiredo M, Saez M, Cabrera León A. Impact of lockdown on COVID-19 incidence and mortality in China: an interrupted time series study. Bull World Health Organ. 2020;. 6. 6.Chatterjee S, Sarkar A, Chatterjee S, Karmakar M, Paul R. Studying the progress of COVID-19 outbreak in India using SIRD model. medRxiv. 2020;. 7. 7.Lau H, Khosrawipour V, Kocbach P, Mikolajczyk A, Schubert J, Bania J, et al. The positive impact of lockdown in Wuhan on containing the COVID-19 outbreak in China. Journal of travel medicine. 2020;27(3):taaa037. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/jtm/taaa037&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32523080&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F23%2F2020.10.14.20212878.atom) 8. 8.Start TD. Dealing with foreign exchange reserves;. 9. 9.Worldometer CC. Worldometer.(2020) 1–22;. 10. 10.Kermack WO, McKendrick AG. A contribution to the mathematical theory of epidemics. Proceedings of the royal society of london Series A, Containing papers of a mathematical and physical character. 1927;115(772):700–721. 11. 11.Standard TB. Educational institutions to remain closed till March 31;. 12. 12.Kermack WO, McKendrick AG. Contributions to the mathematical theory of epidemics. II.—The problem of endemicity. Proceedings of the Royal Society of London Series A, containing papers of a mathematical and physical character. 1932;138(834):55–83. 13. 13.Wang Y, Wang Y, Chen Y, Qin Q. Unique epidemiological and clinical features of the emerging 2019 novel coronavirus pneumonia (COVID-19) implicate special control measures. Journal of medical virology. 2020;92(6):568–576. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F23%2F2020.10.14.20212878.atom) 14. 14.Kupferschmidt K. Why do some COVID-19 patients infect many others, whereas most don’t spread the virus at all. Science. 2020;10. 15. 15.Carletti T, Fanelli D, Piazza F. COVID-19: The unreasonable effectiveness of simple models. arXiv preprint arxiv:200511085. 2020;. 16. 16.Gatto M, Bertuzzo E, Mari L, Miccoli S, Carraro L, Casagrandi R, et al. Spread and dynamics of the COVID-19 epidemic in Italy: Effects of emergency containment measures. Proceedings of the National Academy of Sciences. 2020;117(19):10484–10491. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTE3LzE5LzEwNDg0IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMTIvMjMvMjAyMC4xMC4xNC4yMDIxMjg3OC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 17. 17.Giordano G, Blanchini F, Bruno R, Colaneri P, Di Filippo A, Di Matteo A, et al. Modelling the COVID-19 epidemic and implementation of population-wide interventions in Italy. Nature Medicine. 2020; p. 1–6. 18. 18.Chinazzi M, Davis JT, Ajelli M, Gioannini C, Litvinova M, Merler S, et al. The effect of travel restrictions on the spread of the 2019 novel coronavirus (COVID-19) outbreak. Science. 2020;368(6489):395–400. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIzNjgvNjQ4OS8zOTUiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8xMi8yMy8yMDIwLjEwLjE0LjIwMjEyODc4LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 19. 19.Balcan D, Colizza V, Gonçalves B, Hu H, Ramasco JJ, Vespignani A. Multiscale mobility networks and the spatial spreading of infectious diseases. Proceedings of the National Academy of Sciences. 2009;106(51):21484–21489. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTA2LzUxLzIxNDg0IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMTIvMjMvMjAyMC4xMC4xNC4yMDIxMjg3OC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 20. 20.Julier SJ, Uhlmann JK. Unscented filtering and nonlinear estimation. Proceedings of the IEEE. 2004;92(3):401–422. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1109/JPROC.2003.823141&link_type=DOI) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000220005300002&link_type=ISI) 21. 21.Zheng B, Fu P, Li B, Yuan X. A robust adaptive unscented Kalman filter for nonlinear estimation with uncertain noise covariance. Sensors. 2018;18(3):808. 22. 22.Brauer F, Castillo-Chavez C, Feng Z. Challenges, Opportunities and Theoretical Epidemiology. In: Mathematical Models in Epidemiology. Springer; 2019. p. 507–531. 23. 23.Daley DJ, Gani J. Epidemic modelling: an introduction. vol. 15. Cambridge University Press; 2001. 24. 24.Anderson RM, May RM. Population biology of infectious diseases: Part I. Nature. 1979;280(5721):361–367. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/280361a0&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=460412&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F23%2F2020.10.14.20212878.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1979HE95900033&link_type=ISI) 25. 25.May RM, Anderson RM. Population biology of infectious diseases: Part II. Nature. 1979;280(5722):455–461. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/280455a0&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=460424&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F23%2F2020.10.14.20212878.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1979HF85100031&link_type=ISI) 26. 26.k sys. covid-19: Estimating COVID-19’s Rt in Real-Time with PYMC3; 2020. Available from: [https://github.com/k-sys/covid-19](https://github.com/k-sys/covid-19). 27. 27.Hethcote HW. The mathematics of infectious diseases. SIAM review. 2000;42(4):599–653. 28. 28.Perez AM, Ward MP, Charmandarián A, Ritacco V. Simulation model of within-herd transmission of bovine tuberculosis in Argentine dairy herds. Preventive veterinary medicine. 2002;54(4):361–372. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0167-5877(02)00043-0&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12163252&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F23%2F2020.10.14.20212878.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000177663100006&link_type=ISI) 29. 29.Smith G, Grenfell B, et al. Population biology of pseudorabies in swine. American journal of veterinary research. 1990;51(1):148–155. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=2154146&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F23%2F2020.10.14.20212878.atom) 30. 30.Hoque E, Islam MS, Amin MR, Das SK, Mitra DK. Adjusted Dynamics of COVID-19 Pandemic due to Herd Immunity in Bangladesh. medRxiv. 2020;. 31. 31.Sargent RG. Verification and validation of simulation models. In: Proceedings of the 2010 winter simulation conference. IEEE; 2010. p. 166–183. 32. 32.Data NS, (SEDAC) AC. Population Statistics;. 33. 33.Start TD. Dhaka’s economic activities unplanned: analysts;. 34. 34.Standard TB. Enforcing zone-based lockdown up in the air;. 35. 35.Post TW. The United States might have a secret weapon against coronavirus;. 36. 36.Times TNY. Density Is New York City’s Big ‘Enemy’ in the Coronavirus Fight;. 37. 37.Glaeser EL, Gorback CS, Redding SJ. How much does covid-19 increase with mobility? evidence from new york and four other us cities. National Bureau of Economic Research; 2020. [1]: /embed/graphic-2.gif [2]: /embed/graphic-3.gif [3]: /embed/graphic-4.gif [4]: /embed/graphic-5.gif [5]: /embed/graphic-6.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/inline-graphic-1.gif [13]: /embed/graphic-13.gif [14]: /embed/inline-graphic-2.gif [15]: /embed/graphic-14.gif [16]: /embed/graphic-15.gif [17]: /embed/graphic-16.gif [18]: /embed/graphic-17.gif [19]: /embed/graphic-18.gif [20]: /embed/graphic-19.gif [21]: /embed/graphic-20.gif [22]: /embed/graphic-21.gif [23]: /embed/graphic-22.gif [24]: /embed/inline-graphic-3.gif [25]: /embed/graphic-23.gif