polyatest(counts1, counts2, p.adjust.method = "fdr")
- the first paired sample (and the target of the one-sided test)
- the second paired sample (i.e. prewash)
- the method to adjust p.values for multiple testing
a vector of adjusted p values
Uses a Dirichlet multinomial distribution to test for OTUs in the first
vector that are enriched.