@lijiang

Sculpting in time

Do one thing and do it well.
Every story has a beginning and an end.

2-Minute Read

SIR Modeling

The SIR classical model is mainly used in the field of infectious diseases to predict future trends in the number of infections.

S[t] denotes susceptible susceptible population

I[t] denotes infected population already infected

R[t] indicates recovered recovered population

Reference:

MathWorld:

Kermack-McKendrick Model

The original Kermack-McKendrick model was designed to account for changes in the number of people infected over time, like the plague that occurred in 1665-1666 and the cholera that occurred in 1865. The model assumes that the total population is fixed, that the incubation period for infectious diseases is instantaneous, that the duration of infection is the same as the disease cycle, and that the population is assumed to be non-differentiable, without differences by gender or race.

The model uses three nonlinear ordinary differential equations:

SIRModelTraditionalForm SIRmodel

  • beta denotes infection rate infection rate
  • gamma denotes recovery rate recovery rate
  • i indicates infected population
  • r indicates recovered population
  • s indicates susceptible population

Pre-set values: beta = 20%, gamma = 15%, i(0) = 0.1, r(0) = 0, s[0] = 10

Modelica source code

model model SIRModel parameter Real beta = 0.2; parameter Real gamma = 0.15; Real i; Real r; Real s; initial equation s = 10; i = 0.1; r = 0; equation der(s) = (-1) * beta * i * s; der(i) = (-1) * gamma * i + beta * i * s. der(i) = (-1) * gamma * i + beta * i * s; der(r) = gamma * i; end SIRModel;

Model prediction graph:

![SIRModelPlot](/images/mathematica/SIRModelPlot.png)

SystemModeler prediction plot

![SIRModelPlot2](/images/mathematica/SIRModelPlot2.png)



## SEIR Model ##

Doing

Recent Posts

Categories

About

Keep thinking, Stay curious
Always be sensitive to new things