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) # }