RでのSheffeの多重比較

RでのScheffeの多重比較
Scheffeの多重比較を行うには?
• 青木繁伸先生作成の関数scheffeを利用する
• aov関数でのデフォルトの対比を変更して分
散分析を実行する.ただし,この方法は
Scheffe の方法とは少し異なる.
– 得られる t 値はほぼ同じ
– 参考:Crawley, M. J. 『統計学:Rを用いた入門書』
(共立出版)第12章.
データ
• ボーンシュテット & ノーキ『社会統計学』表7.1
のデータを用いる.
– 監視の有無が課題遂行に及ぼす影響.1要因3
水準のデザイン.
– 社会統計学の授業ウェブからデータをCSVでダウ
ンロードできる.
> monitor <- read.csv("table7_1.csv")
> head(monitor)
score condition
1 13
na
2 14
na
3 10
na
4 11
na
na:監視なしー隔離条件
5 12
na
nt:監視なしー共作業条件
6 10
na
m:監視条件
>
各条件の平均値と分散
> tapply(monitor$score, monitor$condition,
mean)
m na nt
8.2 11.8 10.0
>
> tapply(monitor$score, monitor$condition, var)
m
na
nt
2.622222 1.733333 1.333333
>
デフォルトの contrasts
> contrasts(monitor$condition)
na nt
m 0 0
na 1 0
nt 0 1
対比の指定
> contrasts(monitor$condition)
<- cbind(c(1,-1/2,-1/2),c(0,-1,1))
> contrasts(monitor$condition)
[,1] [,2]
m 1.0 0
na -0.5 -1
nt -0.5 1
>
分析の実行
> summary.lm(aov(monitor$score ~ monitor$condition))
Call:
aov(formula = monitor$score ~ monitor$condition)
Residuals:
Min
1Q Median
3Q
Max
-2.200e+00 -9.500e-01 3.188e-17 9.500e-01 2.800e+00
続きは次のスライド
分析の実行(つづき)
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept)
10.0000 0.2514 39.775 < 2e-16 ***
monitor$condition1 -1.8000 0.3556 -5.063 2.57e-05 ***
monitor$condition2 -0.9000 0.3079 -2.923 0.00694 **
--Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Residual standard error: 1.377 on 27 degrees of freedom
Multiple R-squared: 0.5586, Adjusted R-squared: 0.5259
F-statistic: 17.09 on 2 and 27 DF, p-value: 1.603e-05