Skip to contents

Convenience wrapper: builds a forecast over the observed data times including "cases" and overlays observations.

Usage

chlaa_plot_ppc(
  fit,
  pars = NULL,
  data = NULL,
  n_draws = 200,
  burnin = 0.5,
  thin = 1,
  seed = 1,
  dt = 0.25,
  n_particles = 1,
  obs_model = c("nbinom", "mean")
)

Arguments

fit

A `chlaa_fit` object.

pars

Baseline parameter list. If NULL uses `attr(fit, "start_pars")`.

data

Observed data. If NULL uses `attr(fit, "data")`.

n_draws

Number of posterior draws.

burnin

Burn-in for posterior draws.

thin

Thinning for posterior draws.

seed

Seed.

dt

Model time step.

n_particles

Particles per draw.

obs_model

Observation model noise: "nbinom" or "mean".

Value

A ggplot object.