# Dynamic Analysis of an SEIR Model with Distinct Incidence for Exposed and Infectives

Authors: Junhong Li; Ning Cui
Journal: The Scientific World Journal
 (2013)

---

## Abstract

An SEIR model with vaccination strategy that incorporates distinct incidence rates for the exposed and the infected populations is studied. By means of Lyapunov function and LaSalle’s invariant set theorem, we proved the global asymptotical stable results of the disease-free equilibrium. The sufficient conditions for the global stability of the endemic equilibrium are obtained using the compound matrix theory. Furthermore, the method of direct numerical simulation of the system shows that there is a periodic solution, when the system has three equilibrium points.


---

## Body

## 1. Introduction

Mathematical models have become important tools in analyzing the spread and the control of infectious diseases. Many infectious diseases in nature, such as measles, HIV/AIDS, SARS, and tuberculosis (see [1–6]), incubate inside the hosts for a period of time before the hosts become infectious. Li and Fang (see [7]) studied the global stability of an age-structured SEIR model with infectivity in latent period. Yi et al. (see [8]) discussed the dynamical behaviors of an SEIR epidemic system with nonlinear transmission rate. Li and Zhou (see [9]) considered the global stability of an SEIR model with vertical transmission and saturating contact rate.In this paper, we will consider an SEIR model that the diseases can be infected in the latent period and the infected period. The population sizeN is divided into four homogeneous classes: the susceptible S(t), the exposed (in the latent period) E(t), the infective I(t), and the recovered R(t). It is assumed that all the offsprings at birth are susceptible to the disease. The inflow rate (including birth and immigration) and outflow rate (including natural death and emigration) are denoted by b and d, respectively. The rate of disease-caused death is taken as α. We assume that susceptible individuals are vaccinated at a constant per capita rate a. Due to the partial efficiency of the vaccine, only σ fraction of the vaccinated susceptibles goes to the recovered class. The remained 1-σ fraction of the vaccinated susceptibles has no immunity at all and goes to the exposed class after infected by contact with the infectives. If σ=0, it means that the vaccine has no effect at all, and if σ=1, the vaccine is perfectly effective. The positive parameter ε is the rate at which the exposed individuals become infectious. γ is the constant rate, at which the infectious individuals recover with acquiring permanent immunity. The transfer mechanism from the class S(t) to the class E(t) is guided by the function β(I+qE)/N, where β is the force of infection. q denotes the relative measure of infectiousness for the asymptomatic class E(t).Based on these considerations, and with reference to [10–12], the SEIR model is given by the following system of differential equations:
(1)S′=bN-F(S,E,I)-(σa+d)S,E′=F(S,E,I)-(d+ε)E,I′=εE-(α+γ+d)I,R′=σaS+γI-dR,0=S+E+I+R-N,
where the derivative d/dt is denoted by ′ and F(S,E,I)=βS(1-a)(I+qE)/N+βaS(1-σ)(I+qE)/N.Thus, the total population sizeN implies N′=(b-d)N-αI. Let s=S/N, e=E/N, i=I/N, and r=R/N.Because the variableR does not appear in the equations of S, E, and I, we only need to consider the following subsystem:
(2)s′=b-βs(1-σa)(i+qe)-(σa+b)s+sαi,e′=βs(1-σa)(i+qe)-(b+ε)e+eαi,i′=εe-(α+γ+b)i+αi2,0=s+e+i+r-1.
The system (2) is equivalent to (1). From biological considerations, we study (2) in the following closed set:
(3)T={(s,e,i)∈R+3∣0≤s+e+i≤1},
where R+3 denotes the nonnegative cone of R3 including its lower dimensional faces.


## 2. Equilibria and Global Stability

