skbio.diversity.alpha.chao1

skbio.diversity.alpha.chao1(counts, bias_corrected=True)[source]

Calculate chao1 richness estimator.

State: Experimental as of 0.4.0.

Uses the bias-corrected version unless bias_corrected is False and there are both singletons and doubletons.

Parameters
  • counts (1-D array_like, int) – Vector of counts.

  • bias_corrected (bool, optional) – Indicates whether or not to use the bias-corrected version of the equation. If False and there are both singletons and doubletons, the uncorrected version will be used. The biased-corrected version will be used otherwise.

Returns

Computed chao1 richness estimator.

Return type

double

See also

chao1_ci()

Notes

The uncorrected version is based on Equation 6 in 1:

\[chao1=S_{obs}+\frac{F_1^2}{2F_2}\]

where \(F_1\) and \(F_2\) are the count of singletons and doubletons, respectively.

The bias-corrected version is defined as

\[chao1=S_{obs}+\frac{F_1(F_1-1)}{2(F_2+1)}\]

References

1

Chao, A. 1984. Non-parametric estimation of the number of classes in a population. Scandinavian Journal of Statistics 11, 265-270.