分散分析(個人差を考える)

放送大学「心理統計法」第11章　複雑な分散分析２

A<-c(5,4,5,3,4,2,5,4,4,3)
B<-c(3,5,1,3,2,1,3,4,2,1)
C<-c(3,4,2,2,4,3,3,2,4,3)
D<-c(5,5,4,5,3,4,4,5,4,5)
d<-data.frame(A,B,C,D)
d1<-stack(d)
d1$S<-c(rep(1:10,4))

アクティブデータセットをd1に変更
d1$S 数値変数を因子に変換

d1$S <- as.factor(d1$S)
library(xtable)
print(xtable(summary(aov( values ~ ind + S ,data= d1 ))),type="html")

	Df	Sum Sq	Mean Sq	F value	Pr(>F)
ind	3	22.10	7.37	8.32	0.0004
S	9	11.90	1.32	1.49	0.2004
Residuals	27	23.90	0.89

多重比較

library(knitr)
kable(pairwise.t.test(d1$values,d1$ind,p.adjust.method="bonferroni",paired=TRUE)$p.value)

	A	B	C
B	0.0767942	NA	NA
C	0.3256465	1.0000000	NA
D	1.0000000	0.0043229	0.0403411

• BD間 CD間に有意差あり（有意水準５％）