# 1. Balanced one-way: the six marketing panels of the test_market
# data, four outlets per panel, so the procedure is exactly Tukey's
# HSD. Panels 5 and 6 separate from panel 1.
fit <- lm(brand_movement ~ panel, data = test_market)
ci_tukey_kramer(fit)
# 2. Same data via vector / group interface:
ci_tukey_kramer(test_market$brand_movement, group = test_market$panel)
Run the code above in your browser using DataLab