untb (version 1.0-19)

untb: Ecological drift simulation under the Unified Neutral Theory of Biodiversity

Description

Simulates ecological drift under the UNTB. Function untb() carries out the simulation; function select() carries out a single generational step.

Usage

untb(start, prob.of.mutate=0, D=1, gens=150, keep=FALSE, meta=NULL)
select(a, D=length(a), prob.of.mutate=0, meta=NULL)

Arguments

a, start

Starting ecosystem; coerced to class census. Usually, pass an object of class count; see examples. To start with a monoculture of size 10, use start=rep(1,10) and to use start=1:10.

prob.of.mutate

Probability of mutation

D

Number of organisms that die in each timestep

gens

Number of generations to simulate

keep

In function untb() Boolean with default FALSE meaning to return the system at the end of the simulation and TRUE meaning to return a matrix whose rows are the ecosystem at successive times.

meta

Metacommunity, as represented by a vector of integers (representing species; thus being a census). Can be a count or census object, in which case it is coerced to a census and thence to an integer vector.

References

S. P. Hubbell. The Unified Neutral Theory of Biodiversity. Princeton University Press, 2001.

Examples

data(butterflies)
untb(start=butterflies, prob=0, gens=100)

a <- untb(start=1:10,prob=0.005, gens=1000,keep=TRUE)
plot(species.count(a),type="b")
matplot(species.table(a),type="l",lty=1)

Run the code above in your browser using DataCamp Workspace