Next Article in Journal
Accuracy of the QIAxcel Automated System for MIRU-VNTR Genotyping of Mycobacterium tuberculosis in Two Limited Resource Settings
Next Article in Special Issue
The Rate of Underascertainment of Novel Coronavirus (2019-nCoV) Infection: Estimation Using Japanese Passengers Data on Evacuation Flights
Previous Article in Journal
Characteristics and Dysbiosis of the Gut Microbiome in Renal Transplant Recipients
Previous Article in Special Issue
The Extent of Transmission of Novel Coronavirus in Wuhan, China, 2020
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Estimating the Unreported Number of Novel Coronavirus (2019-nCoV) Cases in China in the First Half of January 2020: A Data-Driven Modelling Analysis of the Early Outbreak

1
JC School of Public Health and Primary Care, Chinese University of Hong Kong, Hong Kong 999077, China
2
Shenzhen Research Institute of Chinese University of Hong Kong, Shenzhen 518060, China
3
Department of Applied Mathematics, Hong Kong Polytechnic University, Hong Kong 999077, China
4
Michigan Institute for Data Science, University of Michigan, Ann Arbor, Michigan, MI 48104, USA
5
School of Public Health, Li Ka Shing Faculty of Medicine, University of Hong Kong, Hong Kong 999077, China
6
Department of Orthopaedics and Traumatology, Chinese University of Hong Kong, Hong Kong 999077, China
7
SH Ho Scoliosis Research Lab, Joint Scoliosis Research Center of Chinese University of Hong Kong and Nanjing University, Hong Kong 999077, China
8
School of Mathematics and Statistics, Huaiyin Normal University, Huaian 223300, China
9
School of Nursing, Hong Kong Polytechnic University, Hong Kong 999077, China
10
Department of Mathematics, Shanghai Normal University, Shanghai 200234, China
*
Authors to whom correspondence should be addressed.
J. Clin. Med. 2020, 9(2), 388; https://doi.org/10.3390/jcm9020388
Submission received: 27 January 2020 / Revised: 29 January 2020 / Accepted: 31 January 2020 / Published: 1 February 2020

Abstract

:
Background: In December 2019, an outbreak of respiratory illness caused by a novel coronavirus (2019-nCoV) emerged in Wuhan, China and has swiftly spread to other parts of China and a number of foreign countries. The 2019-nCoV cases might have been under-reported roughly from 1 to 15 January 2020, and thus we estimated the number of unreported cases and the basic reproduction number, R0, of 2019-nCoV. Methods: We modelled the epidemic curve of 2019-nCoV cases, in mainland China from 1 December 2019 to 24 January 2020 through the exponential growth. The number of unreported cases was determined by the maximum likelihood estimation. We used the serial intervals (SI) of infection caused by two other well-known coronaviruses (CoV), Severe Acute Respiratory Syndrome (SARS) and Middle East Respiratory Syndrome (MERS) CoVs, as approximations of the unknown SI for 2019-nCoV to estimate R0. Results: We confirmed that the initial growth phase followed an exponential growth pattern. The under-reporting was likely to have resulted in 469 (95% CI: 403–540) unreported cases from 1 to 15 January 2020. The reporting rate after 17 January 2020 was likely to have increased 21-fold (95% CI: 18–25) in comparison to the situation from 1 to 17 January 2020 on average. We estimated the R0 of 2019-nCoV at 2.56 (95% CI: 2.49–2.63). Conclusion: The under-reporting was likely to have occurred during the first half of January 2020 and should be considered in future investigation.

1. Introduction

