metaComp {EcoVirtual} | R Documentation |
Metapopulation Competition Model
Description
Simulate a metapopulation dynamics with two competing species, a superior and an inferior competitor. Includes the possibility of habitat destruction in the model.
Usage
metaComp(tmax, rw, cl, f01, f02, i1, i2, pe, D = 0, anima = TRUE)
Arguments
tmax |
maximum simulation time. |
rw |
number of rowns for the simulated landscape. |
cl |
number of columns for the simulated landscape. |
f01 |
initial fraction of patches occupied by the superior competitor. |
f02 |
initial fraction of patches occupied by the inferior competitor. |
i1 |
colonization coefficient for the superior competitor. |
i2 |
colonization coefficient for the inferior competitor. |
pe |
probability of extinction (equal for both species). |
D |
proportion of habitat destroyed. |
anima |
logical; if TRUE, show simulation frames. |
Details
This function uses the metapopulationa model with internal colonization
(see function metaCi in metapopulation
) for the superior
competitor. The inferior competitor can only occupy empty patches and is
displaced by the superior competitor if it occupies the same patch.
The argument 'D' inserts the influences of habitat destruction in the model.
The number of patches in the simulated landscape is defined by rw*cl.
Value
'metaComp' returns a graphic with the simulated landscapes and the results of the proportion of patch occupied by both species.
This function also return an invisible array with the simulation results.
Author(s)
Alexandre Adalardo de Oliveira and Paulo Inacio Prado ecovirtualpackage@gmail.com
References
Stevens, M.H.H. 2009. A primer in ecology with R. New York, Springer.
Gotelli, N.J. 1991. Metapopulation models: the rescue effect, the propagule rain, and the core-satellite hypothesis. The American Naturalist 138:768-776.
See Also
comCompete
, http://ecovirtual.ib.usp.br
Examples
## Not run:
metaComp(tmax=100,cl=20,rw=20,f01=0.1,f02=0.4,i1=0.4,i2=0.5,pe=0.25)
metaComp(tmax=100,cl=20,rw=20,f01=0.1,f02=0.4,i1=0.4,i2=0.5,pe=0.25, D=0.1)
## End(Not run)