dataset <- data.frame(
outcome = c(1, 2, 3, 2, 3, 4, 5, 6, 7),
group = factor(c("a", "a", "a", "b", "b", "b", "c", "c", "c"))
)
anova1 <- aov(outcome ~ group, data = dataset)
summary(anova1)
# post-hoc pairwise comparisons with Holm correction (the default)
posthocPairwiseT(anova1)
# Bonferroni correction instead
posthocPairwiseT(anova1, p.adjust.method = "bonferroni")
Run the code above in your browser using DataLab