rand.neutral {untb} | R Documentation |
Random neutral ecosystem
Description
Given the size of the metacommunity J
, and the fundamental
biodiversity number \theta
, generate an object of class
count
using a stochastic mechanism consistent with the
neutral theory.
Usage
rand.neutral(J, theta=NULL, prob.of.mutate=NULL, string = NULL, pad = FALSE)
Arguments
J |
Size of metacommunity |
theta |
Fundamental biodiversity number |
prob.of.mutate |
Probability of mutation |
string |
String to add to species names. By default (ie
|
pad |
Boolean, with default |
Details
Uses the simulation method on page 289 of Hubbell (2001).
Note
If pad
is TRUE
, and you set string
to
“extinct
”, things will break.
Author(s)
Robin K. S. Hankin
References
S. P. Hubbell 2001. “The Unified Neutral Theory of Biodiversity”. Princeton University Press.
See Also
Examples
rand.neutral(1000, 9)
rand.neutral(1000, 9, string="spp.")
data(butterflies)
rand.neutral(no.of.ind(butterflies), optimal.theta(butterflies),string="spp.")
# what is the distribution of abundance of the second ranked species if
# J=10, theta=0.7?
plot(table(replicate(100,rand.neutral(10,theta=0.7,pad=TRUE)[2])))