get_threshold: Get the contact rate fraction threshold for increases

In seananderson/covidseir: Bayesian SEIR Modelling for Multivariate COVID-19 Case Data

Description

Get the contact rate fraction threshold for increases

Usage

get_threshold(
  obj,
  iter = seq_along(obj$post$R0),
  forecast_days = 25,
  fs = seq(0.3, 0.8, length.out = 4),
  show_plot = TRUE,
  window_check = 25,
  ...
)

Arguments

obj	Output from fit_seir().
iter	Vector of MCMC iterations to work with. More iterations will be slower and not necessarily render a different result. ~200 total iterations may be plenty.
forecast_days	Days to use in forecast.
fs	Contact rate fractions to test.
show_plot	Make a diagnostic plot?
window_check	The window of days to use from the last day forecasted.
...	Other arguments for project_seir().

Value

The threshold value.

Examples

# See ?project_seir

seananderson/covidseir documentation built on June 3, 2021, 6:49 a.m.