OJMSiOpen Journal of Modelling and Simulation2327-4018Scientific Research Publishing10.4236/ojmsi.2017.51008OJMSi-73541ArticlesPhysics&Mathematics Global Stability of a SVEIR Epidemic Model: Application to Poliomyelitis Transmission Dynamics L.N. Nkamba1J.M. Ntaganda2H.Abboubakar3J.C. Kamgang4LorenzoCastelli5DIA-University of Trieste, Trieste, ItalyDepartment of Mathematics, Higher Teacher Training College, University of Yaound&amp;eacute; I, Yaound&amp;eacute;, CameroonLaboratoire d’Analyse, Simulation et Essai (LASE), Ngaound&amp;eacute;r&amp;eacute;, CameroonDepartment of Mathematics, School of Science, College of Science and Technology, University of Rwanda, Butare, RwandaDepartment of Mathematics and Computer Science, ENSAI-University of Ngaoundere, Ngaound&amp;eacute;r&amp;eacute;, Cameroon09122016050198112November 12, 2016Accepted: January 15, January 18, 2017© Copyright 2014 by authors and Scientific Research Publishing Inc. 2014This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/

The lack of treatment for poliomyelitis doing that only means of preventing is immunization with live oral polio vaccine (OPV) or/and inactivated polio vaccine (IPV). Poliomyelitis is a very contagious viral infection caused by poliovirus. Children are principally attacked. In this paper, we assess the impact of vaccination in the control of spread of poliomyelitis via a deterministic SVEIR (Susceptible-Vaccinated-Latent-Infectious-Removed) model of infectious disease transmission, where vaccinated individuals are also susceptible, although to a lesser degree. Using Lyapunov-Lasalle methods, we prove the global asymptotic stability of the unique endemic equilibrium whenever . Numerical simulations, using poliomyelitis data from Cameroon, are conducted to approve analytic results and to show the importance of vaccinate coverage in the control of disease spread.

Deterministic SVEIR Model Poliomyelitis Imperfect Vaccine Direct Lyapunov Method Equilibrium States Global Stability
1. Introduction

In the 70s, having noticed that five million children died every year further to an avoidable disease by the vaccination like poliomyelitis, the WHO introduced the Global Immunization Vision and Strategy (GIVS). Poliomyelitis has been eliminated in the most of countries, but recently we observe the upsurge of infectious in some countries  . Since October 2013, Cameroon is classified by the WHO as the exporting country of the poliovirus  . Poliomyelitis is an acute and sometimes devastating viral disease very contagious caused by poliovirus. Human is the only natural host for poliovirus  . Children are principally attacked. Poliovirus is predominantly transmitted via mother and food contaminated. In the most of case, infection is asymptomatic but the persons infected can transmit disease via their feces  . When a susceptible is exposed to infection by a virulent poliovirus, we can observe few days or few weeks three types of responses (minor illness, aseptic meningitis, and paralytic poliovirus). In case of minor illness, after 3 - 5 days, symptoms can be slight, fever, tiredness, headache, sore throat and vomiting. In the minor illness, the patient recovers in a few days 24 to 72 hours. In the case of non paralytic poliomyelitis in addition in some of minor illness signs and symptoms includes stiffness and pain in the back of neck. In the past days of illness, healing will rapid and complete. In the paralytic poliomyelitis, the predominant damage is flaccid paralysis resulting from lower motor neurons damage. The maximal recovery usually occurs after 6 months, but residuals paralysis lasts much longer. There does not exist a specific treatment for poliomyelitis although improved sanitation and hygiene help to limit the spread of poliovirus. The only specific means of preventing polio is immunization with live polio vaccine (OPV) or/and inactivated polio vaccine (IPV)     .

As part of the necessary multi-disciplinary research approach, mathematical models have been extensively used to provide a framework for understanding of poliomyelitis transmission dynamics and the best strategies to control the spread of infection in the human population. In the literature, considerable work can be found on the mathema- tical modeling of poliomyelitis  -  . Some of these works refer to vaccination as polio control mechanism     , using a standard SEIR model  .

