key: cord-258102-7q854ppl authors: Mandal, S.; Kumar, M.; Sarkar, D. title: LOCKDOWN AS A PANDEMIC MITIGATING POLICY INTERVENTION IN INDIA date: 2020-06-20 journal: nan DOI: 10.1101/2020.06.19.20134437 sha: doc_id: 258102 cord_uid: 7q854ppl Abstract. We use publicly available timeline data on the Covid-19 outbreak for nine indian states to calculate the important quantifier of the outbreak, the sought after Rt or the time varying reproduction number of the outbreak. This quantity can be measured in in several ways, e.g. by application of Stochastic compartmentalised SIR (DCM) model, Poissonian likelihood based (ML) model & the exponential growth rate (EGR) model. The third one is known as the effective reproduction number of an outbreak. Here we use, mostly, the second one. It is known as the instantaneous reproduction number for an outbreak. This number can faithfully tell us the success of lockdown measures inside indian states, as containment policy for the spread of Covid-19 viral disease. This can also, indirectly yield notional value of the generation time inteval in different states. In doing this work we employ, pan India serial interval of the outbreak estimated directly from data from January 30th to April 19th, 2020. Simultaneously, in conjunction with the serial interval data, our result is derived from incidences data between March 14th, 2020 to June 1st, 2020, for the said states. We find the lockdown had marked positive effect on the nature of time dependent reproduction number in most of the Indian states, barring a couple. The possible reason for such failures have been investigated. Global pandemic outbreaks are very common nowadays. India is no exception. The Severe Acute Respiratory Syndrome (2003), [1] Avian Inunza (2006) , [2] Swine Flu (2015) [3] are to name a few. There were others that did not touch upon India but were recent events, such as MERS (2012) [4] & EVD (Ebola virus disease) (2000, 2003, 2004 & most recently in 2012) [5] . None of the above, however, touched the global pandemic scale, of what has been attained by Novel Coronavirus, aka Covid-19, in a short span of time, starting at the end of past year [6, 7, 8] . The global community responded to this unprecedented situation by various policy interventions. Wearing masks & face shields [9, 10] in public, social distancing norms [11] were amongst them. More drastic & perhaps draconian step of lockdown [12] was taken by governments across the world, as a containment policy measure [13] . We analyse the eect of lockdown on the propagation of Covid-19 viral disease. The instantaneous version of basic reproduction number [14] of the infection is plotted against time to gauge the success [15] (or lack thereof) [16] of this policy intervention in nine dierent states of India. In the following, it is shown that this pervasive containment policy has borne fruit in most of the considered provinces. Time dependent or instantaneous reproduction number [17] is an accurate projection is an in-situ description of virility or virulence of epidemic diseases. The basic reproduction number [14] gives us the average number of infectee cases per infector from the previous generation, over a given period of time, in a fully susceptible population. Various policy implimentation and containment measures appreciably reduce the number of contacts, in turn reducing, the eective number [18] of susceptible contacts per potential infector. Epidemiologist have devised a time dependent parameter, eective reproduction number to assimilate the eect of policy intervention into the basic reproduction number during an ongoing epidemic. This quantity is dened as follows: Consider an individual, who turns infectious on day t. We denote by R e (t) the expected number of secondary cases this infectious individual causes, in future [19] . The instantaneous reproduction number, on the other hand, compares the number of new infections on day t with the infection pressure (force of infection) [20] from the days prior to t. It can be interpreted as the average number of secondary cases that each symptomatic individual at time t would infect, if the conditions remained as they were at time t. Hence, the stepwise, undulations, crests, troughs & spikes of this estimate is termed as instantaneous or real time measures. There are various ways to calculate this eective instantanous & other time varying reproduction number. They are such to be: 2.1. Stochastic dynamic contact model-based method [21] . A stochastic Susceptible-Infected-Removed (SIR) model is considered in this case, in place of a deterministic one. Stochastic dynamic model has advantages over the standard deterministic one, in that, it accomodates improved variabilities and allows for better quantication of uncertainies of that number as compared to the standard deterministic model. Here, S(t), I(t) & R(t) denote the number of susceptible, infectious and recovered population at time respectively, and that N = S(t) + I(t) + R(t) is the total population. The infectious period of an infected individual is a random variable T ∼ exp (γ) & the reproduction rate is R(t) ≈ βE(t) = β γ , where β & γ are the transmission rate and recovery rate. The mathematical essence of the model can captured by a set of four coupled rst order linear homogenous dierential equations given such to be Here, ∼ signies deterministic (average) counterparts. We set S (0) equals the population of the region, R(0) = 0, I (0) is 10 to 14 times the average number of conrmed cases from Day 0 to Day 7, and γ the inverse of mean infectious period, obtained from the parametrization of serial interval distribution collected directly from data described in section (3) . The main diculty with this time varying reproduction number is that it assumes a constant transmissibility, where it may vary & often peak, during the generation time interval and just before the onset of symptoms [20] . This model also can not accomodate various disease traits like asymptomaticity (non-detection), or human interferences like isolation measures or migration etc. Hence we do not look at this method any further here. [21] . Here it is assumed that the total number of secondary infectees that were infected by a single praimary infector follows a Poisson distribution. The number of individuals infected on (discrete) Date t is usually replaced by the number of daily incidences reported, on the same Date t. Also, the generation time interval is suitably replaced by the corresponding serial interval interval for all practical purposes. Let N t be the number of reported incidences on Day t. Assuming that the serial interval has a maximum of k days and the number of new cases generated by an infected individual is assumed to follow a Poisson distribution with parameter [22] R. The probability that the serial interval of an individual lies in j days is w j , which can be estimated from the empirical distribution of serial interval or by setting up a discretized Gamma prior on it. Note only the nonnegative values of serial interval are used here. Thus, the likelihood function can be reduced into a thinned Poisson distribution as such All rights reserved. No reuse allowed without permission. (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. The copyright holder for this preprint this version posted June 20, 2020. . https://doi.org/10.1101/2020.06.19.20134437 doi: medRxiv preprint The instantaneous reproduction number can then be estimated by maximising the likelihood function as follows 3. Exponential growth rate-based (EGR) method [21] . In the early days of the epidemic the number of infected cases rise exponentially. The growth rate (Malthusian coecient) r can be estimated by tting a non linear least square tting into the daily incidence curve. The probability density function of serial interval of the outbreak is denoted by f λ (t), then the eective reproduction number is given by the Euler Lotka (type) Equation in case we have a non parametric serial interval distribution then we can dene our eective reproduction number as where λ i are the observed serial intervals. Using the publically available data on github [23] , to create a contact list between infector infectee pairs in the pan Indian context, between 30 th January-2020 to 19 th April-2020. The data is then tted with a log normal / gamma distribution to parametrize the values of mean and standard deviation. From the available data on github [23] , the daily conrmed case incidences were collected for nine states, for the duration of 14 th March-2020 to 1 st June-2020. Applying the Poissonian ML method, the instantaneous R i (t)was plotted for each one of them, as given in gure number two to ten. All rights reserved. No reuse allowed without permission. (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. The copyright holder for this preprint this version posted June 20, 2020. . https://doi.org/10.1101/2020.06.19.20134437 doi: medRxiv preprint (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. The copyright holder for this preprint this version posted June 20, 2020. . https://doi.org/10.1101/2020.06.19.20134437 doi: medRxiv preprint (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. The copyright holder for this preprint this version posted June 20, 2020. . https://doi.org/10.1101/2020.06.19.20134437 doi: medRxiv preprint (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. The copyright holder for this preprint this version posted June 20, 2020. . https://doi.org/10.1101/2020.06.19.20134437 doi: medRxiv preprint It has been seen from above that, in all the state except Gujrat & Karnataka, the lockdown as a containment measure has been quite successful. In India the lockdown started from from 25 th March-2020. After 11 to 13 days of the commencement of the same, seven provinces of India has shown us a signicant downtrend of instantaneous reproduction number. This lag between cause & its eect corresponds to serial time or the generation time interval, that varies from state to state. However, these two states show a opposite trend. After the passage of about one generational time interval, the instantaneous reproduction number peaks sharply. This might correspond to migration [24] at the beginning of the lockdown. In Karnataka, however the value of instantaneous reproduction number uctuates moderately. This is perhaps due to clustering [25] or inadequate testing policies [26] , which may be true for both of the states. In what follows next are two province specic case studies. To be doubly sure, that our time dependent reproduction numbers [27] are calculated correctly over time, we shall t the daily incidence graph of the two provinces which showed contrarian nature in reproduction numbers, during the initiation phase of lockdown. The time Dependent reproduction numbers shall be used to t the daily incidences. The result is given in g. (11) . We show that, the lockdown in India was fairly successful barring a couple of places, due to migration or superspreading etc. We note here that a similar study, with a bigger scope has been reported elsewhere [28] . But it assumes the parametric serial interval, which is dierent from ours. We have deduced our own serial interval (cf. sec. 3) by scraping the pan India raw data and by building our own line list & contact list. Hence Our result is presumed to be signicantly dierent from theirs and more representative of the actual scenarios [29] . The eect partial lifting of the lockdown (unlock) is also seen in the results, in terms of increment in R i (t). 8. Acknowledgement SM wishes to thank K. Bhattacharya, K. Samanta for useful discussions. He also thanks I. Mukhopadhyay for useful help with references. The analysis was was performed in R [30] statistical programming language environment. Public Health Interventions and SARS Spread Avian Inuenza Virus (H5N1): A Threat to Human Health The 2015 inuenza A (H1N1) pdm09 outbreak in India Middle East respiratory syndrome coronavirus (MERS-CoV): A review Ebola virus disease epidemic in West Africa: lessons learned and issues arising from West African countries Incubation period of 2019 novel coronavirus (2019-nCoV) infections among travellers from Wuhan, China Severe Acute Respiratory Syndrome Coronavirus 2: Biology and Therapeutic Options Epidemiology, causes, clinical manifestation and diagnosis, prevention and control of coronavirus disease (COVID-19) during the early outbreak period: a scoping review Should We All Be Wearing Face Masks? Here's Why Experts Are So Conicted. Newspaper What We Know About Face Shields and Coronavirus. Newspaper Lockdown guidelines: Govt's Standard Operating Procedure for social distancing at workplace. Newspaper ETBFSI: Corona impact -Modi announces Janata Curfew -urges citizens to stay at home on Sunday. Newspaper Guidelines for demarcation of containment zones to control COVID-19. Newspaper Notes On R 0 TOI: How eective has India's lockdown been in controlling Covid. Newspaper Why India's Lockdown has been a spectacular failure. Newspaper Teunis: Dierent epidemic curves for severe acute respiratory syndrome reveal similar impacts of control measures How will country-based mitigation measures inuence the course of the COVID-19 epidemic? How Generation Intervals Shape the Relationship Between Growth Rates and Reproductive Numbers Force of infection is key to understanding the epidemiology of Plasmodium falciparum malaria in Papua New Guinean children Estimation of the time-varying reproduction number of COVID-19 outbreak in China Improved inference of time-varying reproduction numbers during infectious disease outbreaks Git Hub [24] PTI: 1,800 people stranded in Uttarakhand to return to Gujarat in 28 buses. Newspaper Ballari's Jindal steel plant emerges COVID-19 cluster. Newspaper Even as COVID-19 Spikes in Rural Karnataka, State Tests 50Less. Newspaper BoÃlle: The R0 package: a toolbox to estimate reproduction numbers for epidemic outbreaks Lockdown Eect on COVID-19Spreadin India: National Data Masking State-Level Trends How Data Became One of the Most Powerful Tools to Fight an Epidemic. Newspaper R: A Language and Environment for Statistical Computing. Vienna, Austria: R Foundation for Statistical Computing