readsdr (version 0.3.0)

sd_data_generator_fun: Function factory for SBC

Description

Function factory for SBC

Usage

sd_data_generator_fun(
  filepath,
  estimated_params,
  meas_mdl,
  start_time = NULL,
  stop_time = NULL,
  timestep = NULL,
  integ_method = "euler"
)

Value

A function.

Arguments

filepath

A string that indicates a path to a file with extension .stmx or .xmile. Vensim files (.mdl) are not xmile files. They must be exported from Vensim with extension .xmile

estimated_params

A list of lists. Each sublist describes each parameter that will be estimated in the inference stage. To construct this description, the user can avail of the function `sd_prior`.

meas_mdl

A list of strings. Each string corresponds to a sampling statement written in Stan language.

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")
  meas_mdl <- list("y ~ poisson(net_flow(C))")
  estimated_params <- list(
    sd_prior("par_beta", "lognormal", c(0, 1)),
    sd_prior("par_rho", "beta", c(2, 2)),
    sd_prior("I0", "lognormal", c(0, 1), "init"))
  sd_data_generator_fun(filepath, estimated_params, meas_mdl)