- 先讀入資料檔。
ex<-read.table("ex16_1.csv",header=T,sep=",")
- 進行獨立樣本二因子變異數分析。
aovex<-aov(score~type*class,data=ex)
- 列出分析結果。
summary(aovex)
- 由ANOVA表可得知有交互作用(interaction)存在,因此必須再去探討單純主要效果(simple main effects)。
- 先將各資料分組(先以b因子分組;亦即是探討同在同樣的a因子之下會對不同的b因子造成什麼顯響。)
ex.b1<-subset(ex,class=="b1")
ex.b2<-subset(ex,class=="b2")
ex.b3<-subset(ex,class=="b3")
- 之後分別進行單因子變異數分析。
summary(aov(score~type,data=ex.b1))
summary(aov(score~type,data=ex.b2))
summary(aov(score~type,data=ex.b3))
No comments:
Post a Comment