Skip to contents

Estimate Rt via a random walk

Source: R/model_infections.R
R_estimate_rw.Rd

This option estimates the effective reproduction number over time using a random walk.

Usage

R_estimate_rw(
  intercept_prior_mu = 1,
  intercept_prior_sigma = 0.8,
  sd_prior_mu = 0,
  sd_prior_sigma = 0.1,
  link = "inv_softplus",
  R_max = 6,
  differenced = FALSE,
  noncentered = TRUE,
  modeldata = modeldata_init()
)

Arguments

intercept_prior_mu

Prior (mean) on the intercept of the random walk.

intercept_prior_sigma

Prior (standard deviation) on the intercept of the random walk.

sd_prior_mu

Prior (mean) on the standard deviation of the random walk.

sd_prior_sigma

Prior (standard deviation) on the standard deviation of the random walk.

Link function. Currently supported are inv_softplus (default) and scaled_logit. Both of these links are configured to behave approximately like the identity function around R=1, but become increasingly non-linear below (and in the case of scaled_logit also above) R=1.

R_max

If link=scaled_logit is used, a maximum reproduction number must be assumed. This should be higher than any realistic R value for the modeled pathogen. Default is 6.

differenced

If FALSE (default), the random walk is applied to the absolute Rt time series. If TRUE, it is instead applied to the differenced time series, i.e. now the trend is modeled as a random walk.

noncentered

If TRUE (default), a non-centered parameterization is used to model the innovations of the random walk (for better sampling efficiency).

modeldata

A modeldata object to which the above model specifications should be added. Default is an empty model given by modeldata_init(). Can also be an already partly specified model returned by other EpiSewer modeling functions.

Value

A modeldata object containing data and specifications of the model to be fitted. Can be passed on to other EpiSewer modeling functions to add further data and model specifications.

The modeldata object also includes information about parameter initialization (init), meta data (.metainfo), and checks to be performed before model fitting (.checks).

Details

The smoothness of Rt estimates is influenced by the prior on the standard deviation of the random walk. It also influences the uncertainty of Rt estimates towards the present / date of estimation, when limited data signal is available. The prior on the intercept of the random walk should reflect your expectation of Rt at the beginning of the time series. If estimating from the start of an epidemic, you might want to use a prior with mean > 1 for the intercept.

The priors of this component have the following functional form:

  • intercept of the random walk: Normal

  • standard deviation of the random walk: Truncated normal

See also

Other Rt models: R_estimate_approx(), R_estimate_ets(), R_estimate_splines()