Volume 26, Number 6—June 2020
Research
Estimating Risk for Death from Coronavirus Disease, China, January–February 2020
Abstract
Since December 2019, when the first case of coronavirus disease (COVID-19) was identified in the city of Wuhan in the Hubei Province of China, the epidemic has generated tens of thousands of cases throughout China. As of February 28, 2020, the cumulative number of reported deaths in China was 2,858. We estimated the time-delay adjusted risk for death from COVID-19 in Wuhan, as well as for China excluding Wuhan, to assess the severity of the epidemic in the country. Our estimates of the risk for death in Wuhan reached values as high as 12% in the epicenter of the epidemic and ≈1% in other, more mildly affected areas. The elevated death risk estimates are probably associated with a breakdown of the healthcare system, indicating that enhanced public health interventions, including social distancing and movement restrictions, should be implemented to bring the COVID-19 epidemic under control.
Since the first case of coronavirus disease (COVID-19) was identified in December 2019 in the city of Wuhan in the Hubei Province of China, the novel virus (severe acute respiratory syndrome coronavirus 2 [SARS-CoV-2]) has continued to spread around the world, resulting in several thousand reported cases in multiple countries. In China, the cumulative number of reported deaths was 2,858 as of February 28, 2020, a figure that already dwarfed the number of persons that succumbed to severe acute respiratory syndrome during 2002–2003 (1).
In the context of an emerging infectious disease with pandemic potential, assessing its efficiency at spreading between humans is critical, as is determining the associated risk for death from the disease. In particular, the type and intensity of public health interventions are often set as a function of these epidemiologic metrics. In the absence of vaccines against SARS-CoV-2 or antiviral drugs for the treatment of COVID-19, the implementation of handwashing and other hygiene-related interventions, as well as nonpharmaceutical interventions such as social distancing and movement restrictions (all of which are the basic strategies available to mitigate disease spread in the population), also generate considerable pressure on the global economy (2).
As interventions are gradually implemented and calibrated during the course of an outbreak, early estimates of the case-fatality ratio (CFR) provide crucial information for policymakers to decide the intensity, timing, and duration of interventions. However, the assessment of epidemiologic characteristics, including the CFR, during the course of an outbreak tends to be affected by right censoring and ascertainment bias (3–5). The phenomenon of right censoring is caused by the gap in illness onset to death between the vulnerable population and the healthy population, resulting in underestimation, whereas ascertainment bias is attributable to the unreported bulk of infected persons who have mild symptoms or asymptomatic infections, potentially leading to overestimation. Assuming that ascertainment bias is consistent, we can adjust for right censoring by using established methods and available data (6,7). To assess the current severity of the epidemic in China, we derived estimates (and quantified uncertainty) of the time-delay adjusted CFR for COVID-19 for the city of Wuhan and for China excluding Wuhan, with quantified uncertainty.
Data Sources
We used 2 different types of epidemiologic data in our analysis. First, we extracted the daily series of confirmed cases and deaths in China from daily reports published by the respective governments of China, Hubei Province, and the city of Wuhan (8–11). Because >50% of the deaths are occurring in Hubei Province, and most of these have occurred in Wuhan, we categorized the data by geographic area: Wuhan City, Hubei Province excluding Wuhan, or China excluding Hubei Province. Diagnosis of COVID-19 relies solely on PCR testing because rapid diagnostic tests for this novel coronavirus are not widely available. Our analysis relies on epidemiologic data reported through February 11, 2020, because of the change in case definition that was announced by the government of China on February 12 (12).
We then obtained from several sources a total of 50 epidemiologic descriptions of patients who died from COVID-19 (9–11). After we checked for duplication and missing data, the sample size with data available was 39 patients for observed delays from illness onset to death and 33 for observed delays from hospitalization to death. We fitted a gamma distribution, an exponential distribution, and a lognormal distribution to these distributions and selected the best model based on the Akaike information criterion (AIC) (Appendix 1). The gamma distribution yielded the best fit for the distribution of delays from hospitalization to death (AIC 202.0), whereas the log-normal distribution gave the best fit for the distribution of delays from illness onset to death (AIC 263.3). On the basis of these 2 delay distributions, we incorporated the distribution of delays from hospitalization to death into the model.
Case-Fatality Ratio
We defined crude CFR as the number of cumulative deaths divided by the number of cumulative cases at a specific point in time. To estimate CFR in real time, we used the delay from hospitalization to death, hs, which is assumed to be given by hs = H(s) – H(s-1) for s>0 where H(s) is a cumulative density function of the delay from hospitalization to death and follows a gamma distribution with mean 10.1 days and SD 5.4 days, obtained from the available observed data. If πa,ti is the time-delay adjusted CFR on reported day ti in area a, the likelihood function of the estimate πa,ti is
where ca,t represents the number of new cases with reported day t in area a, and Da,ti is the cumulative number of deaths until reported day ti in area a (6,7). Among the cumulative cases with reported day t in area a, Da,ti have died, and the remainder have survived the infection. The contribution of those who have died with biased death risk is shown in the middle parenthetical term, and the contribution of survivors is shown in the right parenthetical term. We assume that Da,ti is the result of the binomial sampling process with probability πa,ti.
We estimated model parameters by using a Markov chain Monte Carlo method in a Bayesian framework. We estimated posterior distributions of the model parameters by sampling from the 3 Markov chains. For each chain, we drew 100,000 samples from the posterior distribution after a burn-in of 20,000 iterations. We evaluated convergence of Markov chain Monte Carlo chains by using the potential scale reduction statistic (13,14). Estimates and 95% credibility intervals (CrIs) for these estimates are based on the posterior probability distribution of each parameter and based on the samples drawn from the posterior distributions. All statistical analyses were conducted in R version 3.6.1 (R Foundation for Statistical Computing, https://www.r-project.org) using the rstan package.
As of February 11, 2020, a total of 44,795 cases of COVID-19 had been reported in China, 1,117 of which had resulted in death (9–11; Appendix 2). Of the 44,795 cases reported in China, 19,559 cases (43.7%) occurred in Wuhan, 13,894 cases (31.0%) occurred in Hubei Province excluding Wuhan, and 11,342 cases (25.3%) occurred in China excluding Hubei Province. Of the 1,117 deaths in China, 820 (73.4%) occurred in Wuhan, 248 (22.2%) occurred in Hubei Province excluding Wuhan, and 49 (4.4%) occurred in China excluding Hubei Province.
We charted the cumulative cases and deaths in Wuhan, Hubei Province excluding Wuhan, and China excluding Hubei Province (Figure 1). The curve of the cumulative number of deaths grows after that of the cumulative number of cases. Moreover, the increase in the number of deaths in Wuhan occurred more rapidly and the associated mortality rate was much higher than for the rest of China, whereas the cumulative case counts for the 3 areas in China are relatively similar.
We also charted the observed and model-based posterior estimates of crude CFR and the model-based posterior estimates of the time-delay adjusted CFR for Wuhan, Hubei Province excluding Wuhan, and China excluding Hubei Province (Figure 2). Our model-based crude CFR fitted the observed data well throughout the course of the epidemic except for the very early stage. During the course of the outbreak, our model-based posterior estimates of time-delay adjusted CFR have much higher values than the observed crude CFR, except for the early stage in Wuhan and the later stage in China excluding Hubei Province. Our estimates of the time-delay adjusted CFR appear to be decreasing almost consistently in Hubei Province excluding Wuhan and in China excluding Hubei Province, whereas in Wuhan, estimates were low at the early stage and then increased and peaked in the middle of the study period; the Wuhan estimates then followed a decreasing trend similar to the other 2 areas, reaching ≈12%.
As of February 11, estimates of the time-delay adjusted CFR were 12.2% (95% CrI 11.3%–13.1%) in Wuhan, 4.2% (95% CrI 3.7%–4.7%) in Hubei Province excluding Wuhan, and 0.9% (95% CrI 0.7%–1.1%) in China excluding Hubei Province. The observed crude CFR was 4.2% (95% CI 3.9%–4.5%) in Wuhan, 1.8% (95% CI 1.6%–2.0%) in Hubei Province excluding Wuhan, and 0.43% (95% CI 0.32%–0.57%) in China excluding Hubei Province (Table; Figure 3).
We have derived estimates of the CFR for the ongoing COVID-19 epidemic in China. We have estimated time-delay adjusted CFR in 3 different geographic areas in China and found that the most severely affected areas were Wuhan as well as Hubei Province excluding Wuhan, whereas the rest of China (China excluding Hubei Province) experienced a less severe impact.
Our latest estimates (as of February 11, 2020) of the delay-adjusted CFR in Wuhan reach values as high as 12.2% (95% CrI 11.3%–13.1%), an estimate that is 3-fold higher than our estimate for Hubei Province excluding Wuhan and ≈14-fold higher than our estimate for China excluding Hubei Province. These findings suggest that the situation in Wuhan has been particularly dire compared with the other affected areas in China. We note that the upward trend of CFR during the early phase generally indicates increasing ascertainment bias.
An upward trend in the CFR should be interpreted with caution. Diagnosing cases of COVID-19 is difficult because the associated symptoms are not specific. Further, the fraction of asymptomatic patients with SAR-CoV-2 infection and COVID-19 patients who have mild symptoms is not minor; this fact complicates detection and diagnosis early after illness onset, leading to ascertainment bias (15,16). Indeed, out of a total of 566 residents of Japan who evacuated Wuhan by government-chartered plane during January 29–31, a total of 5 asymptomatic and 4 symptomatic COVID-19 patients were detected after undergoing detailed medical examinations (17). However, considering that this underestimation occurred during the course of outbreak and the number of deaths is reported fairly accurately, the upward trend indicates that the temporal disease burden exceeded the capacity of healthcare facilities and the surveillance system probably missed many cases during the early phase. In addition, hospital-based transmission has occurred, potentially affecting healthcare workers, inpatients, and visitors at healthcare facilities, which might explain an increasing trend and the elevated CFR estimates. Indeed, thousands of healthcare workers have succumbed to the disease in China (18), a pattern that resembles past nosocomial outbreaks of Middle East respiratory syndrome (MERS) and severe acute respiratory syndrome (19,20). During past MERS outbreaks, inpatients with underlying disease or elderly persons infected in the hospital setting have raised the CFR to values as high as 20% (21,22). A growing body of evidence indicates that COVID-19 transmission is facilitated in confined settings; for example, a large cluster (634 confirmed cases) of COVID-19 secondary infections occurred aboard a cruise ship in Japan, representing about one fifth of the persons aboard who were tested for the virus. This finding indicates the high transmissibility of COVID-19 in enclosed spaces (23,24).
A downward trend in CFR is suggestive of the extent of improvements in epidemiologic surveillance. In addition, this pattern indirectly indicates a substantial number of mild or asymptomatic cases in Wuhan and that the underlying transmission might prolong the end of the outbreak or further transmission to other areas unless effective social distancing measures are implemented until a vaccine becomes available. Furthermore, given that the delay-adjusted CFR and crude CFR estimates in Wuhan are ≈14-fold higher than our estimates for China excluding Hubei Province, a breakdown in healthcare delivery probably occurred, underscoring the critical need for urgent medical support in the epicenter of the epidemic.
We also found that the estimates of the delay-adjusted CFR for Hubei Province excluding Wuhan and for China excluding Hubei Province showed a declining trend as the epidemic progressed. A similar trend was previously reported for the 2015 MERS outbreak in South Korea, where a substantial fraction of the case-patients were elderly or had underlying conditions (19,20). The high proportion of vulnerable case-patients at the early phase of the outbreak and the smaller number in the later stage could partly explain the observed decline. However, because the epidemic had yet to peak, this time-dependent decrease was probably caused by ascertainment bias. Moreover, the latest estimates of the delay-adjusted CFR and crude CFR in Hubei Province are ≈5-fold higher than our estimate for China excluding Hubei Province, where the healthcare system has not been overwhelmed. These findings also indicate the need to anticipate additional medical support to deliver medical care to the most vulnerable patients, including those with preexisting health conditions, who are at the highest risk for succumbing to the disease. For comparison, the crude CFR has been estimated at 0.9% in Beijing (25), 1.4% among 1,099 patients across China (26), and 4.3% in a meta-analysis among 50,466 hospitalized patients (27).
Our study has limitations. First, our CFR estimate is influenced by ascertainment bias, which might influence estimates upward. For those infectious diseases characterized by a large fraction of patients with mild illness or asymptomatic infections, the infection-fatality risk (e.g., the number of deaths divided by the total number of persons infected) is a more appropriate index of disease burden (28,29). Therefore, mass serologic surveillance and surveys to assess the presence or absence of symptoms is strongly recommended to disentangle the threat of emerging infectious diseases, including COVID-19. In addition, because our estimates of CFR are based on the number of confirmed cases reported before the February 12 change in the case definition, caution will be needed when comparing our estimates with other CFR estimates that include epidemiologic data from on or after February 12, which would be lower. Second, in our estimation we employed a distribution of delays from illness onset to death (n = 39 patients), which was obtained from secondary sources, but the available epidemiologic data does not include either the date of illness onset or the date of confirmation. For this reason, we used the time delay from hospitalization to death (n = 33 patients).
In conclusion, our estimates of the risk for death from COVID-19 in China as of February 11, 2020, were as high as 12% in the epicenter of the epidemic and as low as ≈1% in the less severely affected areas in China. Because the risk for death from COVID-19 is probably associated with a breakdown of the healthcare system in the absence of pharmaceutical interventions (i.e., vaccination and antiviral drugs), enhanced public health interventions (including social distancing measures, quarantine, enhanced infection control in healthcare settings, and movement restrictions), as well as enhanced hygienic measures in the general population and an increase in healthcare system capacity, should be implemented to rapidly contain the epidemic.
Dr. Mizumoto works as an assistant professor at the Graduate School of Advanced Integrated Studies in Human Survivability, Kyoto University, Japan. His research interests include mathematical and statistical epidemiology and modeling of infectious diseases.
Dr. Chowell is professor of Epidemiology and Biostatistics and chair of the Department of Population Health Sciences at the Georgia State University School of Public Health. He is also a senior research fellow at the Fogarty International Center of the National Institutes of Health. His research interests include mathematical modeling of infectious disease dynamics and control.
Acknowledgment
K.M. acknowledges support from the Japan Society for the Promotion of Science (KAKENHI grant no. 15K20936), the Program for Advancing Strategic International Networks to Accelerate the Circulation of Talented Researchers (grant no. G2801), and the Leading Initiative for Excellent Young Researchers from the Ministry of Education, Culture, Sport, Science, and Technology of Japan. G.C. acknowledges support from the National Science Foundation (grant no. 1414374) as part of the joint National Science Foundation–National Institutes of Health–US Department of Agriculture Ecology and Evolution of Infectious Diseases Program, and support from the United Kingdom Biotechnology and Biological Sciences Research Council (grant no. BB/M008894/1).
References
- World Health Organization. Summary of probable SARS cases with onset of illness from 1 November 2002 to 31 July 2003 [cited 2020 Feb 7]. https://www.who.int/csr/sars/country/table2004_04_21
- Chowell G, Fenimore PW, Castillo-Garsow MA, Castillo-Chavez C. SARS outbreaks in Ontario, Hong Kong and Singapore: the role of diagnosis and isolation as a control mechanism. J Theor Biol. 2003;224:1–8. DOIPubMedGoogle Scholar
- Ghani AC, Donnelly CA, Cox DR, Griffin JT, Fraser C, Lam TH, et al. Methods for estimating the case fatality ratio for a novel, emerging infectious disease. Am J Epidemiol. 2005;162:479–86. DOIPubMedGoogle Scholar
- Kucharski AJ, Edmunds WJ. Case fatality rate for Ebola virus disease in west Africa. Lancet. 2014;384:1260. DOIPubMedGoogle Scholar
- Jewell NP, Lei X, Ghani AC, Donnelly CA, Leung GM, Ho LM, et al. Non-parametric estimation of the case fatality ratio with competing risks data: an application to Severe Acute Respiratory Syndrome (SARS). Stat Med. 2007;26:1982–98. DOIPubMedGoogle Scholar
- Nishiura H, Klinkenberg D, Roberts M, Heesterbeek JA. Early epidemiological assessment of the virulence of emerging infectious diseases: a case study of an influenza pandemic. PLoS One. 2009;4:
e6852 . DOIPubMedGoogle Scholar - Tsuzuki S, Lee H, Miura F, Chan YH, Jung SM, Akhmetzhanov AR, et al. Dynamics of the pneumonic plague epidemic in Madagascar, August to October 2017. Euro Surveill. 2017;22:22. DOIPubMedGoogle Scholar
- Garske T, Legrand J, Donnelly CA, Ward H, Cauchemez S, Fraser C, et al. Assessing the severity of the novel influenza A/H1N1 pandemic. BMJ. 2009;339(jul14 3):b2840.
- The State Council of the People’s Republic of China. Update on new coronavirus pneumonia [in Chinese] [cited 2020 Feb 7]. http://www.nhc.gov.cn/yjb/pzhgli/new_list.shtml
- Health Commission of Hubei Province. China. Pneumonia epidemic prevention and control of new coronavirus infection [in Chinese] [cited 2020 Feb 7]. http://wjw.hubei.gov.cn/fbjd/dtyw
- Health Commission of Wuhan City. Hubei Province, China. Wuhan Municipal Commission of Health and Health on pneumonia of new coronavirus infection [in Chinese] [cited 2020 Feb 7]. http://wjw.hubei.gov.cn/fbjd/dtyw
- The State Council of the People’s Republic of China. Clinical guideline for COVID-19, version 5 [in Chinese]. 2020 [cited 2020 Feb 29]. http://www.gov.cn/zhengce/zhengceku/2020-02/05/5474791/files/de44557832ad4be1929091dcbcfca891.pdf
- Gamerman D, Lopes HF. Markov chain Monte Carlo: stochastic simulation for Bayesian inference. 2nd edition. London: Chapman and Hall/CRC Press; 2006.
- Gelman A, Rubin DB. Inference from iterative simulation using multiple sequences. Stat Sci. 1992;7:457–72. DOIGoogle Scholar
- Rothe C, Schunk M, Sothmann P, Bretzel G, Froeschl G, Wallrauch C, et al. Transmission of 2019-nCoV infection from an asymptomatic contact in Germany. N Engl J Med. 2020;382:970–1. DOIPubMedGoogle Scholar
- Kupferschmidt K. Study claiming new coronavirus can be transmitted by people without symptoms was flawed [cited 2020 Feb 4]. https://www.sciencemag.org/news/2020/02/paper-non-symptomatic-patient-transmitting-coronavirus-wrong
- Ministory of Health, Labour and Welfare, Japan. List of press releases (new coronavirus) [in Japanese]. 2020 [cited 2020 Feb 28]. https://www.mhlw.go.jp/stf/seisakunitsuite/bunya/0000121431_00086.html
- Gan N, Thomas N, Culver D. CNN. Over 1,700 frontline medics infected with coronavirus in China, presenting new crisis for the government. 2020 [cited 2020 Mar 2]. https://edition.cnn.com/2020/02/13/asia/coronavirus-health-care-workers-infected-intl-hnk/index.html
- Chowell G, Abdirizak F, Lee S, Lee J, Jung E, Nishiura H, et al. Transmission characteristics of MERS and SARS in the healthcare setting: a comparative study. BMC Med. 2015;13:210. DOIPubMedGoogle Scholar
- Abdirizak F, Lewis R, Chowell G. Evaluating the potential impact of targeted vaccination strategies against severe acute respiratory syndrome coronavirus (SARS-CoV) and Middle East respiratory syndrome coronavirus (MERS-CoV) outbreaks in the healthcare setting. Theor Biol Med Model. 2019;16:16. DOIPubMedGoogle Scholar
- Mizumoto K, Endo A, Chowell G, Miyamatsu Y, Saitoh M, Nishiura H. Real-time characterization of risks of death associated with the Middle East respiratory syndrome (MERS) in the Republic of Korea, 2015. BMC Med. 2015;13:228. DOIPubMedGoogle Scholar
- Mizumoto K, Saitoh M, Chowell G, Miyamatsu Y, Nishiura H. Estimating the risk of Middle East respiratory syndrome (MERS) death during the course of the outbreak in the Republic of Korea, 2015. Int J Infect Dis. 2015;39:7–9. DOIPubMedGoogle Scholar
- Mizumoto K, Kagaya K, Zarebski A, Chowell G. Estimating the asymptomatic proportion of coronavirus disease 2019 (COVID-19) cases on board the Diamond Princess cruise ship, Yokohama, Japan, 2020. Euro Surveill. 2020;25:25. DOIGoogle Scholar
- Mizumoto K, Chowell G. Transmission potential of the novel coronavirus (COVID-19) onboard the Diamond Princess Cruises ship, 2020. Infect Dis Model. 2020;5:264–70. DOIGoogle Scholar
- Tian S, Hu N, Lou J, Chen K, Kang X, Xiang Z, et al. Characteristics of COVID-19 infection in Beijing. J Infect. 2020;
S0163-4453(20)30101-8 ; [Epub ahead of print].PubMedGoogle Scholar - Guan WJ, Ni ZY, Hu Y, Liang WH, Ou CQ, He JX, et al. China Medical Treatment Expert Group for Covid-19. Clinical characteristics of coronavirus disease 2019 in China. N Engl J Med. 2020;
NEJMoa2002032 ; [Epub ahead of print]. DOIGoogle Scholar - Sun P, Qie S, Liu Z, Ren J, Li K, Xi J. Clinical characteristics of 50466 hospitalized patients with 2019-nCoV infection. J Med Virol. 2020;
jmv.25735 ; [Epub ahead of print]. DOIGoogle Scholar - Wong JY, Wu P, Nishiura H, Goldstein E, Lau EH, Yang L, et al. Infection fatality risk of the pandemic A(H1N1)2009 virus in Hong Kong. Am J Epidemiol. 2013;177:834–40. DOIPubMedGoogle Scholar
- Presanis AM, De Angelis D, Hagy A, Reed C, Riley S, Cooper BS, et al.; New York City Swine Flu Investigation Team. The severity of pandemic H1N1 influenza in the United States, from April to July 2009: a Bayesian analysis. PLoS Med. 2009;6:
e1000207 . DOIPubMedGoogle Scholar
Figures
Table
Cite This ArticleOriginal Publication Date: March 13, 2020
Table of Contents – Volume 26, Number 6—June 2020
EID Search Options |
---|
Advanced Article Search – Search articles by author and/or keyword. |
Articles by Country Search – Search articles by the topic country. |
Article Type Search – Search articles by article type and issue. |
Please use the form below to submit correspondence to the authors or contact them at the following address:
Kenji Mizumoto, Graduate School of Advanced Integrated Studies in Human Survivability, Kyoto University, 1 Yoshida-Nakaadachi-cho, Sakyo-ku, Kyoto 606-8306, Japan
Top