serosv (version 1.1.0)

sir_basic_model: Basic SIR model

Description

Refers to section 3.1.3.

Usage

sir_basic_model(times, state, parameters)

Value

list of class sir_basic_model with the following items

parameters

list of parameters used for fitting the model

output

matrix of population for each compartment over time

Arguments

times

time sequence.

state

the initial state of the model.

parameters

the parameters of the model.

Details

In state:

- S: number of susceptible

- I: number of infected

- R: number of recovered

In parameters:

- alpha: disease-related death rate

- mu: natural death rate (= 1/life expectancy)

- beta: transmission rate

- nu: recovery rate

- p: percent of population vaccinated at birth

Examples

state <- c(S=4999, I=1, R=0)
parameters <- c(
  mu=1/75, # 1 divided by life expectancy (75 years old)
  alpha=0, # no disease-related death
  beta=0.0005, # transmission rate
  nu=1, # 1 year for infected to recover
  p=0 # no vaccination at birth
)
times <- seq(0, 250, by=0.1)
model <- sir_basic_model(times, state, parameters)
model