StepFRM {smfsb} | R Documentation |
This function creates a function for advancing the state of an SPN model using Gillespie's first reaction method. The resulting function (closure) can be used in conjunction with other functions (such as simTs
) for simulating realisations of SPN models.
StepFRM(N)
N |
An R list with named components representing a stochastic
Petri net. Should contain |
An R function which can be used to advance the state of the SPN model N
by using Gillespie's first reaction method. The function closure has interface function(x0,t0,deltat,...)
, where x0
and t0
represent the initial state and time, and deltat
represents the amount of time by which the process should be advanced. The function closure returns a vector representing the simulated state of the system at the new time.
StepEulerSPN
, StepGillespie
, simTs
, simSample
# load the LV model data(spnModels) # create a stepping function stepLV = StepFRM(LV) # step the function print(stepLV(c(x1=50,x2=100),0,1)) # simulate a realisation of the process and plot it out = simTs(c(x1=50,x2=100),0,100,0.1,stepLV) plot(out,plot.type="single",lty=1:2)