rhubbell: Simulates a community from Hubbell's species generator.

In hubbell: Functions for Hubbell's Grand Unified Theory

Description

Simulates a community using Hubbell's species generator.

Usage

rhubbell(theta, J)

Arguments

theta	Hubbell's ultimate diversity number θ.
J	Size of generated community.

Details

Hubbell's species generator P = θ/(θ + j -1) gives the probability that jth collected individual belongs to a new species. If, by chance, this happens, a new species is added to the community with abundance 1. If the jth individual belongs to an old species, abundance 1 is added to some of the old species. The probability of getting this new indivdiual is proportional to species abundances in community at the moment.

Value

Returns a non-ranked vector of species abundances (integers).

Note

The R code is crystal clear – at least much clearer than Hubbell's description (p. 290–291, incomprehensible spaghetti diagram in Fig. 9.1). Use the force, read the source.

Author(s)

Jari Oksanen

References

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

See Also

radfit for displaying results, rbstick for an alternative model (but giving very similar results).

Examples

## A nemoral forest
comm <- rhubbell(20, 200)
require(vegan)
plot(radfit(comm))

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