hubbell: Plays a Hubbell game.

In hubbell: Functions for Hubbell's Grand Unified Theory

Description

Hubbell game is a zero-sum game, where a number of individuals is killed during a timestep, and these are replaced either from within the remaining community or from the surrounding metacommunity. In either case, the replacement probabilities are given by species frequencies.

Usage

hubbell(comm, D, m=0, P=NULL)
hubb2vect(vec)
## S3 method for class 'hubbell'
plot(x, sp.max=64, ...)

Arguments

comm, x	A hubbell object giving the counts of each species in a community.
D	Number of individuals killed during time step.
m	Probability of immigration.
P	Relative frequencies of species in the metacommunity.
vec	A vector giving counts of each species.
sp.max	Species number parameter for plots.
...	Other parameters passed to functions.

Details

Hubbell game is zero-sum game for a local community. In each time step, D individuals are killed. These may be replaced either by the remaining species in the local community, or if m > 0, they may be replaced by the surrounding metacommunity. Metacommunity frequencies P must be given by the user – a Hubbellian way is to use rhubbell.

The community comm is an object of class hubbell. This is simply a vector giving the species id. number for each individual in the community. Function hubb2vect can be used to produce the initial community. Function plot tries to displays the canopy using colour codes for species, together with a ranked abundance – diversity (RAD) curve. The parameter sp.max determines the number of colours in plot and the scaling of horizontal axis in RAD plot.

Value

A hubbell object.

Author(s)

Jari Oksanen

References

Hubbell, S.P. (2001). The Unified Neutral Theory of Biodiversity and Biogeography. Princeton Univ. Press.

See Also

radfit for alternative visualization and analysis.

Examples

## Island Biogeography: completely isolated community 
## starting with 40 species, each with 30 individuals (J=1200)
## and going for 100 time steps. Kill 240 individuals each year.
comm <- hubb2vect(rep(40,30))
comm
plot(comm, sp.max=30)
comm <- hubbell(comm, 240)
plot(comm, sp.max=30)
## Look for 100 time steps ahead
for (t in 1:100) comm <- hubbell(comm, D=240)
plot(comm, sp.max=30)
## Make a Hubbell metacommunity, and start immigration
metacomm <- rhubbell(20, 16000)
for (t in 1:100) comm <- hubbell(comm, D=240, m=0.1, P=metacomm)
comm  ## More species now
plot(comm)

hubbell documentation built on May 2, 2019, 5:40 p.m.