Some SVEIR models are used to assess the potential impact of an imperfect SARS vaccine like SARS vaccine  , Hepatitis B vaccine  , Tuberculosis vaccine  , HIV vaccine   , to mention only these four diseases. From a mathematical point of view, to show the global asymptotic stability of equilibrium points in general, and especially, the global asymptotic stability of the endemic equilibrium, is not an easy task. This requires, in most cases, the use of several different techniques, such as the theory of compound matrix   , the comparison theorem  , or the use of Lyapunov functions associated with the Lassalle invariance principle  , to name a few techniques commonly used by authors. For example, in  , the authors used compound matrix techniques to show the global stability of the endemic equilibrium under some constraints on the parameters of the system. Huiming Wei et al.  proposed an SVEIR model with time delay, and analyzed the dynamic behavior under pulse vaccination. Using comparison theorem, they showed that the infection-free periodic solution is globally attractive. Yu Jiang et al.  modified that model by adding saturation incidence, and used too the comparison theorem to show the global stability of “infection-free” periodic solution.

In this paper, we study the impact of vaccination in the control of poliomyelitis spread via an SVEIR model of infectious disease transmission. Individuals are classified as one of susceptible, vaccinated, exposed, infectious, or recovered. The model is based on a standard SEIR model  , but allows that susceptible individuals may be given an imperfect vaccine that reduces their susceptibility to the disease. Since we consider a leaky vaccine, the V-compartment of vaccinated individuals is considered as a susceptible compartment, and thus we are dealing with a differential susceptibility system with bilinear mass action as in Hyman and Li  . However, we include one-way flow between these two compartments due to vaccination making the model studied here distinct from the model in  . For the case where the basic reproduction number is less than one, the global stability of the disease-free equilibrium has been shown by Gumel et al. in 2006  . However, the global dynamics when the basic reproduction number is greater than one have not been resolved before. By allowing different death rates for each of the compartments, the model studied in this paper is slight generalization of the model studied in  . Using Lyapunov-LaSalle methods, we fully resolve the global dynamics of the model for the full parameter space. We demonstrate that the model exhibits threshold behavior with a globally stable disease-free equilibrium if the basic reproduction number is less than unity and a globally stable endemic equilibrium if the basic reproduction number is greater than unity. Thus, we also fully resolve the global dynamics for the model studied in  .

In order to study the stability of a positive endemic equilibrium state, we use Lyapunov’s direct method and LaSalle’s Invariance Principle with a Lyapunov function of the form:

where are constants, is the population of ith compartment and is the equilibrium level. Lyapunov functions of this type have also proven to be useful for Lotka-Voltera predator-preys systems  , and it appears that they can be useful for more complex compartmental epidemic models as well   .

The main aim of the present paper is to show that our model has a unique endemic equilibrium which is globally asymptotically stable.

This SVEIR model could be used to assess the potential impact of an extended vaccination program (such as for the monovalent serogroup A conjugate MenVacAfric, an anti-meningococcal vaccine introduced in 2011 in Sub-saharan Africa), in order to compare with the impact of a pulse vaccination program.

In the next section, we present our SVEIR epidemic model. Section 3 presents some basic properties like the computation of the basic reproduction ration, , and such as the existence of the equilibrium points. In Section 4, we study the stability properties of the model and in Section 7, numerical simulations will be done with Cameroon data which deal with the vaccination campaign against polio. An conclusion round the paper.

2. Model Description

We divide the entire population into 5 sub-populations of epidemiological significance: susceptible, vaccinated, exposed, infective, and removed compartments with respective sizes, , , and. The latent compartment, , takes into account the delay between the moment of the infection and the moment when an infected individual becomes infectious. The per capita death rates for susceptible, vaccinated, exposed, infective and recovered individuals are, , , and, respectively. The recruitment rate into the susceptible class is assumed to be constant and denoted by. The per capita vaccination rate is.

We assume mass action incidence for susceptible. Vaccination reduces the risk of infection by a factor. Thus, we have mass action incidence for vaccinated individuals and the efficacy of the vaccine is. The case corresponds to a perfect vaccine and the case corresponds to a vaccine with no effect. Each of these cases can be dealt with more simply and directly by studying the basic SEIR model.

The average duration of latency in class before progressing to class is, and the average time spent in class before recovery is. All parameters of the system are assumed to be positive.

Our model consists of the following system of ordinary differential equations:

with initial conditions which satisfy. A schematic of the model is shown in Figure 1.

