我经常使用ezANOVA
包(或afex
)来进行重复测量方差分析 - 我学会了如何在SPSS中运行自定义对比,但在R中找不到简单的解决方案。
以下是一个例子:
require(plyr)
require(ez)
require(afex)
data(ANT)
#Aggregate
x<-ddply(ANT, .(subnum, cue, block), summarise, rt=mean(rt))
x$block<-factor(x$block)
#Anova
ezANOVA(x
,dv=rt
,wid=subnum
,within=.(cue, block)
,type=3
,detailed=T)
#Alternative Anova Method
a<-ez.glm("subnum","rt",x, within=c("cue","block"), type=3, return="Anova")
nice.anova(a, correction = c("none"))
这将返回ANOVA。我想做的是对数据进行手动对比。例如,对于提示类型,提示B不同于提示C&amp; D(0,-1,.5,.5)?