A novel coronavirus (2019-nCoV) infected pneumonia infection, which is deadly [1], was first identified in Wuhan, China in December 2019 [2]. The virus causes a range of symptoms including fever, cough, and shortness of breath [3]. The cumulative number of reported cases slowly increased to cumulative 41 cases by 1 January 2020, and rapidly increased after 16 January 2020. As of 26 January 2020, the still ongoing outbreak had resulted in 2066 (618 of them are in Wuhan) confirmed cases and 56 (45 of them were in Wuhan) deaths in mainland China [4], and sporadic cases exported from Wuhan were reported in Thailand, Japan, Republic of Korea, Hong Kong, Taiwan, Australia, and the United States, please see the World Health Organization (WHO) news release via https://www.who.int/csr/don/en/ from 14 to 21 January 2020. Using the number of cases exported from Wuhan to other countries, a research group at Imperial College London estimated that there had been 4000 (95%CI: 1000–9700) cases in Wuhan with symptoms onset by 18 January 2020, and the basic reproduction number (R0) was estimated at 2.6 (95%CI: 1.5–3.5) [5]. Leung et al. drew a similar conclusion and estimated the number of cases exported from Wuhan to other major cities in China [6], and the potentials of travel related risks of disease spreading was also indicated by [7].

2. Objectives and Methods

Due to an unknown reason, the cumulative number of cases remained at 41 from 1 to 15 January 2020 according to the official report, i.e., no new case was reported during these 15 days, which appears inconsistent with the following rapid growth of the epidemic curve since 16 January 2020. We suspect that the 2019-nCoV cases were under-reported roughly from 1 to 15 January 2020. In this study, we estimated the number of unreported cases and the basic reproduction number, R0, of 2019-nCoV in Wuhan from 1 to 15 January 2020 based on the limited data in the early outbreak.
The time series data of 2019-nCoV cases in mainland China were initially released by the Wuhan Municipal Health Commission from 10 to 20 January 2020 [8], and later by the National Health Commission of China after 21 January 2020 [9]. The case time series data in December 2019 were obtained from a published study [3]. All cases were laboratory confirmed following the case definition by the national health commission of China [10]. We chose the data up to 24 January 2020 instead of to the present study completion date. Given the lag between timings of case confirmation and news release of new cases, the data of the most recent few days were most likely to be tentative, and thus they were excluded from the analysis to be consistent.
We suspected that there was a number of cases, denoted by ξ, under-reported from 1 to 15 January 2020. The cumulative total number of cases, denoted by Ci, of the i-th day since 1 December 2019 is the summation of the cumulative reported, ci, and cumulative unreported cases, Ξi. We have Ci = ci + Ξi, where ci is observed from the data, and Ξi is 0 for i before 1 January and ξ for i after 15 January 2020. Following previous studies [11,12], we modelled the epidemic curve, i.e., the Ci series, as an exponential growing Poisson process. Since the data from 1 to 15 January 2020 appeared constant due to unclear reason(s), we removed these data from the fitting of exponential growth. The ξ and the intrinsic growth rate (γ) of the exponential growth were to be estimated based on the log-likelihood, denoted by , from the Poisson priors. The 95% confidence interval (95% CI) of ξ was estimated by the profile likelihood estimation framework with cutoff threshold determined by a Chi-square quantile [13], χ2pr = 0.95, df = 1. With γ estimated, the basic reproduction number could be obtained by R0 = 1/M(−γ) with 100% susceptibility for 2019-nCoV presumed at this early stage. Here, the function M(∙) was the Laplace transform, i.e., the moment generating function, of the probability distribution for the serial interval (SI) of the disease [11,14], denoted by h(k) and k is the mean SI. Since the transmission chain of 2019-nCoV remained unclear, we adopted the SI information from Severe Acute Respiratory Syndrome (SARS) and Middle East Respiratory Syndrome (MERS), which share the similar pathogen as 2019-nCoV [15,16,17]. We modelled h(k) as Gamma distributions with mean of 8.0 days and standard deviation (SD) of 3.6 days by averaging the SI mean and SD of SARS, mean of 7.6 days and SD of 3.4 days [18], and MERS, mean of 8.4 days and SD of 3.8 days [19].
We were also interested in inferring the patterns of the daily number of cases, denoted by εi for the i-th day, and thus it is obviously that Ci = Ci−1 + εi. A simulation framework was developed for the iterative Poisson process such that E[εi] = Ci−1 × [exp(γ) − 1], where function E[∙] denoted the expectation. The simulation was implemented starting from 1 January 2020 with a cumulative number of cases seed of 40, the same as reported on 31 December 2019. We conducted 1000 samples and calculated the median and 95% CI.

