# Global Stability Results in a SVIR Epidemic Model with Immunity Loss Rate Depending on the Vaccine-Age

Authors: Raúl Peralta; Cruz Vargas-De-León; Pedro Miramontes
Journal: Abstract and Applied Analysis
 (2015)

---

## Abstract

We formulate a susceptible-vaccinated-infected-recovered (SVIR) model by incorporating the vaccination of newborns, vaccine-age, and mortality induced by the disease into the SIR epidemic model. It is assumed that the period of immunity induced by vaccines varies depending on the vaccine-age. Using the direct Lyapunov method with Volterra-type Lyapunov function, we show the global asymptotic stability of the infection-free and endemic steady states.


---

## Body

## 1. Introduction

The persistence of protection induced by vaccines plays a central role in the implementation of human vaccination policies. Two types of evidence document the duration of protection induced by a vaccine: on the one hand, immunological parameters such as antibody kinetics over time and on the other hand the disease incidence for which the vaccine was implemented [1]. In this context, measuring the persistence of circulating antibodies is widely used in human vaccines. The specific antibodies induced by the vaccine are the main parameters used for monitoring the duration of the vaccine-induced immune response. This monitor uses population-based studies that compare antibody levels in a protected group with those of a susceptible group, leading to an estimation of disease risk as a function of postvaccination antibody titers. For several vaccines, data resulting from long-term follow-up vaccine studies have been used to create mathematical models to predict the kinetics of antibody persistence over longer periods of time [2]. An example of such modeling results is the recent publication of the duration of protection induced by a human papillomavirus vaccine, while the trials provided data over a 6.4 year follow-up period, the models predicted a long protection [2]. Further, the waning of protection can also be demonstrated by an increased of disease incidence among subjects previously vaccinated. For example, the implementation of a universal varicella vaccination program in the USA in 1995 has resulted in a substantial reduction in morbidity. However, despite this reduction, new cases of varicella continue to occur, mostly in highly vaccinated school communities. Several studies of these cases have suggested that the time since vaccination may be associated with the risk of varicella. Waning of immunity after varicella vaccination in terms of measurable antibodies has been demonstrated to occur in children previously vaccinated. In this regard, the time since vaccination and the age at vaccination were identified as predictors of this loss of immunity [3–5]. The number of vaccine-preventable diseases targeted by routine vaccination schemes is growing, at time, at least 17 vaccine-preventable diseases were targeted by routine vaccination schemes in USA, but potentially all diseases of infectious etiology can be controlled by vaccination schemes and the evaluation of the persistence of vaccine-induced protection with human vaccines is essential for the susceptible population to determine cases of loss protection. The prevention of infection correlates with the induction of specific antibodies. Although loss of antibody after vaccination may render vaccines again susceptible to some infections. The time of vaccine-induced immunity may vary according to individual characteristics and the vaccine-age [6]. The vaccination of newborns could have a central role in the evolution of vaccine-induced immunity [7].In most infectious disease models, it has been assumed that all vaccinated individuals have equal immunity provided by the vaccine during their vaccination period. However, the duration of protection provided by vaccination varies and is influenced by many factors, particularly the individual characteristics and the vaccine-age [6, 7]. Therefore, there are enough reasons to study the possible effects of loss of vaccine-induced immunity on the dynamics of infectious diseases. Recently, the vaccine-age-dependent epidemic models have been considered [8–10], more precisely, the vaccine-age is used to track the vaccine-protection period and the vaccine wane rate.Iannelli et al. in [9] developed an SIS epidemic model with super-infection, perfect vaccination and vaccine-age and discussed the asymptotical behavior of the equilibria. Li et al. in [10] modeled an SIS disease with vaccine-age and treatment rates and exhibited backward bifurcations under some conditions. Duan et al. in [8] formulated an SIR epidemic model with age of vaccination and proved the global stability of the equilibria. In these models the susceptible population is subjected to a vaccination scheme. Here, we will consider that a proportion of individuals get vaccinated immediately after birth or in newborn.In this paper, by incorporating the vaccination of newborns in vaccine-age and disease related death rate into the SIR model, we formulate a SVIR epidemiological model with a waning rate of vaccine-induced immunity that depends on the vaccine-age. Here, we will analyze the global stability of the steady states in our age-structured model by means of suitable Volterra-type Lyapunov functions of the form:(1)GX=X-1-ln⁡X.Recently, McCluskey and coauthors in [11, 12] proved the global stability of the endemic steady state of SIR models with infection-age-dependent by means of Volterra-type Lyapunov functions. Melnik and Korobeinikov [13] used them to establish global stability for infectious diseases models with age-dependent susceptibility. The global-stability problem for basic age-structured models for the transmission of a vector-borne infectious disease is solved using the Volterra-type functions by Vargas-De-León and colleagues [14].The remainder of this paper is organized as follows. Section2 is devoted to the model derivation. In Section 3 the existence of steady states is proved and an explicit formula for the reproductive number is derived. In Section 4, we will give a very simple and elegant proof of the global stability of steady states. Section 5 contains the discussion and concluding remarks.