Since does not appear in the equations for the other variables, we will consider the following system (model system (3) without the compartment):

Schematic of the compartmental model

with initial conditions which satisfy.

3. Basic Properties and Equilibriums3.1. A Compact Positively Invariant Absorbing Set

In order that the model be well-posed, it is necessary that the state variables, , and remain nonnegative for all. That is, the nonnegative orthant must be positively invariant. Let

where, and.

Lemma 1. The compact set is a positively invariant and attracting.

Proof. For each of the variables, , and, when the variable is equal to zero, the derivative of that variable is non-negative in. It then follows from (  , Proposition 2.1) that is positively invariant.

Let. Then

Consequently,

Similarly, when and so

Let. Then for a given initial condition, there exists such that for all. Then,

for. Thus,

This holds for all and so

W

Since is a positively invariant absorbing set is sufficient to consider the dynamics of the flow generated by system (3) in.

It is easy to see that the model system (3) has a disease-free equilibrium given by

Additionally, an endemic equilibrium may also exist.

3.2. Basic Reproduction Ratio and Equilibrium

Using the method of the references   , the basic reproduction number is

Replacing and by their values in (9), we obtain:

When there is no vaccination (), system (3) is the standard model with

From Equation (10), we claim the following result.

Proposition 1. if and only if.

Proof. It follows from (11) that

Thus, is equivalent to

from which the result follows. W

The value of determines whether or not there exists an endemic equilibrium (  , Theorem 2.3).

Theorem 1. If, then there are no endemic equilibria. If, then there exists a unique endemic equilibrium).

(See Appendix for proof).

4. Stability Analysis of Equilibriums4.1. Stability Analysis of the DFE

For local stability of the disease-free equilibrium, we claim the following:

Theorem 2. If, then the disease-free equilibrium is locally asymptotically stable and unstable if.

Proof. The Jacobian matrix of model (3) evaluate at the disease-free equilibrium is given by

The eigenvalues of are, , and those of the following sub-matrices

The characteristic polynomial of is given by

It clear that the roots of have negative real parts if and only if. It follows that the disease-free equilibrium is locally asymptotically stable whenever and unstable when. This end the proof. W

The following result is proven in (  , Theorem 4.1).

Theorem 3. If, then the disease-free equilibrium is globally asymptotically stable.

If, then the disease-free equilibrium is unstable.

4.2. Stability Analysis of the Endemic Equilibrium

Our main result is the following theorem.

Theorem 4. If, then the endemic equilibrium point is globally asymptotically stable in.

Proof. Consider the following candidate Lyapunov function

Differentiating along solutions to (3) gives:

Since

and,

Since arithmetical mean is greater than geometrical mean, we have the following inequalities

Therefore. Thank’s to the direct Lyapunov theorem of stability, we conclude that is stable.

It remain to prove that is asymptotically stable using the Lasalle invariance principle.

set

it’s clear that;

Backing to the above relations, we have the following implications.

If we set

Finally we have,

At the endemic equilibrium, we have

Replacing by their values given by (24) in the second equation of system (25) yields

If we compare relation (26) with the last equation of (25), then we have:

Consequently:

Finally

Thus, the largest invariant set contained in is .

Then the global stability of follows according to the Lasalle invariance principle  . W

5. Numerical Simulations

In this section we show via numerical simulations that when is lower than one (minor illness), disease will be eliminated from the community, and when is greater than one (meningitis and paralytic form of polio), and epidemics will occurs or the disease will persists in the community. We explore also the impact of vaccination coverage in the spread of poliomyelitis.

Parameters Description and Values

Most of parameters values are from Cameroon, like natural rate of mortality. We assume that the natural rates of mortality of susceptible, recovered, exposed are the same. Value of vaccine efficacy, recovery rate and rate of apparition of clinical symptoms are coming from WHO. For vaccination coverage, we take different values in order to explore different situations. The recruitment rate of susceptible humans, , likely is actually the birth rate, and are taken in   . See Table 1 for the description of parameters and their based line or range value.

6. Numerical Results and Interpretations

Figure 2 illustrate the minor illness form of polio. We assume that, so , and we have showed analytically that If, then the disease-free equilibrium is globally asymptotically stable. We see that in this case, healthy carriers and infectious tend toward horizontal axis, and the infection is eradicated after around 6 months.

