richness_calc: Calculates Richness

Description Usage Arguments Details Value Author(s) References Examples

Description

This function calculates the unique number of features seen in each sample.

Usage

1
richness_calc(omicsData, index = c("observed", "chao1", "ace", "break"))

Arguments

omicsData

an object of the class 'seqData' created by as.seqData.

index

a character vector stating which of the calculations to perform - "observed" for the observed richness, "chao1" the bias-corrected chao1 richness estimator, and/or "ace" for the abundance-based coverage richness estimator. Default is to perform all 3 calculations.

Details

Calculates richness of count data

Value

An object of class richRes (also a data.frame) containing the richness value for every sample in the data object.

Author(s)

Allison Thompson

References

Chao, Anne. Species Richness Estimation.

Examples

1
2
3
4
5
6
7
8
9
## Not run: 
library(mintJansson)
data(rRNA_data)
rRNA_richness <- richness_calc(omicsData = rRNA_data)
rRNA_richness
summary(rRNA_richness)
plot(rRNA_richness)

## End(Not run)

pmartR/pmartRseq documentation built on May 25, 2019, 9:20 a.m.