我试图理解这个问题(我的老师正在度假),我希望得到一些帮助。
使用“contr.sum”对比选项,进行双向方差分析(ANOVA),其中包括时间和面积作为主效应以及两种主效应之间的相互作用。
我不确定如何使用对比来进行方差分析。我的答案似乎是分开的。我创建了一个像这样的模型
modelCO1 = aov(CO~Time+Area+(Time*Area), data = WorkplaceCO)
但这与线性对比无关。每当我尝试使用此代码时
modelCO1$contrasts$Time
我得到的输出只是说
"contr.sum"
哪个不能告诉我任何事情。或者,我已经完成了这个
options(contrasts=c("contr.sum", "contr.poly"))
contrasts(WorkplaceCO$Area)
contrasts(WorkplaceCO$Time)
这给了我一个输出
> contrasts(WorkplaceCO$Area)
[,1]
Nonsmoking 1
Smoking -1
> contrasts(WorkplaceCO$Time)
[,1] [,2] [,3] [,4] [,5]
7:00am 1 0 0 0 0
10:00am 0 1 0 0 0
11:00am 0 0 1 0 0
1:20pm 0 0 0 1 0
4:20pm 0 0 0 0 1
7:00pm -1 -1 -1 -1 -1
但是,ANOVA又适合哪里?很有责任。
答案 0 :(得分:1)
在适合模型之前,需要进行设置对比。因此,如果您在调用options(contrasts=c("contr.sum", "contr.poly"))
之前运行aov()
,那么您将获得所需的模型。请注意,options()
调用会将以后调用的默认值更改为contrasts()
;如果你不想要你可以设置一个因素的对比而不改变这样的默认值:
contrasts(WorkplaceCO$Area) <- contr.sum
contrasts(WorkplaceCO$Time) <- contr.sum
无论您采用哪种方式,都需要在aov()
之前完成。