It is easy to visualize that (2) always has a disease-free equilibrium P0(b/(σa+b),0,0). The Jacobian matrix of (2) at an arbitrary point P(s,e,i) takes the following form:
(4)J(P)=[J11-qβ(1-σa)sJ13J21J22J230εJ33],
where
(5)J11=-β(1-σa)(i+qe)-(σa+b)+αi,J13=αs-β(1-σa)s,J23=αe+β(1-σa)s,J33=2αi-(α+γ+b),J22=βq(1-σa)s-(b+ε)+αi,J21=β(1-σa)(i+qe).Theorem 1.
IfR0<1, the disease-free equilibrium P0 is locally asymptotically stable, where
(6)R0=bβ(1-σa)(qα+qγ+qb+ε)(b+ε)(α+b+γ)(σa+b).Proof.
Let(7)λ1=(ε+qα+qγ+qb)(α+γ+2b+ε)(α+γ+b)(b+ε)qR0-1.
We calculate the characteristic equation ofJ(P0) as follows:
(8)(λ+σa+b)×(λ2+bqβ(1-σa)λ1λ(σa+b)+(b+ε)(α+γ+b)(1-R0)bqβ(1-σa)λ1λ(σa+b))=0.
The stability ofP0 is equivalent to all eigenvalues of (8) being with negative real parts, which can be guaranteed by R0<1. Consequently, the disease-free equilibrium is local asymptotical stability. This proves the theorem.Theorem 2.
IfR01≤1, the disease-free equilibrium P0 is globally asymptotically stable, where
(9)R01=β(1-σa)(qα+qγ+qb+ε)(b+ε)(α+b+γ).Proof.
Consider the following function:(10)L=(α+γ+b)e+β(1-σa)i.
Its derivative along the solutions to the system (2) is as follows:
(11)L′=i[(α+γ+b)β(1-σa)s-(1-σa)(α+γ+b-αi)]+e[(α+γ+b)qβ(1-σa)s-(α+γ+b)(b+ε+αi)+(1-σa)ε]≤e(α+γ+b)(b+ε)(R01-1)≤0.
Furthermore,L′=0 only if e=0. The maximum invariant set in {(s,e,i)∈T:L′=0} is the singleton {P0}. When R01≤1, the global stability of P0 follows from LaSalle’s invariance principle (see [13]). This completes the proof.Theorem 3.
Equation (2) has a unique endemic equilibrium P*(s*,e*,i*) if b≥α and Δ>1, where
(12)Δ=(b+ε-α)(γb+1)×(1ε-(σa+b-α)β(1-σa)(qb+qγ+ε)).Proof.
Let the right side of each of the first three differential equations equal to zero in (2); we obtain the following:
(13)βs*(1-σa)(i*+qe*)-(b+ε)e*+αe*i*=0,εe*-(α+γ+b)i*+αi*2=0,b-βs*(1-σa)(i*+qe*)-(σa+b)s*+αs*i*=0,
with s*>0, e*>0 and i*>0. So we get
(14)bs*=β(1-σa)(i*+qe*)-(σa+b)+αi*,βs*(1-σa)(i*+qe*)e*=(b+ε)+αi*,εe*=(α+γ+b)i*-αi*2.
When the three equations of (14) are multiplied together, we obtain the following:
(15)bβε(1-σa)[(α+γ+b)qε+1-αqi*ε]=(b+ε-αi*)(α+γ+b-αi*)×[-αqβ(1-σa)i*2ε+(β(1-σa)(αq+γq+bq+ε)ε+α)i*-(σa+b)-αqβ(1-σa)i*2ε].
Define the following:
(16)f(i)=f1(i)(b+ε-αi)(α+γ+b-αi),g(i)=bβε(1-σa)[(α+γ+b)qε+1-αqiε],
where
(17)f1(i)=[-αqβ(1-σa)i2ε+(β(1-σa)(αq+γq+bq+ε)ε+α)i-(σa+b)-αqβ(1-σa)i2ε],
and the roots of f(i) are i1=(b+ε)/α>1, i2=(α+γ+b)/α>1, and the other two are i3, i4 which satisfy f1(i)=0. i0=(α+γ+b)/α+ε/(αq) is the root of g(i). Direct calculations show the following:
(18)g(1)=bβε(1-σa)(qb+qγ+ε)>0,f(0)=-(σa+b)(b+ε)(α+γ+b)<0,f(1)=Δg(1),i0>i1,i0>i2,i3i4=ε(σa+b)αqβ(1-σa)>0,i3+i4=1+b+γα+ε(1αq+1qβ(1-σa))>2.
Since Δ>1, the linear function g(i) has exactly one intersection with the function f(i) where i lies in the interval (0,1). Furthermore, s* and e* can be uniquely determined from i* by the following:
(19)e*=i*(α+γ+b-αi*)ε,s*=b[β(1-σa)(i*+e*)+(σa+b)+αi*].
From this, we can easily see that (2) has a unique endemic equilibrium. This completes the proof.Denote the interior ofT by T∘. In this paper, we obtain sufficient conditions that the equilibrium is globally asymptotically stable using the geometrical approach of Li and Muldowney in [14].Theorem 4.
The unique endemic equilibriumP* is globally asymptotically stable in T∘, when
(20)b≥α,Δ>1,R01>1,α<min{σa,b,ε}.Proof.
SinceR01>1, namely, Δ<R0 and P0 is unstable, we can easy see that (2) satisfies the assumptions (H1) and (H2) (see [14]) in the interior of its feasible region T. The unique equilibrium is locally asymptotically stable using simple calculation.
Letx=(s,e,i) and f(x) denotes the vector field of (2) and
(21)J[2](P)=(a11qβs(1-σa)+αe-αs+βs(1-σa)εa22-qβs(1-σa)0β(1-σa)(i+qe)a33),
where
(22)a11=-β(1-σa)(i+qe-qs)-(σa+2b+ε)+2αi,a22=-β(1-σa)(i+qe)-(2b-d+α+β)+3αi,a33=-qβs(1-σa)-(2b+ε+α+β)+3αi.
Set the the following function:
(23)P(x)=P(s,e,i)=diag(1,ei,ei).
Then the matrix B=PfP-1+PJ[2]P-1 can be written in block form as follow:
(24)B=(B11B12B21B22),
where
(25)B11=a11,B21=(εei0)′,B12=(β(1-σa)sie+αi(β(1-σa)-α)sie),B221=e′e-i′i+a22,B222=e′e-i′i+a33,B22=(B221-qβs(1-σa)β(1-σa)(i+qe)B222).
Let(u,v,w) denote the vectors in R3≅R(32), we select a norm in R3 as |(u,v,w)|=max{|u|,|v|+|w|} and let ρ denote the Lozinskii measure with respect to this norm. Using the method of estimating ρ in [15], we have ρ(B)≤sup(g1,g2), where
(26)g1=ρ1(B11)+|B12|,g2=|B21|+ρ1(B22),|B12|, |B21| are matrix norms with respect to the l1 vector norm, and ρ1 denotes the Lozinskii measure with respect to l1 norm. More specifically, ρ1(B11)=B11, |B21|=εe/i, and |B12|=β(1-σa)si/e+αi.
Rewriting the system (2), we have the following:
(27)εei=i′i+(α+γ+b)-αi.
Therefore,
(28)g1=e′e-β(1-σa)(i+qe)-(σa+b)+2αi,g2=e′e+2αi-b-min{σa,ε}.
Since b≥α and α<min{σa,b,ε}, there is m=max{σa+b-2α,b+ε-2α}>0 and
(29)ρ(B)≤sup(g1,g2)≤e′e-mfort>t-.
Along each solution x(t,x0) to (2) such that x0∈K, where K is the compact absorbing set, we thus have the following:
(30)1t∫0tρ(B)ds≤1tloge(t)e(t-)+1t∫0t-ρ(B)ds-mt-t-t,
which implies
(31)q-2≤-m2<0.
This completes the proof.


