MATHEMATICAL MODELING FOR TRANSMISSIBILITY OF COVID-19 VIA MOTORCYCLES ====================================================================== * Benard Okelo ## Abstract In this paper, we present a mathematical model of trigonometric type for transmissibility and deaths as a result of COVID-19. In the model, we analyze the spread of COVID-19 by considering a new parameter, the motor-cycle as a means of public transport, which has not been considered in several other models for COVID-19. We use the mathematical model to predict the spread and deaths and we suggest strategies that can be put in place to prevent the spread caused by motorcycle as a means of public transport. Keywords * Mathematical model * COVID-19 * Motorcycle ## 1. Introduction The outbreak of COVID-19 was detected in December 2019 in mainland China with the city of Wuhan as the recognized epicenter(see [1]-[29] and the references therein). COVID-19 has since then been exported to other countries all over the world[7]. To describe and determine the dynamics of COVID-19, several mathematical models have been constructed by several researchers[2]-[6]. We realize that the insights that can be drawn from these models are discussed, especially as inputs for designing strategies to control the epidemics. Proposed model-based strategies on how to prevent the spread of the disease in local setting, such as during large social gatherings, are also discussed in details[8]-[13]. One of the parameters considered in most of the models is the exposure time which is an instrumental factor in spreading the disease[14]. With a basic reproduction number equal to 2, and 14-day infectious period, an infected person staying more than 9 hours in the event could infect other people [15]. As at now, there is an urgent need to develop a mathematical model to estimate the transmissibility and dynamic of the transmission of the virus[16]-[19]. Several researches [17], [18], [21], [22], [22] and [25] have focused on calculating the basic reproduction number *R* by using the serial intervals and intrinsic growth rate or using ordinary differential equations and Markov Chain Monte Carlo methods[23] and [24]. A lot of these researches have certain similar characteristics: They show that social distancing can significantly slow down the spread of the COVID-19; They stipulated that the effectiveness depends very much on the initial and boundary conditions of social distancing and the results strongly vary by orders of magnitude if the process is initiated a few days earlier or later or the effectiveness of the measure is only slightly higher or lower; The details of the studies are difficult to understand for many readers and therefore a widespread trust on the efficiency on the actual measures regarding the COVID-19 pandemic can not easily be derived from scientific literature. Results from [26] show that the effectiveness of social distancing is particularly great for pandemic with a basic reproductive number of *R* = 1.5 - 2.5. For the COVID-19 pandemic, a basic reproduction number *R* = 2.6 is estimated, with an uncertainty range from 1.5 to 3.5 as seen in [27]. *R* specifies how many new cases are caused by an existing case on average. With *R*¿ 1, the virus spreads exponentially at times, but for *R* ¡1 the spread stops on its own [30]. A large number of studies have now calculated the number of basic reproductions, the incubation period and the duration of the infectiousness based on the number of cases from China [30]. However, the transmission through use of motorcycles as a means of public transport particularly in developing economies have not considered in the published models. The COVID-19 therefore is still a world pandemic and is currently a challenge for several countries all over the the world but worse for the developing countries. The goal of this paper is to present a mathematical model of trigonometric form for the spread of the COVID-19, where we concentrate on a new parameter called the use of motorcycles as a means of public transport to predict spread and number of deaths due to the COVID-19. ## 2. Materials and methods The reported cases of COVID-19, were collected for the modelling study from the World Health Organization official website. The epidemic curve from a sample of 30 days is presented. Simulation methods and statistical analysis was done using MATLAB version 9.7. The fourth-order RungeKutta method, with tolerance set at 0.001, was used to perform curve fitting. Selected preliminary models of COVID-19, published in official webpages of academic/research institutions, as preprints, or as journal articles since the outbreak are reviewed. Important insights from the results of these models are discussed. Moreover, a model using a Susceptible-Exposed-Infected (SEI) framework is has been carefully reviewed to propose measures to prevent epidemics during large events, e.g., during parties or concerts with huge crowds. The parameter values used in this model are based on the known information about COVID-19. The model is represented by a system of instantaneous and noninstantaneous differential equations, and simulated using MATLAB version 9.7. The Differential equations are then reduced into a trigonometric form and later analysed. Motorcycles are used in most developing countries for public transport. In Kenya, they are called “Boda Boda”, while in Nigeria they are commonly known as “Okada.” Most people in the cities, towns and even in rural areas prefer motorcycle because they are faster, relatively cheap and can be used to avoid traffic jam in roads. Despite their usefulness, motorcycles contributes immensely to the spread of COVID-19 because the operators do not follow the traffic rules. For instance, they don’t put masks and helmets, they carry excess passengers, they don’t wash their hands as required. Hence, social distancing is not taken care of. Due to these factors, we study the transmissibility of COVID-19 as a result of the risk factors from the motorcycle operators. ## 3. Derivation of the mathematical model We consider the use of motorcycles as a parameter presented in form of differential equations with instantaneous impulses and we compare it with the case of noninstantaneous impulses. ### Case I. Instantaneous impulses Let an increasing sequence of points ![Graphic][1] be given such that 0 ≤ *t**i* *< t**i*+1, *i* = 1, 2, …, *p*, lim*k*→∞ *t**k* = ∞ Let *t* ∈ *R**n* be a given arbitrary point. Without loss of generality we assume that *t* [0, *t*1). Consider the Initial Value Problem (IVP) for the nonlinear Instantaneous Impulsive Differential equation (IDE) ![Formula][2] where *x, x* ∈ ℝ*n*, *f* : [*t*, ∞) *×*ℝ*n* → ℝ*n*, *I**k* : ℝ*n* → ℝ*n*, (*k* = 1, 2, 3, …). The points *t**k*, *k* = 1, 2, … are called points of instantaneous impulses and the functions *I**k*(*x*), *k* = 1, 2, … are called instantaneous impulsive functions. The solution *x* (*t*; *t*, *x*) in the general case is a piecewise continuous function which is satisfying the integral equation ![Formula][3] ### Case II. Noninstantaneous impulses Let two increasing finite sequences of points ![Graphic][4] and ![Graphic][5] be given such that *s* = 0 *< t**i* ≤ *s**i* *< t**i*+1, *i* = 1, 2, 3, …, *p* and points *t*, *T* ∈ ℝ+ are given such that *s* = 0 *< t* *< t*1, *t**p* *< T t**p*+1, *p* is a natural number. Consider the Initial Value Problem(IVP) for the nonlinear Noninstantaneous Impulsive Differential Equation(NIDE) ![Formula][6] where *x, x* ∈ ℝ, *f**k* : [*s**k*, *t**k*+1] ∩ [*t*, *T*] *×*ℝ → ℝ, *k* = 0, 1, 2, …, *p, ϕ**k* : [*t**k*, *s**k*] *×*ℝ *×* ℝ → ℝ, *k* = 1, 2, 3, … The intervals (*s**k*, *t**k*+1], *k* = 1, 2, … are called intervals of noninstantaneous impulse and the function *ϕ**k*(*t, x, y*), *k* = 1, 2, …, are called noninstantaneous impulsive function. In the special case *s**k* = *t**k*+1, *k* = 0, 1, 2, … each interval of noninstantaneous impulses is reduced to a point, and the problem in Equation 3.3 is reduced to an IVP for an IDE in Equation 3.1 with points of impulses *t**k* and impulsive functions *x* (*t**k* + 0) = *I**k* (*x* (*t**k* − 0)) ≡ *ϕ**k* (*t**k*, *x* (*t**k* − 0), *x* (*t**k* − 0)) *x* (*t**k* − 0) Let *k* ≥ 0 be a given integer, *τ* ∈ [*t**k*, *s**k*) be a given point and consider the corresponding IVP for ODE ![Formula][7] We give a brief description of the solution of IVP for NIDE in Equation 3.4 The solution *x* (*t*; *t*, *x*) of Equation 3.3 is given by ![Formula][8] where; on the interval [*t*, *s*] the solution coincides with *X*(*t*) which is the solution of IVP for ODE in Equation 3.4 for *τ* = *t*, *k* = 0 and ![Graphic][9] on the interval (*s*, *t*1] the solution *x* (*t*; *t*, *x*) satisfies the equation *x* (*t*; *t*, *x*) = *ϕ* (*t, x* (*t*; *t*, *x*), *X* (*s* − 0)); on the interval (*t*1, *s*1] the solution coincides with *X*1(*t*) which is the solution of IVP for ODE in Equation 3.4 for *τ* = *t*1, *k* = 1 and ![Graphic][10]; on the interval (*s*1, *t*2] the solution *x* (*t*; *t*, *x*) satisfies the equation *x* (*t*; *t*, *x*) = *ϕ*1 (*t, x* (*t*; *t*, *x*), *X*1 (*s*1 0)); on the interval (*t*2, *s*2] the solution coincides with *X*2(*t*) which is the solution of IVP for ODE in Equation 3.4 for *τ* = *t*2, *k* = 2 and ![Graphic][11]; and so on. Also, the solution *x* (*t*; *t*, *x*), *t ≥ t* of Equation 3.3 satisfies the following system of integral and algebraic equations ![Formula][12] **Example 3.1**. Consider the IVP for the scalar NIDE ![Formula][13] where *x, x* ∈ ℝ, *A**k*, *k* = 0, 1, 2, …, are constants. The solution of Equation 3.7 is given by ![Formula][14] Now, we consider a number of cases: **Case 1**. Let *A**k* = *A, ϕ**k*(*t, x, y*) = *a**k*(*t*)*x*2*y, a**k* : [*s**k*, *t**k*+1] → ℝ*/{**}, k* = 0, 1, 2, … and *x* ≠ 0. since the nontrivial solution of *x* = *a**k*(*t*)*x*2*y* for any *y*∈ ℝ, *y* ≠ 0 is ![Graphic][15], then the solution of NIDE in Equation 3.7 is given by ![Formula][16] A particular case where *A* = 1, *s**k* = 2*k*− 1, *t**k* = 2*k, a**k*(*t*) = *t, k* = 0, 1, 2 …, 5 has solutions for two different initial values *x* = 0.5 and *x* = 0.3. **Case2**. Let *A**k* = *A, ϕ**k*(*t, x, y*) = *a**k*(*t*)*y, a**k* : [*s**k*, *t**k*+1] → ℝ, *k* = 0, 1, 2, 3, … The solution NIDE in Equation 3.7 is given by ![Formula][17] In the particular case *A* = −1, and *t**k* = 2*k, s**k* = 2*k*+1, *a**k*(*t*) = *t* for *k* = 0, 1, 2, … we have the solutions for the different initial values *x* = 1, 0.5 and 0.2. **Case 3**. Let *A* = 0 and *ϕ**k*(*t, x, y*) = *a**k*(*t*)*y, a**k* : [*t**k*, *s**k*] → ℝ, *k* = 1, 2, …, *p* The solution of NIDE in Equation 3.7 is given by ![Formula][18] **Example 3.2**. Comparisons between the behavior motorcycle operators that leads to the spread of COVID-19 with the corresponding impulsive differential equation and the equation with non-instantaneous impulses are given in the following cases. **I**. Ordinary Differential Equations. Consider the IVP for the scalar ODE ![Formula][19] where *x* ∈ ℝ. The solution of Equation 3.10 is ![Graphic][20] and it approaches 0 as *t* →∞. **II**. Impulsive differential equations. Consider the IVP for the scalar Impulsive Differential Equation (IDE) ![Formula][21] where *x* ∈ ℝ, *a**k*, *k* = 1, 2,… *· k* are constants. The solution of Equation 3.11 is given by ![Graphic][22] for *t* ∈ (*t**k*, *t**k*+1]. The behavior of the solution depends significantly on the amount of the impulse, i.e., the value of the constants *a**k*. **Case II.1**. If *a**k* = 1, *k* = 1, 2, …, then the problem in Equation 3.11 coincides with the IVP for ODE in Equation 3.10 and the solution approaches 0 as *t* → ∞. **Case II.2**. If there exists a natural number *k* : *a**k* = 0, then the solution of Equation 3.11 is zero for *t > t**k*, i.e., all solutions in spite of the initial value *x* coincide for *t > t**k*. **Case II.3**. If *a**k* = *e*−*k*, *t**k* = *k, k* = 1, 2, …, then the solution of the problem in Equation 3.11 is *x*(*t*) = *x**e**t*−*k* for *t* ∈ (*k, k* + 1] and the solution is a periodic function. **Case II.4**. If *a**k* = *e*−2*k*, *t**k* = *k, k* = 1, 2, …, then the solution of the problem in Equation 3.11 is *x*(*t*) = *x**e**t*−2*k* for *t* ∈ (*k, k* + 1] and the solution being a piecewise continuous function approaches 0. **Case II.5**. If *a**k* = *e*−2*k*, *t**k* = *k, k* = 1, 2, …, then the solution of the problem in Equation 3.11 is *x*(*t*) = *x**e**t*+0.5*k* for *t* ∈ (*k, k* + 1] and the solution is an unbounded function. **III**. Differential equations with non-instantaneous impulses. Consider the IVP for the scalar NIDE ![Formula][23] where *x* ∈ ℝ, *a**k*, *k* = 0, 1, 2, … *k* are constants, *t**k* = *k, s**k* : *k < s**k* ≤ *k* + 1 for *k* = 0, 1, … Note in the case *t**k*+1 = *s**k* the second equation of 3.12 is reduced to *x* (*t**k* + 0) = *a**k**x* (*t**k* − 0). We consider the case *t**k*+1 ≠ *s**k*. For the spontaneous transmission of COVID-19 by the motorcycles, the solution is given by ![Formula][24] If *a**k* = 0, *k* = 1, 2, …, then the solution of (3.1.12) is given by ![Formula][25] The solutions do not depend on the initial value *x* for *t > t*1 At this point, we consider the mathematical models with motorcycles as a parameter trigonometrically. Consider Equation 3.13. With some complex and tedious computation, we can trigonometrically represent it as ![Formula][26] where *k* is a positive real number. We note that the days in the period at which ![Graphic][27] and ![Graphic][28] are points when COVID-19 is neither being transmitted nor spread and no deaths are realized. Now taking the following in account: At ![Graphic][29] and ![Graphic][30] (as the only solution). The day considered to be the point (0, 0) is the maximal day point which is stable as per the fourth-order RungeKutta method, with tolerance set at 0.001, was used to perform curve fitting. Since Equation 3.14 is positive definite while ![Graphic][31], the origin is Lyapunov stable with a trajectory curve ![Graphic][32]. Considering the maximal point and the initial conditions of Equation 3.14 we get ![Formula][33] From Figure 1 below, considering different transmission levels i.e. *C* = 0, *C* = 5 and *C* = 10 and putting *k* = 1, the curve in the phase plane has symmetry in the *x* axis only. ![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/24/2020.04.18.20070797/F1.medium.gif) [Figure 1.](http://medrxiv.org/content/early/2020/04/24/2020.04.18.20070797/F1) Figure 1. COVID-19 transmissibility at the maximal points From ![Formula][34] It can be shown that a Lyapunov function *V* (*x, y*) for the Equation 3.15 is given by ![Formula][35] and the origin (0, 0) is a stable equilibrium point of the system of differential Equation 3.15. In fact Equation 3.14 is but a particular case of Equation 3.15. One easily deduces that though the origin (0, 0) is stable, it is not asymptotically stable. Data was collected over 30 days with regard to positivity indices from COVID-19 cases and questions were asked from the respondents with regard to use of motorcycles a public means of transport in different countries all over the world. This leads to the following results: View this table: [Table1](http://medrxiv.org/content/early/2020/04/24/2020.04.18.20070797/T1) The curve in Figure 2 shows that the COVID-19 is spread through the use of motorcycle and the rate of infection is continuous. The curve reaches the peak but the rate of reduction in terms of death and infections is not rapid. Therefore, the pandemic might take longer than expected for it to be contained and to be declared not to be a pandemic anymore. ![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/24/2020.04.18.20070797/F2.medium.gif) [Figure 2.](http://medrxiv.org/content/early/2020/04/24/2020.04.18.20070797/F2) Figure 2. The curve of COVID-19 transmissibility over thirty day period through motorcycle operations. ## 4. Conclusion The results from this study is in agreement with the recommendation of World Health Organization [22] that it is still possible to interrupt virus spread, provided that countries put in place strong measures to detect disease early, for example, development of rapid diagnostic tests, and increasing effectiveness of passenger screening in airports in which thermal screening for COVID-19 infection is estimated using simulation to be 46 percent. Moreover, motorcycle operators are not indispensable when it comes to observing recommendation of World Health Organization and therefore if they adhere strictly to these recommendations, there will be a tremendous decrease in new infections and consequently reduces deaths due to COVID-19. ## Data Availability Data was simulated from sotwares * Received April 18, 2020. * Revision received April 18, 2020. * Accepted April 24, 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/) ## 1.1 References 1. 1.Bogoch II, Watts A, Thomas-Bachli A, Huber C, Kraemer MUG, Khan K. Pneumonia of unknown etiology in Wuhan, China: Potential for International Spread Via Commercial Air Travel. J Travel Med. 2020. [https://doi.org/10.1093/jtm/taaa008](https://doi.org/10.1093/jtm/taaa008). 2. 2.Chan JF, Yuan S, Kok KH, To KK, Chu H, Yang J, et al. A familial cluster of pneumonia associated with the 2019 novel coronavirus indicating personto-person transmission: a study of a family cluster. Lancet. 2020. [https://doi.org/10.1016/S0140-6736(20)30154-9](https://doi.org/10.1016/S0140-6736(20)30154-9). 3. 3.Chen T, Rui J, Wang Q, Zhao Z, Cui J-A, Yin L. A mathematical model for simulating the transmission of Wuhan novel Coronavirus. bioRxiv. 2020: 2020.2001.2019.911669. Accessed 13 Feb 2020. 4. 4.Chen T, Ka-Kit Leung R, Liu R, Chen F, Zhang X, Zhao J, et al. Risk of imported Ebola virus disease in China. Travel Med Infect Dis. 2014;12:6508. 5. 5.Chen T, Leung RK, Zhou Z, Liu R, Zhang X, Zhang L. Investigation of key interventions for shigellosis outbreak control in China. PLoS One. 2014;9: e95006. 6. 6.Chen S, Yang D, Liu R, Zhao J, Yang K, Chen T. Estimating the transmissibility of hand, foot, and mouth disease by a dynamic model. Public Health. 2019;174:428. 7. 7.Cui J-A, Zhao S, Guo S, Bai Y, Wang X, Chen T. Global dynamics of an epidemiological model with acute and chronic HCV infections. Appl Math Lett. 2020;103:106203. 8. 8.Huang C, Wang Y, Li X, Ren L, Zhao J, Hu Y, et al. Clinical features of patients infected with 2019 novel coronavirus in Wuhan, China. Lancet. 2020. [https://doi.org/10.1016/S0140-6736(20)30183-5](https://doi.org/10.1016/S0140-6736(20)30183-5). 9. 9.Ivorra B, Ngom D and Ramos AM, Be-CoDiS: A Mathematical Model to Predict the Risk of Human Diseases Spread Between CountriesValidation and Application to the 20142015 Ebola Virus Disease Epidemic. Bulletin of Mathematical Biology. 2015. Volume 77, Issue 9, pp 16681704. DOI: [http://dx.doi.org/10.1007/s11538-015-0100-x](http://dx.doi.org/10.1007/s11538-015-0100-x) -Preprint: [https://eprints.ucm.es/28809/](https://eprints.ucm.es/28809/) 10. 10. Japan MoH. Development of patients with pneumonia associated with new coronavirus. Ministry of Health Japan. cited 16 Jan 2020. 11. 11.Jonathan M, Jessica RE, Bridgen AT, Cummings AH, Chris PJ. Novel coronavirus 2019-nCoV: early estimation of epidemiological parameters and epidemic predictions. 2020. DOI: [http://dx.doi.org/10.1101/2020.01.23.20018549](http://dx.doi.org/10.1101/2020.01.23.20018549) 12. 12.Kelso JK, Milne GJ, Kelly H, Simulation suggests that rapid activation of social distancing can arrest epidemic development due to a novel strain of influenza, BMC Public Health, 2009, 9, 117. 13. 13.Li 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. [https://doi.org/10.1056/NEJMoa2001316](https://doi.org/10.1056/NEJMoa2001316). 14. 14.Pueyo, T. (2020): Coronavirus: Why You Must Act Now. ([https://medium.com/@tomaspueyo/coronavirus-act-today-or-people-will-die-f4d3d9cd99ca](https://medium.com/@tomaspueyo/coronavirus-act-today-or-people-will-die-f4d3d9cd99ca)) updated March 19, 2020. 15. 15.Roosa, K., Lee, Y., Luo, R., Kirpich, A., Rothenberg, R., Hyman, J. M., Yan, P., Chowell, G. (2020): Real-time forecasts of the COVID-19 epidemic in China from February 5th to February 24th, 2020. Infect. Dis. Model 5, 256263. 16. 16.Shim E, Optimal strategies of social distancing and vaccination against seasonal influenza. Mathematical Biosciences and Engineering, 2013, 10, 5 and 6 doi: 10.3934/mbe.2013.10.1615 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3934/mbe.2013.10.1615&link_type=DOI) 17. 17.Tang, B., Bragazzi, N. L., Li, Q., Tang, S., Xiao, Y., Wu, J. (2020): An updated estimation of the risk of transmission of the novel coronavirus (2019-nCov). Infect. Dis. Model 5, 248255. 18. 18.Tang, B., Wang, X., Li, Q., Bragazzi, N. L., Tang, S., Xiao, Y., Wu, J. (2020): Estimation of the transmission risk of the 2019-nCoV and its implication for public health interventions. J. Clin. Med. 9(2), 462. 19. 19.Tao L, Transmission dynamics of 2019 novel coronavirus (2019-nCoV). 2020. DOI: [http://dx.doi.org/10.1101/2020.01.25.919787](http://dx.doi.org/10.1101/2020.01.25.919787). 20. 20.Thompson, R. N. (2020): Novel coronavirus outbreak in Wuhan, China, 2020: Intense surveillance Is vital for preventing sustained transmission in new locations. J. Clin. Med. 9(2), 498. 21. 21.World Health Organization. Coronavirus. World Health Organization, cited January 19, 2020. Available: [https://www.who.int/health-topics/coronavirus](https://www.who.int/health-topics/coronavirus). 22. 22.World Health Organization. Novel Coronavirus China. World Health Organization. cited January 20, 2020. Available: [https://www.who.int/csr/don/12-january-2020-novel-coronavirus-china/en/](https://www.who.int/csr/don/12-january-2020-novel-coronavirus-china/en/). Accessed 20 Jan 2020. 23. 23.World Health Organization. Novel Coronavirus Thailand (ex-China). World Health Organization. cited January 20, 2020. Available: [https://www.who.int/csr/don/14-january-2020-novel-coronavirus-thailand/en/](https://www.who.int/csr/don/14-january-2020-novel-coronavirus-thailand/en/). 24. 24.World Health Organization. Novel Coronavirus Japan (ex-China). World Health Organization. cited January 20, 2020. Available: [https://www.who.int/csr/don/17-january-2020-novel-coronavirus-japan-ex-china/en/](https://www.who.int/csr/don/17-january-2020-novel-coronavirus-japan-ex-china/en/). 25. 25.Worldometer ([https://www.worldometers.info/coronavirus/](https://www.worldometers.info/coronavirus/)), Accessed March 20, 2020. 26. 26.Wu JT, Leung K, Leung GM, (2020): Nowcasting and forecasting the potential domestic and international spread of the 2019-nCoV outbreak originating inWuhan, China: a modelling study. The Lancet 395(10225), 689697. 27. 27.Yi B, Chen Y, Ma X, Rui J, Cui JA, Wang H, et al. Incidence dynamics and investigation of key interventions in a dengue outbreak in Ningbo City, China. PLoS Negl Trop Dis. 2019;13:e0007659. 28. 28.Zhou P, Yang XL, Wang XG, Hu B, Zhang L, Zhang W, et al. A pneumonia outbreak associated with a new coronavirus of probable bat origin. Nature. 2020. [https://doi.org/10.1038/s41586-020-2012-7](https://doi.org/10.1038/s41586-020-2012-7). 29. 29.Zhu N, Zhang D, Wang W, Li X, Yang B, Song J, et al. A novel coronavirus from patients with pneumonia in China, 2019. N Engl J Med. 2020. [https://doi.org/10.1056/NEJMoa2001017](https://doi.org/10.1056/NEJMoa2001017). 30. 30.Zhao S, Lin Q, Ran J, Musa SS, 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. [https://doi.org/10.1016/j.ijid.2020.01.050](https://doi.org/10.1016/j.ijid.2020.01.050). 31. 31.Zhao S, Musa SS, 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. [https://doi.org/10.3390/jcm9020388](https://doi.org/10.3390/jcm9020388). 32. 32.Zhang S, Hu Q, Deng Z, Hu S, Liu F, Yu S, et al. Transmissibility of acute haemorrhagic conjunctivitis in small-scale outbreaks in Hunan Province, China. Sci Rep. 2020;10:119. [1]: /embed/inline-graphic-1.gif [2]: /embed/graphic-1.gif [3]: /embed/graphic-2.gif [4]: /embed/inline-graphic-2.gif [5]: /embed/inline-graphic-3.gif [6]: /embed/graphic-3.gif [7]: /embed/graphic-4.gif [8]: /embed/graphic-5.gif [9]: /embed/inline-graphic-4.gif [10]: /embed/inline-graphic-5.gif [11]: /embed/inline-graphic-6.gif [12]: /embed/graphic-6.gif [13]: /embed/graphic-7.gif [14]: /embed/graphic-8.gif [15]: /embed/inline-graphic-7.gif [16]: /embed/graphic-9.gif [17]: /embed/graphic-10.gif [18]: /embed/graphic-11.gif [19]: /embed/graphic-12.gif [20]: /embed/inline-graphic-8.gif [21]: /embed/graphic-13.gif [22]: /embed/inline-graphic-9.gif [23]: /embed/graphic-14.gif [24]: /embed/graphic-15.gif [25]: /embed/graphic-16.gif [26]: /embed/graphic-17.gif [27]: /embed/inline-graphic-10.gif [28]: /embed/inline-graphic-11.gif [29]: /embed/inline-graphic-12.gif [30]: /embed/inline-graphic-13.gif [31]: /embed/inline-graphic-14.gif [32]: /embed/inline-graphic-15.gif [33]: /embed/graphic-18.gif [34]: /embed/graphic-20.gif [35]: /embed/graphic-21.gif