alphaDiversity: Calculate the alpha diversity indices

In camel315/otuSummary: Summarizing OTU table regarding the composition, abundance and beta diversity of abundant and rare biospheres

Description

This function will calculate the alpha diversity indices for the total, abundant and rare biospheres.

Usage

alphaDiversity(otutab, siteInCol = FALSE, taxhead = NULL, threshold = 1,
    percent = TRUE, write = TRUE, ...)

Arguments

otutab	A OTU table of microbial community, which can contain a taxonomic column (if siteInCol) or row (if site in rows). The OTU table can be given in numeric counts or in relative abundance.
siteInCol	Logical, if "TRUE", the OTU table contains samples in columns and taxa in rows. By default in this function, the siteInCol is FALSE, meaning the samples in rows.
taxhead	Character, specify the header of taxonomy. By default this argument is NULL.
threshold	Numeric, the threshold of relative abundance upon which the rare biosphere will be subset.
percent	Logical, whether the input OTU table are given in relative abundance. FALSE means that the input OTU table is in numeric counts.
write	Logical, if TRUE, the result will be written out in a Tab separated data frame.
...	arguments to be passed to write.table().

Details

The rare biosphere is defined by the relative abundance cutoffs (which is the "threshold" argument in this function) (Lynch and Neufeld, 2015). This function call the functions "specnumber", "diversity" from package vegan (Oksanen et al, 2013), "chao1" and "chao2" from package fossil (Vavrek, 2011) and the "gini" function from package reldist (http://www.stat.ucla.edu/~handcock/RelDist).

Value

The function will return a list of length 3, including indices of observed, shannon, simpson, invsimpson, chao1, chao2, evenness and Gini.

allBio	The alpha diversity indices for the whole community
abundBio	The alpha diversity indices for the abundant population
rareBio	The alpha diversity indices for the rare biosphere

Author(s)

Sizhong Yang <[e-mail]>

References

Lynch MDJ, Neufeld JD (2015). Ecology and exploration of the rare biosphere. Nature Reviews Microbiology 13: 217-229.

Oksanen J, Blanchet FG, Kindt R, Legendre P, Minchin PR, O'Hara RB et al (2013). vegan: Community Ecology Package. R package version 2.0-7. http://CRAN.R-project.org/package=vegan.

Vavrek MJ (2011). fossil: Palaeoecological and palaeogeographical analysis tools. Palaeontol Electron 14:1T. http://palaeo-electronica.org/2011_1/238/index.html.

Examples

data(otumothur)

test1 <- alphaDiversity(otutab = otumothur, siteInCol = TRUE,
    taxhead = "taxonomy", threshold = 1, percent = FALSE, write = FALSE)

test2 <- alphaDiversity(otutab = otumothur[,-454], siteInCol = TRUE,
    taxhead = NULL, threshold = 1, percent = FALSE, write = FALSE)

camel315/otuSummary documentation built on May 5, 2019, 12:29 p.m.