Simulate the model using the Gillespie direct method algorithm.

RunSimulation(model, theta, tMax = 100, numObservations = 5)

Arguments

model

initialised CompModel object.

theta

vector of model parameters.

tMax

the maximum time that events will be simulated for.

numObservations

the number of observations to be drawn.

Examples

# NOT RUN {
RunSimulation(model, theta)
RunSimulation(model, c(0.003, 0.1), 100.0, 5)
# }