R中的自定义对比(类似于SPSS中的特殊'对比)

时间:2014-10-10 16:26:47

标签: r statistics anova

我经常使用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)?

0 个答案:

没有答案