Tukey의 HSD 게시물 - 특별 3 개 그룹에 대한 R

Question

나는 prop.test을()하고있어에서 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.041 매우 낮다 e-13). 어떤 그룹의 쌍이 상당한 비율의 차이를 보이는지 알고 싶습니다.

나는 다음과 같은 코드를 실행할 수 있습니다 생각

....

library(agricolae)  
HSD.test(prop.test(t(data))) 

...하지만 HSD.test은 AOV 또는 LM 개체를 처리 할 수 ​​있기 때문에 작동하지 않습니다.

남은 해결책은 무엇입니까? 솔루션이 Tukey 수정 대신 Bonferroni (또는 다른) 수정을 의미하는지 여부와 상관없이.

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