## 3. Conclusion

In this paper, we discuss an SEIR model that the diseases can be infected in the latent period and the infected period. The vaccine effectiveness is also taken into account. We investigate the global dynamics of the reduced proportional system. IfR01≤1, the disease-free equilibrium P0 is globally asymptotically stable. The unique equilibrium P* of the system (2) is globally asymptotically stable in T∘, when b≥α, Δ>1, R01>1, and α<min{σa,b,ε}. When q=a=γ=0, (2) becomes the SEIR model without infectivity in latent and disease-caused death (see [8]). When q=0, (2) becomes the SEIR model without infectious in latent (see [12]).The parameters are considered in the following cases:(32)ε=0.05,γ=0.003,α=0.002,β=0.05,
(see [16]), and
(33)q=1,b=0.00001,σa=0.000001.
At this case, there are three fixed points:
(34)E1(0.9091,0,0),E2(0.0911,0.00018673,0.0019),E3(0.0000083926,0.025,24.7973),
and E1, E3 are unstable, E2 stable. The method of direct numerical simulation of (2) shows that there is a periodic solution. The phase portraits of the system (2) in Figure 1(a), and time series of s,e,i are given in Figures 1(b), 1(c), and 1(d).(a) Phase portraits of (2). (b) Time series of s. (c) Time series of e. (d) Time series of i.
(a)
(b)
(c)
(d)



---

*Source: 871393-2013-05-26.xml*