EcoVirtual (version 1.1)

archip: Species Colonization and Species-Area Relationship in Archipelagos

Description

Simulate species colonization from mainland to islands with different sizes.

Usage

archip(n.isl, ar.min, ar.max, S, seed.rain, abund, tmax = 100, anima = TRUE)

Arguments

n.isl

numeric, number of islands.

ar.min

numeric, area of the smallest island.

ar.max

numeric, area of the biggest island.

S

numeric, number of species (species richness from mainland).

seed.rain

numeric, seed rain. Number of seeds colonizing islands on each time.

abund

numeric, abundance of each species in the seed rain.

tmax

numeric, maximum time for the simulations.

anima

logical; if TRUE, show simulation frames.

Value

'archip' returns 3 graphics:

• The species-area relationship: number of species x island area at the end of the simulation. It also returns the coefficients c and z from species-area relationship \(S=cA^z\).

• Colonization rate curves: colonization (number of species per cycle) x number of species for each island.

• Passive colonization: number of species x time for each island.

  'archip' also returns an invisible array with the simulation results.

Details

The mainland has richness (S) and the evenness can be controled argument abund. The 'abund' argument can be one of these 3 options:

1. a vector with the same length of the species richness, meaning the proportion of each species population;

2. a single value more than 1 representing equal abundance of each species (maximum evenness);

3. a single value between 0 and 1, meaning the model of geometric species rank-abundance distribution. The model is: abund*(1-abund)*((1:S)-1), where S is the number of species.

References

Gotelli, N.J. 2008. A primer of Ecology. 4th ed. Sinauer Associates, 291pp.

See Also

animaColExt, bioGeoIsl, http://ecovirtual.ib.usp.br

Examples

# NOT RUN {

# }
# NOT RUN {
archip(n.isl=10,ar.min=10, ar.max=100, S=1000, seed.rain=100, abund=10, tmax=100, anima=TRUE)

archip(n.isl=10,ar.min=10, ar.max=100, S=1000, seed.rain=100, abund=0.5, tmax=100, anima=TRUE)
# }
# NOT RUN {

# }