Introduction To SEIR Models
Introduction To SEIR Models
Introduction To SEIR Models
Nakul Chitnis
SI Model
SIS Model
SIR Model
SEIR Model
2017-05-08 2
Mathematical Models of Infectious Diseases
Population-based models
I Can be deterministic or stochastic
I Continuous time
• Ordinary differential equations
• Partial differential equations
• Delay differential equations
• Integro-differential equations
I Discrete time
• Difference equations
Agent-based/individual-based models
I Usually stochastic
I Usually discrete time
2017-05-08 3
Mathematical Models of Infectious Diseases
Population-based models
I Can be deterministic or stochastic
I Continuous time
• Ordinary differential equations
• Partial differential equations
• Delay differential equations
• Integro-differential equations
I Discrete time
• Difference equations
Agent-based/individual-based models
I Usually stochastic
I Usually discrete time
2017-05-08 3
Outline
SI Model
SIS Model
SIR Model
SEIR Model
2017-05-08 4
SI Model
rβI/N
S I
dS I
= −rβS
dt N
dI I
= rβS
dt N
S: Susceptible humans
I: Infectious humans
r: Number of contacts per unit time
β: Probability of disease transmission per contact
N: Total population size: N = S + I.
2017-05-08 5
Analyzing the SI Model
Idfe = 0
Iee = N
2017-05-08 6
Analyzing the SI Model
Idfe = 0
Iee = N
2017-05-08 6
Numerical Solution of SI Model
1000
800
Infectious Humans
600
400
200
0
0 5 10 15 20
Time (Years)
With r = 365/3 years−1 , β = 0.005, N = 1000, and I(0) = 1.
2017-05-08 7
Definition of Transmission Parameters
2017-05-08 8
Outline
SI Model
SIS Model
SIR Model
SEIR Model
2017-05-08 9
SIS Model
rβI/N
S I
dS I
= −rβS + γI
dt N
dI I
= rβS − γI
dt N
γ: Per-capita recovery rate
2017-05-08 10
Analyzing the SIS Model
for I(0) = I0 .
Equilibrium Points:
Idfe = 0
(rβ − γ)N
Iee =
rβ
2017-05-08 11
Analyzing the SIS Model
for I(0) = I0 .
Equilibrium Points:
Idfe = 0
(rβ − γ)N
Iee =
rβ
2017-05-08 11
Numerical Solution of SIS Model
1000
800
Infectious Humans
600
400
200
0
0 10 20 30 40 50
Time (Days)
With rβ = 0.5 days−1 , γ = 0.1 days−1 , N = 1000, and I(0) = 1.
2017-05-08 12
Outline
SI Model
SIS Model
SIR Model
SEIR Model
2017-05-08 13
of uncertainty that may be prioritized for urgent res
The Basic Reproductive Number (R0 )
Generation
0 1 2
2017-05-08 16
Reproductive Numbers
2017-05-08 17
Reproductive Numbers
2017-05-08 17
Reproductive Numbers
2017-05-08 17
Reproductive Numbers
2017-05-08 17
Reproductive Numbers
2017-05-08 17
Evaluating Re
Number of Probability of
Duration of
Re (t) = contacts transmission
infection
per unit time per contact
Proportion of
× susceptible
population
2017-05-08 18
The Basic Reproductive Number (R0 )
http://www.cameroonweb.com/CameroonHomePage/NewsArchive/Ebola-How-does-it-compare-316932
2017-05-08 19
Outline
SI Model
SIS Model
SIR Model
SEIR Model
2017-05-08 20
SIR Model
rβI/N γ
S I R
dS I
= −rβS
dt N
dI I
= rβS − γI
dt N
dR
= γI
dt
R: Recovered humans
with N = S + I + R.
2017-05-08 21
Analyzing the SIR Model
2017-05-08 22
R0 for the SIR Model
Number of Probability of
Duration of
R0 = contacts transmission
infection
per unit time per contact
1
R0 = r × β ×
γ
rβ
=
γ
If R0 < 1, introduced cases do not lead to an epidemic (the
number of infectious individuals decreases towards 0).
If R0 > 1, introduced cases can lead to an epidemic
(temporary increase in the number of infectious individuals).
rβ S(t)
Re (t) =
γ N
2017-05-08 23
Phase Portrait of SIR Model
THE MATHEMATICS OF INFECTIOUS DISEASES 605
0.8 σ=3
infective fraction, i
0.6
0.4
0.2
0 ↑
0 0.2 1 0.4
smax = σ 0.6 0.8 1
susceptible fraction, s
Hethcote (2000)
2017-05-08 Fig. 2 Phase plane portrait for the classic SIR epidemic model with contact number σ = 3. 24
Numerical Solution of SIR Model
1000
800
Susceptible
Infectious
Humans
600 Recovered
400
200
0
0 20 40 60 80 100
Time (Days)
With rβ = 0.3 days−1 , γ = 0.1 days−1 , N = 1000, and
S(0) = 999, I(0) = 1 and R(0) = 0.
2017-05-08 25
Numerical Solution of SIR Model
1000
800
Susceptible
Infectious
Humans
600 Recovered
400
200
0
0 20 40 60 80 100
Time (Days)
With rβ = 0.3 days−1 , γ = 0.1 days−1 , N = 1000, and
S(0) = 580, I(0) = 20 and R(0) = 400.
2017-05-08 26
Human Demography
2017-05-08 27
Endemic SIR Model
Λ rβI/N γ
Birth S I R
µ µ µ
dS I
= Λ − rβS − µS
dt N
dI I
= rβS − γI − µI
dt N
dR
= γI − µR
dt
N =S+I +R
Λ Λ(γ + µ)
Sdfe = See =
µ rβµ
Λ(rβ − (γ + µ))
Idfe = 0 Iee =
rβ(γ + µ)
γΛ(rβ − (γ + µ))
Rdfe = 0 Ree =
rβµ(γ + µ)
2017-05-08 29
R0 for the Endemic SIR Model
Number of Probability of
Duration of
R0 = contacts transmission
infection
per unit time per contact
1
R0 = r × β ×
γ+µ
rβ
=
γ+µ
1000
800
Susceptible
Infectious
Humans
600 Recovered
400
200
0
0 20 40 60 80 100
Time (Days)
With rβ = 0.3 days−1 , γ = 0.1 days−1 , µ = 1/60 years−1 ,
Λ = 1000/60 years−1 , and S(0) = 999, I(0) = 1 and R(0) = 0.
2017-05-08 31
Numerical Solution of Endemic SIR Model
1000
Susceptible
Infectious
Recovered
800
Humans
600
400
200
0
0 20 40 60 80 100 120
Time (Years)
With rβ = 0.3 days−1 , γ = 0.1 days−1 , µ = 1/60 years−1 ,
Λ = 1000/60 years−1 , and S(0) = 999, I(0) = 1 and R(0) = 0.
2017-05-08 32
Outline
SI Model
SIS Model
SIR Model
SEIR Model
2017-05-08 33
SEIR Model
Λ rβI/N ε γ
Birth S E I R
µ µ µ µ
2017-05-08 34
SEIR Model
dS I
= Λ − rβS − µS
dt N
dE I
= rβS − εE
dt N
dI
= εE − γI − µI
dt
dR
= γI − µR
dt
with
N = S + E + I + R.
2017-05-08 35
R0 for the Endemic SEIR Model
Number of Probability of
Duration of
R0 = contacts transmission
infection
per unit time per contact
Probabililty of
× surviving
exposed stage
1 ε
R0 = r × β × ×
γ+µ ε+µ
rβε
=
(γ + µ)(ε + µ)
If R0 < 1, the disease-free equilibrium point is globally
asymptotically stable and there is no endemic equilibrium
point (the disease dies out).
If R0 > 1, the disease-free equilibrium point is unstable and a
globally asymptotically stable endemic equilibrium point exists.
2017-05-08 36
Extensions to Compartmental Models
2017-05-08 38