In Figure 3, we are in the presence of the meningitis form of polio. Assuming that

Description and values of parameters of model (3)
ParameterDescriptionBased line value or range
Recruitment rate of susceptible2.5
Effective contact rate0.1
Vaccine efficacy0.8
Rate of development of clinical symptoms0.05 - 0.5
Recovery rate0.05
Vaccination coverage rate
Natural mortality rate of susceptible0.0551
Natural mortality rate of vaccinated0.0551
Natural mortality rate of exposed0.0551
Mortality rate of infectious0.08
Natural mortality rate of recovered0.0551

and vaccine coverage, to have. It is clear that infection is a little more severe and the disease reaches at endemic equilibrium point and does not disappear.

In Figure 4, we are in the presence of the most severe form of polio: the paralytic form with, so. As in the case of meningitis form, the patient takes long time to heal and thus continue to transmit the infection during that time. It is important to note that remark is that the infection takes longer to reach the endemic equilibrium point and remains in the population despite vaccination.

We are in front of paralytic polio. We assume, and explore the effect of immunization on the dynamic of the disease. Figure 5 show that more vaccine coverage is high, the number of healthy carriers and infectious is low at equilibrium point. But it is noted that the infection remains in the population.

Figure 6, we explored three cases:

1) even if the vaccine is perfect and nobody is vaccinated; the infection is and remains high in the population and;

2) The vaccination is made; even if the coverage is low infection decreases and reaches a an equilibrium point and;

3) The last and not realistic situation is that infection is eradicated after one year, and when we have perfect vaccine and maximal vaccination coverage and.

Paralytic form of polio<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/8-2860107x178.png" xlink:type="simple"/></inline-formula> Impact of vaccine coverage Impact of vaccine efficacy
7. Conclusions

We highlighted in this article the importance of vaccination in the control of the propagation of the poliomyelitis. We relied on the compartmentalized SVEIR model that characterizes the infectious diseases. We computed, key parameter related to the Reproduction, which governs the asymptotic behavior of the model. We then constructed a Lyapunov function to prove the global asymptotic stability of the endemic equilibrium whenever.

Using data from AHALA (district of Yaound in Cameroon), we simulated the three different forms of polio namely the minor illness, the meningitis form and the paralytic form. In the case of minor illness of polio, we assumed that. The model also allowed an endemic equilibrium point when is greater than 1. In that case, we simulated both meningitis and paralytic form of polio, respectively with and. We found that, the more the vaccine coverage is high, the more the healthy Carriers and Infectious are low. The simulations show that, to eradicate polio in the population means to have simultaneously a perfect vaccine and maximal vaccine coverage. Therefore, other control strategies are to be issued to finally reach that goal.

Acknowledgements

The first author acknowledges with thanks the High teacher Training College of Yaounde. H. A. would like to thank the Direction of the University Institute of Technology of Ngaoundere for their financial assistance in the context of research missions of September 13, 2016.

Cite this paper

Nkamba, L.N., Nta- ganda, J.M., Abboubakar, H., Kamgang, J.C. and Castelli, L. (2017) Global Stability of a SVEIR Epidemic Model: Application to Po- liomyelitis Transmission Dynamics. Open Journal of Modelling and Simulation, 5, 98- 112. http://dx.doi.org/10.4236/ojmsi.2017.51008

AppendixProof of Theorem 1

Proof. In order to determine the existence of possible endemic equilibrium, that is, equilibrium with all positive components which we denote by

we have to look for the solution of the algebraic system of equations obtained by equating the right hand sides of system (3) to zero. In this way we obtain the implicit system of equations,

where is solution of the following equation

with,

and

Note that coefficient is always negative and coefficient is positive (resp. negative) if and only if is greater (less) than unity. Thus, model system (3) admits only one endemic equilibrium whenever the basic reproduction number is greater than unity. When, we have negative. It follows that the model system (3) does not have any endemic equilibrium point whenever. W

Submit or recommend next manuscript to SCIRP and we will provide best service for you:

A wide selection of journals (inclusive of 9 subjects, more than 200 journals)

Providing 24-hour high-quality service

User-friendly online submission system

Fair and swift peer-review system