Simulate stochastic model
simstoch.RdSimulate a stochastic model generated with
uqsa::generateGillespieModel(), using the solver in this package.
Arguments
- experiments
list of expperiments, same as for the deterministic solvers.
- model.so
compiled C code for the model, this has to be a path with at least one slash in it, e.g.: ./model.so
- parameters
a numeric vector of appropriate size
Details
This will simulate all experimental conditions included in the list of experiments, including applying the inputs:
u <- experiments[[i]]$input - the input will be copied to the end of the model's internal parameter vector.
Like for deterministic models, we assume that there is a vector of
unknown parameter (a Markov chain variable, a vector of
optimization variables) and also known parameters (aka the input
parameters). The model itself does not distinguish between the two,
but one is the same between the experiments and one is different
between different experiments: modelParam <- c(mcmcParam, inputParam)