run: Runs 'SimBIID_model' object

In tjmckinley/SimBIID: Simulation-Based Inference Methods for Infectious Disease Models

Runs SimBIID_model object

Description

Wrapper function that compiles (if necessary) and runs a SimBIID_model object. Returns results in a user-friendly manner as a SimBIID_run object, for which print() and plot() generics are provided.

Usage

run(
  model,
  pars,
  tstart,
  tstop,
  u,
  tspan,
  nrep = 1,
  parallel = FALSE,
  mc.cores = NA
)

Arguments

model	An object of class SimBIID_model.
pars	A named vector of parameters.
tstart	The time at which to start the simulation.
tstop	The time at which to stop the simulation.
u	A named vector of initial states.
tspan	A numeric vector containing the times at which to save the states of the system.
nrep	Specifies the number of simulations to run.
parallel	A logical determining whether to use parallel processing or not.
mc.cores	Number of cores to use if using parallel processing.

Value

An object of class SimBIID_run, essentially a list containing elements:

• sums: a data.frame() with summaries of the model runs. This includes columns run, completed, t, u* (see help file for SimBIID_model for more details);

• runs: a data.frame() object, containing columns: run, t, u* (see help file for SimBIID_model for more details). These contain time series counts for the simulations. Note that this will only be returned if tspan = TRUE in the original SimBIID_model object.

• bootEnd: a time point denoting when bootstrapped estimates end and predictions begin (for predict.PMCMC() method).

See Also

mparseRcpp, print.SimBIID_runs, plot.SimBIID_runs

Examples

## set up SIR simulation model
transitions <- c(
    "S -> beta * S * I -> I", 
    "I -> gamma * I -> R"
)
compartments <- c("S", "I", "R")
pars <- c("beta", "gamma")
model <- mparseRcpp(
    transitions = transitions, 
    compartments = compartments,
    pars = pars
)

## compile and run model
sims <- run(
    model = model,
    pars = c(beta = 0.001, gamma = 0.1),
    tstart = 0,
    tstop = 100,
    u = c(S = 119, I = 1, R = 0)
)
sims

## add tspan option to return
## time series counts at different
## time points
model <- mparseRcpp(
    transitions = transitions, 
    compartments = compartments,
    pars = pars,
    tspan = TRUE
)
sims <- run(
    model = model,
    pars = c(beta = 0.001, gamma = 0.1),
    tstart = 0,
    tstop = 100,
    u = c(S = 119, I = 1, R = 0),
    tspan = seq(1, 100, length.out = 10)
)
sims

## run 100 replicate simulations and
## plot outputs
sims <- run(
    model = model,
    pars = c(beta = 0.001, gamma = 0.1),
    tstart = 0,
    tstop = 100,
    u = c(S = 119, I = 1, R = 0),
    tspan = seq(1, 100, length.out = 10),
    nrep = 100
)
sims
plot(sims)

tjmckinley/SimBIID documentation built on Sept. 11, 2022, 11:58 a.m.