plot.deSolve {deSolve} | R Documentation |
Plot the output of numeric integration routines.
## S3 method for class 'deSolve': plot(x, which = 1:(ncol(x)-1), ask = NULL, ...) ## S3 method for class 'deSolve': hist(x, which = 1:(ncol(x)-1), ask = NULL, ...)
x |
an object of class deSolve , as returned by the
integrators, and to be plotted. |
which |
the name(s) or the index to the variables that should be plotted. Default = all variables. |
ask |
logical; if TRUE , the user is asked before
each plot, if NULL the user is only asked if more than one
page of plots is necessary and the current graphics device is set
interactive, see par(ask=.) and
dev.interactive . |
... |
additional graphics arguments passed to
plot.default or hist |
The number of panels per page is automatically determined up to 3 x 3
(par(mfrow=c(3, 3))
). This default can be overwritten by
specifying user-defined settings for mfrow
or mfcol
.
Other graphical parameters can be passed as well. Parameters
xlab
and ylab
are vectorized, so it is possible to
assign specific axis labels to individual plots.
## A Predator-Prey model with 4 species in matrix formulation LVmatrix <- function(t, n, parms) { with(parms, { dn <- r * n + n * (A %*% n) return(list(c(dn))) }) } parms <- list( r = c(r1 = 0.1, r2 = 0.1, r3 = -0.1, r4 = -0.1), A = matrix(c(0.0, 0.0, -0.2, 0.01, # prey 1 0.0, 0.0, 0.02, -0.1, # prey 2 0.2, 0.02, 0.0, 0.0, # predator 1; prefers prey 1 0.01, 0.1, 0.0, 0.0), # predator 2; prefers prey 2 nrow = 4, ncol = 4, byrow=TRUE) ) times <- seq(from = 0, to = 500, by = 0.1) y <- c(prey1 = 1, prey2 = 1, pred1 = 2, pred2 = 2) out <- ode(y, times, LVmatrix, parms) ## Basic line plot plot(out, type = "l") ## User-specified axis labels plot(out, type = "l", ylab = c("Prey 1", "Prey 2", "Pred 1", "Pred 2"), xlab = "Time (d)", main = "Time Series") hist(out, col="darkblue", breaks = 50)