在 R 中对 prop.test(..) 进行 Tukey HSD 事后测试答案

【问题标题】：Tukey HSD post-hoc test on a prop.test(..) in R在 R 中对 prop.test(..) 进行 Tukey HSD 事后测试
【发布时间】：2013-05-30 14:05:03
【问题描述】：

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

【问题讨论】：

标签： r posthoc

【解决方案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

【讨论】：