Uses a Dirichlet multinomial distribution to test for OTUs in the first vector that are enriched.

Usage

polyatest(counts1, counts2, p.adjust.method = "fdr")

Arguments

counts1
the first paired sample (and the target of the one-sided test)
counts2
the second paired sample (i.e. prewash)