Description Usage Arguments Examples
This function creates density plots for the simulation of two-dimensional stochastic differential equations from TSTraj
1 2 3 4 5 6 7 8 9 10 11 12 13 14 |
mat |
a matrix output from |
dim |
dimensions of the plot; |
xlim |
numeric vectors of length 2, giving the x coordinate range. Default |
ylim |
numeric vectors of length 2, giving the y coordinate range. Default |
contour.levels |
the number of contour levels for the two-dimensional plots (i.e., when |
col2d |
vector of colors to be used in the plot. |
contour.lwd |
line width of contour lines if |
contour.lines |
if |
kde2d.n |
number of grid points in each direction. Can be scalar or a length-2 integer vector. Passes to argument |
xlab |
label for x axis when dim = 2. Default is 'X' |
ylab |
label for y axis when dim = 2. Default is 'Y' |
... |
passes arguments to |
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 | ## Not run:
# First, the parameter values, as found in 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 TSTraj
#Example 1 from article
equationx <- "1.54*x*(1.0-(x/10.14)) - (y*x*x)/(1.0 + x*x)"
equationy <- "((0.476*x*x*y)/(1 + x*x)) - 0.112590*y*y"
# Third, run it, as found in TSTraj
ModelOut <- TSTraj(y0 = model.state, time = model.time, deltat = model.deltat,
x.rhs = equationx, y.rhs = equationy, sigma = model.sigma)
# Fourth, plot it:
# in 1D
TSDensity(ModelOut, dim = 1)
# in 2D
TSDensity(ModelOut, dim = 2, kde2d.n = 20, xlab = "")
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.