Stability Analysis of an Extended SEIR COVID-19 Fractional Model with Vaccination Efficiency

Wali, Mubashara; Arshad, Sadia; Huang, Jianfei

doi:https://doi.org/10.1155/2022/3754051

Computational and Mathematical Methods in Medicine

On this page

Abstract Introduction Preliminaries Discussion Data Availability Conflicts of Interest References Copyright Related Articles

Research Article | Open Access

Volume 2022 | Article ID 3754051 | https://doi.org/10.1155/2022/3754051

Stability Analysis of an Extended SEIR COVID-19 Fractional Model with Vaccination Efficiency

Mubashara Wali,¹Sadia Arshad,¹and Jianfei Huang²

Academic Editor: Khalid Hattaf

Received11 Jun 2022

Revised17 Aug 2022

Accepted01 Sept 2022

Published20 Sept 2022

Abstract

This work is aimed at presenting a new numerical scheme for COVID-19 epidemic model based on Atangana-Baleanu fractional order derivative in Caputo sense (ABC) to investigate the vaccine efficiency. Our construction of the model is based on the classical SEIR, four compartmental models with an additional compartment V of vaccinated people extending it SEIRV model, for the transmission as well as an effort to cure this infectious disease. The point of disease-free equilibrium is calculated, and the stability analysis of the equilibrium point using the reproduction number is performed. The endemic equilibrium’s existence and uniqueness are investigated. For the solution of the nonlinear system presented in the model at different fractional orders, a new numerical scheme based on modified Simpson’s 1/3 method is developed. Convergence and stability of the numerical scheme are thoroughly analyzed. We attempted to develop an epidemiological model presenting the COVID-19 dynamics in Italy. The proposed model’s dynamics are graphically interpreted to observe the effect of vaccination by altering the vaccination rate.

1. Introduction

Emerging in December , COVID-19 became the biggest global threat. The whole world faced a critical situation, and a panic was created throughout the globe. The World Health Organization (WHO) declared it as an international health emergency. Not much was known about this virus as it was a new virus that has not been identified earlier in humans. So, relying on the observations and symptoms shown by the people who caught the virus, strategies of isolation, lockdown, and social distancing were adopted to prevent the disease from spreading. With the passage of time, the virus becomes more threatening because of its new variants. To investigate the dynamics and transmission of the virus, researchers utilized the concept of mathematical modeling [1–4]. These models actually do not provide the cure for the infectious diseases but can contribute in predicting the dynamics and behaviour of these infections through simulations.

Researchers rely more on fractional models because of their nonlocal behaviour, hereditary properties, and memory effects rather than integer order models. Besides, a number of experimental facts exhibit that the natural dynamics go along with fractional calculus. Ahmad et al. [5] utilized the modified Euler’s method to study the dynamics of the model for COVID-19 based on fractional differential equations and simulate the results by using the available data of initial days of the spread in Wuhan City. Zakary et al. [6] investigated the situation of the spread of the Corona Virus under the effect of two classes, quarantined and other who did not respect the quarantine strategy using the mathematical model, and an optimal control strategy to reduce the infections in Morocco. Ahmad et al. [7] presented a brief analysis on the WHO reported data for Pakistan and investigated it using the SEIR fractional model. Askar et al. [8] considered the SITR fractional model to forecast the transmission of deadly virus along with imposed lockdown in India. Bushnaq et al. [9] investigated the effect of physical distance to control the spread of the virus through fractional order model. Abdulwasaa et al. [10] proposed a fractal fractional model for forecasting deaths and new cases of COVID-19 outbreak in India. Baba and Rihan [11] used the Caputo-Fabrizio fractional order model to analyze the dynamics of COVID-19 variants. A fractional order delay differential model for COVID-19 infection was provided by Rihan and Gandhi [12] to understand what triggers the severity of symptoms and disease of polluted lung and respiratory system. A stochastic SIAQR epidemic model was proposed by Rihan and Alsakaji [13] to understand the dynamics of COVID-19. Furati et al. [14] tested the data taken from China, Saudi Arabia, Brazil, and Italy by proposing a numerical scheme using fractional order model performing the two-step generalized exponential time-differencing method. Alqahtani [15] numerically analyzed and investigated the stability of the SIR epidemic model considering the health system. Recently another fractional order derivative with nonsingular kernel named as Generalized Hattaf Fractional (GHF) derivative [16] is introduced. We can see the use of this derivative for epidemiological models in recent literature. By creating a new numerical method, Hattaf [17] investigated the qualitative properties of solutions to fractional differential equations with the new GHFs derivative, including stability, asymptotic stability, and Mittag-Leffler stability and later applied the obtained analytical and numerical results to a biological nonlinear system derived from epidemiology. Cheneke et al. [18] developed and examined a coinfection of HIV and cholera model using GHF derivative and the solution’s behaviour is deciphered. In present work we have used ABC derivative for our proposed model. However, we can extend our work using the GHF derivative for fractional order mathematical models in future.