3. Results and Discussion

The number of 2019-nCoV unreported cases was estimated at 469 (95% CI: 403–540), see Figure 1a, which was significantly larger than 0. This finding implied the occurrence of under-reporting between 1 and 15 January 2020. After accounting for the effect of under-reporting, the R0 was estimated at 2.56 (95% CI: 2.49–2.63), see Figure 1b, which is consistent with many existing online preprints with range from 2 to 4 [5,20,21,22]. With the R0 of 2.56 and ξ of 469, the exponential growing framework fitted the cumulative total number of cases (Ci) remarkably well, see Figure 1c, referring to McFadden’s pseudo-R-squared of 0.99.
Our estimation of R0 rely on the SI of 2019-nCoV, which remains unknown as of 26 January 2020. In this work, we employed the SIs of SARS and MERS as approximations to that of 2019-nCoV. The determination of SI requires the knowledge of the chain of disease transmission that needs a sufficient number of patient samples and periods of time for follow-up [23], and thus this is unlikely to be achieved shortly. However, using SIs of SARS and MERS as approximation could provide an insight into the transmission potential of 2019-nCoV at the early outbreak. We note that slightly varying the mean and SD of SI would not affect our main conclusions. The R0 of 2019-nCoV was estimated at 2.56 (95% CI: 2.49–2.63), and it is generally in line with those of SARS, i.e., 2–5 [19,24,25], and MERS, i.e., 2.7–3.9 [26].
For the simulated daily number of cases (εi), see Figure 1d, we found that εi matched the observed daily number after 17 January 2020, but was significantly larger than the observations from 1 to 17 January 2020. This finding implied that under-reporting was likely to have occurred in the first half of January 2020. We estimated that the reporting rate after 17 January 2020 increased 21-fold (95% CI: 18–25) compared to the situation from 1 to 17 January 2020 on average. One of the possible reasons was that the official diagnostic protocol was released by WHO on 17 January 2020 [27], and the diagnosis and reporting efforts of 2019-nCoV infections probably increased. Thereafter, the daily number of newly reported cases started increasing rapidly after 17 January 2020, see Figure 1d. We conducted additional sensitivity analysis by varying the starting date of the under-reporting time window, e.g., 1 January 2020 in the main results, from 2 December 2019 to 3 January 2020, and we report our estimates largely hold. The exact value of the reporting rate was difficult to determine due to lack of serological surveillance data. The reporting rate can be determined if serological surveillance data are available for a population; we would know who was infected (seropositive) and who was not (seronegative), with high confidence. The reporting rate is the ratio of reported cases over the number of seropositive individuals. It was statistically evident that increasing in reporting was likely, and thus it should be considered in the future investigation of this outbreak.
Previous preprint suggested cumulative cases of 1723 (95% CI: 427–4471) as of 12 January 2020, and 4000 (95% CI: 1000–9700) as of 18 January 2020 based on the aggregated international export cases [5]. Our analysis yielded cumulative cases of 280 (95% CI: 128–613) as of 12 January 2020, and 609 (95% CI: 278–1333) as of 18 January 2020 based on the exponential growing mechanistic in the early outbreak. Although our estimate case number appeared to have a lower mean than those estimated by Imai et al. [5], they are not statistically different. This study applied a different screening effort to detect the 2019-nCoV cases from that in Imai et al. [5]. Imai et al. assumed the average screening effort at overseas airports that covered travelers arriving from Wuhan. Whereas we assumed a constant screening effort applied in Wuhan at the same point of time, and then a number of cases (i.e., ξ) should have been reported yet failed to be reported in the first half of January 2020 due to all sorts of reasons. It is not surprising that different assumptions yielded different results, and this difference in screening effort also partly explained why the detected cases out of China mainly presented mild symptoms. Thus, it was reasonable that our estimates appeared lower than those estimated by Imai et al. [5]. It must be emphasized that such a gap in the knowledge would be resolved by serological survey study (for a large population to approximate the actual positive rate) or an explicit estimation of the actual reporting rate.

