function for species richness estimation

Description

This function implements the species richness estimation procedure outlined in Chao and Bunge (2002).

Usage

1
chao_bunge(data, cutoff=10, print=TRUE, answers=FALSE)

Arguments

data

The sample frequency count table for the population of interest. See breakaway for details.

cutoff

The maximum frequency count to use for prediction. Defaults to 10.

print

Logical: whether the results should be printed to screen.

answers

Logical: whether the function should return an argument.

Value

name

The “name” of the model.

est

The estimate of total (observed plus unobserved) diversity.

seest

The standard error in the diversity estimate.

ci

An asymmetric 95% confidence interval for diversity.

Note

This estimator is based on the negative binomial model and for that reason generally produces poor fits to microbial data. The result is artificially low standard errors. Caution is advised.

Author(s)

Amy Willis

References

Chao, A. and Bunge, J. (2002). Estimating the number of species in a stochastic abundance model. Biometrics, 58.

See Also

breakaway; apples; wlrm_untransformed; chao1

Examples

1