Prediction of 2019-nCov in Italy based on PSO and inversion analysis ==================================================================== * Sun Siyi * Zheng Yangping ## Abstract Although China achieved an early victory of controlling the novel coronavirus (2019-nCov), the overseas situation is overwhelming negative, especially in Italy. Up to March 11, 2020, 2019-nCov thoroughly broke out in Italy with over 10,000 confirmed cases notwithstanding the gradually block of the country since March 9, 2020. Estimation of possible infection population and prospective suggestion of handling spread based on exist data are of crucial importance. Considering of the biology parameters obtained based on Chinese clinical data in Wuhan, other scholars’ work and real spread feature of 2019-nCov in Italy, we built a more applicable model called SEIJR with log-normal distributed time delay to forecast the trend of spreading. Adopting Particle Swarm Optimization (PSO), we estimated the early period average spreading velocity (***α***) and conducted inversion analysis of time point (***T***) when the virus first hit the Italy. Based on fixed *α***** and ***T***, we then obtained the average spreading velocity ***α*1** after the lock by PSO. For the aim of offering expeditious advice, we generated the prediction trends with different ***α*** which we considered would be helpful in addressing the infection. Not only solved the complex, nondifferentiable equation of epidemic model, our research also performs well in inversion analysis based on PSO which conveys informative outcomes for further discussion on precatious action. To conclude, the first day of spread is around February 1, 2020 with the early period average spreading velocity *α*=0.330 which is higher than most cities in China except Wuhan. After locking the country and attaching great attention to public precaution, the *α*1 sharply descended to 0.278, indicting the effectiveness of these measures. Furthermore, in order to cope the disease before mid-April, take actions to control the under 0.25 is necessary. Code can be freely downloaded from [https://github.com/Summerwork/2019-nCov-Prediction](https://github.com/Summerwork/2019-nCov-Prediction). ## 1 Introduction A global epidemic disease known as the novel coronavirus (2019-nCov) had seriously hit the most area around the whole world causing unpredictable loss of manpower and finance during the first quarter of 2020 [21] [1]. The first confirmed 2019-nCov case was reported in Wuhan, China on December 31, 2019 (World Health Organization, 2020a). As of March 19 (21:00 GMT),2020, 2019-nCov has resulted in 81262 confirmed cases and 3250 dead cased in China cumulatively (National Health Commission of the People’s Republic of China, 2020). China as the first country faced by the outbreak of the severe disease, it took strict but effective action to contain the spread of 2019-nCov and attained apparent success till now. [21] Many related works have been done in prediction and precaution via constructing proper model and analyzing parameters. [21] [20] [19] While in the other parts worldwide, the menacing disease just became to spread [7], especially in countries with no preparation and experienced measures for suppressing the possible large-scale infection. In this article, we take Italy which is now experiencing severe situation of 2019-nCov as example to conduct analysis with the aim of offering utilizable suggestion. In the first period (before March 9, 2020), we attempt to inverse the virus spread timeline in Italy and the early period average spreading velocity by adopting PSO to optimize the parameters based on our SEIJR model and existing data (European Centre for Disease Prevention and Control). In the second period (only use data from March 9-16, 2020), we optimize the average spreading velocity in order to show the effect of country blockade. In the last period of our research, we demonstrate latent trends of confirmed cases with various average spreading velocity which of vital importance in controlling the disease. The whole flow chart is illustrated in Figure 1. ![Figure 1](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/13/2020.05.08.20095869/F1.medium.gif) [Figure 1](http://medrxiv.org/content/early/2020/05/13/2020.05.08.20095869/F1) Figure 1 Flow chart ## 2 Data and Model ### 2.1 Data Collection and Processing We collected the daily reported confirmed diagnosed data from the website of European Centre for Disease Prevention and Control (European CDC). All these data is public for everyone. Based on the need of our analysis, we preprocessed the data by adding the daily reported confirmed diagnosed cases to obtain the accumulative amount for following inversion of parameters and prediction. Here are details about the processed data we used in this article (Table 1).Table 1. View this table: [Table 1](http://medrxiv.org/content/early/2020/05/13/2020.05.08.20095869/T1) Table 1 New Confirmed and Accumulative Confirmed Data ### 2.2 SEIJR Model Disease transmission is a complex process with multiple variables and uncertainties making it unable to be accurately solved and predicted [16], [11] [6]. However, models are feasible in forecasting for infectious diseases when different characteristics parameters [18] [17] like transmission mode, immunization mode, mortality and average spreading velocity are offered. Classical models for infectious diseases include SIR model [2] and SEIR model, etc. Considering the actual situation in Italy and the transmission characteristics of 2019-nCov obtained in China, this paper built the SEIJR model with log-normal distributed time-delay terms [15] [3]based on the SEIR model [14] [12]. Figure 2 illustrates the SEIJR model. The model describes the problem by assuming population consists of six types (accumulative value): susceptible population S, exposed population *E*, infectious population *I*, confirmed *J*, recovered population *R* and dead population *D*. *α* is average spreading velocity, *β* is diagnose rate, *γ*1, *γ*2 are die rate and *μ* is cure rate. *t*1(*t*) is the time of incubation period [5] E need to become *I*, *t*2(*t*) is the time of waiting period I need to become *J* and *t*3(*t*) is the duration of hospitalization [5] *J* need to become *R*. *S*(*t*), *E*(*t*), *I*(*t*), *J*(*t*), *R*(*t*) and *D*(*t*) are dependent variables of time *t*, respectively. *α*, *β*, *γ*, *μ* are constants only related to the actual situation. Time-delay *t*1, *t*2, *t*3 are only depend on time *t*, respectively, which can be written as *t*1 (*t*), *t*2 (*t*) and *t*3 (*t*). ![Figure 2](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/13/2020.05.08.20095869/F2.medium.gif) [Figure 2](http://medrxiv.org/content/early/2020/05/13/2020.05.08.20095869/F2) Figure 2 SEIJR model ### 2.3 Interpretation In the model, only *E* and *I* have ability to infect *S* which is a process of contact infection with transient time. E indeed infected but shows no symptoms of 2019-nCov and then transfers to *I* after an incubation period *t*1. *I* has symptoms like fever, cough and shortness of breath. Because the pre-virus symptoms are not obvious [5] and the uneven medical facilities in Italy, *I* will be confirmed as *J* after a period *t*2. Due to the seriousness and infectivity of the virus, it can be considered that when it becomes *J*, *J* will be immediately isolation and lose its ability of infection. Treatment will be started immediately after confirmed. *J* will recovery and become *R* after a duration of hospitalization *t*3. Because of the 2019-nCov has certain lethality, it needs some more assumptions: #### 2.3.1 Assumption 1 The number of deaths during the period of *E* is too tiny to be considered, only need to consider the mortality during the period of *I* and *J*; #### 2.3.2 Assumption 2 The data during the period of I is unavailable. We consider that *I* to *D* and *I* to *J* have the same delay time *t*2 which means the only difference between them is proportion; #### 2.3.3 Assumption 3 The official organization and medical institutions do not give any information on how long for patients in the treatment stage will die, but there exists clinical information that for a confirmed patient how long the patient is needed to be cured [5]. As the same way *J* to *R*, *J* to *D* have the same delay time *t*3 only with different proportion. People in *I* can be diagnosed and receive treatment with diagnose rate. Otherwise, they will die with the die rate *γ*1 *=* 1−. The recover rate for *J* is *μ* and the dead rate is *γ*2 *=* 1 − *μ*. The total mortality proportion is *γ = γ*1 *+ γ*2. #### 2.3.4 Assumption 4 After recovering, people in R will not go out because they are in a frail state and will be considered as isolation. ### 2.4 Time-delay Function Combining the basic principles of epidemiology and etiology, it can get that the time *t*2 from infection to diagnosis, which approximately follows the log-normal distribution [8]. The assumption can be applied to *t*1 and *t*2 without loss of generality. For *t*1, due to the lack of Italian clinical information statistics, it’s more accurate and reliable to use the clinical data of Zhong NS et.al: the median incubation period is 4 days [5], the quartile is 2 and 7 days. With the log-normal distribution’s formula: ![Formula][1] Get *μ = ln*4,![Graphic][2] where *Z*0.75 is normal distribution quartile. Then, the log-normal density function corresponding to *t*1 is: ![Formula][3] For *t*2, the data in Italy is still unavailable, but the relevant distribution function is given in the Chinese research report. In the paper of Yang ZW et.al, the relevant statistics of “the time interval between the most recent stay in Hubei Province and the confirmed diagnosis” are given [4]. The end point of this period corresponds to the start point of *J* in the SEIJR model, but there is no corresponding start point. It can only be known that it’s in the middle of the time between *t*1 and *t*2, it is convincing to regard it as the point when *e* enters *I*, which means the time given in that article corresponds to the *t*2 in this article. As the same way, the log-normal density function corresponding to *t*2 is: ![Formula][4] For *t*3, the Italian official organizations and medical institutions still lack clinical information and official statistics, so we also assumes that *t3* follows the log-normal distribution. However, Zhong NS et.al gave some relevant data: the duration of hospitalization which is *t3*, the median is 12 days, and the quartiles are 10 and 14 days [5]. Adopting the same method used in *t*1 estimation, we can get the distribution density function corresponding to *t*3 is: ![Formula][5] Now give the definition of log-distributed normal time-delay function in SEIJR model: ### Definition 1 (Time-delay Function in SEIJR Model) *The time-delay function describes how many people change this stage from the previous stage at time t.* The general formula is written as: ![Formula][6] where *x*(*t*) is the term which has time-delay, and *si*(*t*) is the distribution density mentioned before. Here we finally get the precise differential equation of SEIJR model as: ![Formula][7] ## 3 Methods ### 3.1 Runge-Kutta Methods The previous differential Eq. (6) corresponding to the SEIJR model is a form with integrals and independent variables on the integral limit which from the time-delay function term, so there is no analytical solution. For this case, numerical methods are useful. We combine iteration and degree four Runge-Kutta to generate the numerical solution with the initial condition. Here is the main principles of Runge–Kutta fourth-order method [8], let the differential equation have the form as follow: ![Formula][8] Then its iterative formula is: ![Formula][9] where ![Formula][10] ### 3.2 Particle Swarm Optimization Usually the epidemic model’s descriptive function was derivative-based without considering time delay or just assume fixed linear time delay, such as SIR model. In this article, the SEIJR model solved by combining iteration and degree four Runge-Kutta mentioned before. Accounting for this, the function of least square method (LSE) employed in addressing this problem for different periods are denoted as follow: ![Formula][11] where *f*1 is LSE of the first period. *f*2 is LSE of the second period.![Graphic][12] are predicted and actual value from February 22 to March 9, 2020.![Graphic][13] are predicted and actual value from March 9 to March 25, 2020. Both of them are nonlinear, complex, discontinuous and nondifferentiable [13]. Traditional method based on derivation is infeasible for minimization. While approaches such as annealing algorithm [10] to search for parameters is more calculative expensive than the particle swarm optimization (PSO). PSO [14] is a population-based search algorithm sparked by the forage behavior of birds within a flock. Individuals gain the ability of searching for better solution areas by learning the fitness information from the environment. In PSO algorithm, the velocity of individual is dynamically changed considering its previous flying experience. The algorithm consists of three main parts: individual best, global best and individual optimization based on the best particle of whole population [9]. In this article, *f*1 and *f*2 are the fitness function during the first and second period. The main process of conducting this algorithm and basic parameters are respectively shown in Figure 3 and Table 2. ![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/13/2020.05.08.20095869/F3.medium.gif) [Figure 3.](http://medrxiv.org/content/early/2020/05/13/2020.05.08.20095869/F3) Figure 3. PSO flow chart View this table: [Table 2](http://medrxiv.org/content/early/2020/05/13/2020.05.08.20095869/T2) Table 2 Parameters for PSO ## 4 Results ### 4.1 The first period During this period of our work, we adopt PSO using *f*1 as its fitness function to generate the optimal *α*,*T*. As shown in the Figure 4, the best result is *T*=21 (after rounding) and *α*=0.33. After obtained the greatest *α*,*T*, we draw the prediction curves for further comparison in the following section. Illustrated in Figure 5a is the prediction value![Graphic][14] with the optimal parameter and the actually value *Jact*1; (b) is the same curve with whole actual value *Jact* if the average spreading velocity remain unchanged. ![Figure 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/13/2020.05.08.20095869/F4.medium.gif) [Figure 4.](http://medrxiv.org/content/early/2020/05/13/2020.05.08.20095869/F4) Figure 4. Vistualization of PSO searched results ### 4.2 The second period Considering of the precautious action like blockade taken by Italy government and the information Figure 5b conveyed, we assumed the average spreading velocity changed after these actions. Based on fixed *α*,*T*, we applied PSO to the optimization problem and got the result of *α*1=0.278 which is cut down from *α*. The Figure 6 shows the predicted population![Graphic][15] with *α*1 from March 9 to 25, 2020 and *Jact*2. ![Figure 5.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/13/2020.05.08.20095869/F5.medium.gif) [Figure 5.](http://medrxiv.org/content/early/2020/05/13/2020.05.08.20095869/F5) Figure 5. ![Graphic][16] and *Jact*1*,Jact* with *α*=0.330 ![Figure 6.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/13/2020.05.08.20095869/F6.medium.gif) [Figure 6.](http://medrxiv.org/content/early/2020/05/13/2020.05.08.20095869/F6) Figure 6. ![Graphic][17] and *Jact*2 with *α*=0.278 ### 4.3 Predicted J trends with different *α* Figure 7 gives prediction curves with different *α* with the aim for further discussion about controlling the disease before mid-April. ![Figure 7](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/13/2020.05.08.20095869/F7.medium.gif) [Figure 7](http://medrxiv.org/content/early/2020/05/13/2020.05.08.20095869/F7) Figure 7 *Jpre* with different *α* ## 5 Discussion ### 5.1 Reliability of SEIJR model It is clearly seen that the fitting curve of *J* in the SEIJR model fits well with the confirmed diagnosis data of the Italian official statistics. *T* *=* 21 represents the initial value of the model which indicates the first E appeared 21 days before February 22 that is February 1. According to the information reported by the Italian government, the first case in Italy, when J appeared, was January 31. Two patients from Wuhan, China arrived in Italy by air at January 23 and visited other cities in Italy. They finally arrived in Rome, feeling physical discomfort at January 30 and was confirmed diagnosis then isolated at January 31. It is reasonable to speculate that they had carried the 2019-nCov in Wuhan before arriving in Italy in January 23, and continued to spread in Italy for eight days after January 23 until January 31. The first case of *E* reported in official data appeared at January 23, and the model’s conclusion appeared at February 1. In fact, because the two senior travelers are 67 and 66 years old with poor mobility, unfamiliar surroundings and the simple factors of interpersonal relationships, they have a lower *α* which also leads to real initial time point is earlier than the initial time point of the theoretical model. To sum up, the model has higher accuracy and stronger credibility. [11] ### 5.2 Effect of blockade Prime Minister Giuseppe Conte extended the quarantine lockdown to cover all the region of Lombardy and 14 other northern provinces on March 8, and all region on March 10. At the same time, the Italian government further banned rallies and sports activities nationwide, announced a national blockade, and unnecessarily stopped going out. Compared the average spread velocity *α*1 = 0.278 with *α* *=* 0.330 at the early stage(before March 9, 2020), a decrease existed after the series of actions which verifies the effectiveness of this measure. ### 5.3 Advice to Italy government According to the model’s prediction, when *α* is in a suitable range, the number of confirmed patients will gradually decrease. It is obvious that the smaller the value of *α* is, the faster the number of confirmed diagnoses decreases. Since Italian government hopes to end the 2019-nCov by mid-April, it can control *α* around 0.25 according to the prediction of the model. Based on this conclusion, Italy government is supposed to strengthen the isolation, reduce gathering activities and make people understand the importance of precautionary measures to constrain *α* under 0.25. ## 6 Conclusion Applying PSO to our SEIJR model with log-normal distributed time delay, we obtained the convincing start time (around February 1, 2020) of 2019-nCov and the average spreading velocity (*α*=0.330) at the early stage. We compared the average spreading velocity during the early period and following period, a conspicuous decrease attributed to the effective measures was found. Based on the prediction interval of possible infected population of different *α*, we strongly recommend Italy to keep *α* under 0.25 if they want the situation take a turn for the better even ended before mid-April. ## Data Availability We collected the daily reported confirmed diagnosed data from the website of European Centre for 51 Disease Prevention and Control . All these data is public for everyone. 52 Based on the need of our analysis, we preprocessed the data by adding the daily reported confirmed 53 diagnosed cases to obtain the accumulative amount for following inversion of parameters and prediction. 54 Here are details about the processed data we used in this article. [https://www.ecdc.europa.eu/en](https://www.ecdc.europa.eu/en) ## Acknowledgment The authors would like to thank Dongmei Ai at University of Science and Technology Beijing for her helpful suggestions and discussion. * Received May 8, 2020. * Revision received May 11, 2020. * Accepted May 13, 2020. * © 2020, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution 4.0 International), CC BY 4.0, as described at [http://creativecommons.org/licenses/by/4.0/](http://creativecommons.org/licenses/by/4.0/) ## References 1. 1. A. Al-Mandhari, D. Samhouri, A. Abubakar, and R. Brennan. Coronavirus disease 2019 outbreak: preparedness and readiness of countries in the eastern mediterranean region. Eastern Mediterranean health journal= La revue de sante de la Mediterranee orientale= al-Majallah al-sihhiyah li-sharq al-mutawassit, 26(2):136, 2020. 2. 2. B. Cantó, C. Coll, and E. Sánchez. Estimation of parameters in a structured sir model. Advances in Difference Equations, 2017(1):33, 2017. 3. 3. Y. Chen, J. Cheng, Y. Jiang, and K. Liu. A time delay dynamical model for outbreak of 2019-ncov and the parameter identification. *arXiv preprint arXiv:2002.00418*, 2020. 4. 4. K. J. e. a. Ding Z, Liu Y. A probability model for estimating the expected number of the newly infected and predicting the trend of the diagnosed. Operations Research Transactions, 24(1):1–12, 2020. 5. 5. W.-j. Guan, Z.-y. Ni, Y. Hu, W.-h. Liang, C.-q. Ou, J.-x. He, L. Liu, H. Shan, C.-l. Lei, D. S. Hui, et al. Clinical characteristics of coronavirus disease 2019 in china. New England Journal of Medicine, 2020. 6. 6. X.-N. Han, S. J. De Vlas, L.-Q. Fang, D. Feng, W.-C. Cao, and J. D. F. Habbema. Mathematical modelling of sars and other infectious diseases in china: a review. Tropical Medicine & International Health, 14:92–100, 2009. 7. 7. B. Haynes, N. E. Messonnier, and M. S. Cetron. First travel-related case of 2019 novel coronavirus detected in united states: press release, tuesday, january 21, 2020. 2020. 8. 8. J. Kalbfleisch and J. Lawless. Estimating the incubation period for aids patients. Nature, 333(6173):504–505, 1988. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=2836736&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F13%2F2020.05.08.20095869.atom) 9. 9. J. Kennedy and R. Eberhart. Particle swarm optimization. 4:1942–1948, 1995. 10. 10. A. Khosravi, S. Nahavandi, D. Creighton, and A. F. Atiya. Lower upper bound estimation method for construction of neural network-based prediction intervals. IEEE transactions on neural networks, 22(3):337–346, 2010. 11. 11. J. Koopman. Modeling infection transmission. Annu. Rev. Public Health, 25:303–326, 2004. [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%2F05%2F13%2F2020.05.08.20095869.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000221130600015&link_type=ISI) 12. 12. P. E. Lekone and B. F. Finkenstadt. Statistical inference in a stochastic epidemic seir model with control intervention: Ebola as a case study. Biometrics, 62(4):1170–1177, 2006. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1541-0420.2006.00609.x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17156292&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F13%2F2020.05.08.20095869.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000242771800024&link_type=ISI) 13. 13. S. A. Levin, B. Grenfell, A. Hastings, and A. S. Perelson. Mathematical and computational challenges in population biology and ecosystems science. Science, 275(5298):334–343, 1997. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIyNzUvNTI5OC8zMzQiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8wNS8xMy8yMDIwLjA1LjA4LjIwMDk1ODY5LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 14. 14. M. Y. Li and J. S. Muldowney. Global stability for the seir model in epidemiology. Mathematical biosciences, 125(2):155–164, 1995. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/0025-5564(95)92756-5&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=7881192&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F13%2F2020.05.08.20095869.atom) 15. 15. G. Medley, L. Billard, D. R. Cox, and R. M. Anderson. The distribution of the incubation period for the acquired immunodeficiency syndrome (aids). Proceedings of the Royal society of London. Series B. Biological sciences, 233(1272):367–377, 1988. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1098/rspb.1988.0026&link_type=DOI) 16. 16. T. W. Ng, G. Turinici, and A. Danchin. A double epidemic model for the sars propagation. BMC Infectious Diseases, 3(1):19, 2003. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12964944&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F13%2F2020.05.08.20095869.atom) 17. 17. J. T. Wu, K. Leung, and G. M. Leung. Nowcasting and forecasting the potential domestic and international spread of the 2019-ncov outbreak originating in wuhan, china: a modelling study. The Lancet, 395(10225):689–697, 2020. 18. 18. P. S. Yip, K. Lam, Y. Xu, P. Chau, J. Xu, W. Chang, Y. Peng, Z. Liu, X. Xie, and H. Lau. Reconstruction of the infection curve for sars epidemic in beijing, china using a back-projection method. Communications in Statistics—Simulation and Computation®, 37(2):425–433, 2008. 19. 19. S. Zhang, M. Diao, W. Yu, L. Pei, Z. Lin, and D. Chen. Estimation of the reproductive number of novel coronavirus (covid-19) and the probable outbreak size on the diamond princess cruise ship: A data-driven analysis. International Journal of Infectious Diseases, 2020. 20. 20. S. Zhao, Q. Lin, J. Ran, S. S. Musa, G. Yang, W. Wang, Y. Lou, D. Gao, L. Yang, D. He, 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, 92:214–217, 2020. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ijid.2020.01.050&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32007643&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F13%2F2020.05.08.20095869.atom) 21. 21. J. Zheng. Sars-cov-2: an emerging coronavirus that causes a global threat. Int J Biol Sci, 16(10):1678–1685, 2020. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.7150/ijbs.45053&link_type=DOI) [1]: /embed/graphic-4.gif [2]: /embed/inline-graphic-1.gif [3]: /embed/graphic-5.gif [4]: /embed/graphic-6.gif [5]: /embed/graphic-7.gif [6]: /embed/graphic-8.gif [7]: /embed/graphic-9.gif [8]: /embed/graphic-10.gif [9]: /embed/graphic-11.gif [10]: /embed/graphic-12.gif [11]: /embed/graphic-13.gif [12]: /embed/inline-graphic-2.gif [13]: /embed/inline-graphic-3.gif [14]: /embed/inline-graphic-4.gif [15]: /embed/inline-graphic-5.gif [16]: F5/embed/inline-graphic-6.gif [17]: F6/embed/inline-graphic-7.gif