COVID-19 Epidemic Outside China: 34 Founders and Exponential Growth =================================================================== * Yi Li * Meng Liang * Xianhong Yin * Xiaoyu Liu * Meng Hao * Zixin Hu * Yi Wang * Li Jin ## Summary **Background** In December 2019, pneumonia infected with a novel coronavirus burst in Wuhan, China. Now the situation is almost controlled in China but is worse outside China. We aimed to build a mathematical model to capture the global trend of epidemics outside China. **Methods** In this retrospective, outside-China diagnosis number reported from Jan 21 to Feb 28, 2020 was downloaded from WHO website. We develop a simple regression model on these numbers: ![Formula][1] where Nt is the total diagnosed patient till the *i*th day, t=1 at Feb 1. **Findings** Based on this model, we estimate that there have been about 34 unobserved founder patients at the beginning of spread outside China. The global trend is approximately exponential, with the rate of 10 folds every 19 days. **Evidence before this study** In December 2019, pneumonia infected with a novel coronavirus burst in Wuhan, China. Now the situation is almost controlled in China but is worse outside China. Now there are 4,691 patients across 51 countries and territories outside China. We searched PubMed and the China National Knowledge Infrastructure database for articles published up to Feb 28, 2020, using the keywords “COVID”, “novel coronavirus”, “2019-nCoV” or “2019 novel coronavirus”. No published work about the global trend of epidemics outside China could be identified. **Added value of this study** We built a simple “log-plus” linear model to capture the global trend of epidemics outside China. We estimate that there have been about 34 unobserved founder patients at the beginning of spread outside China. The global trend is approximately exponential, with the rate of 10 folds every 19 days. **Implications of all the available evidence** With the limited number of data points and the complexity of the real situation, a straightforward model is expected to work better. Our model suggests that the COVID-19 disease follows an approximate exponential growth model stably at the very beginning. We predict that the number of confirmed patients outside China will increase ten folds in every 19 days without strong intervention by applying our model. Powerful actions on public health should be taken to combat this epidemic all over the world. ## Introduction In early December of 2019, pneumonia cases of unknown cause emerged in Wuhan, the capital of Hubei province, China.1 A novel coronavirus (now named as COVID-19) was verified and identified as the seventh member of enveloped RNA coronavirus (subgenus sarbecovirus, Orthocoronavirinae subfamily) through high-throughput sequencing.2-4 Human-to-human transmission in hospital and family settings had been accumulating 5-7 and has occurred among close contacts since the middle of December 2019 based on evidence from early transmission dynamics.8 According to WHO statistics, the accumulating number of diagnosed patients in China was 78,961 in Feb 28, 2020.9 COVID-19 raised intension not only within China but internationally. Since the first case of pneumonia caused by COVID-19 was reported from Wuhan, COVID-19 patients had rapidly been diagnosed in other cities of China and the neighboring countries including Thailand, South Korea, Japan, and even to a few Western countries.10-12 On January13 2020, the Ministry of Public Health of Thailand reported the first imported case of lab-confirmed novel coronavirus (2019-nCoV). 13 As of Feb28 2020, cases had been reported in 52 countries and territories. The number of confirmed COVID-19 cases had reached 4,691 outside China, 9 and the number is still increasing. South Korea, the worst affected country outside China, has reported another jump. The 571 new cases brought South Korea’s total confirmed cases to 2,337. A surge in cases of COVID-19 in Italy, Japan and Iran also heightened fears that the world is on the brink of a pandemic.14 In Feb 28, WHO increased the assessment of the risk of spread and risk of impact of COVID-19 to very high at the global level. Recently, many researchers devoted to more details about the analysis of the spread of epidemic.14,15 Some paralleled results have shown that the reproductive number R0 estimated of COVID-19 is bigger than that of SARS based on different models. 16-18 With the limited number of data points and the complexity of the real situation, a straightforward model is expected to work better (see Discussion). In this study, we propose a “log-plus” model for the situation prediction which only requires daily number of total diagnosis outside China. This model assumes that there were some unobserved founder patients at the beginning of spread outside China and exponential growth later. Despite the simplicity of our model, it fits the data well (R2=0.991). The prediction is expected to provide practical significance on social application and evidence for enhancing public health interventions to avoid severe outbreaks. ## Methods ### Data Daily diagnosis number of countries outside China is download from WHO situation reports ([https://www.who.int/emergencies/diseases/novel-coronavirus-2019/situation-reports](https://www.who.int/emergencies/diseases/novel-coronavirus-2019/situation-reports)). The data in this analysis starts from Jan 21 and currently it ends at Feb 28. ### Model We first explored the data by plotting log transformed daily number. We observed a good linear trend recently, but a relatively poor fit at the beginning of spread. We therefore consider that there are some unobserved founder patients at the beginning. Based on exploratory data analysis and mathematical intuition we proposed the following model: ![Formula][2] where Nt is the number of diagnosed patients outside China according to WHO at the *t*th day, t=1 at Feb 1; u is the number of unobserved founder patients at the beginning of spread outside China; a and b are simple linear regression parameters. We enumerated u from 0 to 100 with step size 1. For each u, we calculated Pearson’s correlation R2 between t and *log*10 (*N**t* +*u*). We selected the *û* that maximize R2 and estimate corresponding *â* and ![Graphic][3] using a simple linear regression between t and *log*10 (*N**t* +*û*). ### Availability of source code The source code of model estimation is at: [https://github.com/wangyi-fudan/COVID-19\_Global\_Model](https://github.com/wangyi-fudan/COVID-19_Global_Model) ## Results ### Data Table The WHO daily count of diagnosis number outside China and its “log-plus” transform as well as model fit were shown in table 1. View this table: [Table 1:](http://medrxiv.org/content/early/2020/03/05/2020.03.01.20029819/T1) Table 1: WHO Daily Diagnosis Number Outside China ### Parameter Estimation *û, â* and ![Graphic][4] were estimated as 34, 0·0515 and 2·075 respectively according to Feb 28 data (figure 1) ![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/03/05/2020.03.01.20029819/F1.medium.gif) [Figure 1:](http://medrxiv.org/content/early/2020/03/05/2020.03.01.20029819/F1) Figure 1: Estimation of u Parameter by Enumeration ### Global Trend Model We plot *log*10 (*N**t* + 34) against time to visualize the model fitting (figure 2). The model R2 is 0·991 indicating an excellent model fit. ![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/03/05/2020.03.01.20029819/F2.medium.gif) [Figure 2:](http://medrxiv.org/content/early/2020/03/05/2020.03.01.20029819/F2) Figure 2: The Exponential Growth of COVID-19 Infection Outside China ### Future Number Prediction Currently the number of diagnoses at Febuary 28 is 4,691. We predict that the number of diagnoses outside China will expand exponentially with the rate of 10 folds every 19 days without strong intervention. ### Conclusion Based on this model, we estimated that there have been about 34 unobserved founder patients at the beginning of spread outside China. Our model suggests that the COVID-19 disease follows an approximate exponential growth model stably at the very beginning. The situation is dangerous as we expect 10 folds increases of patients every 19 days without strong intervention. We call on world-wide strong actions on public health referring to experiences learnt from China and Singapore. ## Discussion In this report, we only analyzed the total number of diagnoses outside China. The country scale data is available but is not as complete as total number. Here, we limited this analysis to capture the global trend. This model is a minimal extension of the ‘default’ exponential growth model with the estimation of 34 unobserved founder patients outside China. An almost perfect model fitting (R2=0·991) indicates that the spread of disease does follow our model. A simple and straightforward linear model inherits with some advantages: (1) it works on small sample size due to limited observation or somewhat unperfect data. (2) it is relatively robust in complex situations. The virus spreading pattern is complex and yet varies across the world. A simple model simplified the situation and provides coarse-grained trend estimation. (3) A linear model has nice extrapolation ability comparing to complex models (e.g. Neural network). The existent of 34 unobserved founder patients is not surprising. They may have mild symptoms thus did not went to the hospitals. However, we did not preclude the possibility that they already existed before or parallel with Wuhan’s burst. ## Data Availability Daily diagnosis number of countries outside China is download from WHO situation reports ([https://www.who.int/emergencies/diseases/novel-coronavirus-2019/situation-reports](https://www.who.int/emergencies/diseases/novel-coronavirus-2019/situation-reports)). [https://www.who.int/emergencies/diseases/novel-coronavirus-2019/situation-reports](https://www.who.int/emergencies/diseases/novel-coronavirus-2019/situation-reports) ## Contributors YW conceived the idea and wrote the source code. YW, YL, ML and LJ contributed data analysis, generating tables and figures, and manuscript writing. YL, ML, XY, LX, MH, ZH, YW, LJ contributed the theoretical analysis and manuscript revision. All authors contributed to final revision of the manuscript. ## Declaration of interests The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest. ## Acknowledgments We thank the Fudan University High-End Computing Center for supporting computations involved in this study. * Received March 1, 2020. * Revision received March 3, 2020. * Accepted March 5, 2020. * © 2020, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/) ## References 1. 1.Huang C, Wang Y, Li X, et al. Clinical features of patients with 2019 novel coronavirus in Wuhan, China. Lancet. 2020; doi: 10.1016/S0140-6736(20)30183-5 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(20)30183-5&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31986264&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F05%2F2020.03.01.20029819.atom) 2. 2.Lu R, Zhao X, Li J, et al. Genomic characterization and epidemiology of 2019 novel coronavirus: implications of virus origins and receptor binding. Lancet. 2020; doi: 10.1016/S0140-6736(20)30251-8 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(20)30251-8&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32007145&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F05%2F2020.03.01.20029819.atom) 3. 3.[https://www.who.int/dg/speeches/detail/who-director-general-s-remarks-at-the-media-briefing-on-2019-ncov-on-11-february-2020](https://www.who.int/dg/speeches/detail/who-director-general-s-remarks-at-the-media-briefing-on-2019-ncov-on-11-february-2020) 4. 4.Zhu N, Zhang D, Wang W, et al. A novel coronavirus from patients with pneumonia in China, 2019. N Engl J Med. 2020; doi: 10.1056/NEJMoa2001017 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa2001017&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31978945&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F05%2F2020.03.01.20029819.atom) 5. 5.Chan JF, Yuan S, Kok KH, et al. A familial cluster of pneumonia associated with the 2019 novel coronavirus indicating person-to-person transmission: a study of a family cluster. Lancet. 2020; doi: 10.1016/S0140-6736(20)30154-9 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(20)30154-9&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31986261&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F05%2F2020.03.01.20029819.atom) 6. 6.Phan LT, Nguyen TV, Luong QC, et al. Importation and human-to-human transmission of a novel coronavirus in Vietnam. N Engl J Med. 2020; doi: 10.1056/NEJMc2001272 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMc2001272&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31991079&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F05%2F2020.03.01.20029819.atom) 7. 7.Wu JT, Leung K, Leung GM. Nowcasting and forecasting the potential domestic and international spread of the 2019-nCoV outbreak originating in Wuhan, China: A modeling study. Lancet. 2020; 8. 8.Li Q, Guan X, Wu P, et al. Early Transmission Dynamics in Wuhan, China, of Novel Coronavirus-Infected Pneumonia. N Engl J Med 2020; published online. doi: 10.1056/NEJMoa2001316 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa2001316&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31995857&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F03%2F05%2F2020.03.01.20029819.atom) 9. 9.[https://www.who.int/docs/default-source/coronaviruse/situation-reports/20200226-sitrep-37-covid-19.pdf?sfvrsn=2146841e\_2](https://www.who.int/docs/default-source/coronaviruse/situation-reports/20200226-sitrep-37-covid-19.pdf?sfvrsn=2146841e_2) 10. 10.Kim H. Outbreak of novel coronavirus (COVID-19): What is the role of radiologists? 2020. 11. 11.Shao N, Pan H, Li X, et al. CoVID-19 in Japan: What could happen in the future? medRxiv, 2020. 12. 12.Ryu S, Ali S T, Lim J, et al. The estimate of infected individuals of the 2019-Novel Coronavirus in South Korea by incoming international students from the countries of risk of 2019-Novel Coronavirus: a simulation study. medRxiv, 2020. 13. 13.[https://www.who.int/docs/default-source/coronaviruse/situation-reports/20200227-sitrep-38-covid-19.pdf?sfvrsn=9f98940c\_2](https://www.who.int/docs/default-source/coronaviruse/situation-reports/20200227-sitrep-38-covid-19.pdf?sfvrsn=9f98940c_2) 14. 14.Shao N, Cheng J, Chen W. The reproductive number R0 of COVID-19 based on estimate of a statistical time delay dynamical system. medRxiv 2020; doi:10.1101/2020.02.17.20023747 (preprint). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMC4wMi4xNy4yMDAyMzc0N3YyIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDMvMDUvMjAyMC4wMy4wMS4yMDAyOTgxOS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 15. 15.Zhao Q, Chen Y, Small D S. Analysis of the epidemic growth of the early 2019-nCoV outbreak using internationally confirmed cases. medRxiv, 2020. 16. 16.Xiong H, Yan H. Simulating the infected population and spread trend of 2019-nCov under different policy by EIR model. Available at SSRN 3537083, 2020. 17. 17.Mizumoto K, Kagaya K, Chowell G. Early epidemiological assessment of the transmission potential and virulence of 2019 Novel Coronavirus in Wuhan City: China, 2019-2020. medRxiv, 2020. 18. 18.Zhao S. 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. International Journal of Infectious Diseases, 2020. [1]: /embed/graphic-1.gif [2]: /embed/graphic-2.gif [3]: /embed/inline-graphic-1.gif [4]: /embed/inline-graphic-2.gif