Despite all the measures taken to control the virus, COVID-19 remained a threat because of its high propagation rate leaving the majority of the population susceptible to the disease and escalating the need to develop a suitable and effective vaccine. The process to gain immunity against COVID-19 is still under discussion. Most of the people infected with COVID-19 developed the immunity against the virus, but still it is unknown how long this immunity will last as there are some reported cases of getting the infection for a second time. In this situation, the only solution to gain herd immunity is mass vaccination. Hence, to achieve this milestone, vaccination is required for a huge portion of the global population. In this regard, WHO helped governments to launch and refine their vaccination programmes along with UNICEF, Gavi, and partners. According to reference [19], 43.5% of the world population has been given at least one shot of COVID-19 vaccine till the 21^st of September 2021. Modified mathematical models with an additional compartment of the vaccinated class are discussed in recent literature to foresee the efficacy, effect, and dynamics of the vaccination. Meng et al. [20] established an SEIRV model to differentiate between the effects of mandatory and voluntary vaccination methods on heterogeneous networks. Ramos et al. [21] investigated the impact of variants and vaccines on the spread of COVID-19 applying the model to Italy as a test case. Tetteh et al. [22] analyzed the mass and ring vaccination strategies to control the transmission of SARS-CoV-2 by varying vaccine efficacy rate. Harizi et al. [23] used the SIRV model to check the future dynamics of the pandemic by varying vaccination rates in Canada. Karabay et al. [24] presented an SEIR simulator to model effective immunization and sterilizing outcomes in order to assess the vaccination effect. Tuteja [25] numerically solved the modified SEIR model including the vaccinated class to investigate the viral dynamics of the epidemic. Ghostine [26] proposed an SEIR model with an extension of the vaccination compartment to investigate the effect of vaccination in Saudi Arabia. Wintachaia and Prathomb [27] presented a detailed stability analysis of the COVID-19 epidemic model to investigate the vaccine efficiency. Alsakaji et al. [28] examined the dynamics of COVID-19 in the UAE using an exptended SEIR epidemic model with vaccination, temporal delays, and random noise. For transmission of the infectious disease, Kou et al. [29] established a multiscale agent-based model to investigate the impact of nonpharmaceutical interventions and vaccination. Acuña-Zegarra. [30] formulated a COVID-19 optimal control problem to study the efficacy and responses of the vaccine induced immunity. Albani et al. [31] estimated the parameters of SEIR like epidemic model using daily reports of COVID-19 of New York and Chicago to analyze the impact caused by vaccination delay.

Motivated from the above studies, we have presented in this work a modified SEIRV epidemic model with an additional compartment of vaccinated individuals to examine the effectiveness and efficacy of the vaccine by varying its rate considering the ABC derivative. We have used ABC fractional order derivative because a number of experimental facts exhibit that the natural dynamics go along with fractional calculus as the fractional order derivatives are nonlocal in behaviour and possess the hereditary properties and memory effects. Also, we have constructed a new, more efficient, and effective numerical scheme for fractional differential equation with Atangana-Baleanu derivative in Caputo sense using Simpson’s 1/3 method for the approximation of the integral function involved in the problem.

The paper comprises different sections. In section 2, some basic definitions are recalled which have been used in the work later. SEIRV epidemic model for COVID-19 is presented in section 3. The equilibrium points and the basic reproduction number using next generation method are calculated, and the stability of these points is also analyzed in this section. For the solution of the proposed model modified Simpson’s 1/3 method alongside complete stability and convergence analysis, it is presented in section 4. Conversion of mathematical model to fractional order model is given in section 5. Numerical simulations are also presented in the same section.

