sir: SIR models.

Description

euler.sir is a pomp object encoding a simple seasonal SIR model. Simulation is performed using an Euler multinomial approximation. gillespie.sir has the same model implemented using Gillespie's algorithm. bbs is a nonseasonal SIR model together with data from a 1978 outbreak of influenza in a British boarding school.

Usage

data(euler.sir)
data(gillespie.sir)
data(bbs)

Arguments

Details

This example is discussed extensively in the Introduction to pomp and Advanced topics in pomp vignettes.

The codes that construct these pomp objects can be found in the data-R directory in the installed package. Do file.show(system.file("data-R/sir.R",package="pomp")) to view these codes. For the basic rprocess, dmeasure, rmeasure, and skeleton functions, these codes use compiled native routines built into the package's library. View src/sir.c in the package source or file.show("examples/sir.c") from an Rsession to view these codes.

The boarding school influenza outbreak is described in Anonymous (1978).

References

Anonymous (1978). Influenza in a boarding school. British Medical Journal 1:587

Examples

Run this code

data(euler.sir)
plot(euler.sir)
plot(simulate(euler.sir,seed=20348585))
coef(euler.sir)

data(gillespie.sir)
plot(gillespie.sir)
plot(simulate(gillespie.sir,seed=20348585))
coef(gillespie.sir)

data(bbs)
plot(bbs)
coef(bbs)