## 2. SIR Model with Vaccination of Newborns

In this section, we introduce an SIR epidemic model and its reformulations. We can assign each vaccinated individual its class age, called vaccine-age and denote it byθ, which is the amount of time that has passed since the moment that initiates the vaccination.The model divides the total population into the following subgroups that are unvaccinated susceptible individualsS(t), infectious individuals I(t), recovered individuals R(t), and vaccinated individuals V(t). We stratify the vaccinated part of the population V(t) according to vaccine-age(2)V(t)=∫0+∞v(θ,t)dθ,
where v(θ,t) denotes the distribution of the vaccinated individuals at time t. In this context, then, ∫θ1θ2v(θ,t)dθ is simply the number of vaccinated individuals between vaccine-ages θ1 and θ2 at time t.We assume a simple demographic process in which newly recruited individuals (corresponding to births) enter the population at a rateΛ and exit (due to natural death) at a rate μ. A proportion 0<p<1 of the newly born individuals are vaccinated at birth and appear in the vaccinated class with class age 0, giving rise to a boundary condition v(0,t)=pΛ. The parameter β is the successful contact rate between infected and noninfected individuals, γ is the recovery rate and η is the disease related death rate. The function 0<σ(θ)≤1 denotes the vaccine efficacy, which we assume is a decreasing function of vaccine-age, θ. The function α(θ) denotes the rate at which the vaccine wanes that depends on the vaccine-age (and such vaccinated individuals move to the unvaccinated susceptible class). We assume that α(θ) and is nonnegative and bounded integrable function.Combining the above formulations and assumptions, it follows that the model for the transmission dynamics of the disease in the presence of loss of vaccine-induced immunity is given by the following system of nonlinear partial differential equations:(3)dS(t)dt=(1-p)Λ-βS(t)I(t)-μS(t)+∫0+∞α(θ)v(θ,t)dθ,dI(t)dt=βS(t)I(t)+∫0+∞1-σθβvθ,tItdθ-(μ+γ+η)I(t),∂v(θ,t)∂t+∂v(θ,t)∂θ=-μ+αθv(θ,t)-1-σθβv(θ,t)I(t),
with the following boundary condition v(0,t)=pΛ and initial conditions S(0)=S0>0, I(0)=I0>0, and v(θ,0)=v0(θ)>0, where S0 and I0 are initial populations of the susceptible individuals and infectious individuals, respectively, and an initial distribution of the vaccinated individuals with respect to their vaccine-age v0(θ) is a Lebesgue integrable function. The equation for R(t) is omitted here.In the following section, we will discuss existence of steady states, namely, infection-free steady state and endemic steady state. Also, we will derive an explicit formula for the reproductive number.


## 3. Steady States and Effective Reproductive Number

