R中的线性对比和Anova

时间:2015-03-31 01:40:14

标签: r statistics models linear anova

我试图理解这个问题(我的老师正在度假),我希望得到一些帮助。

  

使用“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又适合哪里?很有责任。

1 个答案:

答案 0 :(得分:1)

在适合模型之前,需要进行设置对比。因此,如果您在调用options(contrasts=c("contr.sum", "contr.poly"))之前运行aov(),那么您将获得所需的模型。请注意,options()调用会将以后调用的默认值更改为contrasts();如果你不想要你可以设置一个因素的对比而不改变这样的默认值:

contrasts(WorkplaceCO$Area) <- contr.sum
contrasts(WorkplaceCO$Time) <- contr.sum

无论您采用哪种方式,都需要在aov()之前完成。