smfsb (version 1.5)

discretise: Discretise output from a discrete event simulation algorithm

Description

This function discretise output from a discrete event simulation algorithm such as gillespie onto a regular time grid, and returns the results as an R ts object.

Usage

discretise(out, dt=1, start=0)

Value

An R ts object containing the discretised output.

Arguments

out

A list containing discrete event simulation output in the form of that produced by gillespie.

dt

The time step required for the output of the discretisation process. Defaults to one time unit.

start

The start time for the output. Defaults to zero.

See Also

simpleEuler, rdiff, gillespie, gillespied, ts

Examples

# load LV model
data(spnModels)
# simulate a realisation of the process and plot it
out = gillespie(LV,10000)
op=par(mfrow=c(2,2))
plot(stepfun(out$t,out$x[,1]),pch="")
plot(stepfun(out$t,out$x[,2]),pch="")
plot(out$x,type="l")

# use the "discretise" function to map it to an R "ts" object
plot(discretise(out,dt=0.01),plot.type="single",lty=c(1,2))
par(op)