Example of SIRS model generated from the pomp package.

This is a stochastic (SIR)-> (SIR) model (two sets of infections in sequence) created from EpiCompare::simulate_agents(). The model is S1->I1->R1->S2->I2->R2. It is two SIR models pasted together with a transition between first recovery and second suscepbiility. The parameters are \(beta_1 = .5\), \(\beta_2 = .08\), \(\gamma_1 = .4\), \(\gamma_2 = .06\), \(\eta = .08\) where \(\eta\) is the transition from R1 to S2.

sirs_data

Format

The object is a data frame of dimension 7525x8. The columns are

sim

simulation number between 1 and 25

t

time between 0 and 300

S1

susceptible for the first time

I1

infectious for first time

R1

recovered for first time

S2

susceptible for second time

I2

infectious for second time

R2

recovered for second and final time

Examples

library(ggplot2)
library(tidyr)
head(sirs_data)
#> # A tibble: 6 x 8
#> # Groups:   sim [1]
#>     sim     t    S1    I1    R1    S2    I2    R2
#>   <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl>
#> 1     1     0  1150     0     0     0     0     0
#> 2     1     1   950    50     0   100    50     0
#> 3     1     2   906    78    16   100    50     0
#> 4     1     3   857    96    44   102    49     2
#> 5     1     4   798   109    87   102    51     3
#> 6     1     5   744   122   121   107    49     7