diagnostics.deSolve {deSolve} | R Documentation |
Prints several diagnostics of the simulation to the screen, e.g. number of steps taken, the last step size, ...
## S3 method for class 'deSolve' diagnostics(obj, Full = FALSE, ...)
obj |
is the output matrix as produced by one of the integration routines. |
Full |
when |
... |
optional arguments allowing to extend |
When the integration output is saved as a data.frame
, then the required
attributes are lost and method diagnostics
will not work anymore.
The integer and real vector with diagnostic values; for function lsodar
also the root information.
See tables 2 and 3 in vignette("deSolve") for what these vectors contain.
Note: the number of function evaluations are *without* the extra calls performed to generate the ordinary output variables (if present).
## The famous Lorenz equations: chaos in the earth's atmosphere ## Lorenz 1963. J. Atmos. Sci. 20, 130-141. chaos <- function(t, state, parameters) { with(as.list(c(state)), { dx <- -8/3 * x + y * z dy <- -10 * (y - z) dz <- -x * y + 28 * y - z list(c(dx, dy, dz)) }) } state <- c(x = 1, y = 1, z = 1) times <- seq(0, 50, 0.01) out <- vode(state, times, chaos, 0) pairs(out, pch = ".") diagnostics(out)