Li
Li
Li
Abstract. This article concerns the SIRC epidemiological model for influenza
A, which efficiently describes the mechanism of disease spreading, including
the susceptible (S), the infected (I) and the recovered (R), along with a cross-
immune class (C) that recovers after being inflected by different strains of the
same viral subtype. The dynamics of the model is completely determined by
the basic reproduction number R0 . If R0 ≤ 1, the disease-free equilibrium
of the SIRC model is globally asymptotically stable, which means influenza
A will die out. Otherwise, the SIRC model may have exactly one endemic
equilibrium which is globally asymptotically stable under certain parametric
conditions. Also, numerical simulations are given to support our analytical
results.
1. Introduction
Influenza [7], a member of the family Orthomyxoviridae, is a RNA virus which
can give rise to epidemic disease between mankind and animals. In general, in-
fluenza is primarily divided into A, B and C types and every type contains a wide
variety of subtypes according to hemagglutinin (HA) and neuraminidase (NA) dif-
ferences. From an epidemiologic standpoint, influenza A is the most common and
the most terrible virus among three types, and can result in the highest pathogenic-
ity because of the easiest way to generate variation. What’s worse, the virus
has brought about more than century pandemic influenza in the past years. The
pathogenicity of influenza B virus is the same as type A, but performed studies are
shown that type B virus does not contribute to century pandemic. Last, influenza C
virus just leads to unconspicuous or feeble respiratory infection and almost doesn’t
incur a pandemic either.
In this article, we pay attention to influenza A virus whose surface often has a
tiny variation, which is referred to as drift. That is, the virus camouflages itself by a
subtle change, and it thereby can elude identificaton by the human immune system.
As a result of the drift, all kinds of flu strains appear every year. Therefore, human
annually need to take a vaccine against influenza for prevention. On the other hand,
shift means influenza A takes place genetic mutations causing a new “subtype”. In
this case, it can lead to outbreaking a global pandemic influenza. For example,
Spanish outbroke H1N1 influenza virus in 1918-1919 about the death toll up to
20-40 million [24]. In 1957-1958, Asia influenza erupted due to the result of H2N2
virus. Furthermore, influenza A virus which threats the lives of even the healthiest
individuals brings serious economic loss for many countries. Some research data
indicate the direct economic loss up to 10-30 billion dollars in America, and the
potential economic losses at 10-15 billion dollars [23]. Consequently, devoting to
studying the spread of influenza A is extremely indispensable.
Mathematically modelling, as a bridge, plays a crucial role in understanding of
the spread of the epidemic. Researchers transform their focus from considering
threshold value of pathophoresis to prevent epidemic propagation by taking some
strategies for control [2]. The model of influenza is parallel to the traditional stan-
dard SIR model in which the crowd are segmented into three compartments: sus-
ceptible (S), infectious (I), recovered (R). A number of compartmental models have
been established based on this idea [26]. Nonetheless, some scholars indicate that
the traditional SIR model is not adequate to describe actual situation of influenza
spreading because every type can evolve different subtypes. Hence, Andreasen
et al. [1] initially considered the dynamics of system with multiple virus strains
consisting of partial cross-immunity [22], concluding that the system exists stable
equilibrium when the number of strains are less than or equal to three. However,
in the work of Lin et al. [13], they demonstrated that oscillations can be sustained
under a linear chain of three cocirculating influenza A strains. These researchers
only consider a special case where cross-protection is symmetric to analyze complex
system because the analysis and computation of multiple viral strains are much in-
tractable. Recently, Minayev and Ferguson [18] proposed a deterministic model
of multi-strain pathogens with symmetric equilibrium, self-organized strain struc-
tures, regular periodic and chaotic regimes, which is determined by cross-immune
response function. In addition to these properties, Kooi et al. [11] revealed bifurca-
tion analysis, Lyapunov exponent calculation as well as quantitative and qualitative
results by numerical simulations in three multi-strain compartment models. Nuño
et al. [21] developed a general multiple strains pathogens model with all kinds of
cross-immunity structures. They illustrated the weaker cross-immunity structures
are more likely to appear instability in the strain coexistence mode. For seasonal
influenza, there is another case where strong and weak cross-immunity can result
in coexistence with the following pandemic. On the contrary, the intermediate level
one may take the place of the seasonal subtype. Chung and Liu [6] extended the
result of Nuño et al. [20] and elucidated the local asymptotic stability of a two-
strain influenza model. Also, the stability may be lost when two strains are far
apart. All the above published works didn’t analyze global asymptotical stability
of models with cross-immune class. In this paper, we investigate a SIRC model
with cross-immune class proposed by Casagrandi et al. [5], which is represented by
EJDE-2017/121 DYNAMICS OF A SIRC EPIDEMIOLOGICAL MODEL 3
the literature due to the high dimensionality for the model. Our proof is based on
a theoretical approach developed in [14, 15].
From a mathematical viewpoint, it is interesting to notice that, in the absence
of cross-immunity (µ = 1), the two classes S and C are immunologically indistin-
guishable, since
d(S + C)
= γ[1 − (S + C)] − β(S + C)I + δR.
dt
Consequently, in the limit of µ → 1, the SIRC model reduces to the following
classical SIRS model
dS(t)
= γ(1 − S) − βSI + δR,
dt
dI(t)
= βSI − (γ + α)I, (1.2)
dt
dR(t)
= αI − (γ + δ)R.
dt
Thus, our results generalize some earlier results on the SIRS model. In this paper,
we shall prove the global stability of a unique endemic equilibrium of (1.2) when
the basic reproduction number R0 is greater than 1.
This article is organized as follows. The positively invariant set is shown in
section 2. Section 3 is devoted to the existence of equilibrium including a disease-
free equilibrium and a unique endemic equilibrium. The stability of the disease-
free equilibrium and the endemic equilibrium is presented in sections 4 and 6,
respectively. Section 5 is devoted to the persistence of the model and the work of
numerical simulation reveals in section 7. Finally, this paper ends up with a brief
conclusion in section 8.
Note that
dS(t) dI(t)
|S=0 = r + ηC > 0, |I=0 = 0,
dt dt
dR(t) dC(t)
|R=0 = (1 − µ)βCI + αI ≥ 0, |C=0 = δR ≥ 0.
dt dt
Then it follows from Lemma 2.1 that R4+ is invariant set.
3. Existence of equilibria
In this section, we consider the existence, type, and number of the equilib-
ria. Obviously, model (1.1) has two equilibria: one is disease-free equilibrium
E0 = (1, 0, 0, 0) which exists for all parameter values; and the other is the endemic
equilibrium E ∗ = (S ∗ , I ∗ , R∗ , C ∗ ), which is a positive solution of the following
equation
γ(1 − S ∗ ) − βS ∗ I ∗ + ηC ∗ = 0,
βS ∗ I ∗ + µβC ∗ I ∗ − (γ + α)I ∗ = 0,
(3.1)
(1 − µ)βC ∗ I ∗ + αI ∗ − (γ + δ)R∗ = 0,
δR∗ − βC ∗ I ∗ − (γ + η)C ∗ = 0.
It follows from the third and forth equations of (3.1) that
δαI ∗
C∗ = . (3.2)
(δµ + γ)βI ∗ + (γ + δ)(γ + η)
Next, substituting (3.2) in the second equation of (3.1), we obtain
γ+α µδα
S∗ = − ∗
· I ∗. (3.3)
β (δµ + γ)βI + (γ + δ)(γ + η)
Combining (3.2) with (3.3) and then substituting them in the first equation of (3.1),
we have
aI ∗2 + bI ∗ + d = 0, (3.4)
6 H. LI, S. GUO EJDE-2017/121
where
a = β 2 (γ + α)(δµ + γ) − β 2 µδα = β 2 (γδµ + γ 2 + αγ) > 0,
b = βγ(γ + α)(δµ + γ) + β(γ + α)(γ + δ)(γ + η)
− γβµδα − γβ 2 (δµ + δ) − βηδα,
d = γ(γ + α)(γ + δ)(γ + η) − γβ(γ + δ)(γ + η).
Therefore, we have the following result.
Theorem 3.1. If R0 ≤ 1, model (1.1) always has a disease-free equilibrium E0 =
(1, 0, 0, 0). If R0 > 1, the model has exactly one endemic equilibrium E ∗ =
(S ∗ , I ∗ , R∗ , C ∗ ).
Proof. It is not difficult to observe that model (1.1) has a disease-free equilib-
rium E0 (1, 0, 0, 0). Now, we consider the existence of the endemic equilibrium
E ∗ (S ∗ , I ∗ , R∗ , C ∗ ) when R0 > 1. It follows from R0 > 1 that d < 0 and hence
that (3.4) has exactly one positive solution I ∗ , which, together with (3.2), implies
that C ∗ is positive. Similarly, from the first and forth equations of (3.1), we see
that both S ∗ and R∗ are positive. Therefore, model (1.1) has exactly one endemic
equilibrium E ∗ when R0 > 1.
(0, 0, 0). By LaSalle’s invariant principle [12, 16], we conclude that model (4.1) is
global asymptotically stable at (0, 0, 0).
Remark 4.3. From the perspective of epidemiological significance, Theorem 4.2
demonstrates that the disease ultimately dies out regardless of the initial values of
model (1.1) when R0 ≤ 1.
8 H. LI, S. GUO EJDE-2017/121
Take µ = 1, then (1.1) can reduce to the classical SIRS model of the form (1.2).
Hence, we can deduce that the disease-free equilibrium G0 = (1, 0, 0) of (1.2) is
globally asymptotically stable when R0 ≤ 1. Namely, we have the following result.
Corollary 4.4. The disease-free equilibrium for model (1.2) is globally asymptoti-
cally stable when R0 ≤ 1.
In fact, it is straightforward to show that the disease-free equilibrium G0 =
(1, 0, 0) is locally asymptotically stable. We can also illustrate that G0 is globally
asymptotically stable by considering the following Lyapunov function
1 2γ γ
V2 = (S − 1 + I + R)2 + I + R2 .
2 β α
5. Persistence
In the section, we will establish the persistence theorem of disease when R0 > 1.
Persistence implies that the infected individuals will persist in the future. In this
paper, we first use the definition given by Butler and Waltman [4]. That is, model
(1.1) is said to be uniformly persistent if there exists a positive number b such that
n o
min lim inf (S(t)), lim inf (I(t)), lim inf (R(t)), lim inf (C(t)) = b (5.1)
t→∞ t→∞ t→∞ t→∞
for every trajectory with positive initial conditions.
For a region E, denote by ∂E and E̊ the boundary and the interior of E, re-
spectively. Denote by ∂F the restriction of the flow F to ∂E and note that ∂E
is, in general, not positively invariant. Let N be the maximal invariant set of ∂E.
Suppose N is a closed invariant set and there exists a cover {Nα }α∈A of N , where
A is a nonempty index set. Nα ⊂ ∂E, N ⊂ ∪α∈A Nα and Nα (α ∈ A) are pairwise
disjoint closed invariant sets. Furthermore, we propose the following hypotheses:
(H1) All Nα are isolated invariant sets of the flow F;
(H2) Nα (α ∈ A) is acyclic, that is, any finite subset of Nα (α ∈ A) does not
form a cycle;
(H3) Any compact subset of ∂E contains, at most, finitely many sets of Nα
(α ∈ A).
The following lemma plays an important role in analyzing the uniformly persistence.
Lemma 5.1 ([9]). Let E be a closed positively invariant subset of X on which a
continuous flow F is defined. Suppose there is a constant α > 0 such that F is
point dissipative on S[∂E, α] ∩ E̊ and the assumptions (H1)–(H3) hold. Then the
flow F is uniformly persistent if and only if W + (Nα ) ∩ S[∂E, α] ∩ E̊ = ∅ for all
α ∈ A, where W + (Nα ) = {y ∈ X| Λ+ (y) ⊂ Nα }.
Theorem 5.2. Model (1.1) is uniformly persistent in D̊ if and only if R0 > 1.
Proof. It is easy to prove the necessity by means of Theorems 4.1 and 4.2 because
the asymptotical stability of E0 excludes all kinds of persistence. Now, we prove
the sufficiency of this theorem by using Lemma 5.1. Choose X = R4 and E = D.
We only need to prove that model (1.1) satisfies all the conditions of Lemma 5.1.
Note that the maximal invariant set on the boundary ∂D only contains a point
E0 which is isolated. Then the assumptions (H1)–(H3) are satisfied. By Lemma
5.1, we observe that the uniform persistence of model (1.1) is equivalent to the
instability of the disease-free equilibrium. Thereby, the proof is complete.
EJDE-2017/121 DYNAMICS OF A SIRC EPIDEMIOLOGICAL MODEL 9
Remark 5.3. It follows from Theorem 5.2 that the uniform persistence of model
(1.1) in the bounded set D̊ is equivalent to the existence of a compact attractor
K ⊂ D̊.
where
∂f [2] −1
B = Pf P −1 + P P ,
∂x
and the matrix Pf is obtained by replacing each entry pij of P by its derivative in
the direction of f , pijf . The quantity µ(B) is Lozinskiĭ measure of B with respect
to a vector norm | · | in RN , N = (n2 ) × (n2 ), defined by
|I + hB| − 1
µ(B) = lim ,
h→0+ h
see [17]. The following global stability result is [14, Theorem 3.5].
Lemma 6.2. Assume that D is simple connected and that assumptions (H3) and
(H4) hold. Then the unique equilibrium x̄ of (6.2) is global stable in D if q 2 < 0.
Proof. From the discussions of Theorem 3.1 and Remark 5.3, we conclude that
model (1.1) satisfies the assumptions (H3) and (H4) in D. From the proof of
Theorem 2.3, we have limt→∞ N (t)=1. Thus, we just need to consider the following
EJDE-2017/121 DYNAMICS OF A SIRC EPIDEMIOLOGICAL MODEL 11
Therefore,
Ṡ
µ1 (B11 ) = − γ − βI, |B12 | = max{µβI, η}, |B21 | = 2δ + βC,
S
Ṡ
µ1 (B22 ) = − γ + max{−δ − βI − η − βS − µβC + α, βS − δ − βI − η}.
S
Thus, we have
Ṡ Ṡ
g1 = − γ + max{(µ − 1)βI, −βI + η} < − γ + η,
S S
Ṡ
g2 = − γ + max{δ − βI − η − βS + (1 − µ)βC + α, βS + δ − βI − η + βC}
S
Ṡ
< − γ + δ − η + max{α + βC, βS + βC}
S
Ṡ
< − γ + δ − η + α + β.
S
This leads to
Ṡ
µ1 (B) ≤ − γ + max{η, δ − η + α + β},
S
where ω = max{η, δ − η + α + β} < γ. Consequently,
1 t
Z
1 S(t)
µ1 (B)ds ≤ log − (γ − ω),
t 0 t S(0)
which yields q2 < 0. The proof is complete.
Remark 6.4. Theorems 5.2 and 6.3 describe that the disease always persists and
becomes endemic at an endemic level, no matter how small size the initial value of
infections has. To eradicate the disease, what we need to do is to reduce the key of
threshold value R0 to below 1.
In view of Theorem 6.3, we can obtain the global asymptotical stability of en-
demic equilibrium E ∗ of model (1.1) when R0 > 1, η < γ and δ−η+α+β < γ. In the
following theorem, we shall see that the constraint that η < γ and δ−η+α+β < γ is
not necessary if we consider the global asymptotical stability of endemic equilibrium
of model (1.2), which can be regarded as a special case of model (1.1).
Theorem 6.5. If R0 > 1 then the unique endemic equilibrium of (1.2) is globally
asymptotically stable.
Proof. It is easy to see that if R0 > 1 then model (1.2) has a unique endemic
equilibrium G∗ = (S1∗ , I1∗ , R1∗ ), where
γ+α (γ + δ)(β − γ − α) α(β − γ − α)
S1∗ = , I1∗ = , R1∗ = .
β (γ + δ + α)β (γ + δ + α)β
Meanwhile, we further see that G∗ is locally asymptotically stable. Next, we show
that G∗ is globally asymptotically stable by considering the following Lyapunov
function
1 I k5
V3 = (S − S1∗ + I − I1∗ + R − R1∗ )2 + k4 (I − I1∗ + I1∗ log ∗ ) + (S − S1∗ + I − I1∗ )2 ,
2 I1 2
EJDE-2017/121 DYNAMICS OF A SIRC EPIDEMIOLOGICAL MODEL 13
where k4 and k5 are positive constants to be determined later. The time derivative
of V3 along the solutions of (1.2) is
dV3 I − I1∗ ˙
= (S −S1∗ +I −I1∗ +R−R1∗ )(Ṡ + I˙ + Ṙ)+k4 ˙
I +k5 (S −S1∗ +I −I1∗ )(Ṡ + I).
dt I
Note that
γ(1 − S1∗ ) − βS1∗ I1∗ + δR1∗ = 0,
βS1∗ I1∗ − (γ + α)I1∗ = 0,
αI1∗ − (γ + δ)R1∗ = 0.
Then, we have
dV3
= (S − S1∗ + I − I1∗ + R − R1∗ )(−γ)(S − S1∗ + I − I1∗ + R − R1∗ )
dt
+ k4 (I − I1∗ )β(S − S1∗ ) + k5 (S − S1∗ + I − I1∗ )(−γ(S − S1∗ )
+ δ(R − R1∗ ) − (γ + α)(I − I1∗ ))
= −γ(1 + k5 )(S − S1∗ )2 − (γ + k5 (γ + α))(I − I1∗ )2 − γ(R − R1∗ )2
+ (−2γ + k4 β − k5 (2γ + α))(S − S1∗ )(I − I1∗ )
+ (−2γ + k5 δ)(S − S1∗ )(R − R1∗ ) + (−2γ + k5 δ)(I − I1∗ )(R − R1∗ ).
Take k4 = 2γ(δ + 2γ + α)/(δβ) and k5 = 2γ/δ, then we have
dV3 2γ 2γ
= −γ(1 + )(S − S1∗ )2 − (γ + (γ + α))(I − I1∗ )2 − γ(R − R1∗ )2 ≤ 0.
dt δ δ
Therefore, the LaSalle’s invariant principle [12, 16] implies that G∗ is globally
asymptotically stable.
7. Numerical simulation
In this section, we aim to provide a numerical simulation to substantiate the
theoretical results established in the previous sections by using the Runge-Kutta
fourth order iterative method. Consider model (1.1) with the parameters given as
follows: γ = 0.3, η = 0.2, µ = 0.05, α = 0.5, δ = 0.5, β = 0.4.
In Figure 1, we set up two sets of initial values. One case is that S(0) =
0.3, I(0) = 0.5, R(0) = 0, C(0) = 0. Another case is that S(0) = 20, I(0) =
20, R(0) = 0, C(0) = 0. It is observed in Figure 1 that all trajectories of model
(1.1) eventually stay in the positively invariant set D regardless of whether or not
the initial values are in D and that we can obtain the pivotal threshold value
R0 = 0.5 for the choice of parameters. In this case, it follows from Theorem 3.1
that (1.1) has a unique equilibrium E0 = (1, 0, 0, 0). Theorem 4.2 means that this
disease-free equilibrium is globally asymptotically stable.
Figure 1 shows that the infected individuals are eventually eradicated from the
crowd, while the susceptible individuals will ultimately approach the maximum
value. The epidemiological implication of Figure 2(a) is that the infected population
vanish over time. In other words, the disease will die out in the long time. If we
change the value of β into 0.98, then R0 = 1.225 and hence E0 is unstable (see
Figure 2(b)).
If (γ, η, µ, α, δ, β) = (0.3, 0.2, 0.05, 0.5, 0.5, 0.9), then there exists a unique en-
demic equilibrium E ∗ which is locally asymptotically stable (see Theorem 6.1).
Figure 3 implies that the number of infected individuals persist and gradually tend
14 H. LI, S. GUO EJDE-2017/121
(a) (b)
35
S
I
30 R
1
S C
I 25
0.8 R
Populations
Populations
C
20
0.6
15
0.4
10
0.2
5
0 0
0 5 10 15 20 25 0 5 10 15 20 25 30
Time Time
(a) (b)
0.14 0.1
0.12
0.08
0.1
0.06
0.08
R
0.06
0.04
0.04
0.02
0.02
0 0
0.8 0.4
0.6 1 0.3 1
0.8 0.8
0.4 0.2
0.6 0.6
0.2 0.1
0.4 0.4
I 0 0.2 S I 0 0.2 S
6.3 are satisfied, and hence that E ∗ is globally asymptotically stable. The epidemi-
ological implication of Figure 4 is that the infected individuals always exist if its
initial value is non-negative.
0.9
0.8
S
0.7
I
Populations
0.6 R
C
0.5
0.4
0.3
0.2
0.1
0
0 2 4 6 8 10 12 14 16 18 20
Time
(a) (b)
1
0.9
0.03
0.8
S
I 0.025
0.7
R
0.02
Populations
0.6 C
0.015
R
0.5
0.01
0.4
0.005
0.3
0
0.5
0.2
0.4 1
0.3 0.8
0.1
0.2 0.6
0.1 0.4
0
0 10 20 30 40 50 60 I 0 0.2 S
Time
In what follows, we carry out sensitivity analysises for (1.1) by the change of
the recovery rate α as well as the contact rate β. In Figure 5, we study the effect
16 H. LI, S. GUO EJDE-2017/121
0.5 α=0.005
0.8
α=0.08
0.4 α=0.9
0.6
α=0.005
S
0.3
I
0.4 α=0.08
0.2
α=0.9
0.2
0.1
0 0
0 5 10 15 20 25 30 0 1 2 3 4 5 6 7
Time Time
0.2
α=0.005 α=0.005
0.15
α=0.08 α=0.08
0.1
α=0.9 α=0.9
C
R
0.1
0.05
0.05
0 0
0 1 2 3 4 5 6 7 0 2 4 6 8 10
Time
Time
1
β=0.02
0.8 0.5 β=0.5
β=0.9
0.4
0.6 β=0.02
S
β=0.5 0.3
I
0.4 β=0.9 0.2
0.2 0.1
0 0
0 1 2 3 4 5 6 7 0 2 4 6 8
Time Time
0.2 0.1
β=0.02
β=0.02
β=0.5
0.15
0.08 β=0.5
β=0.9
β=0.9
0.06
R
C
0.1
0.04
0.05
0.02
0 0
0 2 4 6 8 0 5 10 15
Time Time
It is noted that the model always has a disease-free equilibrium, which is globally
asymptotically stable when R0 ≤ 1. When R0 > 1, we apply the Routh-Hurwitz
criterion to prove that the model has a unique endemic equilibrium, which is locally
asymptotically stable. In this case, the disease-free equilibrium become unstable.
Based on Li-Muldowney’s global-stability criterion [14], we show that the unique
endemic equilibrium can be globally asymptotically stable in a feasible region, i.e.,
influenza A becomes endemic. Although we have established the global stability of
the unique endemic equilibrium E ∗ when R0 > 1, our results are obtained under
the assumptions that η < γ and δ −η +α+β < γ. From Theorem 6.5, we conjecture
that the condition that η < γ and δ − η + α + β < γ is not necessary. Therefore,
the perspective of our work is to show the assumption that R0 > 1 is a sufficient
and necessary condition ensuring the globally asymptotical stability of the unique
endemic equilibrium E ∗ .
Acknowledgments. This work has been supported by the Natural Science Foun-
dation of China (Grant No. 11671123).
References
[1] V. Andreasen, J. Lin, S. A. Levin; The dynamics of cocirculating influenza strains conferring
partial cross-immunity, Journal of Mathematical Biology, 35 (1997), 825–842.
[2] R. M. Anderson, R. M. Mary; Infectious Diseases of Humans: Dynamics and Control, Oxford
University Press, 1991.
[3] R. M. Anderson, R. M. May; Population biology of infectious diseases: Part I, Nature, 280
(1979), 361–367.
[4] G. J. Butler, P. Waltman; Persisitence in dynamical systems, Journal of Differential Equa-
tions 63 (1986), 255–263.
[5] R. Casagrandi, L. Bolzoni, S. A. Levin, V. Andreasen; The SIRC model and influenza A,
Mathematical Biosciences 200 (2006), 152–169.
18 H. LI, S. GUO EJDE-2017/121
[6] K. W. Chung, R. Lui; Dynamics of two-strain influenza model with cross-immunity and no
quarantine class, Journal of Mathematical Biology, 73 (2016), 1467–1489.
[7] D. J. Earn, J. Dushoff, S. A. Levin; Ecology and evolution of the flu, Trends in Ecology and
Evolution, 17 (2002), 334–340.
[8] M. Fiedler; Additive compound matrices and an inequality for eigenvalues of symmetric
stochastic matrices, Czechoslovak Mathematical Journal, 99 (1974), 392–402.
[9] H. I. Freedman, S. Ruan, M. Tang; Uniform persistence and flows near a closed positively
invariant set, Journal of Dynamics and Differential Equations, 6 (1994), 583–600.
[10] W. O. Kermack, A. G. McKendrick; Contributions to the mathematical theory of epidemics:
II. The problem of endemicity, Proceedings of Royal Social London. Series A, Containing
Paper of a Mathemeatical and Physiacal Character, 138 (1932), 55–83.
[11] B. W. Kooi, M. Aguiar, N. Stollenwerk; Bifurcation analysis of a family of multi-strain
epidemiology models, Journal of Computational and Applied Mathematics, 252 (2013), 148–
158.
[12] J. P. LaSalle; The stability of dynamical systems, Society for Industrial and Applied Mathe-
matics, 27 (1976), 1121–1130.
[13] J. Lin, V. Andreasen, S. A. Levin; Dynamics of influenza A drift: the linear three-strain
model, Mathematical Biosciences, 162 (1999), 33–51.
[14] M. Y. Li, J. S. Muldowney; A geometric approach to global stability problems, SIAM Journal
on Mathematical Analysis, 27 (1996), 1070–1083.
[15] M. Y. Li, J. S. Muldowney; On R. A. Smith’s autonomous convergence theorem, The Rocky
Mountain Journal of Mathematics, 25 (1995), 365–378.
[16] A. M. Lyapunov; The general problem of the stability of motion, International Journal of
Control, 31 (1992), 531–534.
[17] R. H. Martin, Jr.; Logarithmic norms and projections applied to linear differential systems,
Journal of Mathematical Analysis and Applications, 45 (1974), 432–454.
[18] P. Minayev, N. Ferguson; Improving the realism of deterministic multi-strain models: impli-
cations for modelling influenza A, Journal of The Royal Society Interface, 6 (2009), 509–518.
[19] J. S. Muldowney; Compound matrices and ordinary differential equations, Rocky Mountain
Journal of Mathematics, 20 (1990), 857–872.
[20] M. Nuño, Z. Feng, M. Martcheva, C. Castillo-Chavez; Dynamics of two-strain influenza with
isolation and partial cross-immunity. SIAM Journal on Applied Mathematics, 65 (2005),
964–982.
[21] M. Nuño, M. Martcheve, C. Castillo-Chavez; Immune level approach for multiple strain
pathogens, Journal of Biological Systems, 17 (2009), 713–737.
[22] D. J. Smith, S. Forrest, D. H. Ackley, A. S. Perelson; Variable efficacy of repeated annual
influenza vaccination. Proceedings of the National Academy of Science, 96 (1999) ,14001–
14006.
[23] T. Szucs; The socio-economic burden of influenza, Journal of Antimicrobial Chemotherapy,
44 (1999), 11–15.
[24] J. K.Taubenberger, A. H. Reid, A. E. Krafft, K. E. Bijwaard, T. G. Fanning; Initial genetic
characterization of the 1918 Spanish influenza virus, Science 275 (1997), 1793–1796.
[25] X. Yang, L. Chen, J. Chen; Permanence and positive periodic solution for the single-species
nonautonomous delay diffusive models, Computers and Mathematics with Applications, 32
(1996) 109–116.
[26] X. Zhong, S. Guo, M. Peng; Stability of stochastic SIRS epidemic models with saturated
incidence rates and delay, Stochastic Analysis and Applications, 35 (2017), 1–26.
Haijiao Li
School of Business Administration, and College of Mathematics and Econometrics,
Hunan University, Changsha, Hunan 410082, China
E-mail address: [email protected]
Shangjiang Guo
College of Mathematics and Econometrics, Hunan University, Changsha, Hunan 410082,
China
E-mail address: [email protected]