inst/extexamples/hpd.summary-coda.R

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

Try the emmeans package in your browser

Any scripts or data that you put into this service are public.

emmeans documentation built on July 3, 2024, 5:11 p.m.