Get the contact rate fraction threshold for increases

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