Learn R Programming

entropart (version 1.4-8)

Coverage: Sample coverage of a community

Description

Calculates an estimator of the sample coverage of a community described by its abundance vector.

Usage

Coverage(Ns, Estimator = "ZhangHuang", CheckArguments = TRUE)

Arguments

Ns

A numeric vector containing species abundances.

Estimator

A string containing one of the possible estimators: "ZhangHuang" (default), "Chao", "Turing".

CheckArguments

Logical; if TRUE, the function arguments are verified. Should be set to FALSE to save time when the arguments have been checked elsewhere.

Value

A named number equal to the calculated sample coverage. The name is that of the estimator used.

Details

The sample coverage \(C\) of a community is the total probability of occurence of the species observed in the sample. \(1-C\) is the probability for an individual of the whole community to belong to a species that has not been sampled.

The historical estimator is due to Turing (Good, 1953). It only relies on singletons (species observed only once). Chao's (Chao and Shen, 2010) estimator uses doubletons too and Zhang-Huang's (Chao et al., 1988; Zhang and Huang, 2007) uses the whole distribution.

References

Chao, A., Lee, S.-M. and Chen, T.-C. (1988). A generalized Good's nonparametric coverage estimator. Chinese Journal of Mathematics 16: 189-199.

Chao, A. and Shen, T.-J. (2010). Program SPADE: Species Prediction And Diversity Estimation. Program and user's guide. CARE, Hsin-Chu, Taiwan.

Good, I. J. (1953). On the Population Frequency of Species and the Estimation of Population Parameters. Biometrika 40(3/4): 237-264.

Zhang, Z. and Huang, H. (2007). Turing's formula revisited. Journal of Quantitative Linguistics 14(2-3): 222-241.

Examples

Run this code
# NOT RUN {
  # Load Paracou data (number of trees per species in two 1-ha plot of a tropical forest)
  data(Paracou618)
  # Ns is the vector of abundances of the metacommunity
  Ns <- Paracou618.MC$Ns
  # Calculate the sample coverage of the metacommunity
  Coverage(Ns)    # Stored in Paracou618.SampleCoverage
# }

Run the code above in your browser using DataLab