Description Usage Arguments Value See Also Examples

This function creates a function for advancing the state of an SPN model
using a simple Euler-Maruyama discretisation of the CLE on a 2D regular grid. The resulting function (closure) can be
used in conjunction with other functions (such as `simTs2D`

)
for simulating realisations of SPN models in space and time.

1 |

`N` |
An R list with named components representing a stochastic
Petri net (SPN). Should contain |

`d` |
A vector of diffusion coefficients - one coefficient for each reacting species, in order. The coefficient is the reaction rate for a reaction for a molecule moving into an adjacent compartment. The hazard for a given molecule leaving the compartment is therefore four times this value (as it can leave in one of 4 directions). |

`dt` |
Time step for the Euler-Maruyama discretisation. |

An R function which can be used to advance the state of the SPN model
`N`

by using a simple Euler-Maruyama algorithm. The function closure has
interface `function(x0,t0,deltat,...)`

, where `x0`

is a 3D array
with rows corresponding to species and columns corresponding to voxels,
representing the initial condition (with dimensions species, x, and y), `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 matrix representing the simulated state of the system at the new time.

`StepGillespie2D`

,`StepCLE`

,
`simTs1D`

, `StepCLE1D`

1 2 3 4 5 6 7 8 9 10 11 12 13 14 | ```
m=150
n=100
T=15
data(spnModels)
x0=array(0,c(2,m,n))
dimnames(x0)[[1]]=c("x1","x2")
x0[,round(m/2),round(n/2)]=LV$M
stepLV2D = StepCLE2D(LV,c(0.6,0.6),dt=0.05)
xx = simTs2D(x0,0,T,0.5,stepLV2D,verb=TRUE)
N = dim(xx)[4]
op=par(mfrow=c(1,2))
image(xx[1,,,N],main="Prey",xlab="Space",ylab="Time")
image(xx[2,,,N],main="Predator",xlab="Space",ylab="Time")
par(op)
``` |

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.