2. Preliminaries

This part will review some of the fundamental definitions of fractional calculus, which will be useful in the subsequent work.

Definition 1. For and , the th–order ABC derivative and integral of are defined [32] as where represents the normalization function given by with and is the Mittag-Leffler function defined as .

The lemma below presents the Gronwall inequality [33], which will be used for the proof of our main result.

Lemma 2.1. The Gronwall inequality for fractional order differential equation states
Let independent of , and {} satisfy the inequality with Then where the parameter is the Mittage-Leffler function. Particularly when inequality (5) becomes

In this study, we used an SEIR model with vaccination effectiveness to forecast the COVID-19 condition when a vaccine is released.

3. Mathematical Model

Mathematical models are of great importance in predicting the behaviour of the viral epidemic diseases. There are a number of models which can be used to describe and investigate the spread of the disease. Recently COVID-19 emerged as a viral epidemic, and we are not aware enough about the behaviour and spread of this virus. To investigate the COVID-19 situation under the effect of vaccination, we consider the following SEIRV epidemic disease model:

To define the model equations, the total population, represented by , is separated into susceptible, exposed, infected, recovered, and vaccinated classes denoted by , , , , and . The flow chart of the above system is shown in Figure 1.

The above system of equation (7) is subject to the initial conditions

and the following are the parameters:

per-capita natural death rate,

the rate at which infectious people recover,

average fatality rate for virus-infected people,

disease transmission probability per contact times the amount of contacts made per unit of time (dimensionless),

the birth rate per-capita,

progression rate from E to I,

contagious period,

incubation period,

vaccination rate,

vaccine ineffectiveness,

vaccine effectiveness,

having T units per unit of time.

3.1. Stability Analysis of SEIRV COVID-19 Model

In this section, stability analysis will be carried out to find the disease-free equilibrium and endemic equilibrium points. In order to get these points, each of the equations given in system of equation (7) will be equated to zero.

3.2. Disease-Free Equilibrium and the Basic Reproduction Number

Because there is no disease propagation in disease-free equilibrium, so , and hence we get so the disease-free equilibrium point is .

The average number of secondary infections induced by an infectious individual is given by the basic reproduction number . It estimates the growth of the virus outbreak. The basic reproduction number is determined using the next generation method [34]. Using and from system of equation (7), the Jacobian matrices of the disease-free equilibrium are as follows:

Here is the next generation matrix and we can obtain the basic reproduction number as , so

Substituting the values of and , the following value of the reproduction number is obtained for the system of equation (7):

Theorem 3.1. If then the COVID-19 model of equation (7) disease-free equilibrium point is locally asymptotically stable.

Proof. The Jacobian of system of equation (7) is given as Jacobian around disease-free equilibrium point is We can observe from the preceding matrix that three of its eigenvalues , are negative and we are left with Now from above matrix, the last two eigenvalues are where and , and for , we have and .
Hence for the disease-free equilibrium is locally asymptotically stable.

3.3. Existence and Uniqueness of Endemic Equilibrium

In order to investigate the existence and uniqueness of the endemic equilibrium for , we will set all the derivatives of model of equation (7) equals to zero. Let be the endemic equilibrium of model of equation (7). Then by solving the following equations we can obtain the endemic equilibrium.

Now, solving first, third, and fourth equations of equation (17), we get

Also equation (17) gives

The following quadratic equation for is obtained: where

It is evident that is always positive while is negative for . Hence, there exists a positive and unique value of , which results in a unique endemic equilibrium for .

4. Numerical Approximation

Consider the nonlinear fractional order differential equation where The vector function in the differential equation is real valued and continuous and satisfies the Lipschitz condition.

To derive the numerical scheme for the above mentioned nonlinear fractional order differential equation, the given interval will be subdivided into equal intervals with time step , where for For convenience, the approximate solution of will be denoted by with , .

Applying the integral operator of equation (2) on equation (22), we get where the integral in the above equation is the ABC integral. Using equation (2) to discretize the integral , we get

