Forecast analysis of the epidemics trend of COVID-19 in the United States by a generalized fractional-order SEIR model ====================================================================================================================== * Conghui Xu * Yongguang Yu * YangQuan Chen * Zhenzhen Lu ## Abstract In this paper, a generalized fractional-order SEIR model is proposed, denoted by SEIQRP model, which has a basic guiding significance for the prediction of the possible outbreak of infectious diseases like COVID-19 and other insect diseases in the future. Firstly, some qualitative properties of the model are analyzed. The basic reproduction number *R* is derived. When *R* < 1, the disease-free equilibrium point is unique and locally asymptotically stable. When *R* > 1, the endemic equilibrium point is also unique. Furthermore, some conditions are established to ensure the local asymptotic stability of disease-free and endemic equilibrium points. The trend of COVID-19 spread in the United States is predicted. Considering the influence of the individual behavior and government mitigation measurement, a modified SEIQRP model is proposed, defined as SEIQRPD model. According to the real data of the United States, it is found that our improved model has a better prediction ability for the epidemic trend in the next two weeks. Hence, the epidemic trend of the United States in the next two weeks is investigated, and the peak of isolated cases are predicted. The modified SEIQRP model successfully capture the development process of COVID-19, which provides an important reference for understanding the trend of the outbreak. Keywords * COVID-19 * Fractional order * Generalized SEIR model * Epidemic * Peak prediction ## 1 Introduction The outbreak of the coronavirus disease (COVID-19) in 2019 occurred in Wuhan, China, at the end of 2019. This is a severe respiratory syndrome caused by the novel coronavirus of zoonotic origin [1]. The Chinese government has implemented many measures, including the establishment of specialized hospitals and restrictions on travel, to reduce the spread. By April 20, 2020, the outbreak in China has been basically controlled. However, the outbreak is still rampant all over the world. At present, the United States, Italy, S-pain and other countries are still in the rising stage of the outbreak. It has posed a great threat to the public health and safety of the world. At present, many countries have adopted mitigation measures to restrict travel and public gatherings, which have a serious impact on the economy. Therefore, it is very important to predict the development trend of this epidemic and estimate the peak of the isolated cases. The epidemic model is a basic tool to research the dynamic behaviours of disease and predict the spreading trend of disease. Establishing a reasonable epidemic model can effectively characterize the development process of the disease. Therefore, a generalized SEIR epidemic model is proposed in this paper, and is defined as SEIQRP model. Some qualitative properties of this model are first analyzed, including the existence and uniqueness of the disease-free and the endemic equilibrium points. Then, conditions are also established to ensure the local asymptotic stability of both disease-free and endemic equilibrium points. So far, many scholars have researched COVID-19 from different perspectives [1–5, 23]. In [6], the epidemics trend of COVID-19 in China was predicted under public health interventions. In [7], the basic reproduction number of COVID-19 in China was estimated and the data-driven analysis was performed in the early phase of the outbreak. In order to predict the cumulative number of confirmed cases and combine the actual measures taken by the government on the outbreak of COVID-19, we further put forward an improved SEIQRP model, denoted by SEIQRPD model. At present, we can obtain the epidemic data of COVID-19 outbreak in the U-nited States before April 20, 2020. We use these data and the improved model to predict the epidemic trend of the U-nited States in the next two weeks, and estimate the peak of isolated cases. Firstly, the data before April 5 was selected to identify the model parameters, and the prediction ability of the improved model with the epidemiological data from April 6 to April 20 is verified. Then, the cumulative number of confirmed cases and isolated cases are predicted in the next two weeks. The peak of isolated cases is thus predicted. In recent years, with the continuous development of fractional calculus theory, fractional order systems modeling approaches have been applied in various engineering and non-engineering fields [19–21]. Compared with the short memory of the integer derivative, the fractional derivative has the information of the whole time interval or long memory. It is more accurate to describe the biological behavior of population by using fractional differential equation. The rest of this manuscript is structured as below. In Section 2, the fractional SEIQRP model and the modified SEIQRP model is proposed. Some qualitative properties of the SEIQRP model are discussed in Section 3. In Section 4, the prediction ability of the modified SEIQRP model is verified by using real data. The disease development in the United States in the next two weeks after April 21, 2020 is predicted. Finally, conclusion and some future works are discussed in Section 5. ## 2 Preliminaries and Model Derivation ### 2.1 Preliminaries In this section, some useful lemmas and definitions will be given to analyze some results of this paper. ### Definition 1 [8] The Caputo fractional order derivative is given below ![Formula][1] where *n* − 1 ≤ *α* ≤ *n*. In order to simplify the symbolic representation, the Caputo fractional differential operator ![Graphic][2] in this paper are represented by *Dα*. ### Definition 2 [8, 22] The Mittag-Leffler function is given below ![Formula][3] where *n* − 1 ≤ *α* ≤ *n*. Consider the following *n*-dimensional fractional order differential equation system ![Formula][4] where *α* ∈ (0,1), *X*(*t*) = (*x*1(*t*),*x*2(*t*)*,…,xn*(*t*))*T* is an *n* dimensional state vector and *A* is an *n* × *n* constant matrix. ### Lemma 1 *[9] For the corresponding linear time-invariant system (1), the following results are true:* *(i) The zero solution is asymptotically stable, if and only if all eigenvalues sj* (*j* = 1,2,…,*n*) *of A satisfy* ![Graphic][5]. *(ii) The zero solution is stable, if and only if all eigenvalues sj of A satisfy ![Graphic][6] and eigenvalues with ![Graphic][7]![Graphic][8] have the same algebraic multiplicity and geometric multiplicity*. ### 2.2 SEIQRP model The outbreak of COVID-19 has had a great impact on the economic growth of any country and daily life of any human. In order to control and prevent the possible outbreak of infectious diseases like the COVID-19 or other insect diseases in the future, it is very important to establish an appropriate model. The transmission diagram of the generalized SEIR model proposed in this paper is shown in Figure 1. We divide the total population into six distinct epidemic classes: susceptible, exposed, infectious, quarantined, recovered and insusceptible. We will represent the number of individuals at time *t* in the above classes by *S*(*t*), *E*(*t*), *I*(*t*), *Q*(*t*), *R*(*t*) and *P*(*t*), respectively. The specific explanations of the above six categories are as follows: * Susceptible *S*(*t*): the number of uninfected individuals at the time *t*. * Exposed *E*(*t*): the number of infected individuals at the time *t* but still in incubation period (without clinical symptoms and low infectivity). * Infectious *I*(*t*): the number of infected individuals at the time *t* (with obvious clinical symptoms). * Quarantined *Q*(*t*): the number of individuals who have been diagnosed and isolated at the time *t*. * Recovered *R*(*t*): the number of recovered individuals at the time *t*. * Insusceptible *P*(*t*): the number of susceptible individuals who are not exposed to the external environment at the time *t*. ![Fig. 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/03/2020.04.24.20078493/F1.medium.gif) [Fig. 1:](http://medrxiv.org/content/early/2020/05/03/2020.04.24.20078493/F1) Fig. 1: Transmission diagram for the model involving six population classes. The incidence rate plays a very important role in the epidemic, which can describe the evolution of infectious disease. According to the spread of different diseases in different regions, there are many forms of incidence rate. [10–12]. We expect to establish a generalized incidence rate, which can contain most of these specific forms. Therefore, the SEIQRP model with general incidence rate is given by ![Formula][9] where *Λ* is the inflow number of susceptible individuals, *β*1 and *β*2 denote the infection rates of the infected individuals and the exposed individuals, respectively, *ρ* is the protection rate, *ε* represent the incubation rate. *δ* is the rate at which symptomatic infections are diagnosed and quarantined. *λ* and *κ* represent the cure rate of isolated individuals and the death rate caused by the disease, respectively, *µ* is the natural mortality, the incidence rate *β*1 *f*1(*S*)*g*1(*I*), *β*2 *f*2(*S*)*g*2(*E*) are used to describe the transmission of diseases, which satisfy the following conditions [16]: 1. *f*1(0) = *f*2(0) = *g*1(0) = *g*2(0) = 0, 2. *f*1(*S*) *>* 0, *f*2(*S*) *>* 0, *g*1(*I*) *>* 0, *g*2(*E*) *>* 0, for any *S*, *E*, *I >* 0, 3. ![Graphic][10], for any *S*, *E*, *I >* 0, 4. ![Graphic][11], for any *E*, *I >* 0. ### 2.3 Modified SEIQRP model In the following, we will combine the actual situation and government mitigation policy to improve the SEIQRP model (2). Firstly, during the outbreak of COVID-19, we need to make reasonable assumptions according to the actual mitigation policies and circumstances. * During the COVID-19 outbreak, population mobility is strictly controlled by many countries. Especially, the policy of city closure was implemented in Hubei province, China. Therefore, the impact of migration will not be considered in the improved model. * For the prediction of this short-term virus transmission, the impact of natural mortality will not be considered. It is assumed that novel coronavirus is the only cause of death during the outbreak. * In order to predict the trend of cumulative confirmed cases, it is necessary to simulate the number of death cases. Therefore, a new class *D*(*t*) to the SEIQRP model will be added in the modified model, which denotes the number of death cases at time *t*. * The time-varying cure rate, mortality rate and infection rate will be applied to the improved model. This can better simulate the impact of the improvement of medical conditions and government control on individuals in reality. Based on the above assumptions and analysis, the following improved model for COVID-19 is proposed. ![Formula][12] where *β*1(*t*) = *σ*1 exp(*−σ*2*t*), *λ*(*t*) = *λ*1(1 *−* exp(*λ*2*t*)) and *κ*(*t*) = *κ*1 exp(*−κ*2*t*). The parameters *σ*1, *σ*2, *λ*1, *λ*2, *κ*1 and *κ*2 are all positive, where *σ* is used to simulate the intensity of government control. It should be emphasized that the protection rate *ρ* for susceptible individuals also reflects the intensity of government control. ## 3 Qualitative Analysis of the SEIQRP Model ### 3.1 The Existence and Uniqueness of Equilibrium Point Obviously, the right side of system (2) satisfies the local Lipschitz condition, then there exists a unique solution of system (2) [8,15]. ### Theorem 1 *The solutions of system (2) are non-negative, and the closed set ![Graphic][13]is a positive invariant set of system (2)*. ### Proof In order to investigate the non-negativity of solutions of the system (2), we consider the following system ![Formula][14] where the initial conditions are *S*10 = 0, *E*10 = 0, *I*10 = 0, *Q*10 = 0, *R*10 = 0 and *P*10 = 0. (0,0,0,0,0,0) is the unique solution of the above system. According to the fractional-order comparison theorem [17], one can deduce that the solutions of system (2) satisfy *S*(*t*) ≥ 0, *E*(*t*) ≥ 0, *I*(*t*) ≥ 0, *Q*(*t*) ≥ 0, *R*(*t*) ≥ 0 and *P*(*t*) ≥ 0. By adding six equations of system (2), one can deduce ![Formula][15] By applying the fractional-order comparison theorem, one has ![Formula][16] When ![Graphic][17], since *Eα*(*−µtα*) ≥ 0, we have ![Formula][18] Thus we can draw the result of Theorem 1. System (2) has an obvious disease-free equilibrium point *M* = (*S*,0,0,0,0*,P*), where ![Formula][19] In order to obtain the endemic equilibrium point of the system (2), we set: ![Formula][20] which implies ![Formula][21] Combining the above equations and the second equation of (5), one has ![Formula][22] Define ![Formula][23] Note that *φ*(0) = 0 and ![Graphic][24]. In order to show that *φ*(*E*) = 0 has at least one positive root in the interval ![Graphic][25], we need to prove that *φ′*(0) *>* 0. Hence ![Formula][26] Therefore, we have ![Formula][27] where the basic reproduction number is given by ![Formula][28] If *R* *>* 1, the system (2) has at least one endemic equilibrium point *M**= (*S*,E*,I*,Q*,R*,P**), where ![Formula][29] In the following section, we show that endemic equilibrium point *M** is unique. According to the above analysis and hypothesis (i)-(iii), one has ![Formula][30] By the hypothesis (iv), this implies *φ′*(*E**) *<* 0. If there is another endemic equilibrium point *M***, then *φ′*(*E***) ≥ 0 holds, which contradicts the previous discussion. Hence, system (2) has a unique endemic equilibrium point *M** when *R* *>* 1. Based on the above analysis, the following result can be obtained. ### Theorem 2 *System (2) has a unique disease-free equilibrium point M**, if R* *<* 1*. System (2) has a unique endemic equilibrium point M*, if R* *>* 1. ### 3.2 Stability Analysis In this section, the local asymptotic stability of disease-free equilibrium point *M* and endemic equilibrium point *M** for system (2) is discussed. ### Theorem 3 *With regard to system (2), the disease-free equilibrium point M* *is locally asymptotic stability, if R*o < 1*; the disease-free equilibrium point M* *is unstable, if R*o > 1. ### Proof The Jacobian matrix of the system (2) at the disease-free equilibrium point *M* is ![Formula][31] where ![Formula][32] The corresponding characteristic equation is ![Formula][33] where ![Formula][34] The characteristic equation *H*(*s*) = 0 has four obvious negative characteristic roots, which are denoted by *s*1 = *s*2 = *−µ*, *s*3 = *−µ − ρ* and *s*4 = *−µ − λ − κ*, respectively. The discriminant of *H*1(*s*) in quadratic form is ![Formula][35] This implies that the other two eigenvalues *s*5 and *s*6 of characteristic equation *H*(*s*) = 0 are real roots. Hence ![Formula][36] If *R* *<* 1, then one can obtain *s*5 + *s*6 *<* 0 and *s*5*s*6 *>* 0, which imply *s*5 *<* 0 and *s*6 *<* 0. If *R* *>* 1, one has *s*5*s*6 *<* 0, which imply *s*5 *>* 0 or *s*6 *>* 0. It follows from Lemma 1 that the proof is completed. Further, we will show the locally asymptotic stability of the endemic equilibrium point *M** of system (2). Similarly, the corresponding Jacobian matrix of the system (2) at *M** is ![Formula][37] where ![Formula][38] with ![Formula][39] Hence, the corresponding characteristic equation is ![Formula][40] where ![Formula][41] with ![Formula][42] By the corresponding results in [13,14], let ![Formula][43] Then, the following result can be obtained. ### Theorem 4 *With regard to system (2), assume that R* *>* 1, *(i) If D*1(*L*1(*s*)) *>* 0*, a*1 *>* 0*, a*3 *>* 0 *and a*1*a*2 *> a*3*, then the endemic equilibrium point M* is locally asymptotically stable*. *(ii) If D*1(*L*1(*s*)) *<* 0*, a*1 *>* 0*, a*2 *>* 0 *and a*3 *>* 0*, then the endemic equilibrium point M* is locally asymptotically stable for ![Graphic][44]*. *(iii) If D*1(*L*1(*s*)) *<* 0*, a*1 *<* 0 *and a*2 *<* 0*, then the endemic equilibrium point M* is unstable for ![Graphic][45]*. *(iv) If D*1(*L*1(*s*)) *<* 0*, a*1 *>* 0*, a*2 *>* 0 *and a*1*a*2 = *a*3*, then for α* ∈ (0,1)*, the endemic equilibrium point M* is locally asymptotically stable*. ### Proof Based on the previous discussion, the characteristic equation *L*(*s*) = 0 has three obvious negative roots *s*1 = *s*2 = *− µ* and *s*3 = *−µ −λ −κ*. In order to investigate the stability of equilibrium point *M**, we only need to discuss the range of the root of *L*1(*s*) = 0, denoted by *s*4, *s*5 and *s*6. (i) By the results in [14], if *D*1(*L*1(*s*)) *>* 0, then *s*4, *s*5 and *s*6 are real roots. Further, by Routh-Hurwitz criterion, the necessary and sufficient conditions for *si*(*i* = 4,5,6) to lie in the left half plane are ![Formula][46] That is to say, under the above conditions, the roots of *L*1(*s*) = 0 satisfy ![Formula][47] Therefore, *M** is locally asymptotically stable and (i) holds (ii) By the results in [14], if *D*1(*L*1(*s*)) *<* 0, then *L*1(*s*) = 0 has a real root and a pair of conjugate complex roots, denoted by *s*4, *m* + *ni* and *m − ni*, respectively. Thus, one has ![Formula][48] By calculation, ![Formula][49] The conditions *a*1 *>* 0, *a*2 *>* 0 and *a*3 *>* 0 imply that ![Formula][50] Further, one has ![Formula][51] That is to say tan2 *q >* 3, which implies ![Formula][52] Therefore, in order to ensure the establishment of ![Graphic][53]![Graphic][54], we must have![Graphic][55]. Thus (ii) holds. The proof of conclusions (iii) and (iv) is similar to that of conclusion (ii), hence we omit it. ## 4 Numerical simulations ### 4.1 Data Sources The data used in this paper are from the Johns Hopkins University Center for Systems Science and Engineering ([https://github.com/CSSEGISandData/COVID-19](https://github.com/CSSEGISandData/COVID-19)). the Johns Hopkins University publishes data of accumulated and newly confirmed cases, recovered cases and death cases of COVID-19 from January 22, 2020. ### 4.2 Analysis of the SEIQRP Model In the following discussion, the standard incidence rate[18] is used to describe the transmission of COVID-19, and is given by ![Formula][56] where *N* represent the total population of the region at the initial time. Hence, ![Formula][57] The effectiveness of the model (15) in describing the spread of COVID-19 is illustrated by selecting the confirmed, cured and dead cases in the United States. According to the real data provided by the Johns Hopkins University, the outbreak in the United States has not been brought under full control. The data of confirmed cases, cured cases and death cases are selected from January 22, 2020, to April 20, 2020. Assuming that the confirmed individuals will be isolated, then ![Formula][58] This hypothesis is in line with the actual situation. Hence, we can obtain the real data of isolated cases. Through the fractional predictor corrector method and the least squares fitting [24], we can identify the parameters of the model (15) through the real data, which can be found in Table 1. View this table: [Table 1:](http://medrxiv.org/content/early/2020/05/03/2020.04.24.20078493/T1) Table 1: Summary table of the parameter identification for model (15) after using least squares fitting to real data from January 22, 2020, to April 20, 2020 ![Fig. 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/03/2020.04.24.20078493/F2.medium.gif) [Fig. 2:](http://medrxiv.org/content/early/2020/05/03/2020.04.24.20078493/F2) Fig. 2: Number of isolated cases predicted and recovered cases predicted by the model (15) for the United States. Based on the parameters in Table 1, we can make a simple prediction of isolated cases and recovered cases in the U-nited States, which can be found in Figure 1. We need to emphasize that the peak here represents the number of isolated cases rather than the cumulative number of confirmed cases. Under the parameters in Table 1, it can be calculated by using (9) that *R* = 0.2041 *<* 1. By the conclusion in Theorem 3, the disease-free equilibrium point is local asymptotic stable. The SEIQRP model has a basic guiding significance for predicting and fitting spreading dynamics of COVID-19. However, the prediction of this model for COVID-19 is relatively rough, we still need to improve model (15) according to actual mitigation policies and research objectives. According to the analysis in Section 2, we choose the SEIQRPD model to predict the trend of the epidemic in the United States under reasonable assumptions. ### 4.3 The SEIQRPD Model for the Prediction of COVID-19 Similarly, the standard incidence rate is used to describe the transmission of COVID-19, the fractional SEIQRPD model can be expressed as ![Formula][59] When *α* = 1, the fractional-order SEIQRPD model is similar to the integer-order model used in [24]. According to the data provided by Johns Hopkins University, by April 20, 2020, the outbreak in China has been basically controlled. In many provinces of China, the number of new cases per day is increasing in single digits. This means that the data in China contains more information about the spreading dynamics of COVID-19. Therefore, the data in Hubei, Guangdong, Hunan and Zhejiang are selected to research the fitting effect of the model (17). According to the real data of these four regions in China, the parameters of the model (17) are identified, and the results are shown in Table 2. The model (17) successfully capture the trend of the outbreak, which can be seen in Figure 3. View this table: [Table 2:](http://medrxiv.org/content/early/2020/05/03/2020.04.24.20078493/T2) Table 2: Summary of parameter identification of model (17) (data used from January 22, 2020, to April 20, 2020). ![Fig. 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/03/2020.04.24.20078493/F3.medium.gif) [Fig. 3:](http://medrxiv.org/content/early/2020/05/03/2020.04.24.20078493/F3) Fig. 3: The fitting effect of the improved model (17) on the outbreak in Hubei, Guangdong, Hunan and Zhejiang. At present, we can obtain the epidemic data of the U-nited States from January 22, 2020 to April 20, 2020. We need to preprocess the data to remove the data smaller than0.5% of the current highest number of confirmed cases. In order to test the prediction ability of the SEIQRPD model (17) for the development process of the epidemic in the United States, we select the data before April 5 to identify the parameters of the SEIQRPD model (17). Furthermore, in order to illustrate the ability of the SEIQRPD model (17) in predicting the outbreak, we compared the real data and fitted data after April 5, which can be seen in Table 3, Table 4 and Figure 4. According to the results in Table 3 and Table 4, it can be found that the real values of current isolated cases and cumulative confirmed cases fall within the range of 95% - 105% of the predicted values. This shows that the SEIQRPD model (17) can effectively predict the data in the next two weeks. View this table: [Table 3:](http://medrxiv.org/content/early/2020/05/03/2020.04.24.20078493/T3) Table 3: Summary of real and fitted data for the cumulative confirmed cases in the United States from April 6, 2020 to April 20, 2020. View this table: [Table 4:](http://medrxiv.org/content/early/2020/05/03/2020.04.24.20078493/T4) Table 4: Summary of real and fitted data for the isolated cases in the United States from April 6, 2020 to April 20, 2020. ![Fig. 4:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/03/2020.04.24.20078493/F4.medium.gif) [Fig. 4:](http://medrxiv.org/content/early/2020/05/03/2020.04.24.20078493/F4) Fig. 4: Based on the data of the United States from January 22 to April 20, 2020 to verify the accuracy of the forecast for the next 15 days. The data before April 20, 2020 are selected to identify the parameters of the improved model (17), and the results are shown in Table 5. The cumulative number of confirmed cases and the number of isolated cases after two weeks are predicted, which can be seen in Table 6. View this table: [Table 5:](http://medrxiv.org/content/early/2020/05/03/2020.04.24.20078493/T5) Table 5: Summary table of the parameter identification for model (17) (data used from January 22, 2020, to April 20, 2020) View this table: [Table 6:](http://medrxiv.org/content/early/2020/05/03/2020.04.24.20078493/T6) Table 6: Summary of predicted data for the United States from April 21, 2020 to May 5, 2020. The isolated cases in the United States will peak on June 18, 2020, with the peak of 1.0431 × 106, which can be seen in Figure 5. ![Fig. 5:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/03/2020.04.24.20078493/F5.medium.gif) [Fig. 5:](http://medrxiv.org/content/early/2020/05/03/2020.04.24.20078493/F5) Fig. 5: Based on the data of the United States from January 22 to April 20, 2020, to verify the accuracy of the forecast for the next 15 days. ## 5 Conclusion We first propose the fractional SEIQRP model with generalized incidence rates. Some qualitative properties of the SEIQRP model are discussed. In order to predict COVID-19 effectively, we propose an improved SEIQRPD model according to the actual mitigation situation. According to the data of the United States before April 5, 2020, the trend of the outbreak in the United States from April 6 to April 20 is successfully predicted as compared to the real records. Then, using the data before April 20, 2020, we forecast the trend of the outbreak in the United States in the next two weeks, and estimate the peak of isolated cases and the date of the peak. The improved SEIQRP model proposed in this paper successfully captures the trend of COVID-19. The long-term prediction needs to adjust the model appropriately according to the change of policy and medical level. We will discuss in the future work. ## Data Availability our results are reproducible. Our data sets are from JHU COVID-19 database referenced in the paper. ## Conflict of Interest We declare that we have no conflict of interest. ## Acknowledgment The plots in this paper were plotted using the plot code adapted from [24]. ## Footnotes * * E-mail: ygyu{at}bjtu.edu.cn * This work is supported by the National Nature Science Foundation of China under Grant 61772063, Beijing Natural Science Foundation under Grant Z180005. * Received April 24, 2020. * Revision received April 29, 2020. * Accepted May 3, 2020. * © 2020, Posted by Cold Spring Harbor Laboratory The copyright holder for this pre-print is the author. All rights reserved. The material may not be redistributed, re-used or adapted without the author's permission. ## References 1. 1.Lin, Q.Y., Zhao, S. Gao, D.Z., Luo, Y.J., Yang, S., Musa. S.S., Wang, M.H., Cai, Y.L., Wang, W.M., Yang, L., He, D.h.: A conceptual model for the coronavirus disease 2019 (COVID-19) outbreak in Wuhan, China with individual reaction and governmental action. International Journal of Infectious Diseases. 93, 211–216 (2020). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ijid.2020.02.058&link_type=DOI) 2. 2.Peng, L., Yang, W., Zhang, D., Zhuge, C., Hong, L.: Epidemic analysis of COVID-19 in China by dynamical modeling. Cold Spring Harbor Laboratory. (2020). arXiv: 2002.06563. 3. 3.Amjad, S. S., Iqbal, N. S., Kottakkaran, S.N.: A Mathematical model of COVID-19 using fractional derivative: outbreak in India with dynamics of transmission and control. (2020). doi:10.20944/preprints202004.0140.v1. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.20944/preprints202004.0140.v1&link_type=DOI) 4. 4.Chen, Y., Cheng, J., Jiang, X., Xu, X.: The reconstruction and prediction algorithm of the fractional TDD for the local outbreak of COVID-19. (2020). arXiv:2002.10302. 5. 5.Cheng, Z.J., Shan, J.: 2019..Cnovel Coronavirus: Where We are and What We Know. Infection. (2020). DOI 10.1007/s15010-020-01401-y. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s15010-020-01401-y&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32072569&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F03%2F2020.04.24.20078493.atom) 6. 6.Yang, Z.F., Zeng, Z.Q., Wang, K., Wong, S., Liang, W.H., Zanin, M., Liu, P., Cao, X.D., Gao, Z.Q., Mai, Z.T., Liang, J.Y., Liu, X.Q., Li, S.Y., Li, Y.M., Ye, F., Guan, W.J., Yang, Y.F., Li, F., Luo, S.M., Xie, Y.Q., Liu, B., Wang, Z.L., Zhang, S.B., Wang, Y.N., Zhong, N.S., He, J.X.: Modified SEIR and AI prediction of the epidemics trend of COVID-19 in China under public health interventions. Journal of Thoracic Disease 12(2), 165–174 (2020). 7. 7.Zhao, S., Lin, Q, Ran J., Musa, S.S., Yang, G.P., Wang, W.M., Lou, Y.J., Gao, D.Z., Yang, L., He, D.H., Wang, M.H.: 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. 214–217 (2020). 8. 8.Podlubny, I.: Fractional differential equations, Academic Press, New York, (1999). 9. 9.Diethelm, K.: The analysis of fractional differential equations: an application-oriented exposition using differential operators of Caputo type, Springer Science, Business Media, (2010). 10. 10.Kuniya, T.: Hopf bifurcation in an age-structured SIR epidemic model. Applied Mathematics Letters. 92, 22–28 (2019). 11. 11.Zhang, X.B., Huo, H., Xiang, H.F., Xiang, H., Meng, X.Y.: An SIRS epidemic model with pulse vaccination and non-monotonic incidence rate. Nonlinear Analysis: Hybrid Systems. 13–21 (2013). 12. 12.Cai, Y., Kang, Y., Wang, W.: A stochastic SIRS epidemic model with nonlinear incidence rate. Applied Mathematics and Computation. 221–240 (2017). 13. 13.Ahmed, E., Elgazzar, A.S.: On fractional order differential equations model for nonlocal epidemics. Physica A. 379(2), 607–614 (2007). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.physa.2007.01.010&link_type=DOI) 14. 14.Rocca, A., West, B.J.: Fractional calculus and the evolution of fractal phenomena. Physica A. 265(3–4), 535–546 (1999). 15. 15.Jalilian, Y., Jalilian, R.: Existence of solution for delay fractional differential equations. Mediterranean Journal of Mathematics. 10(4), 1731–1747 (2013). 16. 16.Yang, Y., Xu, L.: Stability of a fractional order SEIR model with general incidence. Applied Mathematics Letters, 105, (2020). 17. 17.Wang, H., Yu, Y.G., Wen, G.G., Zhang, S., Yu, J.Z.: Global stability analysis of fractional-order Hopfield neural networks with time delay. Neurocomputing. 154, 15–23 (2015). 18. 18.Hu, Z., Liu, S., Wang, H.: Backward bifurcation of an epidemic model with standard incidence rate and treatment rate. Nonlinear Analysis-real World Applications. 9(5), 2302–2312 (2008). 19. 19.Sun, H.G., Zhang, Y., Baleanu, D., Chen, W., Chen, Y.Q.: A new collection of real world applications of fractional calculus in science and engineering. Communications in Nonlinear Science and Numerical Simulation. 59(5), 213–231 (2018). 20. 20.Cao, K.C., Chen, Y.Q.: Fractional order crowd dynamics: Cyber-Human systems modeling and control. (Invited book project. Volume 4 of the De Gryuter Monograph Series ‘Fractional Calculus in Applied Sciences and Engineering‘). ISBN 978-3-11-047398-8. 21. 21.West, B.J.: Fractional Calculus View of Complexity: Tomorrow’s Science, CRC Press, (2015). 22. 22.Li, Y., Chen, Y.Q, Podlubny, I.: Mittag-Leffler stability of fractional order nonlinear dynamic systems. Automatica. 45(8), 1965–1969 (2009). 23. 23.Lu, Z.Z., Yu, Y.G., Chen, Y.Q., Ren, G.J., Xu, C.H., Yin, Z.: A fractional-order SEIHDR model for COVID-19 with inter-city networked coupling effects. Nonlinear Dynamics (Special Issue on ‘Nonlinear dynamics of COVID-19 pandemic: modeling, control, and future perspectives‘). (April 2020). 24. 24.Cheynet, E.: Generalized SEIR epidemic model (fitting and computation). ([https://www.github.com/ECheynet/SEIR](https://www.github.com/ECheynet/SEIR)), GitHub. [1]: /embed/graphic-1.gif [2]: /embed/inline-graphic-1.gif [3]: /embed/graphic-2.gif [4]: /embed/graphic-3.gif [5]: /embed/inline-graphic-2.gif [6]: /embed/inline-graphic-3.gif [7]: /embed/inline-graphic-4.gif [8]: /embed/inline-graphic-5.gif [9]: /embed/graphic-5.gif [10]: /embed/inline-graphic-6.gif [11]: /embed/inline-graphic-7.gif [12]: /embed/graphic-6.gif [13]: /embed/inline-graphic-8.gif [14]: /embed/graphic-7.gif [15]: /embed/graphic-8.gif [16]: /embed/graphic-9.gif [17]: /embed/inline-graphic-9.gif [18]: /embed/graphic-10.gif [19]: /embed/graphic-11.gif [20]: /embed/graphic-12.gif [21]: /embed/graphic-13.gif [22]: /embed/graphic-14.gif [23]: /embed/graphic-15.gif [24]: /embed/inline-graphic-10.gif [25]: /embed/inline-graphic-11.gif [26]: /embed/graphic-16.gif [27]: /embed/graphic-17.gif [28]: /embed/graphic-18.gif [29]: /embed/graphic-19.gif [30]: /embed/graphic-20.gif [31]: /embed/graphic-21.gif [32]: /embed/graphic-22.gif [33]: /embed/graphic-23.gif [34]: /embed/graphic-24.gif [35]: /embed/graphic-25.gif [36]: /embed/graphic-26.gif [37]: /embed/graphic-27.gif [38]: /embed/graphic-28.gif [39]: /embed/graphic-29.gif [40]: /embed/graphic-30.gif [41]: /embed/graphic-31.gif [42]: /embed/graphic-32.gif [43]: /embed/graphic-33.gif [44]: /embed/inline-graphic-12.gif [45]: /embed/inline-graphic-13.gif [46]: /embed/graphic-34.gif [47]: /embed/graphic-35.gif [48]: /embed/graphic-36.gif [49]: /embed/graphic-37.gif [50]: /embed/graphic-38.gif [51]: /embed/graphic-39.gif [52]: /embed/graphic-40.gif [53]: /embed/inline-graphic-14.gif [54]: /embed/inline-graphic-15.gif [55]: /embed/inline-graphic-16.gif [56]: /embed/graphic-41.gif [57]: /embed/graphic-42.gif [58]: /embed/graphic-43.gif [59]: /embed/graphic-46.gif