The epidemic systems can have one or two nonnegative steady states. Specifically, the epidemic system always has an infection-free steady state; in this steady state the component of infectious individuals is zero. Also can have an endemic steady state with all components positive.The system (3) has an infection-free steady state E0=(S0,I0,v0(θ)), where
(4)S0=Λμ1-p+p∫0+∞α(θ)K0(θ)dθ,I0=0,v0(θ)=pΛK0(θ),
with
(5)K0(θ):=exp⁡-∫0θ(μ+α(τ))dτ.K0(θ) is the survival rate at duration θ in the vaccinated class.A possible formula for the reproductive number can be derived by determination of the condition for the instability of infection-free steady stateE0. Here, we will derive the reproductive number by means of linearization around E0, and its condition of instability.Then we perturb the infection-free steady state by letting(6)xt=St-S0,yt=It,z(θ,t)=v(θ,t)-v0(θ).Substituting the variables above in (3) and neglecting terms of order bigger than two, we obtain the following linear system for x(t), y(t), and z(θ,t):
(7)dx(t)dt=-βS0y(t)-μx(t)+∫0+∞α(θ)z(θ,t)dθ,dytdt=βS0y(t)+βy(t)·∫0+∞1-σθv0(θ)dθ-(μ+γ+η)y(t),∂z(θ,t)∂t+∂z(θ,t)∂θ=-μ+αθzθ,t-1-σθβv0θyt.We consider exponential solutions of the form:(8)xt=x¯exp⁡λt,yt=y¯exp⁡λt,z(θ,t)=z¯(θ)exp⁡(λt),
where x¯ and y¯ are arbitrary constants, z¯(θ) is a function of θ, and λ is a real or a complex number. Substituting them in system (7), we obtain
(9)λx¯=-βS0y¯-μx¯+∫0+∞α(θ)z¯(θ)dθ,λy¯=βS0y¯+βy¯∫0+∞1-σθv0(θ)dθ-(μ+γ+η)y¯,λz¯(θ)+∂z¯θ∂θ=-μ+αθz¯(θ)-1-σθβv0(θ)y¯,z(t,0)=0.
If the system (9) has at least one positive eigenvalue then the infection-free steady state E0 is unstable. So one eigenvalue of (9) is
(10)λ1=βS0+β∫0+∞1-σθv0(θ)dθ-(μ+γ+η),=(μ+γ+η)·βμ+γ+ηS0+∫0+∞1-σθv0(θ)dθ-1.
Due to space constraints, we omit the proof that the other two eigenvalues (λ2 and λ3) of system (9) are negative (or negative real part). Consider
(11)Ifβμ+γ+ηS0+∫0+∞1-σθv0(θ)dθ>1,thentheeigenvalueλ1ispositive.
Setting
(12)Rv:=βS0μ+γ+η+βμ+γ+η∫0+∞1-σθv0(θ)dθ
yields the expression for theeffective reproductive number of the disease. In interpreting this expression, note that 1/(μ+γ+η) is the average residence time in the infectious class. The first term βS0/(μ+γ+η) is the average number of secondary infections produced by a single infectious individual during its infectious period among susceptible population. The second term β∫0+∞(1-σ(θ))v0(θ)dθ/(μ+γ+η) is the average number of secondary infections produced by a single infectious individual during its infectious period among vaccinated population. Hence the total expression in (12) represents the average number of secondary infections caused by a single infectious individual introduced in an entirely susceptible and vaccinated population, over its entire infectious period, 1/(μ+γ+η).We will now turn to study the conditions of existence of the endemic steady state. System (3) has an endemic steady state E*=(S*,I*,v*(θ)), where S*,  I*, and v*(θ) satisfy the equalities
(13)0=(1-p)Λ-βS*I*-μS*+∫0+∞α(θ)v*(θ)dθ,(14)0=βS*+∫0+∞1-σθβv*(θ)dθ-(μ+γ+η),(15)dv*(θ)dθ=-μ+αθv*(θ)-1-σθβv*(θ)I*,(16)v*(0)=pΛ.
From (13), it follows that
(17)S*=1μ+βI*(1-p)Λ+∫0+∞α(θ)v*(θ)dθ=1μ+βI*·1-pΛ∫0+∞αθv0θexp⁡-∫0θ1-σθβI*dτhh+∫0+∞αθv0θexp⁡-∫0θ1-σθβI*dτdθ.Equations (15) and (16) give
(18)v*θ=pΛexp⁡-∫0θ(μ+α(τ)+1-σθβI*)dτ=v0(θ)exp⁡-∫0θ1-σθβI*dτ.Combining this with (14), (17), and (18), we obtain the equality
(19)βμ+γ+ηS*+∫0+∞1-σθv*(θ)dθ=βμ+γ+ημ+βI*·1-pΛ∫0+∞αθv0θexp⁡-∫0θ1-σθβI*dτhh+∫0+∞αθv0θexp⁡-∫0θ1-σθβI*dτdθ+βμ+γ+η·∫0+∞1-σθv0θexp⁡-∫0θ1-σθβI*dτdθ=1.It is easy to see that function(20)HI*=βμ+γ+η·1μ+βI*-βI*∫0θ1-σθdτd·(1-p)Λ∫0+∞αθv0θexp⁡-βI*∫0θ1-σθdτdθ+∫0+∞αθv0θ·exp⁡-βI*∫0θ1-σθdτdθ+∫0+∞1-σθv0θ·exp⁡-βI*∫0θ1-σθdτdθ
is continuous and monotonically decreases with the growth of I* and that equalities
(21)H0=Rv,lim⁡I*→+∞HI*=0hold.
Hence, for all Rv>1 there exists I*∈(0,∞) such that equality H(I*)=1 holds. Corresponding S* and v*(θ) are defined by (17) and (18), respectively.Thus, we just proved the following Theorem.Theorem 1.
System (3) always has the infection-free steady state E0=(S0,I0,v0(θ)). If Rv>1, then system also has a unique endemic steady state E*=(S*,I*,v*(θ)).