4. Conclusions

Under-reporting was likely to have occurred and resulted in 469 (95% CI: 403–540) unreported cases from 1 to 15 January 2020. The reporting rate after 17 January 2020 was likely to have increased 21-fold (95% CI: 18–25) compared with the situation from 1 to 17 January 2020 on average, and it should be considered in future investigation. We estimated the R0 at 2019-nCoV to be 2.56 (95% CI: 2.49–2.63).

Author Contributions

All authors conceived the study, carried out the analysis, discussed the results, drafted the first manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the National Natural Science Foundation of China (Grant number 61672013) and the Huaian Key Laboratory for Infectious Diseases Control and Prevention (Grant number HAP201704).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wu, P.; Hao, X.; Lau, E.H.; Wong, J.Y.; Leung, K.S.; Wu, J.T.; Cowling, B.J.; Leung, G.M. Real-time tentative assessment of the epidemiological characteristics of novel coronavirus infections in Wuhan, China, as at 22 January 2020. Eurosurveillance 2020, 25. [Google Scholar] [CrossRef] [PubMed]
  2. Pneumonia of Unknown Cause—China. Available online: https://www.who.int/csr/don/05-january-2020-pneumonia-of-unkown-cause-china/en/ (accessed on 27 January 2020).
  3. Huang, C.; Wang, Y.; Li, X.; Ren, L.; Zhao, J.; Hu, Y.; Zhang, L.; Fan, G.; Xu, J.; Gu, X.; et al. Clinical features of patients infected with 2019 novel coronavirus in Wuhan, China. Lancet 2020. [Google Scholar] [CrossRef] [Green Version]
  4. Situation Report of the Pneumonia Cases Caused by the Novel Coronavirus. Available online: http://www.nhc.gov.cn/yjb/s3578/202001/a3c8b5144067417889d8760254b1a7ca.shtml (accessed on 27 January 2020).
  5. Imai, N.; Dorigatti, I.; Cori, A.; Riley, S.; Ferguson, N.M. Estimating the Potential Total Number of Novel Coronavirus (2019-nCoV) Cases in Wuhan City, China. Available online: https://www.imperial.ac.uk/mrc-global-infectious-disease-analysis/news--wuhan-coronavirus/ (accessed on 27 January 2020).
  6. Nowcasting and Forecasting the Wuhan 2019-nCoV Outbreak. Available online: https://files.sph.hku.hk/download/wuhan_exportation_preprint.pdf (accessed on 27 January 2020).
  7. Bogoch, I.I.; Watts, A.; Thomas-Bachli, A.; Huber, C.; Kraemer, M.U.; Khan, K. Pneumonia of unknown etiology in Wuhan, China: Potential for international spread via commercial air travel. J. Travel Med. 2020. [Google Scholar] [CrossRef] [PubMed]
  8. News Press and Situation Reports of the Pneumonia Caused by Novel Coronavirus. Available online: http://wjw.wuhan.gov.cn/front/web/list2nd/no/710 (accessed on 27 January 2020).
  9. An Outbreak Situation Update on the Pneumonia Caused by the Novel Coronavirus (2019-nCoV) Infection. Available online: http://www.nhc.gov.cn/xcs/yqtb/list_gzbd.shtml (accessed on 27 January 2020).
  10. Definition of Suspected Cases of Unexplained Pneumonia. Available online: http://www.nhc.gov.cn/ (accessed on 27 January 2020).
  11. Zhao, S.; Musa, S.S.; Fu, H.; He, D.; Qin, J. Simple framework for real-time forecast in a data-limited situation: The Zika virus (ZIKV) outbreaks in Brazil from 2015 to 2016 as an example. Parasit. Vectors 2019, 12, 344. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. De Silva, U.; Warachit, J.; Waicharoen, S.; Chittaganpitch, M. A preliminary analysis of the epidemiology of influenza A(H1N1)v virus infection in Thailand from early outbreak data, June–July 2009. Eurosurveillance 2009, 14, 19292. [Google Scholar] [CrossRef] [PubMed]
  13. Fan, J.; Huang, T. Profile likelihood inferences on semiparametric varying-coefficient partially linear models. Bernoulli 2005, 11, 1031–1057. [Google Scholar] [CrossRef]
  14. Wallinga, J.; Lipsitch, M. How generation intervals shape the relationship between growth rates and reproductive numbers. Proc. R. Soc. B Biol. Sci. 2007, 274, 599–604. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Wu, F.; Zhao, S.; Yu, B.; Chen, Y.-M.; Wang, W.; Hu, Y.; Song, Z.-G.; Tao, Z.-W.; Tian, J.-H.; Pei, Y.-Y.; et al. Complete genome characterisation of a novel coronavirus associated with severe human respiratory disease in Wuhan, China. bioRxiv 2020. [Google Scholar] [CrossRef] [Green Version]
  16. Zhou, P.; Yang, X.-L.; Wang, X.-G.; Hu, B.; Zhang, L.; Zhang, W.; Si, H.-R.; Zhu, Y.; Li, B.; Huang, C.-L.; et al. Discovery of a novel coronavirus associated with the recent pneumonia outbreak in humans and its potential bat origin. bioRxiv 2020. [Google Scholar] [CrossRef] [Green Version]
  17. Zhu, N.; Zhang, D.; Wang, W.; Li, X.; Yang, B.; Song, J.; Zhao, X.; Huang, B.; Shi, W.; Lu, R.; et al. A Novel coronavirus from patients with pneumonia in china, 2019. N. Engl. J. Med. 2020. [Google Scholar] [CrossRef] [PubMed]
  18. Assiri, A.; McGeer, A.; Perl, T.M.; Price, C.S.; Al Rabeeah, A.A.; Cummings, D.A.; Alabdullatif, Z.N.; Assad, M.; Almulhim, A.; Makhdoom, H. Hospital outbreak of Middle East respiratory syndrome coronavirus. N. Engl. J. Med. 2013, 369, 407–416. [Google Scholar] [CrossRef] [PubMed]
  19. Lipsitch, M.; Cohen, T.; Cooper, B.; Robins, J.M.; Ma, S.; James, L.; Gopalakrishna, G.; Chew, S.K.; Tan, C.C.; Samore, M.H. Transmission dynamics and control of severe acute respiratory syndrome. Science 2003, 300, 1966–1970. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Riou, J.; Althaus, C.L. Pattern of early human-to-human transmission of Wuhan 2019-nCoV. bioRxiv 2020. [Google Scholar] [CrossRef] [Green Version]
  21. Read, J.M.; Bridgen, J.R.; Cummings, D.A.; Ho, A.; Jewell, C.P. Novel coronavirus 2019-nCoV: Early estimation of epidemiological parameters and epidemic predictions. medRxiv 2020. [Google Scholar] [CrossRef] [Green Version]
  22. Shen, M.; Peng, Z.; Xiao, Y.; Zhang, L. Modelling the epidemic trend of the 2019 novel coronavirus outbreak in China. bioRxiv 2020. [Google Scholar] [CrossRef]
  23. Cowling, B.J.; Fang, V.J.; Riley, S.; Peiris, J.M.; Leung, G.M. Estimation of the serial interval of influenza. Epidemiology 2009, 20, 344. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Wallinga, J.; Teunis, P. Different epidemic curves for severe acute respiratory syndrome reveal similar impacts of control measures. Am. J. Epidemiol 2004, 160, 509–516. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Bauch, C.T.; Lloyd-Smith, J.O.; Coffee, M.P.; Galvani, A.P. Dynamically modeling SARS and other newly emerging respiratory illnesses: Past, present, and future. Epidemiology 2005, 16, 791–801. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Lin, Q.; Chiu, A.P.; Zhao, S.; He, D. Modeling the spread of Middle East respiratory syndrome coronavirus in Saudi Arabia. Stat. Methods Med. Res. 2018, 27, 1968–1978. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Laboratory Testing for 2019 Novel Coronavirus (2019-nCoV) in Suspected Human Cases. Available online: https://www.who.int/health-topics/coronavirus/laboratory-diagnostics-for-novel-coronavirus (accessed on 27 January 2020).
