corncob (version 0.4.2)

HDIbetabinom: Get highest density interval of beta-binomial

Description

Get highest density interval of beta-binomial

Usage

HDIbetabinom(percent, M, mu, phi)

Value

List where lower represents the lower bound and upper represents the upper bound

Arguments

percent

Numeric. Percent interval desired.

M

Numeric vector of sequencing depth

mu

Numeric vector of abundance parameter

phi

Numeric vector of dispersion parameter

Examples

data(soil_phylum_small_otu1)
mod <- bbdml(formula = cbind(W, M - W) ~ DayAmdmt,
phi.formula = ~ DayAmdmt,
data = soil_phylum_small_otu1)
HDIbetabinom(.95, M = mod$M[1], mu = mod$mu.resp[1], phi = mod$phi.resp[1])