Simple model for Covid-19 epidemics – back-casting in China and forecasting in the US ======================================================================================= * Slav W. Hermanowicz ## Abstract and Findings In our previous work, we analyze, in near-real time, evolution of Covid-19 epidemic in China for the first 22 days of reliable data (up to February 6, 2020). In this work, we used the data for the whole 87 days (up to March 13, 2020) in China and the US data available till March 31 (day 70) for systematic evaluation of the logistic model to predict epidemic growth. We sequentially estimated sets of model parameters (maximum number of cases *K*, growth rate *r*, and half-time *t***) and the epidemic “end time” *t**95* (defined as the time when the number of cases, predicted or actual, reached 95% of the maximum). The estimates of these parameters were done for sequences of reported cases growing daily (back-casting for China and forecasting for the US). In both countries, the estimates of *K* grew very much in time during the exponential and nearly exponential phases making longer term forecasting not reliable. For the US, the current estimate of the maximum number of cases *K* is about 265,000 but it is very likely that it will grow in the future. However, running estimates of the “end time” *t**95* were in a much smaller interval for China (60 – 70 days *vs*. the actual value of 67). For the US, the values estimated from the data sequences going back two weeks from now range from 70 to 80 days. If the behavior of the US epidemic is similar to the previous Chinese development, the **number of reported cases** could reach a **maximum around April 10 to 14**. ## Introduction In our previous recent work (Hermanowicz 2020), we used a simple logistic model to analyze the evolution of data on Covid-19 cases as reported in mainland China by the National Health Commission of the People’s Republic of China (NHC 2020). This initial analysis was done in three phases in near-real time (up to Jan. 30, Feb. 3, and Feb 6). The analysis resulted in a sequence of continually updated forecasts of the epidemic dynamics. Although the predicted maximum numbers of cases increased as new data became available, they systematically underestimated final reported numbers (even without accounting for a major change in reporting criteria in China on February 12 resulting in a large surge of new reported cases). However, the successive **predicted dynamics** of the epidemic were **remarkably close** to the final real-world outcome. We explore this issue further in this work where we use the full reported dataset for mainland China for **systematic back-casting** of the epidemic. At the same time, when this work is being conducted, in the US we are experiencing a fast developing Covid-19 epidemic that in some aspects is similar to China but with a delay. We use the experience from our analysis of the situation in China to **forecast** further evolution of coronavirus cases in the United States. Obviously, our predictions are based on the data currently available. These predictions do not account for any possible other secondary sources of infection, changes of diagnostics or reporting, or virus mutation. Modelling the epidemics and infection dynamics is very important and numerous results have been recently reported for worldwide Covid-19 outbreak (Chen, Rui et al. 2020, Chen, Cheng et al. 2020, Hong, He et al. 2020, Hui, Azhar et al. 2020, Imai, Cori et al. 2020, Imai, Dorigatti et al. 2020, Li, Guan et al. 2020, Linton, Kobayashi et al. 2020, Liu, Hu et al. 2020, Liu, Hewings et al. 2020, Roosa, Lee et al. 2020, Shen, Peng et al. 2020, Wu, Leung et al. 2020, Zhang and Wang 2020, Zhao, Lin et al. 2020, Zhou, Hong et al. 2020, Ziff and Ziff 2020). Many reported models are complicated, incorporate tentative assumptions and need parameters estimates that are not reliable as underscored by Fong and coworkers (Fong, Li et al. 2020). In this work, we present the results of fitting a very simple logistic model to the available data and a forecast of new infections. In contrast to other models, the logistic model does not include any external assumptions and is derived completely from available data. ### Logistic Model and its Application The logistic model is one of the simplest used in population dynamics and has been used specifically for epidemics for a long time (Bailey 1950, Cockburn 1960, Mansfield and Hensley 1960, Jowett, Browning et al. 1974, Waggoner and Aylor 2000, Koopman 2004, Bangert, Molyneux et al. 2017). In our earlier analysis of the developing epidemic in China, we used a discrete version of the model due to its uncomplicated structure and easy calculations. In discrete time, more appropriate to daily reported infection cases, the logistic model becomes ![Formula][1] where *P*(*t*) and *P*(*t*+*1*) are populations (cases) on consecutive days, *R**** is the growth rate (basic reproduction number in epidemiology) at the beginning of the logistic growth, and *K* is the limiting population (maximum cases). However, expressing the growth of population *P* in continuous time *t* allows to formulate the model as an ordinary first-order differential equation describing dynamic evolution of the population *P* (in our case the number of infected individuals) being controlled by the growth rate *r* and maximum cases *K* due to limited growth. In continuous time *t*, the change of *P* is ![Formula][2] Initially, the growth of *P* is close to exponential since the term (1− *P*/*K*)is almost one. When *P* becomes larger (commensurate with *K*) the growth rate slows down with ![Formula][3] becoming an effective instantaneous growth rate. The solution to Eq. (2) is the well-known sigmoidal function (logistic function) ![Formula][4] where *t*** is the time when the population reaches one-half of the maximum value *P*(*t***) *= 1/2 K*. Using a differential version of the model is more convenient since a closed-form solution exists and allows for direct estimation of three model parameters: *K, r*, and *t***. The logistic model may be adequate for the analysis of mainland China and the US as a whole since at this time each country can be treated as a unit where a vast majority of cases occurred without any ***significant*** “import” or “export” of cases. ### Data Analysis For China, we used data reported daily by the National Health Commission of the People’s Republic of China (NHC 2020) up to March 13, 2020 (day 87 from the outbreak) when only 11 new cases were reported – less than 2·10−4 of total case number, effectively ending the epidemic on the national level. There is a considerable controversy as to the exact date of the outbreak with most reports pointing to mid-December (Li, Guan et al. 2020, Wang, Horby et al. 2020) while one analysis suggest multiple sources of original infection (Nishiura, Jung et al. 2020). Initially, the outbreak was not recognized and number of ***confirmed*** cases is not fully known (Wu, Hao et al. 2020). In our analysis, we adopted December 17, 2019 as the best estimate of the outbreak following the work of Zhang and co-workers (Zhang and Wang 2020, Zhang and Wang 2020). In addition, in the initial stages of the epidemic the reported numbers of cases may severely underestimate the actual numbers due to asymptomatic carriers (Zhao, Musa et al. 2020). More accurate numbers can be only estimated after the epidemic (Wu and McGoogan 2020). The data from China are show in Figure 1 and in the Appendix. As seen in this figure, the cumulative number of cases grew in a sigmoidal fashion. On February 12, 2020 (day 57) the reporting criteria were changed resulting in a one-day increase of about 15,000 cases. We analyzed the subsequent data with this jump and also ignoring it. The logistic model turned out to be quite robust. While the numerical estimates of the maximum number of cases *K* was significantly affected by the inclusion or exclusion the jump, the dynamics of the model was much less impacted (see further discussion). In this work, we decided to report the results with the jump as parsing data of new and cumulative cases past day 57 is very unreliable. ![Figure 1](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/03/2020.03.31.20049486/F1.medium.gif) [Figure 1](http://medrxiv.org/content/early/2020/04/03/2020.03.31.20049486/F1) Figure 1 Evolution of reported cases in China and the US. (Solid points show data used to fit exponential functions – dashed lines) In the US, there are numerous agencies tracking the Covid-19 cases at different administrative levels in accordance with multi-tiered structure of the US government. For this reason, we relied on web-based data aggregators (mostly at state level) following daily updates (Covid-19 2020). We used January 21, 2020 as the start of the epidemic in the US (day 1) since on that day the first non-repatriated Covid-19 case was reported. It should be noted that like in China the number of cases in the initial stages is underreported due to the same factors (asymptomatic carriers, initial administrative confusion, limited testing). Despite these problems, the data used in this work are the best currently available. The current evolution of the epidemic in the US is presented in Figure 1 and in the Appendix. Despite smaller population, the number of cases in the US at the time of writing already exceeds that of China. Reasons for such higher numbers are outside the scope of this analysis. Both, China and the US data are also plotted in the semi-log format (??) to underscore exponential growth of the epidemic in the initial stages – straight line on these plots. In China, exponential growth occurred from day 30 till day 42 as shown in ?? and ?? by the dashed lines. After day 42 (January 28), the number of cases still increased but the rate of growth was becoming lower and the line representing the cumulative cases deviated from the exponential curve. This feature was also clearly recognized in a previous report (Zhao, Lin et al. 2020). In the US, exponential growth occurred for a longer period and only very recently (March 27 – 30) **perhaps** starts to deviate from the exponential curve. ## Model Estimates and Results For each dataset of reported cumulative case numbers (China and the US), we estimated three parameters of the logistic model (maximum case number *K*, growth rate *r*, half-time *t***) fitting model predictions to the data. We use a custom nonlinear curve fitting procedure employing the Levenberg-Marquardt method for minimization of the residual sum of squares. Similarly to our previous work (Hermanowicz 2020), we estimated model parameters sequentially from datasets growing day after day. For China, the first dataset contained 5 days from day 38 through 42. The next estimate used 6 days from day 38 through 43. This process was repeated until day 87 when the entire China dataset was used. For the US, the first dataset contained also 5 days from day 43 through 47. The last day of the available sequence for the US was day 70 (March 30, 2020). All resulting estimates are show in the Appendix. In case of China, where the epidemic growth has essentially ended, the sequential estimation process (back-casting) simulated near-real time analysis of the dynamics. In the US, where we are still in the substantial growth phase, the sequential estimates are indeed performed in near-real time. In addition to three model parameters, for each day we also estimated predicted time for the epidemic to end. For this purpose. we chose arbitrarily time when the predicted number of cases reach 95% of the predicted maximum *K*. This time, *t**95* was calculated from Eq. (4) by setting *P*(*t*) *= 0*.*95 K* and is also shown in the Appendix. ### China Figure 3 shows the development of sequential estimates of the maximum predicted cases *K* for China and Figure 4 presents the corresponding “ending times” *t**95*. As reported earlier (Hermanowicz 2020) for actual near-real time analysis, estimates of predicted maximum cases *K* depended very heavily on the length of the dataset used for estimation. As examples, Figure 5 shows a few logistic curves corresponding to selected model parameter estimates on specific days for China. As seen in this figure, the initial estimates (close to the exponential phase) were below 20,000 but they grew to about 100,000 as more data became available and was used for estimate refinement. Obviously, the estimate of *K* obtained from the full dataset (up to day 87) matches closely the actual reported number of cases (80,780 *vs*. 80,807) but it should be noted that the estimates of *K* from day 65 (more than 20 days in advance) converged very closely to the actual maximum. ![Figure 2](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/03/2020.03.31.20049486/F2.medium.gif) [Figure 2](http://medrxiv.org/content/early/2020/04/03/2020.03.31.20049486/F2) Figure 2 Exponential growth of epidemics (see Figure 1 for notes) ![Figure 3](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/03/2020.03.31.20049486/F3.medium.gif) [Figure 3](http://medrxiv.org/content/early/2020/04/03/2020.03.31.20049486/F3) Figure 3 Sequential estimates of K (solid data points show days for which logistic model predictions are shown in Figure 5 for China and Figure 6 for the US) ![Figure 4](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/03/2020.03.31.20049486/F4.medium.gif) [Figure 4](http://medrxiv.org/content/early/2020/04/03/2020.03.31.20049486/F4) Figure 4 Sequential estimates of “end time” t95 (see notes in Figure 3) ![Figure 5](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/03/2020.03.31.20049486/F5.medium.gif) [Figure 5](http://medrxiv.org/content/early/2020/04/03/2020.03.31.20049486/F5) Figure 5 Examples of model predictions for China for selected parameter estimates on specific days Due to a change in reporting criteria in China on day 57 (February 12), there was a major increase in the number of reported cases that could not be incorporated in the model. As mentioned before, we decided to use all data (with the jump) in any subsequent analysis. The immediate result was a very large increase in the estimated maximum cases *K* (Figure 3), approximately doubling it from about 50,000 on day 56 to about 99,400 on day 59. This big increase underscores again the sensitivity of the logistic model to data quality. However, the model is also robust in a longer term as the *K* estimates quickly converged toward their ultimate values. Similarly, there was also a significant jump in the estimates of *t**95* (see Figure 4) due to the jump in the case numbers. However, it is remarkable that the estimates of the “end time” *t**95* were much more constrained and much closer to the actual value. Even as early as three weeks before the end of the epidemic the estimates of the “end time” were between 60 and 70 days, very close to the final value of 67 days. ### United States Unlike China, the epidemic in the US is still at the growth phase, perhaps deviating slightly from the exponential growth. Thus, the available dataset is much smaller and the logistic model estimates are burdened with much larger uncertainty. Currently available estimates of the maximum cases *K* are also shown in Figure 3. They exhibit very large variations increasing from approximately 1,100 on day 47 to approximately 265,000 on day 70 without any sign of leveling off. This large variation is not unexpected since the nearly exponential growth does not contain sufficient information on the actual maximum. In other words, the derivative *dP*/*dt* in Eq. (2) is dominated by the term *r P* while (*1*-*P*/*K*) ≈ *1*. This behavior is also seen in the examples of predicted logistic curves for selected sets of estimated parameters *K, r, t*** (Figure 6). ![Figure 6](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/03/2020.03.31.20049486/F6.medium.gif) [Figure 6](http://medrxiv.org/content/early/2020/04/03/2020.03.31.20049486/F6) Figure 6 Examples of model predictions for the US for selected parameter estimates on specific days The sequential estimates of the “end time” *t**95* are plotted in Figure 4. They also tend to increase in time with the increasing number of reported cases but unlike the estimates of *K* (and similarly to the Chinese case) it varies in a much smaller interval – 70 to 80 days in the past two weeks, If the behavior of the US epidemic is similar to the Chinese case, we could expect further leveling off of *t**95* at slightly above 80 days. If this **bold prediction** will hold, we could see the end of the epidemic growth around 80 – 85 days after its origin, around April 10 to 14. Of course, the end here is defined as the cessation of new cases and not the complete recovery of the infected patients. ## Data Availability data in the paper ## Declarations ### Ethics approval and consent to participate The ethical approval or individual consent was not applicable. ### Availability of data and materials All data and materials used in this work were publicly available. ### Consent for publication Not applicable. ### Funding This work was not funded. ### Disclaimer 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. ### Conflict of Interests The author declared no competing interests. ## Appendix View this table: [Table A- 1](http://medrxiv.org/content/early/2020/04/03/2020.03.31.20049486/T1) Table A- 1 Number of cases in China View this table: [Table A- 2](http://medrxiv.org/content/early/2020/04/03/2020.03.31.20049486/T2) Table A- 2 Number of cases in the US View this table: [Table A- 3](http://medrxiv.org/content/early/2020/04/03/2020.03.31.20049486/T3) Table A- 3 Sequential estimates of logistic model for China View this table: [Table A- 4](http://medrxiv.org/content/early/2020/04/03/2020.03.31.20049486/T4) Table A- 4 Sequential estimates of logistic model for the US * Received March 31, 2020. * Revision received March 31, 2020. * Accepted April 3, 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. Bailey, N. T. J. (1950). “A SIMPLE STOCHASTIC EPIDEMIC.” Biometrika 37(3-4): 193–202. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/biomet/37.3-4.193&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=14801045&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F04%2F03%2F2020.03.31.20049486.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1950UZ95800001&link_type=ISI) 2. Bangert, M., D. H. Molyneux, S. W. Lindsay, C. Fitzpatrick and D. Engels (2017). “The cross-cutting contribution of the end of neglected tropical diseases to the sustainable development goals.” Infectious Diseases of Poverty 6. 3. Chen, T., J. Rui, Q. Wang, Z. Zhao, J.-A. Cui and L. Yin (2020) “A mathematical model for simulating the transmission of Wuhan novel Coronavirus.” bioRxiv, 2020.2001.2019.911669 DOI: 10.1101/2020.01.19.911669. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiYmlvcnhpdiI7czo1OiJyZXNpZCI7czoxOToiMjAyMC4wMS4xOS45MTE2Njl2MSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzA0LzAzLzIwMjAuMDMuMzEuMjAwNDk0ODYuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 4. Chen, Y., J. Cheng, Y. Jiang and K. J. Liu (2020). “A time delay dynamic system with external source for the local outbreak of 2019-nCoV.” Applicable Analysis. 5. Cockburn, T. A. (1960). “EPIDEMIC CRISIS IN EAST PAKISTAN - APRIL-JULY, 1958.” Public Health Reports 75(1): 26–36. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2307/4590706&link_type=DOI) 6. Covid-19. (2020). “Template:2019–20 coronavirus pandemic data/United States medical cases.” Retrieved March 30. 2020, 2020, from [https://en.wikipedia.org/wiki/Template:2019%E2%80%9320\_coronavirus\_pandemic\_data/United\_States\_medical\_cases](https://en.wikipedia.org/wiki/Template:2019%E2%80%9320\_coronavirus\_pandemic\_data/United_States_medical_cases). 7. Fong, S. J., G. Li, N. Dey, R. G. Crespo and E. Herrera-Viedma (2020). “Finding an Accurate Early Forecasting Model from Small Dataset: A Case of 2019-nCoV Novel Coronavirus Outbreak.” International Journal of Interactive Multimedia and Artificial Intelligence 6(1): 132–140. 8. Hermanowicz, S. W. (2020) “Forecasting the Wuhan coronavirus (2019-nCoV) epidemics using a simple (simplistic) model - update (Feb. 8, 2020).” medRxiv, 2020.2002.2004.20020461 DOI: 10.1101/2020.02.04.20020461. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMC4wMi4wNC4yMDAyMDQ2MXYyIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDQvMDMvMjAyMC4wMy4zMS4yMDA0OTQ4Ni5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 9. Hong, N., J. He, Y. Ma, H. Jiang, L. Han, L. Su, W. Zhu and Y. Long (2020) “Evaluating the secondary transmission pattern and epidemic prediction of the COVID-19 in metropolitan areas of China.” medRxiv, 2020.2003.2006.20032177 DOI: 10.1101/2020.03.06.20032177. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMC4wMy4wNi4yMDAzMjE3N3YzIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDQvMDMvMjAyMC4wMy4zMS4yMDA0OTQ4Ni5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 10. Hui, D. S., E. Azhar, T. A. Madani, F. Ntoumi, R. Kock, O. Dar, G. Ippolito, T. D. McHugh, Z.A. Memish, C. Drosten, A. Zumla and E. Petersen (2020). “The continuing 2019-nCoV epidemic threat of novel coronaviruses to global health — The latest 2019 novel coronavirus outbreak in Wuhan, China.” International Journal of Infectious Diseases 91: 264–266. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ijid.2020.01.009&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31953166&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F04%2F03%2F2020.03.31.20049486.atom) 11. Imai, N., A. Cori, I. Dorigatti, M. Baguelin, C. Donnelly, S. Riley and N. M. Ferguson (2020). Report 3: Transmissibility of 2019-nCoV London, UK, Imperial College London. 12. Imai, N., I. Dorigatti, A. Cori, C. Donnelly, S. Riley and N. M. Ferguson (2020). Report 2: Estimating the potential total number of novel Coronavirus cases in Wuhan City, China London, UK, Imperial College London. 13. Jowett, D., J. A. Browning and B. C. Haning (1974). NONLINEAR DISEASE PROGRESS CURVES. 14. Koopman, J. (2004). “Modeling infection transmission.” Annual Review of Public Health 25: 303–326. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1146/annurev.publhealth.25.102802.124353&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15015922&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F04%2F03%2F2020.03.31.20049486.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000221130600015&link_type=ISI) 15. Li, Q., X. Guan, P. Wu, X. Wang, L. Zhou, Y. Tong, R. Ren, K. S. M. Leung, E. H. Y. Lau, J. Y. Wong, X. Xing, N. Xiang, Y. Wu, C. Li, Q. Chen, D. Li, T. Liu, J. Zhao, M. Liu, W. Tu, C. Chen, L. Jin, R. Yang, Q. Wang, S. Zhou, R. Wang, H. Liu, Y. Luo, Y. Liu, G. Shao, H. Li, Z. Tao, Y. Yang, Z. Deng, B. Liu, Z. Ma, Y. Zhang, G. Shi, T. T. Y. Lam, J. T. Wu, G. F. Gao, B.J. Cowling, B. Yang, G. M. Leung and Z. Feng (2020). “Early Transmission Dynamics in Wuhan, China, of Novel Coronavirus–Infected Pneumonia.” New England Journal of Medicine. 16. Linton, N. M., T. Kobayashi, Y. Yang, K. Hayashi, A. R. Akhmetzhanov, S.-m. Jung, B. Yuan, R. Kinoshita and H. Nishiura (2020) “Epidemiological characteristics of novel coronavirus infection: A statistical analysis of publicly available case data.” medRxiv, 2020.2001.2026.20018754 DOI: 10.1101/2020.01.26.20018754. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMC4wMS4yNi4yMDAxODc1NHYyIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDQvMDMvMjAyMC4wMy4zMS4yMDA0OTQ4Ni5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 17. Liu, T., J. Hu, M. Kang, L. Lin, H. Zhong, J. Xiao, G. He, T. Song, Q. Huang, Z. Rong, A. Deng, W. Zeng, X. Tan, S. Zeng, Z. Zhu, J. Li, D. Wan, J. Lu, H. Deng, J. He and W. Ma (2020) “Transmission dynamics of 2019 novel coronavirus (2019-nCoV).” bioRxiv, 2020.2001.2025.919787 DOI: 10.1101/2020.01.25.919787. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiYmlvcnhpdiI7czo1OiJyZXNpZCI7czoxOToiMjAyMC4wMS4yNS45MTk3ODd2MiI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzA0LzAzLzIwMjAuMDMuMzEuMjAwNDk0ODYuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 18. Liu, X., G. J. D. Hewings, S. Wang, M. Qin, X. Xiang, S. Zheng and X. Li (2020) “Modeling the situation of COVID-19 and effects of different containment strategies in China with dynamic differential equations and parameters estimation.” medRxiv, 2020.2003.2009.20033498 DOI: 10.1101/2020.03.09.20033498. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMC4wMy4wOS4yMDAzMzQ5OHYxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDQvMDMvMjAyMC4wMy4zMS4yMDA0OTQ4Ni5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 19. Mansfield, E. and C. Hensley (1960). “THE LOGISTIC PROCESS - TABLES OF THESTOCHASTIC EPIDEMIC CURVE AND APPLICATIONS.” Journal of the Royal Statistical Society Series B-Statistical Methodology 22(2): 332–337. 20. NHC. (2020). “National Health Commission of the People’s Republic of China.” Retrieved Feb. 8, 2020, 2020, from [www.nhc.gov.cn/xcs/xxgzbd/gzbd_index.shtml](http://www.nhc.gov.cn/xcs/xxgzbd/gzbd_index.shtml). 21. Nishiura, H., S.-m. Jung, N. M. Linton, R. Kinoshita, Y. Yang, K. Hayashi, T. Kobayashi, B. Yuan and A. R. Akhmetzhanov (2020). “The Extent of Transmission of Novel Coronavirus in Wuhan, China, 2020.” Journal of Clinical Medicine 9(2): 330. 22. Roosa, K., Y. Lee, R. Luo, A. Kirpich, R. Rothenberg, J. M. Hyman, P. Yan and G. Chowell (2020). “Short-term Forecasts of the COVID-19 Epidemic in Guangdong and Zhejiang, China: February 13-23, 2020.” Journal of Clinical Medicine 9(2). 23. Shen, M., Z. Peng, Y. Xiao and L. Zhang (2020) “Modelling the epidemic trend of the 2019 novel coronavirus outbreak in China.” bioRxiv, 2020.2001.2023.916726 DOI: 10.1101/2020.01.23.916726. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiYmlvcnhpdiI7czo1OiJyZXNpZCI7czoxOToiMjAyMC4wMS4yMy45MTY3MjZ2MSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzA0LzAzLzIwMjAuMDMuMzEuMjAwNDk0ODYuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 24. Waggoner, P. E. and D. E. Aylor (2000). “Epidemiology: A science of patterns.” Annual Review of Phytopathology 38: 71-+. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1146/annurev.phyto.38.1.71&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11701837&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F04%2F03%2F2020.03.31.20049486.atom) 25. Wang, C., P. W. Horby, F. G. Hayden and G. F. Gao (2020). “A novel coronavirus outbreak of global health concern.” The Lancet. 26. Wu, J. T., K. Leung and G. M. Leung (2020). “Nowcasting and forecasting the potential domestic and international spread of the 2019-nCoV outbreak originating in Wuhan, China: a modelling study.” The Lancet. 27. Wu, P., X. Hao, E. H. Y. Lau, J. Y. Wong, K. S. M. Leung, J. T. Wu, B. J. Cowling and G. M. Leung (2020). “Real-time tentative assessment of the epidemiological characteristics of novel coronavirus infections in Wuhan, China, as at 22 January 2020.” Eurosurveillance 25(3): 2000044. 28. Wu, Z. and J. M. McGoogan (2020). “Characteristics of and Important Lessons From the Coronavirus Disease 2019 (COVID-19) Outbreak in China: Summary of a Report of 72?314 Cases From the Chinese Center for Disease Control and Prevention.” JAMA. 29. Zhang, C. and M. Wang (2020) “MRCA time and epidemic dynamics of the 2019 novel coronavirus.” bioRxiv, 2020.2001.2025.919688 DOI: 10.1101/2020.01.25.919688. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiYmlvcnhpdiI7czo1OiJyZXNpZCI7czoxOToiMjAyMC4wMS4yNS45MTk2ODh2MyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzA0LzAzLzIwMjAuMDMuMzEuMjAwNDk0ODYuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 30. Zhang, C. and M. Wang (2020) “Origin time and epidemic dynamics of the 2019 novel coronavirus.” bioRxiv, 2020.2001.2025.919688 DOI: 10.1101/2020.01.25.919688. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiYmlvcnhpdiI7czo1OiJyZXNpZCI7czoxOToiMjAyMC4wMS4yNS45MTk2ODh2MyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzA0LzAzLzIwMjAuMDMuMzEuMjAwNDk0ODYuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 31. Zhao, S., Q. Lin, J. Ran, S. S. Musa, G. Yang, W. Wang, Y. Lou, D. Gao, L. Yang, D. He and M.H. Wang (2020) “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.” bioRxiv, 2020.2001.2023.916395 DOI: 10.1101/2020.01.23.916395. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiYmlvcnhpdiI7czo1OiJyZXNpZCI7czoxOToiMjAyMC4wMS4yMy45MTYzOTV2MiI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzA0LzAzLzIwMjAuMDMuMzEuMjAwNDk0ODYuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 32. Zhao, S., Q. Lin, J. Ran, S. S. Musa, G. Yang, W. Wang, Y. Lou, D. Gao, L. Yang, D. He and M.H. Wang (2020). “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. 33. Zhao, S., S. S. Musa, Q. Y. Lin, J. J. Ran, G. P. Yang, W. M. Wang, Y. J. Lou, L. Yang, D. Z. Gao, D. H. He and M. H. Wang (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(2). 34. Zhou, X., N. Hong, Y. Ma, J. He, H. Jiang, C. Liu, G. Shan, L. Su, W. Zhu and Y. Long (2020) “Forecasting the Worldwide Spread of COVID-19 based on Logistic Model and SEIR Model.” medRxiv, 2020.2003.2026.20044289 DOI: 10.1101/2020.03.26.20044289. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMC4wMy4yNi4yMDA0NDI4OXYyIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDQvMDMvMjAyMC4wMy4zMS4yMDA0OTQ4Ni5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 35. Ziff, A. L. and R. M. Ziff (2020) “Fractal kinetics of COVID-19 pandemic.” medRxiv, 2020.2002.2016.20023820 DOI: 10.1101/2020.02.16.20023820. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMC4wMi4xNi4yMDAyMzgyMHYyIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDQvMDMvMjAyMC4wMy4zMS4yMDA0OTQ4Ni5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) [1]: /embed/graphic-1.gif [2]: /embed/graphic-2.gif [3]: /embed/graphic-3.gif [4]: /embed/graphic-4.gif