EcoVirtual (version 1.1)

metaPop: 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.

Value

Metapopulation functions return graphics with the simulation results. These functions also return an invisible array with the simulation data.

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.

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

http://ecovirtual.ib.usp.br

Examples

# NOT RUN {
# }
# 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)
# }
# NOT RUN {
# }