Simulate the model using the Gillespie direct method algorithm.
RunSimulation(model, theta, tMax = 100, numObservations = 5)
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. |
# NOT RUN { RunSimulation(model, theta) RunSimulation(model, c(0.003, 0.1), 100.0, 5) # }