skbio.math.diversity.alpha.ace

skbio.math.diversity.alpha.ace(counts, rare_threshold=10)[source]

Calculate the ACE metric (Abundance-based Coverage Estimator).

Parameters:

counts : 1-D array_like, int

Vector of counts.

rare_threshold : int, optional

Threshold at which an OTU containing as many or fewer individuals will be considered rare.

Returns:

double

Computed ACE metric.

Raises:

ValueError

If every rare OTU is a singleton.

Notes

ACE was first introduced in [R26] and [R27]. The implementation here is based on the description given in the EstimateS manual [R28].

If no rare OTUs exist, returns the number of abundant OTUs. The default value of 10 for rare_threshold is based on [R29].

If counts contains zeros, indicating OTUs which are known to exist in the environment but did not appear in the sample, they will be ignored for the purpose of calculating the number of rare OTUs.

References

[R26](1, 2) Chao, A. & S.-M Lee. 1992 Estimating the number of classes via sample coverage. Journal of the American Statistical Association 87, 210-217.
[R27](1, 2) Chao, A., M.-C. Ma, & M. C. K. Yang. 1993. Stopping rules and estimation for recapture debugging with unequal failure rates. Biometrika 80, 193-201.
[R28](1, 2) http://viceroy.eeb.uconn.edu/estimates/
[R29](1, 2) Chao, A., W.-H. Hwang, Y.-C. Chen, and C.-Y. Kuo. 2000. Estimating the number of shared species in two communities. Statistica Sinica 10:227-246.