run_code/test_demogcpp.R

In mrc-ide/siR: Individual-based SIR models in R and Rcpp

N <- 20000
t <- 80000
t1 <- demog_test(N, t, 1, age_of_death, equilibrium_age)

plot(tapply(t1$deaths, rep(1:800, each = 100), sum), t = "l")


for(i in 1:10){
  t1 <- demog_test(N, t, 1, age_of_death, equilibrium_age)
  lines(tapply(t1$deaths, rep(1:800, each = 100), sum), t = "l", col = rainbow(10)[i])
}