ecol.death.sim {animation} | R Documentation |
A simulation of the death of two species with certain probabilities.
ecol.death.sim(nr = 10, nc = 10, num.sp = c(50, 50), col.sp = c(1, 2), pch.sp = c(1, 2), col.die = 1, pch.die = 4, cex = 3, ...)
nr, nc |
number of rows and columns of the field (plants grow on a nr x nc grid)
|
num.sp |
number of two plants respectively |
col.sp,pch.sp |
colors and point symbols of the two species respectively |
col.die, pch.die, cex |
the color, point symbol and magnification to annotate the plant which dies (symbol default to be an `X') |
... |
other arguments passed to plot to set up the plot
|
Suppose there are two plant species in a field: A and B. One of them will die at each time and a new plant will grow in the place where the old plant died; the species of the new plant depends on the proportions of two species: the larger the proportion is, the greater the probability for this species to come up will be.
a vector (factor) containing 1's and 2's, denoting the plants finally survived
2 * nmax
image frames will actually be produced, and the option nmax
will be adjusted to 2 * nmax
in the end.
Yihui Xie <http://yihui.name>
This animation is motivated by a question raised from Jing Jiao, a student in biology, to show the evolution of two species.
The original post is in the forum of the ``Capital of Statistics'': http://cos.name/bbs/read.php?tid=14093 (in Chinese)
oopt = ani.options(nmax = 50, interval = 0.3) par(ann = FALSE, mar = rep(0, 4)) ecol.death.sim() ## Not run: ## large scale simulation ani.options(nmax = 1000, interval = 0.02) ecol.death.sim(col.sp = c(8, 2), pch.sp = c(20, 17)) ## End(Not run) ani.options(oopt)