HDIbetabinom | R Documentation |
Get highest density interval of beta-binomial
HDIbetabinom(percent, M, mu, phi)
percent |
Numeric. Percent interval desired. |
M |
Numeric vector of sequencing depth |
mu |
Numeric vector of abundance parameter |
phi |
Numeric vector of dispersion parameter |
List where lower
represents the lower bound and upper
represents the upper bound
data(soil_phylum_small)
mod <- bbdml(formula = OTU.1 ~ DayAmdmt,
phi.formula = ~ DayAmdmt,
data = soil_phylum_small)
HDIbetabinom(.95, M = mod$M[1], mu = mod$mu.resp[1], phi = mod$phi.resp[1])
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.