powered by
GUniFrac is also sensitive to different sequencing depth. To compare microbiomes on an equal basis, rarefaction might be used.
Rarefy(otu.tab, depth = min(rowSums(otu.tab)))
Return a list containing:
rarefied OTU table
IDs of samples that does not reach the specified sequencing depth
OTU count table, row - n sample, column - q OTU
required sequencing depth; If not specified, the lowest sequencing depth is used.
Jun Chen <chen.jun2@mayo.edu>
Chen, J., Bittinger, K., Charlson, E.S., Hoffmann, C., Lewis, J., Wu, G.D., Collman, R.G., Bushman, F.D. and Li, H. (2012). Associating microbiome composition with environmental covariates using generalized UniFrac distances. 28(16): 2106–2113.
data(throat.otu.tab) # Rarefaction otu.tab.rff <- Rarefy(throat.otu.tab, 1024)$otu.tab.rff
Run the code above in your browser using DataLab