割と簡単だったのでメモ
とりあえず multcomp をインストールして、aov で分析して一旦Tukeyで検定してから補正をかけるイメージ。
Shafferが強いと思っていたが、もっと強いのがWestfallだった
res.aov <- aov(outcome ~ treatment, data = df)
res.tukey <- glht(res.aov)
summary(res.tukey, test = adjusted("Shaffer"))
summary(res.tukey, test = adjusted("bonferroni"))
summary(res.tukey, test = adjusted("Westfall"))
summary(res.tukey, test = adjusted("single-step"))
