| metaPop {EcoVirtual} | R Documentation |
Metapopulation Models
Description
Simulate metapopulation dynamics with propagules seed rain, internal colonization and rescue effect.
Usage
metaPop(cl, rw, f0, pi, pe, tmax, anima = TRUE)
metaCi(cl, rw, f0, ci, pe, tmax, anima = TRUE)
metaEr(cl, rw, f0, pi, ce, tmax, anima = TRUE)
metaCiEr(cl, rw, f0, ci, ce, tmax, anima = TRUE)
Arguments
cl |
number of columns for the simulated landscape. |
rw |
number of rows for the simulated landscape. |
f0 |
initial proportion of occupied patches. |
pi |
probability of colonization. |
pe |
probability of extinction. |
tmax |
maximum simulation time. |
anima |
show animation frames. |
ci |
colonization coefficient, represents the maximum probability of colonization (when f=1) and should be a number between 0 and 1. |
ce |
coefficient of extinction, represents the maximum probability of extinction (when f=0) and should be a number between 0 and 1. |
Details
'metaPop' is the seed rain metapopulation model, including only propagules seed rain from a external pool (no extinction).
'metaCi' is the Internal Colonization model, where number of propagules depends on number of occupied patches, there is no external pool.
'metaEr' is the Rescue Effect model, where extinction probability is negatively associated with number of occupied patches.
'metaCiEr' includes both effects: Rescue Effect and Internal Colonization.
The number of patches in the simulated landscape is defined by rw*cl.
Value
Metapopulation functions return graphics with the simulation results. These functions also return an invisible array with the simulation data.
Author(s)
Alexandre Adalardo de Oliveira and Paulo Inacio Prado ecovirtualpackage@gmail.com
References
Gotelli, N.J. 1991. Metapopulation models: the rescue effect, the propagule rain, and the core-satellite hypothesis. The American Naturalist 138:768-776.
Gotelli, N.J. 2008. A primer of Ecology. 4th ed. Sinauer Associates, 291pp.
See Also
Examples
## Not run:
metaPop(cl=10,rw=10,f0=0.5,pi=0.3,pe=0.15, tmax=100)
metaCi(cl=10,rw=10,f0=0.1,ci=1,pe=0.5, tmax=100)
metaEr(cl=10, rw=10, f0=0.2, pi=0.2, ce=0.15, tmax=100)
metaCiEr(cl=10, rw=10, f0=0.2, ci=0.2, ce=0.15, tmax=100)
## End(Not run)