Last data update: 2014.03.03

R: Metapopulation Competition Model
metaCompR 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, fi1, fi2, i1, i2, pe, D=0, anima=TRUE)

Arguments

tmax

maximum simulation time.

cl

number of columns for the simulated landscape.

rw

number of rowns for the simulated landscape.

fi1

initial fraction of patches occupied by the superior competitor.

fi2

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,fi1=0.1,fi2=0.4,i1=0.4,i2=0.5,pe=0.25)
metaComp(tmax=100,cl=20,rw=20,fi1=0.1,fi2=0.4,i1=0.4,i2=0.5,pe=0.25, D=0.1)

## End(Not run)

Results