readsdr (version 0.3.0)

sd_measurements: Generate measurements

Description

Generate measurements

Usage

sd_measurements(
  n_meas,
  meas_model,
  ds_inputs,
  start_time = NULL,
  stop_time = NULL,
  timestep = NULL,
  integ_method = "euler"
)

Value

A data frame.

Arguments

n_meas

Number of measurements. An integer.

meas_model

Measurement model. A list of strings, in which each string corresponds to sampling statement in Stan language.

ds_inputs

A list of deSolve inputs generated by read_xmile

start_time

A number indicating the time at which the simulation begins.

stop_time

A number indicating the time at which the simulation ends.

timestep

A number indicating the time interval for the simulation. Also known as dt.

integ_method

A string indicating the integration method. It can be either "euler" or "rk4"

Examples

  filepath <- system.file("models/", "SEIR.stmx", package = "readsdr")
  mdl      <- read_xmile(filepath)

  mm1        <- "y ~ poisson(C)"
  meas_model <- list(mm1)

  sd_measurements(n_meas       = 2,
                  meas_model   = meas_model,
                  ds_inputs    = mdl$deSolve_components,
                  start_time   = 0,
                  stop_time    = 10,
                  timestep     = 1/16,
                  integ_method = "rk4")