Figure 1. The estimates of the unreported cases between 1 and 15 January 2020, the basic reproduction number (R0), and fitting results of the number of 2019-nCoV cases time series. Panel (a) shows the likelihood profile (, dark green curve) of the estimated number of unreported cases (ξ), and the cutoff threshold (horizontal red dashed line) for the 95% CI. The relationship between the number of unreported cases (ξ) and R0, where the bold curve is the mean estimation, and the dashed curves are the 95% CI of estimated R0. In panels (a,b), the green shading area represents the 95% CI (on the horizontal axis), and the vertical green line represents the maximum likelihood estimate (MLE) of the number of unreported cases. With the MLE of R0 at 2.56, panels (c,d) show the exponential growth fitting results of the cumulative number of cases (Ci) and the daily number of cases (εi) respectively. In panels (c,d), the gold squares are the reported cases, the blue bold curve represents the median of the fitting results, the dashed blue curves are the 95% CI of the fitting results, and the purple shading area represents the time window from 1 to 15 January 2020. In panel (c), the blue dots are the cumulative total, i.e., reported and unreported, number of cases. In panel (d), the grey curves are the 1000 simulation samples.
Figure 1. The estimates of the unreported cases between 1 and 15 January 2020, the basic reproduction number (R0), and fitting results of the number of 2019-nCoV cases time series. Panel (a) shows the likelihood profile (, dark green curve) of the estimated number of unreported cases (ξ), and the cutoff threshold (horizontal red dashed line) for the 95% CI. The relationship between the number of unreported cases (ξ) and R0, where the bold curve is the mean estimation, and the dashed curves are the 95% CI of estimated R0. In panels (a,b), the green shading area represents the 95% CI (on the horizontal axis), and the vertical green line represents the maximum likelihood estimate (MLE) of the number of unreported cases. With the MLE of R0 at 2.56, panels (c,d) show the exponential growth fitting results of the cumulative number of cases (Ci) and the daily number of cases (εi) respectively. In panels (c,d), the gold squares are the reported cases, the blue bold curve represents the median of the fitting results, the dashed blue curves are the 95% CI of the fitting results, and the purple shading area represents the time window from 1 to 15 January 2020. In panel (c), the blue dots are the cumulative total, i.e., reported and unreported, number of cases. In panel (d), the grey curves are the 1000 simulation samples.
Jcm 09 00388 g001

