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>, Lujun Zhang
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