get_doubling_time.Rd
Get the doubling time
get_doubling_time( obj, iter = seq_along(obj$post$R0), forecast_days = 25, window_check = 25, show_plot = TRUE, ... )
obj | Output from |
---|---|
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 |
A vector of doubling times across iterations.
Assumes that prevalence is increasing; negative values imply halving
times. Use show_plot
if you want to check the time series.
# See ?project_seir