data(otumothur)
summaryInfo <- otuReport(otutab = otumothur, siteInCol = TRUE, taxhead = "taxonomy",
platform = "qiime", percent = FALSE, taxlevel = "phylum", collap = ";")
length(summaryInfo)
names(summaryInfo)
summaryInfo[[1]]
bplot(summaryInfo[["taxaFreqs"]])
Run the code above in your browser using DataLab