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)