Plots distribution of reads.
plot_read_distribution(x, groups, plot.type = c("density", "histogram"))
phyloseq-class
object
Metadata variable to check for groups based sequencing effort.
Either density or histogram plot
A ggplot
plot object.
library(microbiome)
data(zackular2014)
ps0 <- zackular2014
p <- plot_read_distribution(ps0, groups = "DiseaseState", plot.type = "density")
print(p)