Calculates the number of species from probability vector. The name is that of the estimator (the bias correction) used.

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 | ```
Richness(NorP, Correction = "Chao1", Alpha = 0.05, JackOver = FALSE,
CheckArguments = TRUE, Ps = NULL, Ns = NULL)
bcRichness(Ns, Correction = "Chao1", Alpha = 0.05, JackOver = FALSE,
CheckArguments = TRUE)
## S3 method for class 'ProbaVector'
Richness(NorP, Correction = "Chao1", Alpha = 0.05, JackOver = FALSE,
CheckArguments = TRUE, Ps = NULL, Ns = NULL)
## S3 method for class 'AbdVector'
Richness(NorP, Correction = "Chao1", Alpha = 0.05, JackOver = FALSE,
CheckArguments = TRUE, Ps = NULL, Ns = NULL)
## S3 method for class 'integer'
Richness(NorP, Correction = "Chao1", Alpha = 0.05, JackOver = FALSE,
CheckArguments = TRUE, Ps = NULL, Ns = NULL)
## S3 method for class 'numeric'
Richness(NorP, Correction = "Chao1", Alpha = 0.05, JackOver = FALSE,
CheckArguments = TRUE, Ps = NULL, Ns = NULL)
``` |

`Ps` |
A probability vector, summing to 1. |

`Ns` |
A numeric vector containing species abundances. |

`NorP` |
A numeric vector, an integer vector, an abundance vector ( |

`Correction` |
A string containing one of the possible corrections: |

`Alpha` |
The risk level, 5% by default, used to optimize the jackknife order. |

`JackOver` |
If |

`CheckArguments` |
Logical; if |

Bias correction requires the number of individuals. Use `bcRichness`

and choose the `Correction`

.

Chao correction techniques are from Chao (1984) and Chiu *et al.* (2015). The Jackknife estimator is calculating by a straight adaptation of the code by Ji-Ping Wang (`jackknife`

in package `SPECIES`

). The optimal order is selected according to Burnham and Overton (1978; 1979). The argument `JackOver`

allows selecting one order over the optimal one.
Many other estimators are available elsewhere, the ones implemented here are necessary for other entropy estimations.

The functions are designed to be used as simply as possible. `Richness`

is a generic method. If its first argument is an abundance vector, an integer vector or a numeric vector which does not sum to 1, the bias corrected function `bcRichness`

is called. Explicit calls to `bcRichness`

(with bias correction) or to `Richness.ProbaVector`

(without correction) are possible to avoid ambiguity. The `.integer`

and `.numeric`

methods accept `Ps`

or `Ns`

arguments instead of `NorP`

for backward compatibility.

A named number equal to the estimated number of species.

Eric Marcon <Eric.Marcon@ecofog.gf>

Burnham, K. P., and Overton,W. S. (1978), Estimation of the Size of a Closed Population When Capture Probabilities Vary Among Animals. *Biometrika*, 65:625-633.

Burnham, K. P., and Overton,W. S. (1979), Robust Estimation of Population Size When Capture Probabilities Vary Among Animals. *Ecology* 60:927-936.

Chao, A. (1984) Nonparametric estimation of the number of classes in a population. *Scandinavian Journal of Statistics* 11:265-270.

Chiu, C.-H., Wang, Y.-T., Walther, B. A., Chao, A. (2014) An Improved Nonparametric Lower Bound of Species Richness via a Modified Good-Turing Frequency Formula. *Biometrics* 70(3):671-682.

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 | ```
# Load Paracou data (number of trees per species in two 1-ha plot of a tropical forest)
data(Paracou618)
# Ns is the total number of trees per species
Ns <- as.AbdVector(Paracou618.MC$Ns)
# Species probabilities
Ps <- as.ProbaVector(Paracou618.MC$Ns)
# Whittaker plot
plot(Ns)
# Number of observed species
Richness(Ps)
# Estimate the actual number of species
bcRichness(Ns, Correction = "Chao1")
bcRichness(Ns, Correction = "iChao1")
bcRichness(Ns, Correction = "Jackknife")
bcRichness(Ns, Correction = "Jackknife", JackOver=TRUE)
``` |

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

All documentation is copyright its authors; we didn't write any of that.