## 4. Global Stability of the Steady States

In this section, our goal is to obtain the global stability results of system (3). We address global asymptotic stability of system (3) by constructing appropriate Volterra-type Lyapunov functions (1). We have the following result.Theorem 2.
(i) WhenRv≤1, then the infection-free steady state E0 is globally asymptotically stable and is the only steady state. (ii) When Rv>1, then the endemic steady state E* is globally asymptotically stable.Proof.
(i) To prove the global asymptotic stability of the infection-free steady state, we consider the following Lyapunov functionL(t)=L1(t)+W(t), where
(22)L1(t)=S0GS(t)S0+It,W(t)=∫0+∞v0(θ)Gv(θ,t)v0(θ)dθ.
Using μS0=(1-p)Λ+∫0+∞α(θ)v0(θ)dθ, we have
(23)dL1dt=(1-p)Λ2-S(t)S0-S0S(t)+∫0+∞α(θ)v0(θ)-S(t)S0+1+ln⁡S(t)S0dθ+∫0+∞α(θ)v0(θ)v(θ,t)v0(θ)-ln⁡S(t)S0-S0v(θ,t)S(t)v0(θ)dθ+βI(t)∫0+∞v0θ1-σθv(θ,t)v0(θ)-1dθ-(μ+γ+η)1-RvI(t).
We have
(24)dWdt=∫0+∞1-v0(θ)v(θ,t)∂v(θ,t)∂tdθ=-∫0+∞1-v0θvθ,thhhhh·∂vθ,t∂θ+μ+αθvθ,thhhhhhhh∂vθ,t∂θ+β1-σ(θ)v(θ,t)I(t)dθ.
Note that
(25)∂∂θGv(θ,t)v0(θ)=1-v0(θ)v(θ,t)∂θv(θ,t)v0(θ)-v(θ,t)·Dθ[v0(θ)]v0(θ)v0(θ),
where ∂θv(θ,t) denotes ∂v(θ,t)/∂θ, and Dθ[v0(θ)] denotes dv0(θ)/dθ.
We knowDθ[v0(θ)]=-(μ+α(θ))v0(θ).
Hence,(26)v0(θ)∂∂θGv(θ,t)v0(θ)=1-v0θvθ,t∂θvθ,t+μ+αθvθ,t.
We have
(27)dWdt=-∫0+∞v0(θ)∂∂θGv(θ,t)v0(θ)dθ-βI(t)∫0+∞v0θ1-σθv(θ,t)v0(θ)-1dθ.
Hence, using integration by parts,
(28)∫0+∞v0(θ)∂∂θGv(θ,t)v0(θ)dθ=v0θGvθ,tv0θθ=0θ=+∞-∫0+∞Dθv0θGvθ,tv0θdθ=v0(θ)Gv(θ,t)v0(θ)θ=+∞-v00Gv0,tv00-∫0+∞Dθ[v0(θ)]Gv(θ,t)v0(θ)dθ
and using the following equalities v0(0)=pΛ, v(0,t)=pΛ, and Dθ[v0(θ)]=-(μ+α(θ))v0(θ), we have
(29)∫0+∞v0(θ)∂∂θGv(θ,t)v0(θ)dθ=v0θGvθ,tv0θθ=+∞+∫0+∞μ+αθv0θGvθ,tv0θdθ.
Substituting the expression (29) in (27)(30)dWdt=-v0(θ)Gv(θ,t)v0(θ)θ=+∞-∫0+∞μv0(θ)Gv(θ,t)v0(θ)dθ-∫0+∞α(θ)v0(θ)v(θ,t)v0(θ)-1-ln⁡v(θ,t)v0(θ)dθ-βI(t)∫0+∞v0θ1-σθv(θ,t)v0(θ)-1dθ=-v0(θ)Gv(θ,t)v0(θ)θ=+∞-∫0+∞μv0(θ)Gv(θ,t)v0(θ)dθ-βIt∫0+∞v0θ1-σθvθ,tv0θ-1dθ+∫0+∞αθv0θhhhhh·-vθ,tv0θ+1+ln⁡S0vθ,tStv0θ+ln⁡StS0dθ,
and combining (23) and (30), we have
(31)dLdt=dL1dt+dWdt=1-pΛ2-StS0-S0St-v0(θ)Gv(θ,t)v0(θ)θ=+∞-∫0+∞μv0(θ)Gv(θ,t)v0(θ)dθ-∫0+∞α(θ)v0(θ)GS(t)S0+GS0v(θ,t)S(t)v0(θ)dθ-(μ+γ+η)1-RvI(t).
Therefore, Rv≤1 ensures that dL/dt≤0 holds in the nonnegative region of the phase space. The strict equality holds only if S(t)=S0 and v(θ,t)=v0(θ) holds simultaneously with either Rv=1 or I(t)=0. It is easy to verify that infection-free steady state E0 is the only invariant set of the system in this set, and hence by LaSalle's invariance principle for semiflows (Corollary 2.3 in [15]), Rv≤1 is sufficient to ensure that this steady state is globally asymptotically stable.
(ii) To prove the global asymptotic stability of the endemic steady state, we consider the following Volterra-type Lyapunov functionU(t)=U1(t)+J(t), where
(32)U1(t)=S*GS(t)S*+I*GI(t)I*,J(t)=∫0+∞v*(θ)Gv(θ,t)v*(θ)dθ.
Using βS*+β∫0+∞(1-σ(θ))v*(θ)dθ=(μ+γ+η) and μS*=(1-p)Λ-βS*I*+∫0+∞α(θ)v*(θ)dθ, we have
(33)dU1dt=(1-p)Λ2-S(t)S*-S*S(t)+∫0+∞α(θ)v*(θ)-S(t)S*+1+ln⁡S(t)S*dθ+∫0+∞αθv*θhhhhhh·v(θ,t)v*(θ)-ln⁡S(t)S*-S*v(θ,t)S(t)v*(θ)dθ+βI*∫0+∞v*θ1-σθhhhhhhhh·v(θ,t)I(t)v*(θ)I*-v(θ,t)v*(θ)-I(t)I*+1dθ.
Similarly, we have(34)dJdt=∫0+∞1-v*(θ)v(θ,t)∂v(θ,t)∂tdθ=-∫0+∞v*(θ)v(θ,t)v*(θ)-1hhhhh·∂θvθ,tvθ,t+μ+αθhhhhhhhh∂θvθ,tvθ,t+β1-σθItdθ=-∫0+∞v*(θ)v(θ,t)v*(θ)-1hhhhh·∂θv(θ,t)v(θ,t)+μ+α(θ)+β1-σ(θ)I*dθ+βI*∫0+∞v*(θ)1-σ(θ)hhhhhhhh·v(θ,t)v*(θ)-1-v(θ,t)I(t)v*(θ)I*+I(t)I*dθ.
Note that
(35)∂∂θGv(θ,t)v*(θ)=1-v*θvθ,t∂θvθ,tv*θ-vθ,t·Dθv*θv*θv*θ.
We know Dθ[v*(θ)]=-(μ+α(θ))v*(θ)-(1-σ(θ))βv*(θ)I*.
Hence,(36)v*(θ)∂∂θGv(θ,t)v*(θ)=1-v*θvθ,t·∂θvθ,t+μ+αθ+1-σθβI*vθ,t.
We have
(37)dJdt=-∫0+∞v*(θ)∂∂θGv(θ,t)v*(θ)dθ+βI*∫0+∞v*(θ)1-σθhhhhhhhh·v(θ,t)v*(θ)-1-v(θ,t)I(t)v*(θ)I*+I(t)I*dθ.
Hence, using integration by parts,
(38)∫0+∞v*(θ)∂∂θGv(θ,t)v*(θ)dθ=v*θGvθ,tv*θθ=0θ=+∞-∫0+∞Dθ[v*(θ)]Gv(θ,t)v*(θ)dθ=v*(θ)Gv(θ,t)v*(θ)θ=+∞-v*(0)Gv(0,t)v*(0)-∫0+∞Dθ[v*(θ)]Gv(θ,t)v*(θ)dθ
and using the following equalities v*(0)=pΛ, v(0,t)=pΛ, and Dθ[v*(θ)]=-μv*(θ)-α(θ)v*(θ)-(1-σ(θ))βv*(θ)I*. we have
(39)∫0+∞v*(θ)∂∂θGv(θ,t)v*(θ)dθ=v*(θ)Gv(θ,t)v*(θ)θ=+∞+∫0+∞μ+αθ+1-σθβI*hhhhhhhh·v*θGvθ,tv*θdθ.
Substituting the expression (39) in (37)(40)dJdt=-v*(θ)Gv(θ,t)v*(θ)θ=+∞-∫0+∞μ+α(θ)v*(θ)Gv(θ,t)v*(θ)dθ+βI*∫0+∞v*(θ)1-σ(θ)hhhhhhhh·ln⁡vθ,tv*θ-vθ,tItv*θI*+ItI*dθ,
and combining (33) and (40), we have
(41)dUdt=dU1dt+dJdt=(1-p)Λ2-S(t)S*-S*S(t)-v*(θ)Gv(θ,t)v*(θ)θ=+∞-∫0+∞μ+(1-σ(θ))βI*v*(θ)Gv(θ,t)v*(θ)dθ-∫0+∞αθv*θhhhhh·GS(t)S*+GS*v(θ,t)S(t)v*(θ)dθ.
That is positive-definite function U(t) has nonpositive derivative dU/dt. Furthermore, the equality dU/dt=0 holds only if S(t)=S* and v(θ,t)=v*(θ), simultaneously. Thus the endemic steady state E* is the only positively invariant subset of the set with dU/dt=0, and hence by LaSalle's invariance principle for semiflows (Corollary 2.3 in [15]), the steady state E* is globally asymptotically stable when it exists.The following Corollary immediately follows from the proof.Corollary 3.
System (3) is uniformly persistent.


