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

但是，ANOVA又适合哪里？很有责任。

Answer 1

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

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

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

R中的线性对比和Anova

1 个答案: