diffusion {simecol} | R Documentation |
simecol example: This is a random walk (basic particle diffusion) model.
data(diffusion)
An S4 object according to the rwalkModel
specification.
The object contains the following slots:
main
parms
ninds
speed
area
times
by
-argument
ignored).init
x
and y
and movement angle a
) of the
particles.To see all details, please have a look into the implementation.
## Not run: data(diffusion) ## (1) basic example plot(sim(diffusion)) ## show "grid of environmental conditions" image(inputs(diffusion)) ## (2) scenario ## with homogeneous environment (no "refuge" in the middle) no_refuge <- diffusion # Cloning of the whole model object inputs(no_refuge) <- matrix(1, 100, 100) plot(sim(no_refuge)) ## (3) example of using an observer function observer(diffusion) <- function(state, ...) { ## numerical output to the screen cat("mean x=", mean(state$x), ", mean y=", mean(state$y), ", sd x=", sd(state$x), ", sd y=", sd(state$y), "\n") # animation par(mfrow=c(2,2)) plot(state$x, state$y, xlab="x", ylab="y", pch=16, col="red", xlim=c(0, 100)) hist(state$y) hist(state$x) # default case: return the state --> iteration stores it in "out" init } sim(diffusion) ## (4) remove the observer observer(diffusion) <- NULL diffusion <- sim(diffusion) ## End(Not run)