run_scenario: Run scenario analysis
In COVID-19-Modelling/BassSIR: Infectious disease modelling with Bass-SIR model

Run scenario analysis

1
2
3

run_scenario(sim, fn)

compare_scenarios(sim, ..., fn_change = c("Y0", "Y1", "Yt"))

`sim`	baseline simulations
`fn`	a function for modifying parameters

library(R2jags)
cases <- as_bass_data(n_covid19$Hubei)
est <- BassSIR::fit(cases, r_rec = 1/22.2, r_death = 1/22.3, type = "BassSIR")

sim <- simulate(est, nsim = 500)

zero_kappa <- function(pars) {
  pars$kappa <- rep(0, length(pars$kappa))
  return(pars)
}

lockdown <- run_scenario(sim, zero_kappa)

cp <- compare_scenarios(sim, Lockdown = lockdown)