## 5. Concluding Remarks

The vaccination of newborn or neonatal immunization may be an effective preventive strategy against some pathogens. The ideal vaccine would be a single dose given at birth and providing immediate and long protection against multiple diseases. In this paper, we introduced and analytically studied the global stability properties of an SIR epidemiological model with vaccination, where it is assumed that the waning rate of vaccine-induced immunity is depending on the vaccine-age. In this model, we introduce the vaccination of newborns or neonatal immunization in the vaccine-age.One of the fundamental questions of mathematical epidemiology is to find the reproductive number, we derived an explicit formula for the effective reproductive number of epidemic system (3) using the condition that the infection-free steady state E0 is unstable if and only if Rv>1.We perform a nonlinear stability analysis, by means of the Lyapunov function techniques and LaSalle's Invariance Principle for semiflows. We show that the classical threshold condition for the effective reproductive number,Rv, holds: Rv>1; then the endemic steady E* is globally asymptotically stable, whereas if Rv≤1, then the infection-free steady state E0 is globally asymptotically stable.It is noteworthy, that we only impose integrability condition on the wane rate function of vaccineα(θ) in the model, and hence, global stability results derived are robust and do not depend on specific forms of this function. Some forms of the function α(θ) have been proposed in [8, 9]. The following function α(θ) is reported in [9]:
(42)α(θ)=0,if0<θ≤θmin⁡,α0,ifθmin⁡<θ<θmax⁡,
where θmin⁡ and θmax⁡ are minimum and maximum vaccine-ages, respectively, and α0 is positive constant. We generalized the specific form given in [8]
(43)α(θ)=0,if0<θ≤θmin⁡,κ0θ-θmin⁡2exp⁡-κ1θ-θmin⁡,ifθmin⁡<θ≤θmed,α0,ifθmed<θ<θmax⁡,
where κ0 and κ1 are positive constants and θmed is intermedium vaccine-age.We notice that the effective reproductive number can be written as(44)Rv=R01-p1-∫0+∞αθ+μ1-σθK0(θ)dθhhhh·1-∫0+∞αθ+μ1-σθK0(θ)dθ,
where R0=βΛ/(μ(μ+γ+η)) is the basic reproductive number of the disease. From this, a measure of vaccine impact (see [16]), ϕ, can be derived as
(45)ϕ=R0-RvR0=p1-∫0+∞αθ+μ1-σθK0(θ)dθ.
It is easy to check that
(46)∂ϕ∂p=1-∫0+∞αθ+μ1-σθK0(θ)dθ>0.
Thus, the vaccine impact ϕ is an increasing function with respect to the argument p.Finally, we note that the approach used in this paper can be applied to a wider variety of the age-dependent epidemic models with waning vaccine-induced immunity.



---

*Source: 341854-2015-02-01.xml*