simTs1D: Simulate a model on a regular grid of times, using a function...

In smfsb: Stochastic Modelling for Systems Biology

Simulate a model on a regular grid of times, using a function (closure) for advancing the state of the model

Description

This function simulates single realisation of a model on a 1D regular spatial grid and regular grid of times using a function (closure) for advancing the state of the model, such as created by StepGillespie1D.

Usage

simTs1D(x0,t0=0,tt=100,dt=0.1,stepFun,verb=FALSE,...)

Arguments

x0	The initial state of the process at time t0, a matrix with rows corresponding to reacting species and columns corresponding to spatial location.
t0	The initial time to be associated with the initial state x0. Defaults to 0.
tt	The terminal time of the simulation.
dt	The time step of the output. Note that this time step relates only to the recorded output, and has no bearing on the accuracy of the simulation process.
stepFun	A function (closure) for advancing the state of the process, such as produced by StepGillespie1D.
verb	Output progress to the console (this function can be very slow).
...	Additional arguments will be passed to stepFun.

Value

An R 3d array representing the simulated process. The dimensions are species, space, and time.

See Also

StepGillespie1D, simTs

Examples

data(spnModels)
N=20; T=30
x0=matrix(0,nrow=2,ncol=N)
rownames(x0)=c("x1","x2")
x0[,round(N/2)]=LV$M
stepLV1D = StepGillespie1D(LV,c(0.6,0.6))
xx = simTs1D(x0,0,T,0.2,stepLV1D,verb=TRUE)
op=par(mfrow=c(1,2))
image(xx[1,,],main="Prey",xlab="Space",ylab="Time")
image(xx[2,,],main="Predator",xlab="Space",ylab="Time")
par(op)

smfsb documentation built on Jan. 13, 2024, 3:02 a.m.