The Basic Reproduction Number of Novel Coronavirus (2019-Ncov) Estimation Based on Exponential Growth in the Early Outbreak in China from 2019 to 2020: A Reply to Dhungana
Shi Zhao,Qianying Lin,Jinjun Ran,Salihu S Musa,Guangpu Yang,Weiming Wang,Yijun Lou,Daozhou Gao,Lin Yang,Daihai He,Maggie H Wang
DOI: https://doi.org/10.1016/j.ijid.2020.02.025
IF: 12.073
2020-01-01
International Journal of Infectious Diseases
Abstract:The ongoing outbreak of the novel coronavirus (2019-nCoV) pneumonia in Wuhan, China and other regions remains a major public health concern. We thank Dhungana for comments to our study, Zhao et al., 2020Zhao S. Lin Q. Ran J. Musa S.S. Yang G. Wang W. et al.Preliminary estimation of the basic reproduction number of novel coronavirus (2019-nCoV) in China, from 2019 to 2020: A data-driven analysis in the early phase of the outbreak.Int J Infect Dis. 2020; 92: 214-217https://doi.org/10.1016/j.ijid.2020.01.050Abstract Full Text Full Text PDF PubMed Scopus (1133) Google Scholar, recently published in the International Journal of Infectious Diseases. The estimates on the basic reproduction number, R0, were carried out early in the outbreak as of January 22, 2020, when the surveillance data and the knowledge of the key epidemiological features of 2019-nCoV were limited. The assumptions of exponential growth as well as other similar growth patterns are commonly accepted and adopted to capture the growth trends during the early phase of an outbreak (Nishiura and Chowell, 2014Nishiura H. Chowell G. Early transmission dynamics of Ebola virus disease (EVD), West Africa, March to August 2014.Eurosurveillance. 2014; 19: 20894Crossref Google Scholar, Chowell et al., 2004Chowell G. Hengartner N.W. Castillo-Chavez C. Fenimore P.W. Hyman J.M. The basic reproductive number of Ebola and the effects of public health measures: the cases of Congo and Uganda.J. Theor Biol. 2004; 229: 119-126Crossref PubMed Scopus (371) Google Scholar, Wearing et al., 2005Wearing H.J. Rohani P. Keeling M.J. Appropriate models for the management of infectious diseases.PLoS Med. 2005; 2: e174Crossref PubMed Scopus (304) Google Scholar). The exponential growth rate (γ), or the intrinsic growth rate, is estimated from the early epidemic curve and used to calculate the R0. We repeat the analysis Zhao et al., 2020Zhao S. Lin Q. Ran J. Musa S.S. Yang G. Wang W. et al.Preliminary estimation of the basic reproduction number of novel coronavirus (2019-nCoV) in China, from 2019 to 2020: A data-driven analysis in the early phase of the outbreak.Int J Infect Dis. 2020; 92: 214-217https://doi.org/10.1016/j.ijid.2020.01.050Abstract Full Text Full Text PDF PubMed Scopus (1133) Google Scholar, γ is estimated at 0.18 (95%CI: 0.14–0.22), 0.15 (95%CI: 0.12–0.18) and 0.11 (95%CI: 0.09–0.13) per day associated with 2-, 4- and 8-fold increase in the reporting rate, respectively. By using the serial interval (SI) estimate (mean ± SD at 7.5 ± 3.4 days) from Li et al., 2020Li Q. Guan X. Wu P. Wang X. Zhou L. Tong Y. et al.Early Transmission Dynamics in Wuhan, China, of Novel Coronavirus–Infected Pneumonia.N Engl J Med. 2020; Crossref Scopus (10144) Google Scholar, we found the R0 at 3.33 (95%CI: 2.17–4.04), 2.69 (95%CI: 2.28–3.17) and 2.13 (95%CI: 1.88–2.42) associated with 2-, 4- and 8-fold increase in the reporting respectively. Our estimates were in line with the WHO estimates in both the early version (2-fold case) and the published version. The key message as we highlighted in the paper is the changes in the reporting rate. This is recently reconfirmed by Tuite and Fishman (Tuite and Fisman, 2020Tuite A.R. Fisman D.N. Reporting, Epidemic Growth, and Reproduction Numbers for the 2019 Novel Coronavirus (2019-nCoV) Epidemic.Ann Internal Med. 2020; https://doi.org/10.7326/M20-0358Crossref Scopus (81) Google Scholar). We thank the editor and Dhungana to give us this opportunity to reclarify our key message that the reporting rate was not constant during the early outbreak and could affect the estimation of R0. There is indeed a large amount of later confirmed cases that were not counted in the early official daily situation reports Li et al., 2020Li Q. Guan X. Wu P. Wang X. Zhou L. Tong Y. et al.Early Transmission Dynamics in Wuhan, China, of Novel Coronavirus–Infected Pneumonia.N Engl J Med. 2020; Crossref Scopus (10144) Google Scholar; Imai et al., 2020Imai 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. Preprint published by the Imperial College London, 2020https://www.imperial.ac.uk/mrc-global-infectious-disease-analysis/news--wuhan-coronavirus/Google Scholar; Riou and Althaus, 2020Riou J. Althaus C.L. Pattern of early human-to-human transmission of Wuhan 2019 novel coronavirus (2019-nCoV), December 2019 to January 2020.Eurosurveillance. 2020; 25Crossref Scopus (801) Google Scholar; Zhao et al., 2020bZhao S. Musa S.S. Lin Q. Ran J. Yang G. Wang W. 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; 9Crossref Scopus (305) Google Scholar). In other words, if the same reporting standard in the second half of January was applied to the first half of January, the number of cases would be much higher. Other teams either used a retrospective dataset which was not publicly available on January 23, 2020 or used overseas reported cases which were not (to a much less extent) affected by the changes in reporting rate. By using the same analysis and dataset as in Zhao et al., 2020Zhao S. Lin Q. Ran J. Musa S.S. Yang G. Wang W. et al.Preliminary estimation of the basic reproduction number of novel coronavirus (2019-nCoV) in China, from 2019 to 2020: A data-driven analysis in the early phase of the outbreak.Int J Infect Dis. 2020; 92: 214-217https://doi.org/10.1016/j.ijid.2020.01.050Abstract Full Text Full Text PDF PubMed Scopus (1133) Google Scholar, an additional sensitivity analysis on the R0 estimates and varying SI and reporting rate was conducted and is shown in Figure 1. We report that R0 estimates increase while the mean SI increases or the SD of SI decreases. By selecting a mean between 7 and 8 days and SD between 3 and 4 days for SI of 2019-nCoV, the R0 estimates are largely consistent within a range from 2 to 4 in many existing reports (Li et al., 2020Li Q. Guan X. Wu P. Wang X. Zhou L. Tong Y. et al.Early Transmission Dynamics in Wuhan, China, of Novel Coronavirus–Infected Pneumonia.N Engl J Med. 2020; Crossref Scopus (10144) Google Scholar, Imai et al., 2020Imai 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. Preprint published by the Imperial College London, 2020https://www.imperial.ac.uk/mrc-global-infectious-disease-analysis/news--wuhan-coronavirus/Google Scholar, Riou and Althaus, 2020Riou J. Althaus C.L. Pattern of early human-to-human transmission of Wuhan 2019 novel coronavirus (2019-nCoV), December 2019 to January 2020.Eurosurveillance. 2020; 25Crossref Scopus (801) Google Scholar, Zhao et al., 2020bZhao S. Musa S.S. Lin Q. Ran J. Yang G. Wang W. 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; 9Crossref Scopus (305) Google Scholar, Wu et al., 2020Wu J.T. Leung K. Leung G.M. Nowcasting and forecasting the potential domestic and international spread of the 2019-nCoV outbreak originating in Wuhan, China: a modelling study.The Lancet. 2020; 395: 689-697https://doi.org/10.1016/S0140-6736(20)30260-9Abstract Full Text Full Text PDF PubMed Scopus (2798) Google Scholar), see panels (f), (g), (j) and (k) Fig. 1. We conclude that our previous estimation and main conclusions in hold based on the reasonable selection of the SI estimates of 2019-nCoV. Not only is our early version (2-fold case) in line with the WHO estimates, but also we pointed out the issue in the reporting rate changes in the official reported cases. The ethical approval or individual consent was not applicable. All data and materials used in this work were publicly available. Not applicable. WW was supported by National Natural Science Foundation of China (Grant Number 61672013) and Huaian Key Laboratory for Infectious Diseases Control and Prevention (Grant Number HAP201704), Huaian, Jiangsu, China.