Як тест HSD Tukey може бути більш значущим, ніж некорекційне P значення t.test?


9

Мені прийшов пост " Пост-пара-паральні порівняння двостороннього ANOVA " (відповідаючи на це повідомлення ), де показано наступне:

dataTwoWayComparisons <- read.csv("http://www.dailyi.org/blogFiles/RTutorialSeries/dataset_ANOVA_TwoWayComparisons.csv")

model1 <- aov(StressReduction~Treatment+Age, data =dataTwoWayComparisons)
summary(model1) # Treatment is signif

pairwise.t.test(dataTwoWayComparisons$StressReduction, dataTwoWayComparisons$Treatment, p.adj = "none")
# no signif pair

TukeyHSD(model1, "Treatment")
# mental-medical   is the signif pair.

(Вихід додається внизу)

Невже хтось може пояснити, чому HSD Tukey здатний знайти значну пару, коли парний (невідрегульований pvalue) тест не справляється?

Дякую.


Ось вихідний код

> model1 <- aov(StressReduction~Treatment+Age, data =dataTwoWayComparisons)
> summary(model1) # Treatment is signif
            Df Sum Sq Mean Sq F value    Pr(>F)    
Treatment    2     18   9.000      11 0.0004883 ***
Age          2    162  81.000      99     1e-11 ***
Residuals   22     18   0.818                      
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1   1 
> 
> pairwise.t.test(dataTwoWayComparisons$StressReduction, dataTwoWayComparisons$Treatment, p.adj = "none")

        Pairwise comparisons using t tests with pooled SD 

data:  dataTwoWayComparisons$StressReduction and dataTwoWayComparisons$Treatment 

         medical mental
mental   0.13    -     
physical 0.45    0.45  

P value adjustment method: none 
> # no signif pair
> 
> TukeyHSD(model1, "Treatment")
  Tukey multiple comparisons of means
    95% family-wise confidence level

Fit: aov(formula = StressReduction ~ Treatment + Age, data = dataTwoWayComparisons)

$Treatment
                 diff         lwr        upr     p adj
mental-medical      2  0.92885267 3.07114733 0.0003172
physical-medical    1 -0.07114733 2.07114733 0.0702309
physical-mental    -1 -2.07114733 0.07114733 0.0702309

> # mental-medical   is the signif pair.

Відповіді:


9

Тому що ваш парний т-тест вище не коригується для віку, і вік пояснює велику дисперсію в StressReduction.


1
Привіт onestop. Я уявив, що це так, але не міг повністю зрозуміти, що на wikipage ( en.wikipedia.org/wiki/Tukey's_range_test ) є будь-які пропозиції щодо посилань, як це робиться?
Тал Галілі
Використовуючи наш веб-сайт, ви визнаєте, що прочитали та зрозуміли наші Політику щодо файлів cookie та Політику конфіденційності.
Licensed under cc by-sa 3.0 with attribution required.