rbstick: Generates a Random Brokenstick Community

In hubbell: Functions for Hubbell's Grand Unified Theory

Description

Generates a random brokenstick community.

Usage

rbstick(n)

Arguments

n	Number of species.

Details

Brokenstick has been a popular null-model for abundance relations since the days of MacArthur. Brokenstick community is generated so that species simultaneously divide the available community (or resources) in random pieces. When these pieces are arranged in rank abundance order, the result looks very similar to log-Normal.

The current functin returns a random community. Expected brokenstick communities can be generated using vegan function bstick.

Value

Returns the ranked abundance proportions (sum = 1).

Author(s)

Jari Oksanen

References

Tokeshi, M. (1993). Species abundance patterns and community structure. Adv. Ecol. Res. 24, 111-186.

Wilson, J.B. (1993). Would we recognize a Broken-Stick community if we found one? Oikos 67: 181-183.

See Also

radfit for displaying results and fitting a brokenstick model, and bstick for the expected brokenstick distribution. Function rhubbell gives a Hubbellian alternative.

Examples

## Generate and plot a broken-stick community
comm <- rbstick(30)
require(vegan)
(mod <- radfit(comm, family=gaussian))
plot(mod) 

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