Get the doubling time

get_doubling_time(
  obj,
  iter = seq_along(obj$post$R0),
  forecast_days = 25,
  window_check = 25,
  show_plot = TRUE,
  ...
)

Arguments

obj

Output from fit_seir().

iter

Vector of MCMC iterations to work with.

forecast_days

Days to use in forecast.

window_check

The window of days to use from the last day forecasted.

show_plot

Make a diagnostic plot? Each line is a posterior sample of the log prevalence time series. Each line should be well approximated by a linear model

...

Other arguments for project_seir().

Value

A vector of doubling times across iterations.

Details

Assumes that prevalence is increasing; negative values imply halving times. Use show_plot if you want to check the time series.

Examples

# See ?project_seir