# NOT RUN {
# Population from Tokars (2018)
nsam <- 1e6L
ndays <- 304L
pop_tok <- sim_reference(
init_pop_size = nsam,
vaccinations = generate_counts(nsam, ndays, 0.55, mean = 100, sd = 50),
cases_novac = generate_counts(nsam, ndays, 0.12, mean = 190, sd = 35),
ve = 0.48,
lag = 14,
deterministic = TRUE
)
head(pop_tok)
sum(pop_tok$avert)
# }
Run the code above in your browser using DataLab