forecast_ensemble: find confidence envelopes by simulation
In bbolker/McMasterPandemic: Pandemic Model

forecast_ensemble

R Documentation

find confidence envelopes by simulation

Description

find confidence envelopes by simulation

Usage

forecast_ensemble(
  fit,
  nsim = 200,
  forecast_args = fit$forecast_args,
  qvec = c(0.05, 0.5, 0.95),
  qnames = c("lwr", "value", "upr"),
  seed = NULL,
  imp_wts = FALSE,
  Sigma = bbmle::vcov(fit$mle2),
  scale_Sigma = 1,
  calc_Rt = FALSE,
  fix_pars_re = "nb_disp",
  raw_ensembles = FALSE,
  parallel = FALSE,
  n_cores = 4,
  .progress = if (interactive()) "text" else "none"
)

Arguments

`fit`	output from `calibrate`
`nsim`	number of simulations
`forecast_args`	arguments to pass to `forecast_sim`
`qvec`	vector of quantiles: NULL to return an array (nt x nvars x nsims) instead of a tibble with date/var+ quantiles
`qnames`	quantile names
`seed`	random-number seed
`imp_wts`	use importance weighting, i.e. weight ensemble based on log-likelihood?
`Sigma`	covariance matrix to pass to `pop_pred_samp`
`scale_Sigma`	multiplier for covariance matrix
`calc_Rt`	calculate and include R(t) in prediction/forecast?
`fix_pars_re`	a regular expression specifying the names of parameters that should be treated as fixed when constructing the parameter ensemble
`raw_ensembles`	(logical) return ensembles (FIXME: should implement return_type=c("array","quantiles","raw") (but not really raw, "dataframe"??))
`parallel`	whether to attempt parallel processing
`n_cores`	number of cores for parallel forecasting (ignored if not parallel) (??)
`.progress`	progress bar?

Other classic_macpan: add_d_log(), add_updated_vaxrate(), aggregate_agecats(), calibrate_comb(), calibrate(), check_age_cat_compatibility(), check_contact_rate_setting(), col_multiply(), condense_age(), condense_params_vax(), condense_state(), condense_vax(), dev_is_tikz(), do_step(), expand_params_age(), expand_params_desc_age(), expand_params_desc_variant(), expand_params_desc_vax(), expand_params_mistry(), expand_params_variant(), expand_params_vax(), expand_state_age(), expand_state_vax(), expand_stateval_testing(), fix_pars(), fix_stored(), forecast_sim(), getData(), get_GI_moments(), get_Gbar(), get_R0(), get_doses_per_day(), get_evec(), get_kernel_moments(), get_opt_pars(), get_r(), invlink_trans(), make_betavec(), make_beta(), make_jac(), make_ratemat(), make_state(), make_test_wtsvec(), make_vaxrate(), mk_Nvec(), mk_agecats(), mk_contact_rate_setting(), mk_mistry_Nvec(), mk_pmat(), mk_vaxcats(), mle_fun(), non_expanded_states, rExp(), read_params(), repair_names_age(), restore(), run_sim_ageify(), run_sim_break(), run_sim_loglin(), run_sim_mobility(), run_sim_range(), run_sim(), show_ratemat(), testify(), texify(), trans_state_vars(), update_contact_rate_setting(), update_foi(), update_params_mistry(), vis_model(), write_params()