key: cord-0880127-drqnrwdl authors: Zhao, Shi; Lin, Qianyin; Ran, Jinjun; Musa, Salihu S.; Yang, Guangpu; Wang, Weiming; Lou, Yijun; Gao, Daozhou; Yang, Lin; He, Daihai; Wang, Maggie H. title: 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 date: 2020-01-30 journal: Int J Infect Dis DOI: 10.1016/j.ijid.2020.01.050 sha: 5a64d8be518b46570d2567a704fc8d6d95bd1008 doc_id: 880127 cord_uid: drqnrwdl BACKGROUNDS: An ongoing outbreak of a novel coronavirus (2019-nCoV) pneumonia hit a major city in China, Wuhan, December 2019 and subsequently reached other provinces/regions of China and other countries. We present estimates of the basic reproduction number, R(0), of 2019-nCoV in the early phase of the outbreak. METHODS: Accounting for the impact of the variations in disease reporting rate, we modelled the epidemic curve of 2019-nCoV cases time series, in mainland China from January 10 to January 24, 2020, through the exponential growth. With the estimated intrinsic growth rate (γ), we estimated R(0) by using the serial intervals (SI) of two other well-known coronavirus diseases, MERS and SARS, as approximations for the true unknown SI. FINDINGS: The early outbreak data largely follows the exponential growth. We estimated that the mean R(0) ranges from 2.24 (95%CI: 1.96–2.55) to 3.58 (95%CI: 2.89–4.39) associated with 8-fold to 2-fold increase in the reporting rate. We demonstrated that changes in reporting rate substantially affect estimates of R(0). CONCLUSION: The mean estimate of R(0) for the 2019-nCoV ranges from 2.24 to 3.58, and is significantly larger than 1. Our findings indicate the potential of 2019-nCoV to cause outbreaks. The atypical pneumonia case, caused by a novel coronavirus (2019-nCoV), was first reported and confirmed in Wuhan, China on December 31, 2019 (World Health Organization, 2020a). As of January 26 (17:00 GMT), 2020, there have been 2033 confirmed cases of 2019-nCoV infections in mainland China, including 56 deaths (National Health Commission of the People's Republic of China, 2020). The 2019-nCoV cases were also reported in Thailand, Japan, Republic of Korea, Hong Kong, Taiwan and the US, and all of these cases were exported from Wuhan; see the World Health Organization (WHO) news release https://www.who.int/csr/don/en/ from January 14-21. The outbreak is still ongoing. A recently published preprint by Imai et al. estimated that a total of 1723 (95% CI: 427-4471) cases of 2019-nCoV infections in Wuhan had onset of symptoms by January 12, 2020 (Imai et al., 2020) . The likelihood of travel related risks of disease spreading is suggested by Bogoch et al. (2020) , which indicates the potentials of regional and global spread (Leung et al., 2020) . To the best of our knowledge, there is no existing peer-reviewed literature quantifying the transmissibility of 2019-nCoV as of January 22, 2020. In this study, we estimated the transmissibility of 2019-nCoV via the basic reproduction number, R 0 , based on the limited data in the early phase of the outbreak. We obtained the number of 2019-nCoV cases time series data in mainland China released by Wuhan Municipal Health Commission, China and National Health Commission of China from January 10 to January 24, 2020 from (Wuhan Municipal Health Commission, China, 2020). All cases were laboratory confirmed following the case definition by the National Health Commission of China (National Health Commission of the People's Republic of China, 2020). Although the date of submission of this study is January 26, we choose to use data up to January 24. Note that the data of the most recent few days contain a number of infections that were infected outside Wuhan due to travel, and thus this part of the infections is excluded from the analysis. Although there were cases confirmed on or before January 16, the official diagnostic protocol was released by WHO on January 17 (World Health Organization, 2020b). To adjust the impact of this event, we considered a time-varying reporting rate that follows a linear increasing trend, motivated by the previous study (Wu et al., 2010) . We assumed that the reporting rate, r(t), started increasing on January 17, and stopped at the maximal level on January 21. The reporting rate increase corresponds to accounts for the announcement on improving the 2019-nCoV surveillance of the Hubei provincial government (Hubei provincial government, 2020). The length of the reporting increasing part roughly equals the average of the incubation periods of two other well-known coronavirus diseases, i.e., the Middle East Respiratory Syndrome (MERS) and the Severe Acute Respiratory Syndrome (SARS), i.e., 5 days (Bauch et al., 2005; Lin et al., 2018; Donnelly et al., 2003) . Denoting the daily reported number of new cases by c(t) for the t-th day, then the adjusted cumulative number of cases, C(t), is CðtÞ ¼ P t t¼0 cðtÞ=rðtÞ. Instead of finding the exact value of r(t), we calculated the fold change in r(t) that is defined by the ratio of r on January 10 over that on January 24 minus 1. We illustrated six scenarios with 0-(no change), 0.5-, 1-, 2-, 4-and 8-fold increase in reporting rate, see Figure 1 (a), (c), (e), (g), (i) and (k). Following previous studies (Zhao et al., 2019; de Silva et al., 2009) , we modelled the epidemic curve obeying the exponential growth. The nonlinear least square (NLS) framework is adopted for data fitting and parameter estimation. The intrinsic growth rate (g) of the exponential growth was estimated, and the basic reproduction number could be obtained by R 0 = 1/M(Àg) with 100% susceptibility for 2019-nCoV at this early stage. The function M(Á) is the Laplace transform, i.e., the moment generating function, of the probability distribution for the serial interval (SI) of the disease (Zhao et al., 2019; Wallinga and Lipsitch, 2007) , denoted by h(k) and k is the mean SI. Since the transmission chain of 2019-nCoV remains unclear, we adopted the SI information from SARS and MERS, which share a similar pathogen as 2019-nCoV. We modelled h(k) as Gamma distributions with a mean of 7.6 days and standard deviation (SD) of 3.4 days for MERS (Assiri et al., 2013) , and mean of 8.4 days and SD of 3.8 days for SARS (Lipsitch et al., 2003) as well as their average, see the row heads in Table 1 for each scenario. The exponential growth fitting results are shown in Figure 1 (b), (d), (f), (h), (j) and (l). The coefficient of determination, R-squared, ranges from 0.91 to 0.92 for all reporting rate changing scenarios, which implies that the early outbreak data were largely following the exponential growth. In Table 1 , we estimated that the R 0 ranges from 2.24 (95%CI: 1.96-2.55) to 5.71 (95%CI: 4.24-7.54) associated with an 8-fold to 0-fold increase in the reporting rate. All R 0 estimates are significantly larger than 1, which indicates the potential of 2019-nCoV to cause outbreaks. Since the official diagnostic protocol was released by WHO on January 17 (World Health Organization, 2020b), an increase in the diagnosis and reporting of 2019-nCoV infections probably occurred. Thereafter, the daily number of newly reported cases started increasing around January 17, see Figure 1 , which implies that more infections were likely being diagnosed and recorded. We suggested that (b) show the scenarios that the reporting rate was unchanged. Panels (c) and (d) show the scenarios that the reporting rate increased by 0.5-fold. Panels (e) and (f) show the scenarios that the reporting rate increased by 1-fold. Panels (g) and (h) show the scenarios that the reporting rate increased by 2-fold. Panels (i) and (j) show the scenarios that the reporting rate increased by 4-fold. Panels (k) and (l) show the scenarios that the reporting rate increased by 8-fold. changes in reporting might exist, and thus it should be considered in the estimation, i.e., 8-, 4-and 2-fold changes are more likely than no change in the reporting efforts. Although six scenarios about the reporting rate were explored in this study, the real situation is difficult to determine given limited data and (almost) equivalent model fitting performance in terms of R-squared. However, with increasing reporting rate, we found the mean R 0 is likely to be between 2 and 3. Our analysis and estimation of R 0 rely on the accuracy of the SI of 2019-nCoV, which remains unknown as of January 25. In this work, we employed the SIs of SARS and MERS as approximations to that of 2019-nCoV. The determination of SI requires knowledge of the chain of disease transmission that needs a sufficient number of patient samples and periods of time for follow-up (Cowling et al., 2009) , and thus this is unlikely to be achieved shortly. However, using SIs of SARS and MERS as approximation could provide an insight to the transmission potential of 2019-nCoV at the early stage of the outbreak. We reported that the mean R 0 of 2019-nCoV is likely to be from 2.24 (8-fold) to 3.58 (2-fold), and it is largely in the range of those of SARS, i.e., 2-5 (Bauch et al., 2005; Lipsitch et al., 2003; Wallinga and Teunis, 2004) , and MERS, i.e., 2.7-3.9 (Lin et al., 2018) . We note that WHO reported the basic reproduction number for the human-to-human (direct) transmission ranged from 1.4 to 2.5 (World Health Organization, 2020b), which is marginally lower than ours. However, many of the existing online preprints estimate the mean R 0 ranging from 2 to 5 (Imai et al., 2020; Riou and Althaus, 2020; Read et al., 2020; Shen et al., 2020) , which is largely consistent with our results. We estimated the mean R 0 of 2019-nCoV ranging from 2.24 (95%CI: 1.96-2.55) to 3.58 (95%CI: 2.89-4.39) if the reporting effort has been increased by a factor of between 8-and 2-fold, respectively, after the diagnostic protocol released on January 17, 2020 and many medical supplies reached Wuhan. The ethical approval or individual consent was not applicable. All data and materials used in this work were publicly available. Not applicable. The funding agencies had no role in the design and conduct of the study; collection, management, analysis, and interpretation of the data; preparation, review, or approval of the manuscript; or decision to submit the manuscript for publication. Table 1 The summary table of the estimated basic reproduction number, R 0 , under different scenarios. The estimated R 0 is shown as in the 'median (95%CI)' format. The 'reporting rate increased' indicates the number of fold increase in the reporting rate from January 17, when WHO released the official diagnostic protocol (World Health Organization, 2020b), to January 20, 2020. Hospital outbreak of Middle East respiratory syndrome coronavirus Dynamically modeling SARS and other newly emerging respiratory illnesses: past, present, and future Pneumonia of unknown etiology in Wuhan, China: potential for international spread via commercial air travel Estimation of the serial interval of influenza A preliminary analysis of the epidemiology of influenza A (H1N1) v virus infection in Thailand from early outbreak data Epidemiological determinants of spread of causal agent of severe acute respiratory syndrome in Hong Kong Strengthening the new coronavirus infection of pneumonia prevention and control. Hubei provincial government Estimating the potential total number of novel Coronavirus (2019-nCoV) cases in Wuhan City. China: Preprint published by the Imperial College London Nowcasting and forecasting the Wuhan 2019-nCoV outbreak. Preprint published by the School of Public Health of the University of Hong Kong Modeling the spread of Middle East respiratory syndrome coronavirus in Saudi Arabia Transmission dynamics and control of severe acute respiratory syndrome Definition of suspected cases of unexplained pneumonia', the National Health Commission of the People's Republic of China Novel coronavirus 2019-nCoV: early estimation of epidemiological parameters and epidemic predictions Pattern of early human-to-human transmission of Wuhan 2019-nCoV Modelling the epidemic trend of the 2019 novel coronavirus outbreak in China How generation intervals shape the relationship between growth rates and reproductive numbers Different epidemic curves for severe acute respiratory syndrome reveal similar impacts of control measures Pneumonia of unknown cause -China', Emergencies preparedness, response, Disease outbreak news, World Health Organization (WHO) World Health Organization. Laboratory testing for 2019 novel coronavirus (2019-nCoV) in suspected human cases, World Health Organization (WHO) School closure and mitigation of pandemic (H1N1) 2009, Hong Kong News press and situation reports of the pneumonia caused by novel coronavirus Simple framework for real-time forecast in a datalimited situation: the Zika virus (ZIKV) outbreaks in Brazil from 2015 to 2016 as an example The authors declare no competing interests. All authors conceived the study, carried out the analysis, discussed the results, drafted the first manuscript, critically read and revised the manuscript, and gave final approval for publication. The authors would like to acknowledge anonymous colleagues for helpful comments.