![MathJax Logo](/templates/jsp/_style2/_tandf/pb2/images/math-jax.gif)
ABSTRACT
Based on game theory, we propose an age-structured model to investigate the imitation dynamics of vaccine uptake. We first obtain the existence and local stability of equilibria. We show that Hopf bifurcation can occur. We also establish the global stability of the boundary equilibria and persistence of the disease. The theoretical results are supported by numerical simulations.
1. Introduction
Infectious diseases are one of the main enemies threatening human's health. Every year millions of people die as a result of various infectious diseases. Mathematical modelling plays a key role in studying and discovering the transmission mechanisms of important diseases. There exists ample theoretical work investigating infectious diseases and very valuable results on the control of epidemic spread have been obtained. In some of the studies [Citation5, Citation16, Citation18, Citation22, Citation29, Citation34] on childhood diseases, one eminent feature is that the transmission rates have salient differences since children have different immunities.
Vaccination has been widely considered as one of the effective methods to reducing the morbidity and mortality from infectious diseases. There are two key vaccination policies: voluntary vaccination and mandatory vaccination. There has been a vigorous debate about voluntary vaccination policy failing to protect population adequately. A rational vaccine decision-making is determined by various factors, such as perceived infection risk, potential side effects from vaccination, and vaccinating behaviours of other individuals. Because of the declining familiarity with the interplay between perceived infection risk and potential side effects from vaccination, parents have various reasons for avoiding the potential side effects for their children, relying instead on enough other children being vaccinated to provide herd immunity. Therefore, rational decision may lead to a reduced number of vaccination intakes. This is a free-riding dilemma between individuals and the public good. Game theory builds a bridge connecting the epidemic models with individual behaviours.
Vaccine decision-making based on game theory has been extensively investigated (see, e.g. [Citation1–3, Citation8, Citation10–15, Citation17, Citation28, Citation35, Citation35–38, Citation40]). In these works, it is usually assumed that individuals have full information for perceived infection risk and potential vaccine health risk. Under this assumption, rational vaccine decision-making will get the highest personal utilities, that is, there exists a Nash equilibrium where no individuals could be better off by unilaterally changing to a different strategy [Citation3, Citation15]. In [Citation35], Xia and Liu employed a computational approach to characterize the impact of social influence on individuals' vaccination decision-making while in [Citation36] they also investigated the impact of the two factors, information of the disease prevalence and the perceived vaccination risk, and fading coefficient of awareness spread. Recently, Xu and Cresmman [Citation38] built a nonlinear epidemic model with the smoothed best response by game-theory based decisions on vaccination. Their results showed that if there is a perceived cost for vaccination, the smoothed best response is very effective in controlling the disease spread, but if vaccination is free, the best response is a good control. Zhang et al. [Citation40] investigated the ‘double-edged sword’ effect on public health conditions for rational decision-making. Shim et al. used game dynamic models to gain insight into the decision-making between vaccine skeptics and vaccine believers [Citation31] and the decision-making with regard to antiviral intervention during an influenza pandemic [Citation30], respectively.
The above-mentioned references focus on vaccine uptake based on the perfect information the parents mastered on perceived infection risk and potential vaccine side effects. In fact, the vaccine uptake behaviours evolve with respect to time. Individuals' decision on vaccination or not is conducted by imitating others who appear to have adopted more successful strategies. The process is called ‘imitation dynamics’ and it has been studied by some researchers. Bauch [Citation2] investigated parents' vaccinating decisions for their children with the assumption that the susceptibles behave strategically in accordance with imitation dynamics and studied the dependence of epidemic prevalence and coverage of vaccination on these strategic decisions. Fu et al. [Citation17] proposed an agent-based model on social network with game theory to shed light on how imitation of peers shapes individual vaccination choices. d'Onofrio et al. [Citation12] assumed that the perceived risk of vaccination is a function with respect to the incidence and studied an susceptible-infected-removed (SIR) transmission model with dynamic vaccine demand based on an imitation mechanism. A common assumption in these existing vaccination models on imitation dynamics is that individuals are homogeneously mixed or heterogeneously mixed on social networks [Citation9, Citation11, Citation37]. This cannot capture the feature of different transmission rates for childhood diseases which we mentioned earlier. In order to address the imitated behaviours on vaccination decision-making, we develop an age-structured epidemic model based on the game theory. We assume that the transmission incidence rate varies according to the infection age and individuals act their behaviours with unperfect information about perceived infection risk and potential vaccine side effects. The payoff of the perceived infection risk is also relative with the infection age. Parents imitate the vaccination decision-making of their neighbours and then adopt a vaccination policy.
The rest of the paper is organized as follows. Based on the game theory, an SIR epidemic model with infection age is introduced in Section 2. We study the existence of equilibria and their local stability in Section 3. A local Hopf bifurcation may occur. Section 4 is devoted to the global stability of boundary equilibria based on the fluctuation lemma and the results for the case without vaccination. We also consider the persistence of disease in Section 4. In Section 5, numerical simulations are provided to demonstrate the theoretical results. The paper concludes with some discussions.
2. The model formulation
In [Citation2], Bauch proposed an epidemic model to predict vaccinating behaviour based on game theory. The total population is divided into three classes: the Susceptible , the Recovered
, and the Infected
at time t. It is assumed that all newborns are susceptible; All the groups decrease due to natural death; Susceptible individuals are infected by infectious individuals and enter the infected group; The infectious individuals exit the group due to natural death and recovery; The recovered group increases due to recovered infectious individuals and the vaccinated newborns and decreases due to natural death; The recovered individuals acquire permanent immunity from the disease or vaccination and never leave this group. The incidence has a bilinear form. In order to understand a strategic interaction between individuals when they are deciding whether or not to vaccinate, a replicator equation for an evolution population game is introduced, and the payoff function varies according to time. If x denotes the proportion of children that are vaccinated at birth, the equation for x formally reads as
where the vaccination cost is
and the non-vaccination cost is
. Here
denotes the perceived probability of suffering significant morbidity upon infection, m quantifies the sensitivity of vaccinating behaviour to changes in prevalence, k denotes the imitation rate. Under these assumptions, Bauch [Citation2] obtained the following model based on game theory,
(1)
(1) where Λ is the input rate (note that we modified the input rate to be Λ instead of μ as in [Citation2]), β is the transmission rate, μ is the natural death rate, γ is the recovery rate.
It is well-known that most childhood diseases have incubation period (the length of time from when a child is exposed to the illness to when the first symptoms appear in that child). For example, the ranges for the incubation period of measles is 7–21 days [Citation7], of mumps is 16–25 days [Citation26], of pertussis is 4–21 days [Citation27]. During the incubation period, whether the exposed child is infectious or not depends on the disease. During the contagious period (the amount of time during which a sick child can give the disease to others), the infectivity of infectious individuals may vary. Moreover, inadequate vaccination coverage may increase congenital rubella syndrome numbers by increasing the average age at infection[Citation25]. To describe these phenomena, we introduce the infection age. Let denote the density of infected individuals of infection age a at time t. Similar to the bilinear form, we assume that the incidence rate depends on the infectious age a in the following way,
where
is the transmission rate with infection age a. In epidemiology,
is called the force of infection at time t.
For convenience, the perceived vaccination cost is ; while the potential risk cost for non-vaccinator mainly depends on the epidemic transmission, so its cost is
, or the payoff benefit is
. Then the payoff gain for a vaccinator compared to a non-vaccinator is
. This is formally analogous with [Citation38]. Therefore, similar to [Citation2], we obtain the replicator equation for the portion x of vaccinated children at birth as follows,
Based on the above assumptions, the model with game theory is given by the following system of ordinary and partial differential equations,
(2)
(2) with the boundary condition
and initial conditions
where
is the set of all integrable functions from
into
.
is the recovered rate with infection age a, δ is the imitation rate which has the familiar biological meaning with k in Equation (Equation1
(1)
(1) ).
Due to the biological background, Λ, μ, δ and r are all positive, and belongs to
with β being not identically zero, where
denotes the set of all bounded and uniformly continuous functions from
into
.
Note that the third equation of (Equation2(2)
(2) ) is decoupled from the others. Then we can ignore it and only consider the following system,
(3)
(3) with the boundary condition
and initial conditions
Due to Iannelli [Citation21] and Magal [Citation23] , if
satisfies the coupling equation
then Equation (Equation3
(3)
(3) ) is well-posed. In the sequel, whenever we say a solution of
Equation (Equation3
(3)
(3) ), the above assumptions on the boundary and initial conditions are satisfied.
Proposition 2.1.
Let be a solution of Equation (Equation3
(3)
(3) ) with the maximal interval of existence
(
is allowed to be
) . Then
and
for
and
.
Proof.
Firstly, from the third equation of (Equation3(3)
(3) ), we have
(4)
(4) Since
, it follows that
for
.
Secondly, by the first equation of (Equation3(3)
(3) ), we see that
for
. Clearly,
for
as
and
for
.
Now, for any , we claim that
for
and
(the idea of the proof is borrowed from Browne and Plyugin [Citation6]). By way of contradiction, assume that there exists
and
such that
. Integrating the second equation of (Equation3
(3)
(3) ) with the boundary condition yields that, for
and
, we have
(5)
(5) where
for
and
for
. Since
, we have
and
. Let
. Then
. By the continuity of
and the definition of
, we know that
for
and
. Note that, with the help of Equation (Equation5
(5)
(5) ),
for
. Since B is continuous on
and
, there exists an
such that
for
. Let
and
. Denote
. Let
and
. Then Y is a complete distance space with the supremum distance and S is a closed subset of Y . Define L on S by
Obviously,
. Note that
for
and
. Then, for
and
, we have
and this shows that L is a mapping from S into itself. Moreover, for
,
, we have
that is, L is a contraction mapping on S. By the Banach Fixed Point Theorem, there exists a unique
such that
for
. By the uniqueness of solutions, we have
for
. This contradicts the definition of
and hence we have proved the claim. Since ξ is arbitrary, we see that
for
and
. This completes the proof.
Proposition 2.1 implies that for
. Setting
for
. With the help of Proposition 2.1, we can obtain
(6)
(6) for
, which implies that
. Therefore, the maximal interval of existence for every solution of Equation (Equation3
(3)
(3) ) is
. Let
be the solution semiflow associated with Equation (Equation3
(3)
(3) ), that is,
Note that, by Equation (Equation6(6)
(6) ),
. Moreover, if
for some
then
for all
. Define
Then we have proved the following result.
Proposition 2.2.
Γ is an attractive and positively invariant set for Equation (Equation3(3)
(3) ).
To end this section, we mention that, by Equation (Equation4(4)
(4) ), both sets
and
are also positively invariant subsets of Equation (Equation3
(3)
(3) ). It is clear that for solutions in
we have
for
. On the other invariant set
, Equation (Equation3
(3)
(3) ) reduces to the case without vaccination, that is,
(7)
(7) This model has been studied by Magal et al. [Citation24] and their main results are as follows.
Theorem 2.3
[Citation24, Theorems 1.2 and 1.3]
If
then the disease free equilibrium
is globally asymptotically stable for the semiflow generated by Equation (Equation7
(7)
(7) ).
Let
and
Assume
. Then every solution of Equation (Equation7
(7)
(7) ) with the initial value in
respectively, in
stays in
respectively, in
. Moreover, each solution with initial value in
converges to
. Furthermore, every solution with an initial value in
converges to the endemic equilibrium
, which is locally asymptotically stable.
The attractivity of the endemic equilibrium in Theorem 2.3(ii) is established by the approach of Lyapunov functionals. We also mention that Theorem 2.3 will play an important role in dealing with the attractivity of the boundary equilibria of Equation (Equation3(3)
(3) ).
3. Existence of equilibria and their local stability
Recall that for
. Define
and
Note that
is the survival probability at infection age a and hence K is the mean value of the infectivity of an exposed child. As the total population size is
,
is the average number of cases that one typical infected individual can produce during the infectious period. In epidemiology,
is called the basic reproduction umber. As we will see soon, the structure of equilibria of Equation (Equation3
(3)
(3) ) depends on the values of
.
is an equilibrium of Equation (Equation3(3)
(3) ) if it satisfies
(8)
(8) By direct computation, we easily see that Equation (Equation8
(8)
(8) ) can have at most four equilibria,
,
,
,
, where
The existence of equilibria of Equation (Equation3
(3)
(3) ) is summarized below.
Theorem 3.1.
If
then Equation (Equation3
(3)
(3) ) has two equilibria
and
.
If
then Equation (Equation3
(3)
(3) ) has three equilibria
and
.
If
then Equation (Equation3
(3)
(3) ) has four equilibria
and
.
From the biological viewpoint, is the disease free and pure vaccinator strategy equilibrium,
is the disease free and non-vaccinator strategy equilibrium,
is the endemic and non-vaccinator strategy equilibrium, and
is the endemic and vaccinator strategy equilibrium.
Let be an equilibrium of Equation (Equation3
(3)
(3) ). Then the linearized system at
is
Letting
,
, and
leads to the characteristic equation at
, which is
where
. Then we can get the local stability of
,
, and
as follows.
Theorem 3.2.
The disease free and pure vaccinator equilibrium
is always unstable.
The disease free and non-vaccinator equilibrium
is locally asymptotically stable if
and it is unstable if
.
The endemic and non-vaccinator equilibrium
is locally asymptotically stable if
and it is unstable if
.
Proof.
(i) The characteristic equation for is
which has a positive root
. This implies that
is unstable.
(ii) The characteristic equation at is
Besides the two negative roots
and
, the other roots are given by the equation
(9)
(9) Since
and
, it follows that Equation (Equation9
(9)
(9) ) has a positive root if
. If
, we claim that all roots of Equation (Equation9
(9)
(9) ) have negative real parts. Otherwise, if
is a root of Equation (Equation9
(9)
(9) ) with non-negative real part, then
a contradiction. This proves the claim. Therefore,
is stable if
and it is unstable if
.
(iii) The characteristic equation at is
First, the root of
is
, which is positive (respectively, negative) if
(respectively,
). Second, we claim that
has no root with non-negative real part. By contradiction, suppose it has a root
with non-negative real part. Then
It follows that
which is a contradiction since
and
and
. This proves the claim. To summarize, we have proved (iii).
Theorem 3.2 tells us that the disease free and pure vaccinator equilibrium is always unstable. This means that if the level of the vaccinated children is very high then the unvaccinated children have no incentives to vaccinate since the herd immunity can protect them and they do not care about the potential risk from vaccination.
The analysis for the stability of is not so easy. In fact, the characteristic equation at
is
(10)
(10) where
. Obviously, 0 is not a root to Equation (Equation10
(10)
(10) ) when
. Then Equation (Equation10
(10)
(10) ) is equivalent to
It follows easily that Equation (Equation10
(10)
(10) ) has no non-negative real roots. However, it is hard to see whether Equation (Equation10
(10)
(10) ) has roots with non-negative real part or not. Actually, it may have roots with non-negative real parts as we will see soon.
In order to have a clear picture of the dynamics of Equation (Equation3(3)
(3) ), we make the following further assumption.
Assumption 3.3.
Assume that the transmission rate and the transfer rate are
Assumption 3.3 means that when the infection age is larger than τ the infectivity of all infectious individuals is the same and constant. We call τ the latent period. Assumption 3.3 also implies that the self-healing only occurs when the infected become infectious. If we add a constant recovered rate for infection age then the analysis is quite similar. However, if we assume that the recovered rates differ at an infection age other than the latent period τ then the analysis will be very complicated. In the following we analyse the role played by τ.
Under Assumption 3.3, . It follows that
Let
and
. Then
if and only if
and
;
if and only if
and
. Moreover,
and Equation (Equation10
(10)
(10) ) becomes
(11)
(11) where
.
Assume . Then Equation (Equation11
(11)
(11) ) reduces to
Note that
. Due to the Routh–Hurwitz criterion,
is locally asymptotically stable if
and unstable otherwise. By Theorem 3.2, we know that
,
, and
all are unstable if
. Therefore, we make one more assumption.
Assumption 3.4.
that is ,
which implies that the total non-vaccinators imitating vaccinators sample rate is small enough.
Next, we consider the possibility of the stability switching for system (Equation3(3)
(3) ). Note that
Equation (Equation11
(11)
(11) ) is a special case of the transcendental equation considered by Beretta and Kuang [Citation4], where they provided practical guidelines that combine graphical information with analytical work to effectively study local stability. The theory was also illustrated by them with first-order and second-order characteristic equations. Here is an application to a third-order characteristic equation. Let
Then one can easily verify assumptions (i)–(iv) in Beretta and Kuang [Citation4, p. 1145]. Under Assumptions 3.3 and 3.4, a stability change at
can only happen when there are eigenvalues crossing the imaginary axis from left to right. Therefore, we seek a pair of purely imaginary eigenvalues
with
for some
since 0 is not an eigenvalue. Substituting
into Equation (Equation11
(11)
(11) ) (we only need to consider one of the roots) and separating the real and imaginary parts yields
(12)
(12) Then
(13)
(13) We square both sides of the equations in (Equation12
(12)
(12) ) and add them up to obtain
or
(14)
(14) with
. According to Descartes' Rule of Signs, Equation (Equation14
(14)
(14) ) has at most one positive root. On the other hand, it is easy to see that Equation (Equation14
(14)
(14) ) has a positive root. Therefore, Equation (Equation14
(14)
(14) ) has a unique positive root, denoted by
with
. As
is a simple root of Equation (Equation14
(14)
(14) ), we have
Let
be the solution of
(obtained by substituting
into the right-hand sides of Equation (Equation13
(13)
(13) )). For any
, define
by
The following result is deduced from Theorem 2.2 of Beretta and Kuang [Citation4].
Proposition 3.5.
Suppose that there exist and
such that
. Then Equation (Equation11
(11)
(11) ) has a pair of simple conjugate pure imaginary roots
at
. Moreover, the pair of simple conjugate imaginary roots crosses the imaginary axis from left to right if
and crosses the imaginary axis from right to left if
where
Denote
The set
is finite though it is a little bit long and difficult to show here. If
, let
for some
with
. It is easy to deduce
as Equation (Equation3
(3)
(3) ) is stable when τ is small enough. The following result follows from Proposition 3.5 and the above discussion.
Theorem 3.6.
Suppose Assumptions 3.3 and 3.4 hold and . Then the following statements are true.
If
then the endemic and vaccinator equilibrium
is locally asymptotically stable for all
.
If
with
then the endemic and vaccinator equilibrium
is locally asymptotically stable for all
. Moreover, if
then a Hopf bifurcation occurs at
.
Under the condition in Theorem 3.6, there may be a stability switch for . A stability switch may occur if there exists
such that
.
4. The attractivity of boundary equilibria and disease persistence
In this section, we first study the attractivity of the boundary equilibria and
by applying Theorem 2.3 and the comparison principle.
To establish the attractivity of , we need the fluctuation lemma. For a function
, we denote
Lemma 4.1
(Fluctuation Lemma [Citation20])
Let be a bounded and continuously differentiable function. Then there exist sequences
and
such that
and
as
.
Theorem 4.2.
Suppose that . Then the equilibrium
attracts all solutions of Equation (Equation3
(3)
(3) ) with
.
Proof.
Let be any solution of Equation (Equation3
(3)
(3) ) with initial condition and boundary conditions from
. Then
for
and hence
satisfies
Let
be the solution of the auxiliary system
(15)
(15) with
and
. With the help of the comparison principle, it is easy to see that
and
for
. Since
, it follows from Theorem 2.3 that
converges to
. Then
which implies that
Next, we show . Choose
. Since
,
implies that
. Then there exists a
such that
for
. By the third equation of (Equation3
(3)
(3) ), we get
Remembering
for
and
, we easily see that
.
Finally, we show . By Lemma 4.1, choose
such that
,
, and
as
. Taking limit in
produces
or
. This, combined with
, yields
.
To summarize, we have shown that converges to
. This completes the proof.
For the attractivity of , we define
where
. With the integrated semigroup approach, one can show (similarly to Magal et al. [Citation24], for example) that
is a positively invariant subset for Equation (Equation3
(3)
(3) ).
Theorem 4.3.
Suppose that . Then the equilibrium
attracts all solutions of Equation (Equation3
(3)
(3) ) with
.
Proof.
Let be a solution of Equation (Equation3
(3)
(3) ) with
.
We first show . From the proof of Theorem 4.2, we know that
(16)
(16) for
. Since
, it follows from Theorem 2.3 that
(17)
(17) Then we can obtain
Choose
such that
, which is possible since
. So there exists
such that
By the third equation of (Equation3
(3)
(3) ), we have
This, together with
and
for
, gives
.
Next, we show that and
in
as
. For any
, it follows from
that there exists
such that
Then, for
,
satisfies
Let
be the solution of the auxiliary system
with
and
. Then applying the comparison principle yields
(18)
(18) for
. Since
, applying Theorem 2.3 again, we see that
(19)
(19) Notice that, for
,
By Equations (Equation16
(16)
(16) )–(Equation19
(19)
(19) ), we have
and
By the arbitrariness of η, we have
. Therefore, we have proved that
and
in
as
. This completes the proof.
Finally, we study the disease persistence of Equation (Equation3(3)
(3) ). Define
, which is the same as
. Let
We distinguish two kinds of persistence.
The disease in Equation (Equation3
(3)
(3) ) is uniformly weakly ρ-persistent if there exists an
, independent of the initial conditions, such that if
then
The disease in Equation (Equation3
(3)
(3) ) is uniformly strongly ρ-persistent if there exists an
, independent of the initial conditions, such that if
then
Lemma 4.4.
Assume . Then Equation (Equation3
(3)
(3) ) is uniformly weakly ρ-persistent.
Proof.
Since , there exists an ε such that
and
(20)
(20) By way of contradiction, we assume that there exists
such that
Then there is a
such that
As before, it follows from
for
,
, and
that
. Therefore, there is a
such that
As a result, we have
which implies that
Then there exists
such that
With the help of Equation (Equation5
(5)
(5) ), we easily see that
(21)
(21) By replacing the initial condition with
, we can assume that Equation (Equation21
(21)
(21) ) holds for all
. Note that both
and
are bounded functions on
and hence their Laplace transforms exist at least on
. It follows from Equation (Equation21
(21)
(21) ) (with
) that
(22)
(22) where
denotes the Laplace transform of a function. As
is not identically zero on
, we know that
for
. It follows from Equation (Equation22
(22)
(22) ) that
for
. In particular,
, which contradicts with Equation (Equation20
(20)
(20) ). This completes the proof.
Next, we establish the uniform strong ρ-persistence. For this purpose, it is crucial to show Φ has a global compact attractor in . A global compact attractor
is a maximal compact invariant set in
such that for any open set that contains
, all solutions of
Equation (Equation3
(3)
(3) ) that start at zero from a bounded set, are contained in that open set, at least for sufficiently large time. To establish the existence of global attractors, we need the following three results.
Lemma 4.5
([Citation19, Theorem 3.4.6])
If
is asymptotically smooth, point dissipative and orbits of bounded sets are bounded, then there exists a global attractor.
A semiflow is called asymptotically smooth if each forward invariant bounded closed set is attracted by a non-empty compact set.
Lemma 4.6
([Citation19, Lemma 3.2.3])
For each suppose
has the property that
is complete continuous and there is a continuous function
such that
as
and
if
. Then
is asymptotically smooth.
The next result will be needed in the discussion of Φ being asymptotically smooth.
Lemma 4.7.
For any there exists
such that
(23)
(23)
Proof.
It is easy to see that for
. Now, let
and
. Then
or
(24)
(24) We now come to estimate
Note that, by Equation (Equation5
(5)
(5) ), we have
Therefore,
Since
and
for
, we can get
This, together with Equation (Equation24
(24)
(24) ) and the fact that β is uniformly continuous, immediately produces Equation (Equation23
(23)
(23) ).
By Lemma 4.4, we know that for
. So it induces a semiflow on
.
Lemma 4.8.
If then there exists a global attractor for the solution semiflow Φ of Equation (Equation3
(3)
(3) ) in
.
Proof.
With the help of Proposition 2.2 and Lemma 4.5, we only need to show that the restricted semiflow on is asymptotically smooth. This is done by applying Lemma 4.6 as follows.
For and
, let
where
(25)
(25) and
(26)
(26) Then
. Clearly, both
and
are non-negative. It follows from Equation (Equation26
(26)
(26) ) that
and hence
satisfies the assumption in Lemma 4.6.
Now, we show that is complete continuous, that is, the set
is precompact for any fixed
and any bounded set
. This is done by applying the Fréchet–Kolmogorov theorem [Citation39]. First, it follows easily from the definitions of
,
, and Γ that
is bounded and this verifies the first condition of the Fréchet–Kolmogrov theorem. Second, by Equation (Equation25
(25)
(25) ), the third condition of Fréchet-Kolmogorov theorem is satisfied. Finally, we verify the second condition of the Fréchet–Kolmogrov theorem. It suffices to show that
(27)
(27) If
then Equation (Equation27
(27)
(27) ) holds automatically since
by Equation (Equation25
(25)
(25) ). Without loss of generality, we assume that
. Since we concern with the limit as h tends to
, we only consider
. Then
as we know
for
and
. This, combined with Lemma 4.7, immediately gives Equation (Equation27
(27)
(27) ) and hence we have completed the proof.
Now, with the assistance of Lemmas 4.4, 4.8, and [Citation33, Theorem 2.3], we can obtain the following result.
Theorem 4.9.
If then Equation (Equation3
(3)
(3) ) is uniformly strongly ρ-persistent.
5. Numerical simulations
In this section, we always assume that Assumption 3.3 holds. We provide some simulations to illustrate the theoretical results obtained in the previous sections. Here for
.
First, set ,
,
,
,
,
, and
. Then
. By Theorem 4.2,
attracts all solutions with initial conditions in
. Figure (a) supports this with the initial condition
. In Figure (b), we plot the phase diagram in the SI-plane with five different initial conditions.
Figure 1. With ,
,
,
,
,
, and
, the equilibrium
attracts all solutions with initial conditions in
. (a) Time series of S, I, and x with five different initial conditions; (b) the phase diagram in the SI-plane with five different initial conditions.
![Figure 1. With Λ=0.4, r=0.9, δ=0.8, τ=1.25, β=0.001, μ=0.01, and γ=0.1, the equilibrium E2 attracts all solutions with initial conditions in {(S,i,x)∈Γ:x∈[0,1)}. (a) Time series of S, I, and x with five different initial conditions; (b) the phase diagram in the SI-plane with five different initial conditions.](/cms/asset/01055054-aab4-4cf2-8b03-6da26fb63ce7/tjbd_a_1099749_f0001_c.jpg)
Second, enlarge the transmission rate to while keeping the other parameters at the same values as in the first case. Then
, which is between 1 and
. It follows from Theorems 3.2 and 4.3 that the endemic and vaccinator equilibrium
is globally asymptotically stable. Figure indicates that Equation (Equation3
(3)
(3) ) evolves towards
.
Figure 2. With ,
,
,
,
,
, and
, the equilibrium
is locally asymptotically stable. (a) Times series of S (blue), I (green), and x (red) with five different initial conditions, (b) the corresponding phase diagram in the SI-plane.
![Figure 2. With Λ=0.4, r=0.9, δ=0.8, τ=1.25, β=0.02, μ=0.01, and γ=0.1, the equilibrium E3 is locally asymptotically stable. (a) Times series of S (blue), I (green), and x (red) with five different initial conditions, (b) the corresponding phase diagram in the SI-plane.](/cms/asset/993648e6-1b78-4aa3-8147-10d3ce00eda7/tjbd_a_1099749_f0002_c.jpg)
Suppose the cost for the vaccination is too small and let . From the third equation of (Equation3
(3)
(3) ) we know that x monotonously increases to 1. Then Equation (Equation3
(3)
(3) ) has a pure vaccinator equilibrium
, which is unstable by Theorem 3.2 (see Figure (a)). if the cost of the vaccination is free or low, then the coverage of the proportion vaccinated is towards high values. However, persons have no incentives to vaccinate and the pure vaccinator equilibrium is unstable.
Figure 3. Representative time series of S (blue), I (green), and x (red). (a) is unstable, here
,
,
,
,
,
, and
; (b)
is locally asymptotically stable with
,
,
,
,
,
, and
.
![Figure 3. Representative time series of S (blue), I (green), and x (red). (a) E1 is unstable, here Λ=0.4, r=0.000001, δ=0.5, τ=1.25, β=0.005, μ=1/50×365, and γ=0.1; (b) E4 is locally asymptotically stable with Λ=0.4, r=0.01, δ=0.5, τ=1.25, β=0.005, μ=1/50×365, and γ=0.1.](/cms/asset/70d97fef-a03f-41fb-bee2-71eb62f3f1d4/tjbd_a_1099749_f0003_c.jpg)
We rise the cost for vaccination to (which is between
(lower cost) and
(higher cost)) and keep the other parameters unchanged as in the second case. Note that
. Theorem 3.2 implies that the endemic and vaccinator equilibrium
is locally asymptotically stable, which is supported in Figure (b).
Theorem 3.2, Proposition 3.5, and Theorem 3.6 imply that δ, τ, and r play an important role in the evolution of (Equation3(3)
(3) ). Recall that δ represents the imitation rate describing the imitated behaviours and the vaccination behaviour for children is determined by the imitation rate. Let
,
,
,
,
, and
. For some larger imitation rate δ, system (Equation3
(3)
(3) ) exhibits limit cycles and hence it can be destabilized by the imitation rate δ. Figure (a) shows the bifurcation diagram by taking δ as a bifurcation parameter, which indicates that the amplitude of the oscillations increases as the imitation rate does; while Figure (b) shows the phase diagram in the Ix-plane with
,
, and
, and the initial condition
.
Figure 4. With ,
,
,
,
, and
, the imitation rate δ can destabilize (Equation3
(3)
(3) ). (a) The bifurcation diagram showing I with respect to the imitation rate δ with δ varying from 0 to 6; (b) representative phase diagram in the Ix-plane with δ being
(blue),
(green), and
(red), respectively, and the same initial condition
.
![Figure 4. With Λ=0.4, r=0.01, τ=1.25, β=0.05, μ=1/50×365, and γ=0.1, the imitation rate δ can destabilize (Equation3(3) dS(t)dt=Λ(1−x)−S(t)∫0∞β(a)i(t,a)da−μS(t),∂i(t,a)∂t+∂i(t,a)∂a=−(μ+γ(a))i(t,a),dx(t)dt=δx(1−x)∫0∞β(a)i(t,a)da−r(3) ). (a) The bifurcation diagram showing I with respect to the imitation rate δ with δ varying from 0 to 6; (b) representative phase diagram in the Ix-plane with δ being 0.5(blue), 1.5(green), and 2.5(red), respectively, and the same initial condition (S0,i0,x0)=(0.5,5e−0.4a,0.1).](/cms/asset/29007d8a-173d-43c2-99b5-66ef0b8e3585/tjbd_a_1099749_f0004_c.jpg)
The perceived vaccine risk r impacts the prevalence of the disease spread. If the cost of the vaccination is high, such as a vaccine scare or a strong side effect, the frequency of the vaccination will gradually decrease and the prevalence of the disease will increase. Figure (a) shows the infection equilibrium structure of (Equation3(3)
(3) ) in terms of the disease prevalence vs. r with
,
,
,
,
, and
. We can see that the disease prevalence increases when the cost of vaccination increases. Moreover, the equilibrium
changes stability with r in the interval
at
from stable to unstable.
Figure 5. With ,
,
,
, and
, (a) infection equilibrium structure of (Equation3
(3)
(3) ) with
; (b) bifurcation diagram with the latent period τ as bifurcation parameter and
.
![Figure 5. With Λ=0.4, δ=0.5, β=0.05, μ=1/50×365, and γ=0.1, (a) infection equilibrium structure of (Equation3(3) dS(t)dt=Λ(1−x)−S(t)∫0∞β(a)i(t,a)da−μS(t),∂i(t,a)∂t+∂i(t,a)∂a=−(μ+γ(a))i(t,a),dx(t)dt=δx(1−x)∫0∞β(a)i(t,a)da−r(3) ) with τ=1.25; (b) bifurcation diagram with the latent period τ as bifurcation parameter and r=0.01.](/cms/asset/5e90fc1f-a4e4-4f5c-b997-5d4ec3fa2018/tjbd_a_1099749_f0005_c.jpg)
The latent period τ is also a key parameter for controlling the disease. For a childhood disease, if some medicine can extend the latent period of the disease, the prevalence of the disease can be lowered. Figure (b) shows that τ can destabilize (Equation3(3)
(3) ) and oscillations occur when τ is from 100 to 500. Here
and the other parameters except τ have the same values as above for the impact of r. Also, the prevalence of the disease gradually increases with respect to τ.
6. Discussion
From the theoretical analysis in the above sections, the pure vaccinator equilibrium is always unstable. Parents have no incentive to vaccinate if the vaccination coverage is high. The disease free and non-vaccinator equilibrium
is globally asymptotically stable if and only if
. If
, then there is an endemic and non-vaccinator equilibrium
, which is locally asymptotically stable if
. The endemic and vaccinator equilibrium
is locally asymptotically stable if
,
and
. If
with
, then a Hopf bifurcation occurs from the equilibrium
, which implies that Equation (Equation3
(3)
(3) ) can be destabilized. From the biological perspective, the childhood diseases are controlled if the transmission rates are small enough; Otherwise, the diseases will spread and become endemic diseases in some region. Moreover, if vaccination payoff is high, parents have a heavy burden from vaccinating their children, then they would avoid vaccinating their children and this leads to the spread of childhood diseases. If the vaccination payoff is not large enough and the total sampling rate for non-vaccinators imitating vaccinators is small enough, even if some of the parents would vaccinate their children, childhood diseases still spread in these areas. If, however, the sampling rate for non-vaccinators imitating vaccinators is large enough, this will lead to imitating phenomena switching, i.e, first, non-vaccinators imitating vaccinators decreases, and then the vaccination side effects will increase; Hence, vaccinators imitating non-vaccinators increases, and then the diseases will spread. Then the phenomenon will continue to cycle.
The imitation rate δ is the main parameter that leads to destabilization of the system and to a Hopf bifurcation, which is proved in [Citation2, Citation12]. The amplitude of the oscillation increases with the increase of the imitated behaviours. As mentioned earlier in Section 5, individuals imitate others more readily as the amount of information for the significant side effects of vaccination increases, which in turn enhances the difficulties in disease control. The differences in infectivity, involved in childhood diseases, is incorporated through , the different transmission ability for the different values of the infection age a. For special cases, the different transmissibility can be described by the parameters τ and β. Suitable higher value of τ will produce a limit cycle.
When there is no cost for vaccination, that is, , the payoff benefit is larger than the cost of vaccination as long as
is not zero. From the third equation of (Equation3
(3)
(3) ), the frequency of vaccination monotonously increases. It follows from the first two equations combined with the third equation of (Equation3
(3)
(3) ) that system (Equation3
(3)
(3) ) evolves to
. On the other hand, if the cost of vaccination is not free, and the initial conditions satisfy
, the frequency of the vaccination monotonously increases. It follows from the first and second equation of (Equation3
(3)
(3) ) that infected individuals decrease until
. System (Equation3
(3)
(3) ) evolves to the endemic and non-vaccinator equilibrium
or a stable endemic and vaccinator equilibrium
which depends on whether the initial vaccinated proportion is equal or larger than 0. When
, it follows from the third equation of (Equation3
(3)
(3) ) that the proportion vaccinated x decreases. Prevalence of the disease increases until
. This implies that Equation (Equation3
(3)
(3) ) also evolves to the endemic and non-vaccinator equilibrium
or the endemic and vaccinator equilibrium
which depends on the initial vaccinated proportion
.
In summary, parents can make rational decisions in favour of disease control if they understand the interplay between the perceived vaccination risk, prevalence of the disease, and the variable transmission abilities of the children. Our investigation implies that high vaccination coverage can not guarantee the elimination of the disease for a voluntary vaccination policy. On the other hand, limited vaccination may be harmful for control of the disease spread. The amount of up-to-date information for the vaccine use has two opposite effects: knowledge about the disease prevalence encourages parents to take the vaccine for their children, on the other hand the potential side effects discourage them to take vaccine for their children. Rational decision-making individuals, depending on updated information about the perceived vaccine risk compared with the prevalence of the disease, decide whether or not to vaccinate their children. Based on the game theory, we conclude that it is important that parents make rational decision whether or not to vaccinate their children.
Disclosure statement
No potential conflict of interest was reported by the authors.
Additional information
Funding
References
- N.E. Basta, D.L. Chao, M.E. Halloran, and I.M. Longini Jr, Strategies for pandemic and seasonal influenza vaccination of schoolchildren in the United States, Am. J. Epidemiol. 170(6) (2009), pp. 679–686. doi: 10.1093/aje/kwp237
- C.T. Bauch, Imitation dynamics predict vaccinating behaviour, Proc R. Soc. B 272(1573) (2005), pp. 1669–1675. doi: 10.1098/rspb.2005.3153
- C.T. Bauch and D.J.D. Earn, Vaccination and the theory of games, Proc. Natl Acad. Sci. USA 101 (2004), pp. 13391–13394. doi: 10.1073/pnas.0403823101
- E. Beretta and Y. Kuang, Geometric stability switch criteria in delay differential systems with delay dependent parameters, SIAM J. Math. Anal. 33(5) (2002), pp. 1144–1165. doi: 10.1137/S0036141000376086
- S.R. Boëlle, Modelling the effects of population structure on childhood disease: The case of varicella, PLoS Comput. Biol. 7(7) (2011), p. e1002105. doi: 10.1371/journal.pcbi.1002105.
- C.J. Browne and S.S. Plyugin, Global analysis of age-structured within-host virus model, Discrete Contin. Dyn. Syst. Ser. B. 18(8) (2013), pp. 1999–2017. doi: 10.3934/dcdsb.2013.18.1999
- Centers for Disease Control and Prevention. CDC Health Information for International Travel 2014. Oxford University Press, New York, 2014.
- F.C. Coelho and C.T. Codeço, Dynamic modeling of vaccinating behavior as a function of individual beliefs, PLoS Comput. Biol. 5(7) (2009), p. e1000425. doi: 10.1371/journal.pcbi.1000425
- D.M. Cornforth, T.C. Reluga, E. Shim, C.T. Bauch, A.P. Galvani, and L.A. Meyers, Erratic flu vaccination emerges from short-sighted behaviour in contact networks, Plos Comput. Biol. 7(1) (2011), p. e1001062. doi: 10.1371/journal.pcbi.1001062
- R. Cressman, Evolutionary Dynamics and Extensive Form Games, MIT Press, Cambridge, MA, 2003.
- A. d'Onofrio and P. Manfredi, Vaccine demand driven by vaccine side effects: Dynamic implications for SIR diseases, J. Theoret. Biol. 264(2) (2010), pp. 237–252. doi: 10.1016/j.jtbi.2010.02.007
- A. d'Onofrio, P. Manfredi, and P. Poletti, The impact of vaccine side effects on the natural history of immunization programmes: An imitation-game approach, J. Theoret. Biol. 273(1) (2011), pp. 63–71. doi: 10.1016/j.jtbi.2010.12.029
- D.J.D. Earn, P. Rohani, B.M. Bolker, and B.T. Grenfell, A simple model for complex dynamical transitions in epidemics, Science 287 (2000), pp. 667–670. doi: 10.1126/science.287.5453.667
- E.P. Fenichel, C. Castillo-Chavez, M.G. Ceddia, G. Chowell, P.A. Parra, G.J. Hickling, G. Holloway, R. Horan, B. Morin, C. Perrings, M. Springborn, L. Velazquez, and C. Villalobos, Adaptive human behavior in epidemiological models, Proc. Natl. Acad. Sci. USA 108 (2011), pp. 6306–6311. doi: 10.1073/pnas.1011250108
- P.E.M. Fine and J.A. Clarkson, Individual versus public priorities in the determination of optimal vaccination policies, Am. J. Epidemiol. 124 (1986), pp. 1012–1020.
- B.F. Finkenstädt and B.T. Grenfell, Time series modelling of childhood diseases: A dynamical systems approach, J. Roy. Statist. Soc. Ser. C 49 (2000), pp. 187–205. doi: 10.1111/1467-9876.00187
- F. Fu, D.I. Rosenbloom, L. Wang, and M.A. Nowak, Imitation dynamics of vaccination behaviour on social networks, Proc. R. Soc. B 278 (2011), pp. 42–49. doi: 10.1098/rspb.2010.1107
- J.T. Griffin, N.M. Ferguson, and A.C. Ghani, Estimates of the changing age-burden of Plasmodium falciparum malaria disease in sub-Saharan Africa, Nat. Commun. 5:3136 (2014). doi:10.1038/ncomms4136.
- J.K. Hale, Asymptotic Behavior of Dissipative Systems, AMS, Providence, 1988.
- W.M. Hirsch, H. Hanisch, and J.-P . Gabriel, Differential equation models of some parasitic infections: Methods for the study of asymptotic behavior, Comm. Pure Appl. Math. 38 (1985), pp. 733–753. doi: 10.1002/cpa.3160380607
- M. Iannelli, Mathematical theory of age-structured population dynamics, in Applied Mathematics Monographs 7, comitato Nazionale per le Scienze Matematiche, Consiglio Nazionale delle Ricerche (C.N.R.), Giardini, Pisa, 1995.
- N.B. Kandala, T.P. Madungu, J.B.O. Emina, K.P.D. Nzita, and F.P. Cappuccio, Malnutrition among children under the age of five in the Democratic Republic of Congo (DRC): Does geographic location matter? BMC Public Health 11 (2011), pp. 261–276. doi: 10.1186/1471-2458-11-261
- P. Magal, Compact attractors for time periodic age-structured population models, Electron. J. Differ. Equ. 65 (2001), pp. 1–35.
- P. Magal, C.C. McCluskey, and G.F. Webb, Lyapunov functional and global asymptotic stability for an infection-age model, Appl. Anal. 89 (2010), pp. 1109–1140. doi: 10.1080/00036810903208122
- C.J.E. Metcalf, J. Lessler, P. Klepac, F. Cutts, and B.T. Grenfell, Impact of birth rate, seasonality and transmission rate on minimum levels of coverage needed for rubella vaccination, Epidemiol. Infect. 140(12) (2012), pp. 2290–2301. doi: 10.1017/S0950268812000131
- National Center for Immunization and Respiratory Diseases (NCIRD). Available at http://www.cdc.gov/mumps/clinical/qa-disease.html.
- Pertussis (Whooping Cough): Questions and Answers. Available at http://www.immunize.org/catg.d/p4212.pdf.
- T.C. Reluga and A.P. Galvani, A general approach for population games with application to vaccination, Math. Biosci. 230(2) (2011), pp. 67–78. doi: 10.1016/j.mbs.2011.01.003
- D. Schenzle, An age-structured model of pre- and post-vaccination measles transmission, Math. Med. Biol. 1 (1984), pp. 169–191. doi: 10.1093/imammb/1.2.169
- E. Shim, G.B. Chapman, and A.P. Galvani, Decision making with regard to antiviral intervention during an influenza pandemic, Med. Decis. Making 30(4) (2010), p. E64-81.
- E. Shim, J.J. Grefenstette, S.M. Albert, B.E. Cakouros, and D.S. Burke, A game dynamic model for vaccine skeptics and vaccine believers: Measles as an example, J. Theoret. Biol. 295 (2012), pp. 194–203. doi: 10.1016/j.jtbi.2011.11.005
- E. Shim, B. Kochin, and A.P. Galvani, Insights from epidemiological game theory into gender-specific vaccination against rubella, Math. Biosci. Eng. 6(4) (2009), pp. 893–854. doi: 10.3934/mbe.2009.6.839
- R.H. Thieme, Uniform persistence and permanence for non-autonomous semiflows in population biology, Math. Biosci. 166 (2000), pp. 173–201. doi: 10.1016/S0025-5564(00)00018-3
- D. Weycker, J. Edelsberg, M.E. Halloran, I.M. Longini, A. Nizam, V. Ciuryla, and G. Oster, Population-wide benefits of routine vaccination of children against influenza, Vaccine 23(10) (2005), pp. 1284–1293. doi: 10.1016/j.vaccine.2004.08.044
- S. Xia and J. Liu, A computational approach to characterizing the impact of social influence on individuals' vaccination decision making, PLoSONE. 8(4) (2013), p. e60373. doi: 10.1371/journal.pone.0060373.
- S. Xia and J. Liu, A belief-based model for characterizing the spread of awareness and its impacts on individuals' vaccination decisions, J. R. Soc. B 11 (2014), p. 20140013.
- S. Xia, J. Liu, and W. Cheung, Identifying the relative priorities of subpopulations for containing infectious disease spread, PLoS ONE. 8(6) (2013), p. e65271. doi: 10.1371/journal.pone.0065271.
- F. Xu and R. Cressman, Disease control through voluntary vaccination decisions based on the smoothed best response, Comput. Math. Methods Med. (2014), Art. ID 825734, 14pp. doi:10.1155/2014/825734.
- K. Yasida, Functional Analysis, 2nd ed., Springer-Verlag, Berlin-Heidelberg, 1968.
- H. Zhang, F. Fu, W. Zhang, and B. Wang, Rational behavior is a ‘double-edged sword’ when considering voluntary vaccination, Phys. A 391 (2012), pp. 4807–4815. doi: 10.1016/j.physa.2012.05.009