HDIbetabinom: Get highest density interval of beta-binomial

In corncob: Count Regression for Correlated Observations with the Beta-Binomial

Get highest density interval of beta-binomial

Description

Get highest density interval of beta-binomial

Usage

HDIbetabinom(percent, M, mu, phi)

Arguments

percent	Numeric. Percent interval desired.
M	Numeric vector of sequencing depth
mu	Numeric vector of abundance parameter
phi	Numeric vector of dispersion parameter

Value

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

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])

corncob documentation built on May 29, 2024, 7:15 a.m.