is approximated on each subinterval by the following piecewise quadratic interpolation polynomial. where is mid point of the interval

Also, the value of will be approximated using the interpolation

Putting equation (27) in equation (26), we get

Using above approximation, equation (25) takes the form

Evaluating all the integrals given in equation (30), we obtain where

Thus we can obtain the following modified fractional Simpson’s 1/3 scheme for equation (22) in ABC sense:

4.1. Convergence and Stability Analysis

The following lemmas are very important to examine the proposed scheme’s convergence and stability.

Lemma 4.1. Assume that , then there exist a constant such that

Proof. Using the Lagrange interpolation error formula, we have From above lemma it is clear that the truncation error of the scheme of equation (33) is .

Lemma 4.2. For , the weights , given in equation (33) have the following order of magnitude:

Proof. According to the definition of it deduces that implies Note that , so the coefficients of power of are less than 1. Thus, the series given above is convergent, and hence the proof is completed for . Similar strategy can be adopted for , so we neglect the proof.

Theorem 4.1. The numerical scheme (33) is conditionally convergent. That is, for and a sufficiently small , there exist a constant such that

Proof. To prove the theorem, considering the difference between actual solution and approximate solution, we get Using Lipschitz condition of equation (23) after applying triangle inequality and using Lemma 4.1, we get Now using Lemma 4.2 under the condition , we get where For any and a sufficiently small , there exist a constant such that Hence by applying Gronwall inequality given in Lemma 2.1 and using the property (43), we conclude

Theorem 4.2. For and a sufficiently small , there exist a constant such that the numerical scheme (33) is stable.

Proof. To prove the theorem, let the perturbation of and be denoted by and , respectively, for . Adding this perturbation to the numerical scheme (33) we get Subtracting equation (33) from equation (45), we get Using the Lipschitz condition followed by triangle inequality, we get Now if then using Lemma 4.2 we get where
Using property (43), we have Hence by applying the Gronwall inequality given in Lemma 2.1, we conclude

5. Fractional Model

Employing Atangana-Baleanu fractional order derivative in Caputo sense of order on system of equation (7), we have the following fractional model, along with initial conditions

Adopting the previously defined procedure, model of equation (51) takes the form

Furthermore applying the numerical scheme(33), system (53) becomes where are already defined.

5.1. The Positivity and Boundness of the Solution

In this subsection we will discuss the positive nature of the epidemiological model of equation (51). We will recall the following lemma for the proof of the main theorem on the positivity of the solution of the fractional model (51).

Lemma 5.1. Generalized Mean Value Theorem: let and Atangana-Baleanu fractional order derivative,for, then with

Remark 1. Let and Atangana-Baleanu fractional order derivative, for . Lemma 5.1 shows that if , then the function is nondecreasing and if then is nonincreasing for all

Theorem 5.1. The solution to the nonnegative initial conditions of equation (52) of the proposed model of equation (51) is unique and limited in .

