私は3つのグループで prop.test() をやっています:
data = data.frame(a = c(85*0.3412,85*(1-0.3412)),b = c(203*0.2217,203*(1-0.2217)), c = c(231,197))
prop.test(t(data))
p.value は非常に低い (1.041e-13)。どのグループのペアが比率の有意差を示しているかを知りたいです。
次のコードを実行できると思いました....
library(agricolae)
HSD.test(prop.test(t(data)))
...しかし、HSD.test は aov または lm オブジェクトしか処理できないため、機能していません。
どのような解決策が残っていますか? ソリューションが Tukey 補正ではなく Bonferroni (またはその他) の補正を意味するかどうかは関係ありません。