ちょっと結果の見かたについて教えてください
例えばRでOne way anovaをかけた結果
> summary(AnovaModel.1)
Df Sum Sq Mean Sq F value Pr(>F)
factor(var1) 2 3235 1617 6.367 0.0266 *
Residuals 7 1778 254
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
でPost hocは
> pairwise.t.test(Dataset$var2, Dataset$var1, p.adj="bonferroni")
Pairwise comparisons using t tests with pooled SD
data: Dataset$var2 and Dataset$var1
1 2
2 1.000 -
3 0.043 0.088
P value adjustment method: bonferroni
> TukeyHSD(AnovaModel.1, "factor(var1)")
Tukey multiple comparisons of means
95% family-wise confidence level
Fit: aov(formula = var2 ~ factor(var1), data = Dataset, na.action = na.omit)
$`factor(var1)`
diff lwr upr p adj
2-1 6.217252 -32.108606 44.54311 0.8837720
3-1 39.490291 3.639733 75.34085 0.0333967
3-2 33.273039 -2.577519 69.12360 0.0668032
Multiple Comparisons of Means: Dunnett Contrasts
Fit: aov(formula = var2 ~ group.factor, data = Dataset)
Linear Hypotheses:
Estimate Std. Error t value Pr(>|t|)
2 - 1 == 0 6.217 13.014 0.478 0.849
3 - 1 == 0 39.490 12.173 3.244 0.025 *
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
(Adjusted p values reported -- single-step method)
質問はDunnettではSignif. code *をつけてくれているんですが、その他の検定方法ではP値が0.05以下でもSignif. code がついていません。これは有意差がないという結論ではないですよね。。。 |
|