Proof. The existence and uniqueness of the model solution of equations (51) and (52) can be achieved in the time interval by [35]. The nonnegative region would seem to be a positive invariant region that must be represented. We observe from the model of equation (51) From Remark 1 and system of equation (56), the solution . Also, in each line bounded the nonnegative octant, the vector field points will remain in . Therefore, the fractional model of equation (51) of a solution ( is nonnegative if the initial condition is set positively to invariant.

6. Numerical Results and Discussion

Considering the case of Italy, we applied the proposed model to analyse the impact of vaccination. The vaccination campaign for COVID-19 started on 27^th December 2020 in Italy managed by the Ministry of Health. In the early months of the vaccination campaign, the Italian government targeted health doctors and administrative staff as well as guests and nursing home personnel followed the vaccination of elderly people and public service personnel in second phase of the campaign. Italian government bought four out of seven, WHO approved vaccines to facilitate its people against the virus. So far, 73.3% of the total population has received at least one dose and 65.5% of the total population is fully vaccinated. As of September 2021, top ten regions of Italy for the total number of administered COVID-19 vaccine doses can be seen in Figure 2. Moreover the number of people who received at least one dose and fully vaccinated according to the data [19] are presented in Figure 3.

Further, graphical presentation of the numerical simulation for the fractional model (51) are given in this section. Table 1 shows the estimated model parameters and initial values used to evaluate the impact of vaccination on the COVID-19 epidemic.

Figure 4 shows the impact by varying the vaccination rates on different classes of the proposed model with representing the model output without vaccination. Although there is no remarkable change in the graph of susceptible class, we can observe a significant decrease in infected and exposed classes by increasing the vaccination rate. We can see the peak of the infected as well as that of exposed individuals gradually decreases with the increase in effective rate of vaccination. From Figure 4, we can observe that for exposed class there is a difference of almost individuals in peak values for no vaccination to a vaccination rate . Similarly, the peak values for no vaccination compared to a vaccination rate of for the infected class differ by almost individuals. So, interpreting graphically the dynamics of the proposed model to observe the effect of vaccination by altering the vaccination rate suggest that, mass vaccination along with the implementation of the control measures until the achievement of complete herd immunity is needed to control the spread of the disease. Hence, majority of susceptible population should be vaccinated to prevent the new waves of COVID-19. The vaccination of the susceptible class will control the spread of the disease resulting in reduction of the mortality rate, which can clearly be seen by recovered class presented in Figure 4. Simulation results show that immunization is an essential strategy for lowering the incidence and severity of COVID-19 infections in the workplace, in communities, and across the globe. Therefore, it is advocated for policymakers to increase public trust and confidence in the efficiency and security of vaccinations as well as the expertise and dependability of the institutions that provide them. Figures 5 and 6 show a comparison of integer order model, fractional model for fractional order , and real data of total number of deaths owing to the COVID-19 outbreak in Italy. Figure 5 presents the comparison for first days starting from February 2020 to May 2020 of the pandemic before vaccination. We can see that after 40 days, real data is better explained by fractional order model in comparison to integer order. Figure 6 is assigned to display the comparison from July 2021 to September 2021 along with administered vaccination effect considering the initial values of the compartments as of July 2021 . The relevant COVID-19 data can be found at https://github.com/pcmdpc/COVID-19. Also we can see from our proposed model, fractional order models are better at fitting the real data than classical models because fractional order derivative and integral operators, due to their nonlocal nature, depend not only on their current state but also on all of their previous states, making them more effective than other classical deterministic operators in predicting the biological model’s future state.

7. Concluding Remarks

Within the Atangana-Baleanu Caputo fractional operator, we studied alternative conditions for qualitative characterization and dynamics of the suggested model. Calculating disease-free and endemic equilibrium points, as well as the reproduction number, are used to examine the prerequisites for local stability. A newly developed numerical scheme is used to obtain the solution of the concerned fractional order model. The proposed scheme’s stability and convergence are briefly examined. To depict the behaviour of the COVID-19 virus, a graphical description of the dynamics of the relevant system is developed within different fractional orders under the effect of vaccination suggesting the need of mass vaccination to gain the herd immunity in order to prevent the further spread of the disease. Due to the limited access of the available vaccine, other preventive measures should also be continued.

Data Availability

The data used to support the findings of this study can be found at https://ourworldindata.org/covid-vaccinations.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

References

Z. Zhang, “A novel covid-19 mathematical model with fractional derivatives: singular and nonsingular kernels,” Chaos, Solitons and Fractals, vol. 139, article 110060, 2020.
View at: Publisher Site | Google Scholar
S. Rezapour, H. Mohammadi, and M. E. Samei, “SEIR epidemic model for COVID-19 transmission by Caputo derivative of fractional order,” Advances in Difference Equations, vol. 2020, no. 1, 2020.
View at: Publisher Site | Google Scholar
A. Atangana and A. S. İğret, “Mathematical model of COVID-19 spread in Turkey and South Africa: theory, methods, and applications,” Advances in Difference Equations, vol. 2020, no. 1, 89 pages, 2020.
View at: Publisher Site | Google Scholar
Z. Zhang, A. Zeb, O. F. Egbelowo, and V. S. Erturk, “Dynamics of a fractional order mathematical model for COVID-19 epidemic,” Adv. Difference Equ., vol. 2020, no. 1, p. 420, 2020.
View at: Publisher Site | Google Scholar
S. Ahmad, Q. M. Aman Ullah, H. Al-Mdallal, K. Khan, and A. K. Shah, “Fractional order mathematical modeling of COVID-19 transmission,” Chaos, Solitons and Fractals, vol. 139, article 110256, 2020.
View at: Publisher Site | Google Scholar
O. Zakary, S. Bidah, M. Rachik, and H. Ferjouchia, “Mathematical model to estimate and predict the covid-19 infections in morocco: optimal control strategy,” Journal of Applied mathematics, vol. 2020, Article ID 9813926, 13 pages, 2020.
View at: Publisher Site | Google Scholar
Z. Ahmad, M. Arif, F. Ali, I. Khan, and K. S. Nisar, “A report on COVIDâ€‘19 epidemic in pakistan using SEIR fractional model,” Scientific Reports, vol. 10, no. 1, 2020.
View at: Publisher Site | Google Scholar
S. S. Askar, D. Ghosh, P. K. Santra, A. Abdelalim, and G. S. M. Elsadany, “A fractional order SITR mathematical model for forecasting of transmission of COVID-19 of India with lockdown effect,” Results in Physics, vol. 24, article 104067, 2021.
View at: Publisher Site | Google Scholar
S. Bushnaq, T. Saeed, D. F. M. Torres, and A. Zeb, “Control of COVID-19 dynamics through a fractional-order model,” Alexandria Engineering Journal, vol. 60, no. 4, pp. 3587–3592, 2021.
View at: Publisher Site | Google Scholar
M. A. Abdulwasaa, M. S. Abdo, K. Shah et al., “Fractal-fractional mathematical modeling and forecasting of new cases and deaths of COVID-19 epidemic outbreaks in India,” Results in Physics, vol. 20, article 103702, 2021.
View at: Publisher Site | Google Scholar
I. A. Baba and A. Rihan, “A fractional-order model with different strains of COVID-19,” Physica A, vol. 603, article 127813, 2022.
View at: Publisher Site | Google Scholar
F. A. Rihan and V. Gandhi, “Dynamics and sensitivity of fractional-oder delay differential model for coronavirus (COVID-19) infection,” Progress in Fractional Differentiation and Applications, vol. 7, no. 1, pp. 43–61, 2021.
View at: Google Scholar
F. A. Rihan and H. J. Alsakaji, “Dynamics of a stochastic delay differential model for COVID-19 infection with asymptomatic infected and interacting people: case study in the UAE,” Results in Physics, vol. 28, article 104658, 2021.
View at: Publisher Site | Google Scholar
K. M. Furati, I. O. Sarumi, and A. Q. M. Khaliq, “Fractional model for the spread of COVID-19 subject to government intervention and public perception,” Applied Mathematical Modelling, vol. 95, pp. 89–105, 2021.
View at: Publisher Site | Google Scholar
R. T. Alqahtani, “Mathematical model of SIR epidemic system (COVID-19) with fractional derivative: stability and numerical analysis,” Advances in Difference Equations, vol. 2021, no. 1, 2021.
View at: Publisher Site | Google Scholar
K. Hattaf, “A new generalized definition of fractional derivative with non-singular kernel,” Computation, vol. 8, no. 2, p. 49, 2020.
View at: Publisher Site | Google Scholar
K. Hattaf, “On the stability and numerical scheme of fractional differential equations with application to biology,” Computation, vol. 10, no. 6, p. 97, 2022.
View at: Publisher Site | Google Scholar
K. Regassa Cheneke, K. Purnachandra Rao, and E. G. Kenassa, “A new generalized fractional-order derivative and bifurcation analysis of cholera and human immunodeficiency co-infection dynamic transmission,” International Journal of Mathematics and Mathematical Sciences, vol. 2022, Article ID 7965145, 15 pages, 2022.
View at: Publisher Site | Google Scholar
https://ourworldindata.org/covid-vaccinations.
X. Meng, Z. Cai, S. Si, and D. Duan, “Analysis of epidemic vaccination strategies on heterogeneous networks: based on SEIRV model and evolutionary game,” Applied Mathematics and Computation, vol. 403, article 126172, 2021.
View at: Publisher Site | Google Scholar
A. M. Ramos, M. Vela-Pérez, M. R. Ferrández, A. B. Kubik, and B. Ivorra, “Modeling the impact of SARS-CoV-2 variants and vaccines on the spread of COVID-19,” Communications in Nonlinear Science and Numerical Simulation, vol. 102, article 105937, 2021.
View at: Publisher Site | Google Scholar
J. N. A. Tetteh, V. K. Nguyen, and E. A. Hernandez-Vargas, “COVID-19 network model to evaluate vaccine strategies towards herd immunity,” 2020, https://www.medrxiv.org/content/10.1101/2020.12.22.20248693v1.
View at: Google Scholar
I. Harizi, S. Berkan, and A. Tayebi, “Modeling the effect of population-wide vaccination on the evolution of COVID-19 epidemic in Canada,” 2020, https://www.medrxiv.org/content/10.1101/2021.02.05.21250572v2.
View at: Google Scholar
A. Karabay, A. Kuzdeuov, M. Lewis, and H. A. Varo, “A vaccination simulator for COVID-19: effective and sterilizing immunization cases,” 2021, https://www.medrxiv.org/content/10.1101/2021.03.28.21254468v2.
View at: Google Scholar
G. S. Tuteja, “Stability and numerical investigation of modified SEIR model with vaccination and life-long immunity,” European Journal of Molecular and Clinical Medicine., vol. 7, pp. 3034–3044, 2020.
View at: Google Scholar
R. Ghostine, M. Gharamti, S. Hassrouny, and I. Hoteit, “An extended SEIR model with vaccination for forecasting the COVID-19 pandemic in Saudi Arabia using an ensemble Kalman filter,” Mathematics, vol. 9, no. 6, p. 636, 2021.
View at: Publisher Site | Google Scholar
P. Wintachaia and K. Prathomb, “Stability analysis of SEIR model related to efficiency of vaccines for COVID-19 situation,” Heliyon, vol. 7, no. 4, article e06812, 2021.
View at: Publisher Site | Google Scholar
H. J. Alsakaji, F. A. Rihan, and A. Hashish, “Dynamics of a stochastic epidemic model with vaccination and multiple time- delays for COVID-19 in the UAE,” Complexity, vol. 2022, Article ID 4247800, 15 pages, 2022.
View at: Publisher Site | Google Scholar
L. Kou, X. Wang, Y. Li, X. Guo, and H. Zhang, “A multi-scale agent-based model of infectious disease transmission to assess the impact of vaccination and non-pharmaceutical interventions: the COVID-19 case,” Journal of Safety Science and Resilience, vol. 2, no. 4, pp. 199–207, 2021.
View at: Publisher Site | Google Scholar
M. A. Acuña-Zegarra, S. Díaz-Infante, D. Baca-Carrasco, and D. Olmos-Liceaga, “COVID-19 optimal vaccination policies: a modeling study on efficacy, natural and vaccine-induced immunity responses,” Mathematical Biosciences, vol. 337, article 108614, 2021.
View at: Publisher Site | Google Scholar
V. V. Albani, J. Loria, E. Massad, and J. P. Zubelli, “The impact of COVID-19 vaccination delay: a data-driven modeling analysis for Chicago and New York City,” Vaccine, vol. 39, no. 41, pp. 6088–6094, 2021.
View at: Publisher Site | Google Scholar
A. Atangana and D. Baleanu, “New fractional derivatives with nonlocal and non-singular kernel: theory and application to heat transfer model,” Thermal Science, vol. 20, no. 2, pp. 763–769, 2016.
View at: Publisher Site | Google Scholar
J. Huang and D. Yang, “A unified difference-spectral method for time-space fractional diffusion equations,” International Journal of Computer Mathematics, vol. 94, no. 6, pp. 1172–1184, 2017.
View at: Publisher Site | Google Scholar
P. Van den Driessche and J. Watmough, “Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission,” Mathematical Biosciences, vol. 180, no. 1-2, pp. 29–48, 2002.
View at: Publisher Site | Google Scholar
W. Lin, “Global existence theory and chaos control of fractional differential equations,” Journal of Mathematical Analysis and Applications, vol. 332, no. 1, pp. 709–726, 2007.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2022 Mubashara Wali et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

262

Downloads

429

Citations