R中的线性对比度和方差分析

Question

我正在努力理解这个问题（我的老师正在休假），希望能得到一些帮助。

使用“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

但同样，方差分析适用于何处？非常感谢。

Answer 1

需要在拟合模型之前设置对比度。因此，如果您在调用aov() 之前运行options(contrasts=c("contr.sum", "contr.poly"))，那么您将获得您想要的模型。请注意，options() 调用会将未来调用的默认值更改为 contrasts()；如果您不希望这样，您可以在不更改默认值的情况下设置一个因子的对比度：

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

无论你怎么做，都需要在aov()之前完成。