QPot (version 1.0)

TSDensity: Density plot from simulation of two-dimensional stochastic differential equations

Description

This function allows you to create density plots for the simulation of two-dimensional stochastic differential equations in TSTraj

Usage

TSDensity(mat, dim = 1, xlim = "NULL", ylim = "NULL",
  contour.levels = 15, col2d = c("blue", "yellow", "orange", "red"),
  contour.lwd = 0.5, contour.lines = TRUE, ...)

Arguments

mat

a matrix output from TSTraj.

dim

dimensions of the plot; dim=1 plots simple density histogram or dim=2 plots the density in state space (i.e., X and Y respectively on the abscissa and ordinate axes).

xlim

numeric vectors of length 2, giving the x coordinate range. Default = 'NULL' automatically sizes plot window.

ylim

numeric vectors of length 2, giving the y coordinate range. Default = 'NULL' automatically sizes plot window.

contour.levels

the number of contour levels for the two-dimensional plots (i.e., when dim=2).

col2d

vector of colors to be used in the plot.

contour.lwd

line width of contour lines if contour.lines=TRUE.

contour.lines

if TRUE, then black countour lines added to the graph.

...

passes arguments to kde2d and plot.

Examples

# First, the parameter values, as found in \\code{\\link{TSTraj}}
	model.state <- c(x = 3, y = 3)
	model.sigma <- 0.2
	model.deltat <- 0.005
	model.time <- 100

# Second, write out the deterministic skeleton of the equations to be simulated,
# as found in \\code{\\link{TSTraj}}
	equationx <- "1.54*x*(1.0-(x/10.14)) - (y*x*x)/(1.0 + x*x)"

# Third, run it, as found in \\code{\\link{TSTraj}}