Simulate vivax model with delay to treatment and MDA, with stochastic model

Simulates a draw of the chosen vivax stochastic model

simulate_vivax_delay_mda_sto(
  parameters,
  STOmodel = model_sto_vivax_delay(),
  STOmodel_mda = model_sto_vivax_delay_mda(),
  maxtime = 1465,
  year = FALSE,
  sto_method = "exact",
  runs = 1,
  seeds = NULL
)

Arguments

parameters

model parameters

STOmodel

stochastic model to be simulated (reactions)

STOmodel_mda

stochastic model with MDA to be simulated (during prophylaxis)

maxtime

maximal time step

year

if TRUE, aggregates outputs per year. if FALSE, returns daily output

sto_method

a scalarindicating which simulation method is used. Default ("exact") is Gillespie's direct method. Other options are "approximate" (tau-leap) or "mixed". cf. the documentation of the TiPS package for more information.

runs

number of draws of the stochastic model

seeds

a vector of the length of runs containing the seeds for each simulation (don't use "0" which has another use in TiPS)

Value

A dataframe containing the simulated model