Share and Cite

MDPI and ACS Style

Zhao, S.; Musa, S.S.; Lin, Q.; Ran, J.; Yang, G.; Wang, W.; Lou, Y.; Yang, L.; Gao, D.; He, D.; et al. Estimating the Unreported Number of Novel Coronavirus (2019-nCoV) Cases in China in the First Half of January 2020: A Data-Driven Modelling Analysis of the Early Outbreak. J. Clin. Med. 2020, 9, 388. https://doi.org/10.3390/jcm9020388

AMA Style

Zhao S, Musa SS, Lin Q, Ran J, Yang G, Wang W, Lou Y, Yang L, Gao D, He D, et al. Estimating the Unreported Number of Novel Coronavirus (2019-nCoV) Cases in China in the First Half of January 2020: A Data-Driven Modelling Analysis of the Early Outbreak. Journal of Clinical Medicine. 2020; 9(2):388. https://doi.org/10.3390/jcm9020388

Chicago/Turabian Style

Zhao, Shi, Salihu S. Musa, Qianying Lin, Jinjun Ran, Guangpu Yang, Weiming Wang, Yijun Lou, Lin Yang, Daozhou Gao, Daihai He, and et al. 2020. "Estimating the Unreported Number of Novel Coronavirus (2019-nCoV) Cases in China in the First Half of January 2020: A Data-Driven Modelling Analysis of the Early Outbreak" Journal of Clinical Medicine 9, no. 2: 388. https://doi.org/10.3390/jcm9020388

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop