Next Article in Journal
On a New Modification of the Erdélyi–Kober Fractional Derivative
Previous Article in Journal
Enhancement in Thermal Energy and Solute Particles Using Hybrid Nanoparticles by Engaging Activation Energy and Chemical Reaction over a Parabolic Surface via Finite Element Approach
Previous Article in Special Issue
Numerical Solution of Nonlinear Fractional Diffusion Equation in Framework of the Yang–Abdel–Cattani Derivative Operator
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multi-Model Selection and Analysis for COVID-19

School of Mathematics and Computer Science, Northwest Minzu University, Lanzhou 730000, China
*
Author to whom correspondence should be addressed.
Fractal Fract. 2021, 5(3), 120; https://doi.org/10.3390/fractalfract5030120
Submission received: 18 August 2021 / Revised: 7 September 2021 / Accepted: 10 September 2021 / Published: 13 September 2021

Abstract

:
In the face of an increasing number of COVID-19 infections, one of the most crucial and challenging problems is to pick out the most reasonable and reliable models. Based on the COVID-19 data of four typical cities/provinces in China, integer-order and fractional SIR, SEIR, SEIR-Q, SEIR-QD, and SEIR-AHQ models are systematically analyzed by the AICc, BIC, RMSE, and R means. Through extensive simulation and comprehensive comparison, we show that the fractional models perform much better than the corresponding integer-order models in representing the epidemiological information contained in the real data. It is further revealed that the inflection point plays a vital role in the prediction. Moreover, the basic reproduction numbers R 0 of all models are highly dependent on the contact rate.

1. Introduction

Since December 2019, a virus named COVID-19 was first confirmed in Wuhan, China and has quickly spread throughout the world. As of 25 June 2021, the WHO reported about 179 million confirmed cases, including more than 3.8 million deaths [1]. In order to control the outbreak, it is crucial to predict the spreading trend, analyze the key epidemiological parameters, and simulate the results of various control measures. All actions of individuals and governments depend heavily on predictions of how many people will be infected. Therefore, many models have been reported to accurately estimate and predict the epidemic trend, such as empirical functions, statistical methods, dynamical and stochastic models [2,3,4]. Among them, integer-order and fractional dynamical models can better describe the dynamic behavior of disease transmission, which has attracted attention widely [5,6,7].
According to the importance of mathematical modelling, integer-order infectious disease models have been widely studied for its simplicity and accessibility. Based on the compartment hypothesis of the population, most models include compartments of susceptible (S), exposed (E), infected (I), recovered (R), death (D), quarantined (Q), asymptomatic (A), insusceptible (P), and hospitalized (H) individuals, such as SIR [8], SEIR [9], SEIR-Q [10], SEIR-QD [11], and SEIR-AHQ [12] models. In addition, many sophisticated models also have been proposed [13,14,15,16].
In the integer-order models, it is usually assumed that the infection rate, conversion rate, and removal rate are constants, which lead to the system states that are independent of history and have no memory. However, it is well known that the evolution, spread, and control of infectious diseases are closely related to the historical states, so the diseases are non-Markov processes. The more distant experience and knowledge have less influence on the present states. Thus, the growth rate of the state is not exponential, but it is more like an inverse power-law function, which can be accurately described by a fractional model. The fractional model is a generalization of the classical integer-order model. On the other hand, fractional calculus proved to be more effective in modeling dynamics with long memory and hereditary properties. The memory property of disease is embodied in integrating more information from the past. In addition, the hereditary property reflects in the genetic profile along with evolution of the disease. All these lead to the fractional model having more advantages in modeling some biological systems than the classical integer-order. Therefore, many fractional-order models are applied to characterize memory-dependent phenomena for COVID-19, such as fractional SEIRD [17], fractional SEIQR [18], fractional SEIQRDP [19,20], fractional SEIHDR [21], and so on.
The main purposes of modeling for COVID-19 are analysis transmission trends and biological characteristics. Many dynamical models, including integer-order and fractional types, have been proposed to predict the epidemic trend. However, different models have given vastly different results of Wuhan [22,23,24,25]: the values of the basic reproduction number from 1.5 to 6, the peak times from mid-February to late March, and the total number of infections ranges from 20,000 to millions. Knowing which model is more reliable is a pending issue.
In the face of so many integer-order and fractional dynamical models, which one would be our best choice to estimate the epidemic trend and epidemiological features for COVID-19? In this paper, this pivotal problem is systematically explored by comparing the forecast ability of five models based on the four classic model selection methods. We consider integer-order and fractional SIR, SEIR, SEIR-Q, SEIR-QD, and SEIR-AHQ models for selection and analysis. Through extensive simulation and comprehensive comparison, we obtain many meaningful results, which really deepen our understanding of disease.
The remainder of the paper is arranged as follows: Section 2 introduces the integer-order and fractional models. In Section 3, four means are given for multi-model selection and analysis. Section 4 carries out comparisons and evaluations. Finally, some concluding remarks are given in Section 5.

2. Integer-Order and Fractional Models

We collect and study some typical integer-order models such as SIR, SEIR, SEIR-Q, SEIR-QD, and SEIR-AHQ models to find out which is the best model for COVID-19. Fractional models are constructed from the corresponding integer-order models. In these models, S represents the susceptible population, I represents the infectious population, E represents the individuals who are infected but in the latent period, P represents the isolated susceptible population, Q represents the individuals who are confirmed and infected, R denotes recovered cases, D denotes death cases, A denotes asymptomatic cases, and H denotes hospitalized cases. In the SIR and SEIR models, the parameters α ,   γ 1 , and β denote the infection rate, the average latent time and the cure rate, separately. In the SEIR-Q model, the parameters α 1 ,   α 2 ,   ε ,   γ ,   κ ,   δ ,   λ , and η denote the asymptomatic transmission rate, the symptomatic transmission rate, the onset rate, the recovery rate of infected, the testing rate, the sensitivity of test, the specificity of test, and the recovery rate of the quarantined population, respectively. In the SEIR-QD model, the parameters δ ,   α ,     γ 1 ,     β 1 ,   λ , and κ denote the protection rate, the infection rate, the average latent time, the average quarantine time, the cure rate, and the mortality rate, separately. In the SEIR-AHQ model, the parameters α ,   ε ,   θ ,   σ ,   μ ,   δ 1 ,     δ 2 ,   λ ,   e ,   ω ,   β 1 ,   β 2 , and β 3 denote the transmission rate of contact, the quarantined rate, the transmissibility ratio between asymptomatic and symptomatic, the transition rate of exposed, the probability of having symptoms among infected, the transition rate of symptomatic infected, the transition rate of quarantined exposed, the rate at which the quarantined uninfected contacts are released into the wider community, the contact rate, the disease-induced death rate, the recovery rate of asymptomatic, the recovery rate of the symptomatic and the recovery rate of quarantined, respectively.
Fractional calculus is a strong tool to characterize the memory effects of the disease. The Caputo fractional derivative contains a power-law function in the convolution kernel, which can reflect the fact that the contribution of the earlier state is significantly smaller than the current state. Based on the integer-order models and the Caputo fractional derivative method [26,27], the fractional models are given as follows:
(1) Fractional SIR model
D 0 , t q C S ( t ) = α S ( t ) I ( t ) / N , D 0 , t q C I ( t ) = α S ( t ) I ( t ) / N β I ( t ) , D 0 , t q C R ( t ) = β I ( t ) .
(2) Fractional SEIR model
D 0 , t q C S ( t ) = α S ( t ) I ( t ) / N , D 0 , t q C E ( t ) = α S ( t ) I ( t ) / N γ E ( t ) , D 0 , t q C I ( t ) = γ E ( t ) β I ( t ) , D 0 , t q C R ( t ) = β I ( t ) .
(3) Fractional SEIR-Q model
D 0 , t q C S ( t ) = ( α 1 E ( t ) + α 2 I ( t ) ) S ( t ) / N κ ( 1 λ ) S ( t ) , D 0 , t q C E ( t ) = ( α 1 E ( t ) + α 2 I ( t ) ) S ( t ) / N ( ε + κ δ ) E ( t ) , D 0 , t q C Q ( t ) = κ ( 1 λ ) S ( t ) + κ δ ( E ( t ) + I ( t ) ) η Q ( t ) , D 0 , t q C I ( t ) = ε E ( t ) ( γ + κ δ ) I ( t ) , D 0 , t q C R ( t ) = γ I ( t ) + η Q ( t ) .
(4) Fractional SEIR-QD model
D 0 , t q C S ( t ) = α S ( t ) I ( t ) / N δ S ( t ) , D 0 , t q C E ( t ) = α S ( t ) I ( t ) / N γ E ( t ) , D 0 , t q C I ( t ) = γ E ( t ) β I ( t ) , D 0 , t q C Q ( t ) = β I ( t ) λ Q ( t ) κ Q ( t ) , D 0 , t q C R ( t ) = λ Q ( t ) , D 0 , t q C D ( t ) = κ Q ( t ) , D 0 , t q C P ( t ) = δ S ( t ) .
(5) Fractional SEIR-AHQ model
D 0 , t q C S ( t ) = ( α e + e ε ( 1 α ) ) S ( t ) ( I ( t ) + θ A ( t ) ) / N + λ P ( t ) , D 0 , t q C E ( t ) = α e ( 1 ε ) S ( t ) ( I ( t ) + θ A ( t ) ) / N σ E ( t ) , D 0 , t q C P ( t ) = e ε ( 1 α ) ( I ( t ) + θ A ( t ) ) λ P ( t ) , D 0 , t q C H ( t ) = δ 1 I ( t ) + δ 2 Q ( t ) ( ω + β 3 ) H ( t ) , D 0 , t q C Q ( t ) = e ε α ( I ( t ) + θ A ( t ) ) δ 2 Q ( t ) , D 0 , t q C I ( t ) = σ μ E ( t ) ( ω + γ 1 + δ 1 ) I ( t ) , D 0 , t q C R ( t ) = β 1 I ( t ) + β 2 A ( t ) + β 3 H ( t ) , D 0 , t q C A ( t ) = σ ( 1 μ ) E ( t ) β 1 I ( t ) .
The D 0 , t q C is the Caputo fractional derivative, which is defined as:
D 0 , t q C x ( t ) = 1 Γ ( n q ) 0 t ( t τ ) n q 1 x ( n ) ( τ ) d τ ,  
where n 1 < q n Z + .
The transfer diagrams for five models are given in Figure 1. When q = 1 , the Caputo fractional derivative reduces to the classic integer-order derivative, so the fractional SIR, SEIR, SEIR-Q, SEIR-QD, and SEIR-AHQ models become corresponding integer-order models.

3. Four Means for Model Evaluation and Analysis

In order to be closer to the disease transmission, more parameters are added, which makes the model structure more complicated. This can easily lead to the models overfitting the real data. It is not trivial to rationally evaluate various dynamical models. Furthermore, many different aspects should be considered at the same time [28], such as complexity and accuracy, fitting and prediction, robustness and sensitivity, and so on. Overemphasizing one aspect and neglecting another, consciously or unconsciously, would lead to unsatisfactory forecasts. We further introduce four means for model selection to solve this problem reasonably.

3.1. The Corrected Akaike Information Criterion (AICc)

The AIC was first proposed by Japanese statistician Akaike and is now a measuring standard of the degree of fitting for models. Based on the concept of entropy, the AIC can balance the model complexity and accuracy of the data fitted by the model. When K > 40 / N , that is too many estimated parameters (K) compare to the size of the data points (N), the following corrected AIC [29] can be used:
A I C c = N l o g ( L ) + 2 K N N K 1 ,  
where N is the sample size, K is the number of parameters, and L is the likelihood function. The AICc is derived under Gaussian assumptions which are weakly dependent, so the use of AICc is highly recommended in practice.

3.2. The Bayesian Information Criterion (BIC)

The BIC was first introduced by Schwarz (1978) as a competitor to AIC. The BIC is an asymptotic approximation to a transformation of the Bayesian posterior probability of the model [30], which is defined as:
B I C = N l n ( L ) + k l n ( N ) .
The smaller value of BIC means better fitting accuracy between model and data. Compared with AIC, the BIC punishes the model parameter number more in the case of large data; therefore, the BIC prefers to choose simple models with fewer parameters.

3.3. The Root Mean Square Error (RMSE)

The RMSE is used to calculate how much the forecasted value differs from the actual value [31]. The smaller value shows the higher prediction accuracy; therefore, the root mean square error is used as an evaluation standard of measurement precision. The RMSE is defined as:
R M S E = i = 1 N ( A i F i ) 2 N ,
where N is the size of data points, and A i and F i represent the actual value and forecasted value, separately. The RMSE can be applied to estimate the accuracy of models.

3.4. The Pearson’s Correlation Coefficient (R)

In order to accurately describe the degree of correlation between the forecasted value and the actual value [31], the Pearson’s correlation coefficient (R) is defined as:
R = N i = 1 N A i F i ( i = 1 N A i ) ( i = 1 N F i ) N ( i = 1 N A i 2 ) ( i = 1 N A i ) 2 N ( i = 1 N F i 2 ) ( i = 1 N F i ) 2 .
The value of R satisfies 1 < R < 1 . The + and show positive and negative linear correlations between the forecasted and actual values, respectively. The R greater than or equal to 0.8 is considered as strong correlation. R 2 demonstrates the strength of the association between A i and F i .

4. Model Evaluation and Analysis Based on the Real Data

To make a quantitative comparison of different models, the outbreak of COVID-19 is studied, which has seriously affected economy and health of people around the world. In this section, we apply four means to evaluate and analyze the most widely used models in the field of epidemiology. Since the data link of COVID-19 in China is intact, we use the data from the China National Health Commission (http://www.chinacdc.cn/, accessed on 12 September 2021) for simulation. The fifty-four-day data from 22 January 2020 to 15 March 2020 of Beijing, Chongqin, Tianjing, and Heilongjiang are used, which are considered as some successful cities during the battle against the disease in China. The unknown parameters involved in integer-order and fractional models are identified by the non-linear least squares method and Simulink Design Optimization of MATLAB. Through a large number of numerical tests, the following conclusions can be drawn.

4.1. The Fractional Models Can Better Fit the Real Data than the Corresponding Integer-Order Models

It is well known that the outbreak of COVID-19 is not only dependent on the number of current infected people, but also relates to those who were infected in the past; that is, the disease has memory. The fractional epidemic model can reflect the memory and heredity properties, which is more compatible with the spread process of disease. The four means are computed to evaluate and analyze the validity of integer-order and fractional models.
The smaller AICc, BIC, RMSE values and the larger R value of model indicate the better fitting accuracy, respectively. The difference of four means between integer-order and fractional models are given as:
D j = I N T j F R A j ,
where j is the AICc, BIC, RMSE, or R, and I N T j and F R A j are the values of the mean j on the integer-order model and the fractional model, respectively. When D AICc > 0 ,   D BIC > 0 ,   D RMSE > 0 , and D R < 0 , the fractional model has better fitting accuracy than the corresponding integer-order model, respectively.
In Figure 2, the D j are calculated from 22 January 2020 to different end times in Beijing, Chongqing, Tianjin, and Heilongjiang. It is found that all of the D j about AICc, BIC, and RMSE are greater than 0, and the D R values are less than 0; that is, the fractional models are better than the corresponding integer-order models for different time windows.
For convenience, all of our subsequent discussions are based on the fractional models.

4.2. No Model Is Reliable for Long-Term Forecasting Based on the Early-Stage Real Data

The rapid spreading of COVID-19 results in a dramatic increase in the number of infections. Therefore, how to predict its spread trend is a tough task. Five fractional models are used to estimate long-term prediction ability. In Figure 3, the predication and fitting results of three cities (Beijing, Chongqing, Tianjin) and Heilongjiang province are given by AICc, BIC, RMSE, and R means. In the early stage of the epidemic with very limited data, all models can fit and predict spread trend in a short-term very well, but no model is reliable for long-term forecasting. All models always underestimate or overestimate the size of the epidemic for long-term prediction.

4.3. The Fractional SEIR-Q and SEIR-QD Models Can More Accurately Describe the COVID-19 Spread Trends

Fractional dynamical models are more in line with the biological significance by specific parameters (fractional order). The fractional dynamical models can be divided into three groups based on their performance. In Figure 4, the fractional SIR and SEIR models look insufficient to describe the spread trends, especially in the final stage, since too few parameters are involved. The fractional SEIR-AHQ model contains redundant parameters, which generates the larger value of AICc. That leads to its fluctuation being the greatest among the models. The fractional SEIR-Q and SEIR-QD systems are more suitable to characterize the infectious characteristics of COVID-19 by appropriately introducing the effect of quarantine.

4.4. The Inflection Point of the Real Data Is Vital for Prediction

The inflection point plays a pivotal role in predicting the epidemic trend, which is in agreement with previous report [28]. In Figure 3, when enough data are available, typically when the number of infected people passes the inflection point, forecast values converge to the true values by fractional models. For example, all models fail to make a convincing forecast in Beijing when using the data before 14 February (early stage and middle stage), which is also the inflection point. However, after the inflection point (late stage), fractional SEIR-Q, SEIR-QD, and SEIR-AHQ models can successfully capture the real data.

4.5. A Single Mean Is Insufficient to Evaluate the Model’s Prediction Capability

When using the models to fit data from the epidemic process, four means are used to select suitable models. The AICc and the BIC consider the trade-off between the fitting effect of model and the principle of least parameters, but the AICc is more sensitive to the parameters. The RMSE measures the degree of deviation between the real data and the fitted data, and the R highlights the correlation between the outcomes of model and the real data. It should be note that a single evaluation criterion is not reliable. For example, in Figure 4a, when the end times are from 4 February to 8 February, the most suitable model is the fractional SIR, but the fractional SEIR-QD model achieves the better trade-off in the rest time period. In Figure 4b, the fractional SEIR-QD model can be picked out by BIC. In Figure 4c, the fractional SEIR-AHQ model is the best model by RMSE. In Figure 4d, the fractional SEIR-Q, SEIR-QD, and SEIR-AHQ models are considered as strong correlations between the forecasted data and the actual data by R. Therefore, the models cannot be evaluated by a single mean; namely, the means should be combined to select suitable models.

4.6. All of the Basic Reproduction Number R 0 Heavily Depend on Contact Rate

The sensitivity analysis for the fractional epidemic models is to study the main intervention measures that affect the COVID-19 dynamics. Using the next-generation matrix method [32], the basic reproduction number of fractional SIR, SEIR, SEIR-Q, SEIR-QD, and SEIR-AHQ models are given as α β , α β , α 1 ( γ + κ δ ) + α 2 ε ( ε + κ δ ) ( γ + κ δ ) , α β and β 1 α μ ( 1 ε ) + α ( β 2 + δ 1 ) ( 1 ε ) ( 1 μ ) β 1 ( β 2 + δ 1 ) , respectively. Latin hypercube sampling (LHS) and partial rank correlation coefficient (PRCC) [33] are used to analyze the influence of the parameters on the model. The higher the PRCC value of a parameter, the greater its impact on the dynamic behavior [34]. The involved parameters are given in Table 1, which includes PRCC and p-values. The bar graph of PRCC results is shown in Figure 5. Furthermore, whatever the model, R 0 are heavily dependent on the contact rate α , except for α 1 and α 2 in the SEIR-Q model. It is obvious that non-medical intervention measures (social distance restriction, self-isolation, staying at home or wearing masks, etc.) have directly impacted on epidemic prevention and control.
Based on the references [6,20,21,35], we could easily obtain unique solutions and a locally asymptotically stable of equilibrium points. We only take fractional SEIR-QD as an example; other models can get similar results.
Remark 4.1.
The initial value problem for the fractional SEIR-QD model has a unique solution in + 7 .
Remark 4.2.
The diseases-free equilibrium (0, 0, 0, 0,   R * ,   D * ,   P * ) of the fractional SEIR-QD system is locally asymptotically stable if R 0 < 1 and unstable if R 0 > 1.

5. Conclusions

In this paper, we systematically studied the prediction capability of five widely used integer-order and fractional-order models with the COVID-19 data in China. Through numerical simulations, many insightful results are obtained by the four means and the PRCC. Firstly, the fractional models have better performance than the corresponding integer-order models. Fractional models can better describe epidemic characteristics in the diseases by an additional parameter (fractional order). Secondly, almost all of the models overestimate or underestimate the infected people of outbreak based on the early-stage data. Thirdly, the fractional SEIR-Q and SEIR-QD models can better predict the scales of epidemic trends, which have the better trade-off between the model complexity and the fitting accuracy. Fourthly, it is found that the inflection points in the data play a very important role. The prediction results of all integer-order and fractional models are closely relate to inflection points. Fifthly, a single evaluation criterion is insufficient to estimate the model’s capability. The best model for modeling COVID-19 can be picked out by AICc and BIC means. Finally, based on the LHS and PRCC methods, we find that non-medical intervention measures play a critical role in COVID-19 control.
This paper provides a new insight about selection and analysis of infectious diseases by fractional modeling. Furthermore, models with vaccine can play an important role in the allocation of vaccines, which are not considered in this paper. We will discuss this in our future work.

Author Contributions

N.M.: Prepared the draft. W.M.: Developed the concept, methodology, and revised the paper finally. Z.L.: Performed the analysis. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Fundamental Research Funds for the Central Universities (Grant No. 31920210018), and the Innovation Team of Intelligent Computing and Dynamical System Analysis and Application of Northwest Minzu University.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. WHO. Coronavirus (COVID-19) Dashboard. Available online: https://covid19.who.int/ (accessed on 12 September 2021).
  2. Giordano, G.; Blanchini, F.; Bruno, R. Modelling the COVID-19 epidemic and implementation of population-wide interventions in Italy. Nat. Med. 2020, 26, 855–860. [Google Scholar] [CrossRef]
  3. Shaikh, A.S.; Shaikh, I.N.; Nisar, K.S. A mathematical model of COVID-19 using fractional derivative outbreak in India with dynamics of transmission and control. Adv. Differ. Equ. 2020, 2020, 373. [Google Scholar] [CrossRef]
  4. Rihan, F.A.; Alsakaji, H.J.; Rajivganthi, C. Stochastic SIRC epidemic model with time-delay for COVID-19. Adv. Differ. Equ. 2020, 2020, 1–20. [Google Scholar] [CrossRef]
  5. Rihan, F.A.; Al-Mdallal, Q.M.; AlSakaji, H.J. A fractional-order epidemic model with time-delay and nonlinear incidence rate. Chaos Solitons Fractals 2019, 126, 97–105. [Google Scholar] [CrossRef]
  6. Ahmad, S.; Ullah, A.; Al-Mdallal, Q.M. Fractional order mathematical modeling of COVID-19 transmission. Chaos Solitons Fractals 2020, 139, 110256. [Google Scholar] [CrossRef]
  7. Ma, W.Y.; Zhao, Y.T.; Guo, L.H.; Chen, Y.Q. Qualitative and quantitative analysis of the COVID-19 pandemic by a two-side fractional-order compartmental model. ISA Trans. 2021. submitted. [Google Scholar]
  8. Comunian, A.; Gaburro, R.; Giudici, M. Inversion of a SIR-based model: A critical analysis about the application to COVID-19 epidemic. Phys. D Nonlinear Phenom. 2020, 413, 132674. [Google Scholar] [CrossRef]
  9. Yang, Z.F.; Zeng, Z.Q.; Wang, K. Modified SEIR and AI prediction of the epidemics trend of COVID-19 in China under public health interventions. J. Thorac. Dis. 2020, 12, 165–174. [Google Scholar] [CrossRef] [PubMed]
  10. Kuniya, T.; Inaba, H. Possible effects of mixed prevention strategy for COVID-19 epidemic: Massive testing, quarantine and social distancing. AIMS Public Health 2020, 7, 490–503. [Google Scholar] [CrossRef] [PubMed]
  11. Peng, L.R.; Yang, W.Y.; Zhang, D.Y.; Zhuge, C.J.; Hong, L. Epidemic analysis of COVID-19 in China by dynamical modeling. arXiv 2020, arXiv:2002.06563. [Google Scholar]
  12. Tang, B.; Bragazzi, N.L.; Li, Q.; Tang, S.; Xiao, Y.; Wu, J. An updated estimation of the risk of transmission of the novel coronavirus (2019-nCov). Infect. Dis. Modelling 2020, 5, 248–255. [Google Scholar] [CrossRef]
  13. Sindhu, T.N.; Shafiq, A.; Al-Mdallal, Q.M. On the analysis of number of deaths due to COVID-19 outbreak data using a new class of distributions. Results Phys. 2021, 21, 103747. [Google Scholar] [CrossRef]
  14. Hedayati, M.; Ezzati, R.; Noeiaghdam, S. New Procedures of a Fractional Order Model of Novel Coronavirus (COVID-19) Outbreak via Wavelets Method. Axioms 2021, 10, 122. [Google Scholar] [CrossRef]
  15. Noeiaghdam, S.; Micula, S.; Nieto, J.J. A Novel Technique to Control the Accuracy of a Nonlinear Fractional Order Model of COVID-19: Application of the CESTAC Method and the CADNA Library. Mathematics 2021, 9, 1321. [Google Scholar] [CrossRef]
  16. Khan, A.; Alshehri, H.M.; Abdeljawad, T. Stability analysis of fractional nabla difference COVID-19 model. Results Phys. 2021, 22, 103888. [Google Scholar] [CrossRef]
  17. Rajagopal, K.; Hasanzadeh, N.; Parastesh, F. A fractional-order model for the novel coronavirus (COVID-19) outbreak. Nonlinear Dyn. 2020, 101, 711–718. [Google Scholar] [CrossRef] [PubMed]
  18. Zhang, Z.Z.; Zeb, A.; Egbelowo, O.F.; Erturk, V.S. Dynamics of a fractional order mathematical model for COVID-19 epidemic. Adv. Differ. Equ. 2020, 2020, 420. [Google Scholar] [CrossRef]
  19. Bahloul, M.A.; Chahid, A.; Laleg-Kirati, T.M. Fractional-Order SEIQRDP Model for Simulating the Dynamics of COVID-19 Epidemic. IEEE Eng. Med. Biol. Mag. 2020, 1, 249–256. [Google Scholar] [CrossRef]
  20. Xu, C.H.; Yu, Y.G.; Chen, Y.Q.; Lu, Z.Z. Forecast analysis of the epidemics trend of COVID-19 in the United States by a generalized fractional-order SEIR model. Nonlinear Dyn. 2020, 101, 1621–1634. [Google Scholar] [CrossRef] [PubMed]
  21. Lu, Z.Z.; Yu, Y.G.; Chen, Y.Q.; Ren, G.J. A fractional-order SEIHDR model for COVID-19 with inter-city networked coupling effects. Nonlinear Dyn. 2020, 101, 1717–1730. [Google Scholar] [CrossRef] [PubMed]
  22. Cyranoski, D. When will the Coronavirus Outbreak Peak? 2020. Available online: http://www.nature.com/articles/d41586-020-00361-5 (accessed on 12 September 2021).
  23. Tang, B.; Wang, X.; Li, Q. Estimation of the transmission risk of the 2019-nCoV and its implication for public health interventions. J. Clin. Med. 2020, 9, 462. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Chen, T.M.; Rui, J.; Wang, Q.P.; Zhao, Z.Y.; Cui, J.A. A mathematical model for simulating the phase-based transmissibility of a novel coronavirus. Infec. Dis. Poverty 2020, 9, 1–8. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Lin, Q.; Zhao, S.; Gao, D.; Lou, Y.; Yang, S.; Musa, S.S. A conceptual model for the Coronavirus disease 2019 (COVID-19) outbreak in Wuhan, China with individual reaction and governmental action. Int. J. Infect. Dis. 2020, 93, 211–216. [Google Scholar] [CrossRef]
  26. Ma, W.Y.; Li, C.P.; Wu, Y.J.; Wu, Y.Q. Synchronization of fractional fuzzy cellular neural networks with interactions. Chaos Interdiscip. J. Nonlinear Sci. 2017, 27, 103106. [Google Scholar] [CrossRef] [PubMed]
  27. Li, Z.M.; Ma, W.Y.; Ma, N.R. Partial topology identification of tempered fractional-order complex networks via synchronization method. Math. Methods Appl. Sci. 2021, 1–14. [Google Scholar] [CrossRef]
  28. Yang, W.Y.; Zhang, D.Y.; Peng, L.R. Rational evaluation of various epidemic models based on the COVID-19 data of China. arXiv 2020, arXiv:2003.05666. [Google Scholar]
  29. Burnham, K.P.; Anderson, D.R. Multimodel Inference: Understanding AIC and BIC in Model Selection. Sociol. Methods Res. 2004, 33, 261–304. [Google Scholar] [CrossRef]
  30. Cremonesi, P.; Koren, Y.; Turrin, R. Performance of recommender algorithms on top-n recommendation tasks. In Proceedings of the Fourth ACM Conference on Recommender Systems, Barcelona, Spain, 26–30 September 2010; pp. 39–46. [Google Scholar]
  31. Singh, P.; Borah, B. An efficient time series forecasting model based on fuzzy time series. Eng. Appl. Artif. Intell. 2013, 26, 2443–2457. [Google Scholar] [CrossRef]
  32. van den Driessche, P.; Watmough, J. Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission. Math. Biosci. 2020, 180, 29–48. [Google Scholar] [CrossRef]
  33. Renardy, M.; Hult, C.; Evans, S. Global sensitivity analysis of biological multiscale models. Curr. Opin. Biomed. Eng. 2019, 11, 109–116. [Google Scholar] [CrossRef]
  34. Ullah, S.; Khan, M.A. Modeling the impact of non-pharmaceutical interventions on the dynamics of novel coronavirus with optimal control analysis with a case study. Chaos Solitons Fractals 2020, 139, 110075. [Google Scholar] [CrossRef] [PubMed]
  35. Owusu-Mensah, I.; Akinyemi, L.; Oduro, B. A fractional order approach to modeling and simulations of the novel COVID-19. Adv. Differ. Equ. 2020, 2020, 1–21. [Google Scholar] [CrossRef] [PubMed]
Figure 1. The transfer diagrams of five models. (a) SIR; (b) SEIR; (c) SELR-Q; (d) SEIR-QD; and (e) SEIR-AHQ.
Figure 1. The transfer diagrams of five models. (a) SIR; (b) SEIR; (c) SELR-Q; (d) SEIR-QD; and (e) SEIR-AHQ.
Fractalfract 05 00120 g001
Figure 2. The difference values D j from 22 January 2020 to different end times in Beijing, Chongqing, Tianjin, and Heilongjiang. (ad) D j of Beijing; (eh) D j of Chongqing; (il) D j of Tianjin; (mp) D j of Heilongjiang.
Figure 2. The difference values D j from 22 January 2020 to different end times in Beijing, Chongqing, Tianjin, and Heilongjiang. (ad) D j of Beijing; (eh) D j of Chongqing; (il) D j of Tianjin; (mp) D j of Heilongjiang.
Fractalfract 05 00120 g002
Figure 3. Prediction of the COVID-19 epidemic trends in Beijing city, Chongqing city, Tianjin city, and Heilongjiang province from 22 January 2020 to 5 March 2020 based on the data of first 10 (early stage), 20 (middle stage), and 30 (late stage) days, respectively. (ac) The number of infected people in Beijing at different stages; (df) The number of infected people in Chongqing at different stages; (gi) The number of infected people in Tianjin at different stages; (jl) The number of infected people in Heilongjiang at different stages.
Figure 3. Prediction of the COVID-19 epidemic trends in Beijing city, Chongqing city, Tianjin city, and Heilongjiang province from 22 January 2020 to 5 March 2020 based on the data of first 10 (early stage), 20 (middle stage), and 30 (late stage) days, respectively. (ac) The number of infected people in Beijing at different stages; (df) The number of infected people in Chongqing at different stages; (gi) The number of infected people in Tianjin at different stages; (jl) The number of infected people in Heilongjiang at different stages.
Fractalfract 05 00120 g003
Figure 4. The values of AICc, BIC, RMSE, and R from 22 January to different end times in Beijing. (a) AICc; (b) BIC; (c) RMSE; (d) R.
Figure 4. The values of AICc, BIC, RMSE, and R from 22 January to different end times in Beijing. (a) AICc; (b) BIC; (c) RMSE; (d) R.
Fractalfract 05 00120 g004
Figure 5. The PRCC results of parameters involved in R 0 . (a) SIR; (b) SEIR; (c) SELR-Q; (d) SEIR-QD; (e) SEIR-AHQ.
Figure 5. The PRCC results of parameters involved in R 0 . (a) SIR; (b) SEIR; (c) SELR-Q; (d) SEIR-QD; (e) SEIR-AHQ.
Fractalfract 05 00120 g005
Table 1. The PRRC values of R 0 with corresponding P-values.
Table 1. The PRRC values of R 0 with corresponding P-values.
ModelsParametersPRCC Valuesp-Values
SIR α 0.88660.0000
β −0.4624 1.5497 × 10 106
SEIR α 0.86670.0000
β −0.4556 4.3039 × 10 103
SEIR-Q α 1 0.73260.0000
α 2 0.2386 2.7566 × 10 27
ε −0.5570 2.1135 × 10 163
κ −0.2734 1.3051 × 10 35
δ −0.1798 5.4878 × 10 16
γ 0.03960.0763
SEIR-QD α 0.88230.0000
β −0.4399 1.9184 × 10 95
SEIR-AHQ α 0.7213 5.2370 × 10 321
μ −0.5645 1.0293 × 10 168
ε −0.3493 1.7440 × 10 58
δ 1 −0.1348 1.4459 × 10 9
β 1 −0.07270.0011
β 2 0.02070.3542
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ma, N.; Ma, W.; Li, Z. Multi-Model Selection and Analysis for COVID-19. Fractal Fract. 2021, 5, 120. https://doi.org/10.3390/fractalfract5030120

AMA Style

Ma N, Ma W, Li Z. Multi-Model Selection and Analysis for COVID-19. Fractal and Fractional. 2021; 5(3):120. https://doi.org/10.3390/fractalfract5030120

Chicago/Turabian Style

Ma, Nuri, Weiyuan Ma, and Zhiming Li. 2021. "Multi-Model Selection and Analysis for COVID-19" Fractal and Fractional 5, no. 3: 120. https://doi.org/10.3390/fractalfract5030120

Article Metrics

Back to TopTop