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)
p.adjust.method
the method to adjust p.values for multiple testing

Value

a vector of adjusted p values

Description

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