Апостериорный тест Тьюки HSD на prop.test(..) в R

Я делаю prop.test() для 3 групп:

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 очень низкое (1,041e-13). Я хочу знать, какая пара групп демонстрирует значительную разницу в пропорциях.

Я думал, что смогу запустить следующий код....

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

... но это не работает, потому что HSD.test может работать только с объектами aov или lm.

Какое решение у меня осталось? Неважно, подразумевает ли решение поправку Бонферрони (или другую) вместо поправки Тьюки.


r posthoc
person Remi.b    schedule 30.05.2013    source источник

Ответы (1)


arrow_upward
3
arrow_downward

По умолчанию выполняется корректировка Бонферрони-Холма; см. документацию для других вариантов.

> 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
person Aaron left Stack Overflow    schedule 30.05.2013