inst/extexamples/hpd.summary-coda.R

In emmeans: Estimated Marginal Means, aka Least-Squares Means

# Create an emmGrid object from a system file
cbpp.rg <- do.call(emmobj, 
                   readRDS(system.file("extdata", "cbpplist", package = "emmeans")))
cbpp.emm <- emmeans(cbpp.rg, "period")

hpd.summary(cbpp.emm)   # or just summary(cbpp.emm) as it gets redirected

# Equivalence test for any two-fold difference
summary(pairs(cbpp.emm), type = "response", delta = log(2))