rSIR: Simulates from a simple stochastic SIR epidemic model.

In jessseliu/epidmod: epidmod

Description

rSIR returns a matrix contains information about time, number of susceptible people, infected people, recover people, total number of people.

Usage

rSIR(
  N0 = 1000,
  S0 = 1000,
  I0 = 0,
  R0 = N0 - I0 - S0,
  days = 100,
  pars = c(1, 2, 1, 0.1, 0.9)
)

Arguments

N0	An integer. The population size at time 0.
S0	An integer. The initial number of susceptible people.
I0	An integer. The initial number of infected people.
R0	An integer. The initial number of recovered people.
days	An integer. The number of days for which to simulate.
pars	A numeric vector: (a.rate, i.rate, r.rate, pI, pIm). Description about pars *a.rate represents rate of arrival of people into the population. *i.rate represents individual effective contact rate. *r.rate represents recovery time for each infected individual is 1/ r.rate. *pI represents probability that an arrival is infected. *pIm represents probability that an infected person is immune after recovery.

Value

A numeric matrix with 5 columns. Row i contains the values of (t, S_t, I_t, R_t, N_t) at time t.

Examples

model1 <- rSIR (N0 = 100, S0 = 100, I0 = 0, R0 = 0,  days = 100, pars = c(1, 4, 2, 0.3, 0.7))
model1
print(model1)
plot(model1)

jessseliu/epidmod documentation built on Sept. 15, 2020, 8:53 p.m.