Abstract
Background India is vulnerable to community infection of COVID-19 due to crowded and poor living condition, high density, slums in urban areas and poor health care system. The number of COVID 19 infection has crossed 300,000 with over 7,500 deaths despite a prolonged period of lock down and restrictions in public spaces. Given the likely scale and magnitude of this pandemic, it is important to understand its impact on the age pattern of mortality under varying scenarios.
Objective The main objective of this paper is to understand the age pattern of mortality under varying scenarios of community infection.
Data and Methods Data from the Sample Registration System (SRS), covidi19india.org and country specific data from worldmeter is used in the analyses. Descriptive statistics, case-fatality ratio, case fatality ratio with 14 days delay, abridged life table,years of potential life lost (YPLL) and disability adjusted life years (DALY) is used.
Results The case fatality ratio (CFR) with 14 days delay for India is at least twice higher (8.0) than CFR of 3.4. Considering 8% mortality rate and varying scenario of community infection by 0.5%, 1% and 2%, India’s life expectancy will reduce by 0.8, 1.5 and 3.0 years and potential life years lost by 12.1 million, 24.3 million and 48.6 million years respectively. A community infection of 0.5% may result in DALY by 6.2 per 1000 population. Major share of PYLL and DALY is accounted by the working ages.
Conclusion COVID-19 has a visible impact on mortality with loss of productive life years in working ages. Sustained effort at containing the transmission at each administrative unit is recommended to arrest mortality owing to COVID-19 pandemic.
What is known?The case fatality rate associated with COVID-19 is low in India compared to many other countries. The mortality level is higher among elderly and people with co-morbidity.
Contribution The case fatality ratio is illusive in the sense that the same with 14 days delay for India is at least twice higher (8.0). The COVID-19 attributable mortality has the potential to reduce the longevity of the population. Unlike developed countries, about half of the COVID-19 attributable mortality would be in the working age group of 45-64 years. With any level of community infection, the years of potential life lost (YPLL) and disability adjusted life years (DALY) world be highest in the working age group (45-64 years).
Introduction
With additional 350 thousand deaths in last five months, the global share of COVID-19 attributable death has increased by 0.05% (https://population.un.org/wpp/DataQuery/). The first case surfacing in Wuhan city of China in December, the estimated number of confirmed COVID-19 cases has crossed 5 million mark by 20th May 2020 (https://www.worldometers.info/). Despite this brief duration of this outbreak, the epidemic has engulfed the entire world. The COVID-19 virus is highly contagious with a longer survival duration affects the vital organs including the respiratory system. It has altered the public health priority globally and nationally and posed a challenge to medical science.
The outbreak of COVID-19 has led to a health emergency and economic crisis worldwide. The global recession following COVID-19 epidemic is projected to be worst ever in human history. In developing countries, the hard hit are the poor, migrant labour, small and petty shopkeeper; who have lost their livelihood (Bhagat et al. 2020). Many adverse economic effects are yet to surface coming days. The adverse impact of COVID-19 on human health includes raised the levels of mortality and morbidity level, treatment complications, increased health spending and reduced spending on non-health goods and services. The case fatality rate from the epidemic is high among elderly and critically ill patients such as cancer patient, hypertensive, diabetic and patients with heart disease and for those with co-morbidity (Chow et al 2020; Bonow et al. 2020; Lippi et al. 2020; Pal & Bhadada 2020; Guan et al. 2020)
The medical know-how on treatment of COVID-19 is limited and there is no effective vaccine to prevent the deadly virus. An effective vaccine for the COVID-19 is set to be in place between six month and two years from now. Besides, the health infrastructure across countries is under severe strain for treatment of COVID-19. Countries such as Italy, United States and Spain with reasonably well-developed health care system in the world could not limit the human loss in wake of COVID-19 epidemic. The Institute for Health Metric has predicted that deaths from COVID-19 in United States would lead to scarcity of medical beds and equipment (Murray et al. 2020).
The first empirical study of Diamond Princes cruise ship of Japan placed the COVID-19 infection fatality rate at 1.3%, corrected case fatality ratio at 2.6% for all age group and 13% among the older aged 70+ (Russell et al. 2020). The mortality rate due to COVID-19 and its bearing on life expectancy is higher in countries with higher proportion of adults and elderly (Lai et al. 2020). The overall crude fatality rate from COVID-19 in China was greater among the elderly aged 60 and above as well as males (Verity et al. 2020; Dudley et al. 2020). The COVID 19 fatality is responsive to population aging and health care resources (Dilcher et al. 2020; Ji et al. 2020). Italy had a larger share of confirmed COVID-19 infection among elderly and those with multiple comorbidities (Onder et al. 2020). Myocardial injuries were significantly associated with the fatal outcome of COVID-19 (Guo et al. 2020).
After 100 days of outbreak of COVID-19 (14th May, 2020), India rank 11th with 2% of the global total COVID-19 confirmed cases and moving at a faster pace to surpass China. The case fatality ratio of India (3.2) is relatively lower but the number of testing per million population is low compared to many other countries. With a population of nearly 1.35 billion, India is at high risk of community transmission of COVID-19. The densely populated cities, with higher proportion of slum population coupled with poor health infrastructure, low income level and greater dependence on informal activities would raise the potential spread of COVID-19 pandemic. Projection suggests that the number of cases in India is reduced due to lockdown measures (Dwivedi et al. 2020) and may rise with ease of movement. The estimated cases are very low owing to low intensity of testing (Goli & James 2020).
Despite nationwide lock down and several effective measures for over two months in India, the number of COVID-19 confirmed cases was close to 170,000 mark with over 4,700 deaths as of 29th May 2020. The disease has reached the stage of community infection in the states of Maharashtra, Gujarat, Delhi and Madhya Pradesh and no attempt has been made to understand the pattern of mortality associated with COVID-19 in India. The aim of this paper is to examine the possible age pattern of mortality under likely event of community infection of COVID-19 in India.
Data and Methods
Data on COVID-19 at national and global level is made available in a seamless manner than any other disease. A number of key variables such as number of infection, deaths and recoveries were available for smallest geographic unit such as district to country level on daily basis. We have used data from the covid19india.org (https://www.covid19india.org/), worldmeter (https://www.worldometers.info/), the Sample Registration System (SRS) (ORGI, 2019) and the population projection of expert committee (MoHFW, 2019).
The covid19india.org provides micro data of patients on various domain such as demographics, death and recovered, testing number, travel history etc. (https://api.covid19india.org/csv). We have used raw data (patients’ data) and death and recovered data, as of 15th May, 2020 (data updated till 9th May). The raw data file provides selected variables such as detected district, state and city, age and sex, current status (hospitalised, recovered and death), nationality and status change date. While data on detected state, nationality and current status are available for most of the cases, data are largely missing for other variables. As of 9th May 2020, there were 62,864 confirmed cases of which 19,301 were recovered and 2,100 were deceased. The age data was available for 7,191 hospitalised and recovered cases and 511 death cases. Besides, data from SRS abridged life table 2013-17 is used for the analyses (ORGI, 2019). We have used the SRS death rate to estimate deaths without COVID-19 since COVID-19 infection started in 2020. The expert committee on population projection of India 2011 has been used in estimation (MoHFW, 2019).
Methods
Descriptive statistics, case-fatality ratio, lag case fatality ratio, estimation of longevity, years of Potential life lost (YPLL) and disability adjusted life years (DALY) are used in analyses. The set of assumptions used in estimation procedure are detailed below.
Assumptions on Mortality without COVID-19 and with COVID-19
All deaths occurring due to COVID-19 are additional deaths in the population that could have been avoided in absence of COVID-19 cases in India
The age pattern of mortality obtained from SRS, 2013-17 was related to deaths without COVID-19 and assumed to remain same for all diseases other than COVID-19 for 2020.
Missing age data on hospitalisation and deaths have similar distribution as that of reported age
Effectiveness of medical treatment of COVID-19 will remain same till the end of the year 2020.
Age data on COVID-19 attributable mortality was available for 24% of all deceased (511) and that of hospitalised and recovery cases were available for 12% of total confirmed cases (7191). We have distributed the total COVID-19 attributed death cases in each of the five-year age group in keeping with the distribution of death for which age data is known. The official population projection of 2020, the age distribution of 2021 (2020 is not available) and the age specific death rate of SRS is used to arrive at expected number of deaths without COVID-19 (MoHFW 2019). The total number of deaths with COVID-19 is sum of deaths with and without COVID-19. The age specific death rate (ASDR) was estimated from all deaths including COVID-19 deaths and the key input in construction of abridge life table. The abridged life table is used to estimate the life expectancy under each of the four scenarios including that of without COVID-19.
Estimation of case fatality ratio and lag case fatality ratio
The case fatality ratio is defined as the cumulative number of confirmed cases and cumulative number of deaths. Though it is a simple indictor and easy to understand, it underestimates the true mortality in the population (Baud et al. 2020). The lag case fatality ratio is the ratio of cumulative number of deaths and cumulative number of confirmed cases 14 days ago.
Estimation of Years of Potential Life Lost (YPLL)
The YPLL is estimated as follows (Gardner & Sanborn, 1990), where,
di = Number of deaths in ith age group
Li = Life expectancy at age i
Estimation of Premature Years of Potential Life Lost (PYPLL)
The PYPLL is estimated as follows (Gardner & Sanborn, 1990), where, i= Age at death
di = Number of deaths at age i
N=Upper age limit (70 years)
Estimation of Disability Adjusted Life Years (DALY)
DALY is a summary measure of population health, combining mortality and disability. The DALY is the sum of years of life lost (YLL) and years lived with disability (YLD).
The basic formula for YLL and YLD are:
YLL= L×N
YLD= I×DW×L
The formula used for estimation of DALY with discounting and age weighting are given below:
YLL = N/r (1- exp (-r*L))
where, N= number of deaths
L= Life expectancy at age of death
r= discount rate (e.g. 3% corresponds to a discount rate of 0.03)
YLD= (I ×DW × L (1- exp (-r*L))/ r
where, I= number of incidence/prevalence cases
DW= disability weight
L= duration of disability
r= discount rate
As COVID-19 is a novel disease, its disability weight is not available. Since COVID-19 is an infectious disease, we have used the weight of Infectious disease (severe) as proxy for COVID-19. Global Burden of Disease and World Health Organization have estimated this disability weight as 0.133 (Salomon et al. 2013). We have used duration of disability as 60 days because the patients of COVID-19 have been hospitalized for about 30 days and after discharged, quarantined for 14-28 days approximately.
Results
Rising Trends in COVID-19 cases
We begin the discussion by presenting the cumulative number of confirmed, deceased and recovered cases of India over a period of over two months (Fig 1). The exponential curve suggests that the average growth rate of 11% for total confirmed cases and 12% for deceased cases. With the current growth rate, the doubling time is estimated at 6.5 days for confirmed cases, 5 days for recovered and 5.5 days for death cases.
Case Fatality Ratio and Lagged Case Fatality Ratio
Fig 2 presents trends in case fatality rate and lagged case fatality rate (LCFR) (with 14 days lag) in India on daily basis for over 47 days. There are at least two limitations of the CFR. First, even if there is no further new infection, there will be additional deaths from hospitalised or infected cases. Second, deaths occurring due to COVID-19 does not relate to infection on same day. The maximum incubation period of COVID-19 infection is 14 days and hence a refined measure that relate the deaths to number of infection 14 days prior to death. The LCFR was fluctuating in early days of infection but remained stable in last one month. From the graph, we infer that 1) the mortality level due to COVID-19 is at least twice higher than CFR 2) The LCFR has declined in initial period but remains stable in recent days.
Table 1 presents the cumulative number of cases, death cases, the CFR and LCFR for selected states of India. Among these states, the LCFR was lowest in Kerala followed by Bihar and highest in West Bengal (29.7) followed by Gujarat (14.4) and Maharashtra (10.3). The LCFR for India is at least twice higher than CFR suggesting that the mortality due to COVID-19 is higher than that is conveyed by CFR.
Table 2 presents the actual number of confirmed cases, adjusted number of cases (moving average) and the case fatality ratio of COVID-19 by age group in India. Total number of adjusted cases in India are distributed as per the distribution of total cases for which age is known and total number of confirmed cases. Findings suggest that about 60% of COVID-19 cases are in the age group of 30-64. About 2.3% deaths are under 15 years, 12.8% are in 15-44 age group, 48.2% deaths are in 45-64 age group and 36.8% are in 65+ age group. A major share of the COVID-19 deaths is in the working age group. The age specific case fatality ratio increases sharply from age 45 years and above.
Fig 3(a) presents the histogram on age of COVID-19 recovery among the hospitalised patients and Fig 3(b) presents that of deaths in India. The median age of hospitalisation cases was 34 years compared to 61 years of deceased. About one-fourth of deceased were below 70 years and that of hospitalised cases were below 46 years.
Probability of Dying with and without COVID-19
Under this hypothetical situation, the age pattern of mortality clearly suggests that the COVID-19 is likely to affect the adults 45-75, majority in working age group. Beyond age 60, the probability of death due to COVID-19 would be at least four times higher among elderly compared with non-elderly. The probability of dying in 15-60 age group would increase from 0.17 without COVID-19 to 0.19 with 0.5% community infection, 0.20% with 1% community infection and 0.23 with 2% community infection.
Life expectancy with and without COVID-19
Table 3 provides life expectancy and reduction in life expectancy at each age under three scenarios. It may be noticed that the life expectancy will be reduced in most of the age group under each scenario. If the community infection reaches at 0.5% by the end of the year, the life expectancy is likely to reduce by 0.8 years from 2013-17 level. If the infection is at 1% level, the life expectancy would reduce by 1.5 years while an infection rate of 2% would reduce the life expectancy by 3.0 years.
Table 4 presents the estimated YPLL in the absence of COVID-19 and at varying level of COVID-19 infection. Without any COVID, YPLL was 399 million from all disease with majority of share among children (55% among 0-1 age group), 15% among 45-64 years and 14% among 65 and above. YPLL was 12 million with 0.5% infection of COVID-19. However, the age pattern of YPLL due to COVID-19 suggests that it is largely affecting the working population, aged 45-64 years. The pattern remains similar even at higher level of community infection.
Table 5 shows the estimated DALY with discounting and age weighting in the three scenarios of COVID-19 infection. The estimated DALY is 6.2 per thousand population with 0.5% infection of COVID-19 which would increase to 12.3 and 24.6 per thousand population with 1% and 2% infection, respectively. Age-group 45-64 makes the highest contribution to the overall DALY. As COVID-19 is a short term but fatal disease, the YLL component of DALY is higher than the YLD.
Discussion and Conclusion
This is the first ever systematic attempt at assessing the likely effect of COVID-19 spread on mortality and longevity in India. It compares a situation with varying levels of community infection of COVID-19. The followings are the salient findings of the study. First, India fares better than global average and many other worst affected countries as regard fatality rate, recovery rate and test detection ratio. However, it the mortality due to COVID-19 is measured using case fatality rate with a lag of 14 days, it is twice as high as conventional CFR. The inter-state variation in LCFR was large with West Bengal, Maharashtra and Gujarat being three worst performing states. Second, the community spread of COVID-19 to be restricted at 0.5% of the population in India would result in 0.8 years compression of the life expectancy at birth. Third, the impact of COVID-19 on premature mortality, particularly among working adults in 45-64 age group is very high. Majority of the deaths being in the ages 0f 45-64, the contribution to YPLL and DALY is likely to be the most by 45-64 years age group. Lastly, this exercise is an instance of effective provisioning of the patient level and aggregated CVID-19 data for research and policy formulation in India. However, care must be taken to sensitise the officials/health workers in providing complete information on key variables to facilitate its scientific treatment.
We provide some explanation and suggestion in support of the findings. The low fatality rate and high test-detection ratio are primarily due to effective measured adopted by Govt. of India in controlling the community spread of the disease. The lockdown measures of about two months is the most proactive step taken by the national government in curbing the epidemic. But the lockdown cannot be continued indefinitely given its economic cost as observed in terms of reassessment of growth. The pace of new infection and deaths due to COVID-19 is soaring every day. It is believed to have reached the community level in the states of Maharashtra, Gujarat and Delhi. The situation is alarming in Maharashtra and particularly in Mumbai city. Many of the COVID-19 patients are not getting hospital in civic run hospitals and families of positive patients are not getting tested anymore. However, it is interesting to note the control of COVID-19 in Kerala and the best practice adopted in the Kerala may well be practised in worst affected states of India. Efforts to minimise the loss of human life needs preparedness of health facilities to cater to severe patients apart from imposing strict adherence of the ‘new normal’ (i.e. social distancing and hand hygiene). Identification of high-risk cases and arrangement of intensive care for such patients would help reducing premature mortality. While community testing in India does not seem viable, it is necessary to provide testing treatment options for this infection in every town and cities of the country.
The study suffers from following limitations. We believe that our estimates for elderly is underestimated as missing data is likely to be higher among elderly compared to younger population. Our assumption of all deaths due to COVID-19 would have been avoided if there were no COVID-19 cases may not always hold true and the assumption of SRS death rate of 2013-17 would prevail is illustrative. Despite these limitations, we believed that this study offers a reasonable understanding of the impact of COVID-19 on longevity in India. It has demonstrated the scope to replicate similar analyses elsewhere and guide the impact assessment of the epidemic in the national and global scene.
Research Ethics Approval
The study does not involve human participation. It used secondary data available in public domain and research ethics approval not required.
Data Availability
Data on COVID-19 at national and global level is made available in a seamless manner than any other disease. A number of key variables such as number of infection, deaths and recoveries were available for smallest geographic unit such as district to country level on daily basis. We have used data from the covid19india.org (https://www.covid19india.org/), worldmeter (https://www.worldometers.info/), the Sample Registration System (SRS) (ORGI, 2019) and the population projection of expert committee (MoHFW, 2019).
https://api.covid19india.org/csv.