SEIR {SimInf} | R Documentation |
Create an SEIR model to be used by the simulation framework.
SEIR(u0, tspan, events = NULL, beta = NULL, epsilon = NULL, gamma = NULL)
u0 |
A |
tspan |
A vector (length >= 2) of increasing time points
where the state of each node is to be returned. Can be either
an |
events |
a |
beta |
The transmission rate from susceptible to exposed. |
epsilon |
The incubation rate from exposed to infected. |
gamma |
The recovery rate from infected to recovered. |
The SEIR model contains four compartments; number of susceptible (S), number of exposed (E) (those who have been infected but are not yet infectious), number of infectious (I), and number of recovered (R). Moreover, it has three state transitions,
S -- beta S I / N --> E
E -- epsilon E --> I
I -- gamma I --> R
where β is the transmission rate, ε is the incubation rate, γ is the recovery rate, and N=S+E+I+R.
The argument u0
must be a data.frame
with one row for
each node with the following columns:
The number of sucsceptible in each node
The number of exposed in each node
The number of infected in each node
The number of recovered in each node
SEIR
## Create a SEIR model object. model <- SEIR(u0 = data.frame(S = 99, E = 0, I = 1, R = 0), tspan = 1:100, beta = 0.16, epsilon = 0.25, gamma = 0.077) ## Run the SEIR model and plot the result. set.seed(3) result <- run(model) plot(result)