对于R中的prop.test（..）的Tukey HSD事后检验

Question

我在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对象。

我还有什么解决方案？ 无论解决方案是否暗示Bonferroni（或其他）校正而不是Tukey校正。

Answer 1

这默认为Bonferroni-Holm调整; 请参阅其他选项的文档。

> pairwise.prop.test(t(data))

    Pairwise comparisons using Pairwise comparison of proportions 

data:  t(data) 

  a      b      
b 0.0490 -      
c 0.0025 3.1e-13

P value adjustment method: holm