Abstract
Novel coronavirus (COVID-19) rapidly spread from China to other parts of the world. Knowledge of incubation period and reproduction number is important in controlling any epidemic. The distribution of these parameters helps estimate the epidemic size and transmission potential of the disease. We estimated incubation period and reproduction number of COVID-19 for India utilizing data reported by Ministry of Health and Family Welfare (MoHFW), Government of India (GOI) and data in public domain. The mean incubation period seems to be larger at 6.93 (SD=±5.87, 95% CI: 6.11-7.75). and 95th percentile estimate for best fit normal distribution is 17.8 days. Weibull distribution, the best fit for the reproduction number estimated pre lockdown reproduction number as 2.6 (95% CI=2.34 - 2.86) and post lockdown reduced to 1.57 (95% CI=1.3 – 1.84) implying effectiveness of the epidemic response strategies. The herd immunity is estimated between 36-61% for R0 of 1.57 and 2.6 respectively.
Introduction
While the novel coronavirus (COVID-19) spread rapidly from China to other developed countries, India saw a steady flow of patients mod early March and by 10 May 2020, it had gripped the country with 63,420 confirmed cases and 2109 deaths1-2. The novel coronavirus (SARS-CoV-2) though related is distinct from severe acute respiratory syndrome (SARS) coronavirus and Middle East respiratory syndrome (MERS) coronavirus3. Many researchers struggled to estimate the magnitude of the epidemic wherein the epidemiological parameters remained uncertain. Knowledge of key epidemiological parameters including incubation period and reproduction number is important in controlling any epidemic. The distribution of these parameters helps estimate the epidemic size and transmission potential3-7 of the disease.
The incubation period is defined as the time from infection to the onset of illness8 and is crucial for epidemiological modelling in predicting the transmission dynamics, infectiousness and quarantine period9. It is also important for several important public health activities like length of active monitoring, surveillance and control.
The reproduction number (R0) is the most fundamental parameter in infectious disease dynamics describing the contagiousness or transmissibility of infectious agents and is defined as the average number of secondary cases caused by a single infectious individual in a entirely susceptible population10. An outbreak is expected to continue if R0 has a value >1 and to end if R0 is <1. R0 fluctuates if the rate of human–to-human or human and vector interactions varies over time or space. There exists scant evidence which supports the applicability of R0 outside the region for which the value was calculated11. Estimation of changes in transmission over time can provide insights into the epidemiological situation and identify whether outbreak control measures are adequate and are having the desired measurable effect, and help in undertaking midcourse corrections.
Herd immunity is defined as the resistance to the spread of a contagious disease within a population that results if a sufficiently high proportion of individuals are immune to the disease, especially through vaccination or immunity post natural infection12. When a high proportion of the population is immune, it is difficult for infectious diseases to spread, because there are not many people who can be infected and the transmission chain gets broken.
Our current understanding of these epidemiological parameters for India is limited. Hence this study was undertaken to address above issue and estimate incubation period and reproduction number of COVID-19 for India utilizing data reported by Ministry of Health and Family Welfare (MoHFW), Government of India (GOI) and the data available in public domain.
Materials and Methods
Data sources
The analysis is based on publicly available data. Data were retrieved from the official website of the MoHFW, GOI1 and covid19indi tracker2. The other sources for data aggregation and triangulation were, World Health Organization (WHO)13, Centre for Disease Control (CDC)14, update reports of Indian Council of Medical Research (ICMR)15, and India Situation Reports on Novel Coronavirus (2019-nCoV)16. Since the analysis is based on publicly available data ethics approval was not required.
Information on demographic characteristics, date on exposure and date of confirmation of disease status/illness onset were extracted for estimating incubation period and adjusted for delay-time. The delay time in getting the test results/diagnosis as confirmed case of COVID-19 was adjusted based on the delay time available from states and MoHFW reports. For estimating reproduction number, details of the source of infection from positive case giving number of probable cases infected by suspect was extracted. The number of people who contracted the disease from primary case.
Statistical distributions for Incubation period and reproduction number
The incubation period data post adjustment of delay-time in test results was subjected to best fit model. Besides normal distribution four other commonly used incubation period distribution (Weibull, Log normal, Gamma and Erlang) were fitted. Estimation of the median incubation period, mean (sd), and quantiles (5th, 25th, 75th, 95th and 97th) was also done.
For reproduction number, the best fit model (Weibull distribution) was used. Besides the best, fit three more distributions based on review of literature for distribution of reproduction number (lognormal, Gamma and Exponential) were considered. The values of R0 were estimated for pre and post lockdown period to evaluate the impact of the epidemic control measures instituted by the Government of India. The herd immunity (HI) estimate was based on R0 value12. (HI= (R0 – 1) or R0 = 1 − 1/R0) The statistical software used were IBM SPSS Statistics for Windows version 23.0 (SPSS Inc., Chicago, Ill., USA) and EASY-FIT, a software system for data fitting in dynamical systems.
Results
Data on 268 lab confirmed cases were extracted, the range of age was from 1.5 years to 89 years with a mean age of 36.45 years (SD=±17.27). Ratio of male to female was 1.5:1 with 60.3% males and 39.7% females.
Incubation Period
Using date of exposure and date of confirmation of disease status/illness onset, the estimates of incubation period were determined. Table 1 gives the incubation period estimates for various distributions. We fitted five distributions to the data: Normal, Weibull, Log normal, Gamma and Erlang. The normal distribution provided best fit for data with median and mean incubation period of 6.93 (SD=±5.87, 95% confidence interval CI: 6.11-7.75). The incubation period ranged from 1 to 19.26 days (5th percentile to 97th percentile) for this best fit. The Weibull distribution was the best next fit with mean incubation period to be higher than normal distribution with mean of 8.17 (SD=8.10). Figure 1 show cumulative distribution functions for best fit. The median incubation period for other distributions ranged 3.46 to 6.06. The probability and cumulative distribution functions for various distributions in order of best fit is provided in the supplementary.
Cumulative distribution function for best fit Normal distribution for Incubation period for India
Reproduction Number and Herd Immunity
The reproduction number was estimated for days before lockdown and post lockdown to assess the impact of various control measures to include social distancing adopted by India. Weibull distribution was the best fit for the reproduction number followed by Log Normal and Gamma distribution. We estimated the initial reproduction number before lockdown by GOI to be 2.6 (95% CI=2.34 - 2.86) and post lockdown the reproduction number reduced to 1.57 (95% CI=1.3 – 1.84). The herd immunity is estimated in the range of 61-62%. Table 2 gives the descriptive statistics of reproduction number before and post lock down.
Discussion
This data driven paper to the best of our knowledge presents the estimates of incubation period and reproduction number of COVID-19 in India for the first time. We characterised the distribution of incubation period and reproduction number for COVID-19 for India. We estimated that the incubation period follows the best fit of normal distribution with around 7 days ranging from 1 to 22 day with 95th percentile of the distribution at 17.8 days. The incubation period seems to be longer for India as compared to 2 to 10 days by WHO13 and 2-12 days by ECDC17. The United States’ CDC18 estimates the incubation period for COVID-19 to be between 2 and 14 days. Most of existing studies are on data originating from Wuhan, China. A study on travellers from Wuhan on 88 cases estimated the mean incubation period to be 6.4 days with a range from 2.1 to 11.1 days. The upper limit of 11.1 days was considered conservative19. Two other studies from Wuhan, China on early transmission dynamics of COVID-19 estimated the mean incubation period as 5.2 days (95% CI, 4.1 to 7.0), with the 95th percentile of the distribution at 12.5 days20 and other study with 2-14 days with Log Normal as best fit21. A study by John Hopkins Bloomberg School of Public Health on regions and countries outside Hubei22 gave median incubation period as 5.1 days (95% CI, 4.5 to 5.8 days), and also stated that 97.5% of those who develop symptoms will do so within 11.5 days (CI, 8.2 to 15.6 days) of infection. Understanding incubation period is crucial so as to introduce and implement effective quarantine duration and preventing the spread of the virus.
Weibull distribution seems to be the best fit for reproduction number distribution in India. The reproduction number for India is estimated to be 2.6 in pre lock down phase and reduced to around 1.6 post lockdown phase with a herd immunity threshold of 61.5%. Though the value of R0 has reduced from 2.6 in pre lock down period to 1.57 in post lock down period. i.e. 60.38% reduction it is still greater than 1. COVID-19 epidemic will increase as long as R0 is greater than 1 and control efforts to bring R0 below 1 needs to be implemented aggressively.
A report based on the impact of the interventions across 11 European countries estimated a posterior mean of 0.97 [0.14-2.14] for Norway and 2.64 [1.40-4.18] for Sweden, with an average of 1.43 across the 11 country posterior means, a 64% reduction compared to the pre-intervention values23. A study from Wuhan, Hubei Province, China20 estimated an R0 of approximately 2.2, however this estimate was limited to early epidemic period up to January 4. The median R0 from Diamond Princess cruise ship24 was estimated to be 2.28 (95% CI=2.06-2.52). Switzerland25 estimated the basic reproduction number R0 of COVID-19 at 2.78 (95% confidence interval, CI: 2.51 - 3.11). Transmission decreased with the strengthening of social distancing measures by 89% (95% CI: 83%-94%) which brought the R0 to below 1 at 0.32. Study25 by Wu et. al. estimated the reproductive number for COVID-19 as 2·68 (95% CI: 2·47–2·86). The wide range of R0 from different studies indicates the challenges in estimating R0. Additionally, R0 is dynamic in the sense that any factor having the potential to influence the contact rate, including population density, social distancing, and seasonality will ultimately affect R0.
Our estimate of herd immunity threshold as 61% (with R0 =2.6) is staggering as 826,873,494 people need to be infected to achieve herd immunity in India which may result in high death rate. Many countries have implemented aggressive lock downs however some countries like Sweden had a different approach to tackle COVID-19 by taking individual responsibility for social distancing and keeping society functioning with no official lock down by the Sweden government and still successful in keeping the reproduction number in control26. Long term lockdowns are not sustainable as its ill effects impacts many other health programmes. As per WHO27 the weekly detection of new tuberculosis cases in India has gone down by nearly 75 per cent during the COVID-19 lockdown. The malaria modelling analysis by the Global Malaria Programme at WHO28 and the HIV modelling study, convened by the WHO and UNAID29, predicts the collateral damage of COVID will be as devastating as pandemic, with millions succumbing to preventable, treatable illness and disease. However, in absence of vaccine and other uncertainties regarding the virus, herd immunity is debatable issue as a preventive measure and achieving herd immunity is likely to be a long-drawn process.
Limitations of the study
The data on complete information for estimating incubation period as well reproduction number could be extracted for limited number of cases. Also, the time delay of test results, with varying turn around time among States, for estimating incubation period was averaged for India based on the public domain information and media reports.
Conclusions
This study provides the estimates for key epidemiological charachteristics of COVID-19 in India. The mean incubation period seems to be larger at 6.93 and 95th percentile estimate for best fit normal distribution to be 17.8 days. Best fit for reproduction number follows weibull distribution. Our estimates for pre lockdown reproduction number was 2.6. and post lockdown the reproduction number reduced to 1.57. This implies and shows that the epidemic response strategies adopted by India are effective. However the herd immunity is estimated in the range of 36-61% for R0 of 1.57 and 2.6 respectively.
Data Availability
The analysis is based on publicly available data