首页
学习
活动
专区
圈层
工具
发布
社区首页 >问答首页 >R中的线性对比与Anova

R中的线性对比与Anova
EN

Stack Overflow用户
提问于 2015-03-31 01:40:13
回答 1查看 675关注 0票数 1

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

利用“contr.sum”对比选项,进行时间和面积作为主效应的方差分析(ANOVA),以及两个主效应之间的相互作用。

我不太清楚如何用对比来做方差分析。我的答案似乎是分开的。我创造了一个这样的模型

代码语言:javascript
复制
modelCO1 = aov(CO~Time+Area+(Time*Area), data = WorkplaceCO)

但这与线性对比无关。每当我尝试使用这段代码

代码语言:javascript
复制
modelCO1$contrasts$Time

我得到的输出是

代码语言:javascript
复制
"contr.sum"

但这并没有告诉我什么。或者,我也这样做了

代码语言:javascript
复制
options(contrasts=c("contr.sum", "contr.poly"))

contrasts(WorkplaceCO$Area)
contrasts(WorkplaceCO$Time)

这给了我一个输出

代码语言:javascript
复制
> 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

但再说一遍,方差分析在哪里适用?非常感激。

EN

回答 1

Stack Overflow用户

发布于 2015-03-31 02:02:03

设置对比需要做之前,你适合的模型。因此,如果在调用options(contrasts=c("contr.sum", "contr.poly"))之前运行aov(),那么您将得到所需的模型。请注意,options()调用将更改对contrasts()的未来调用的默认值;如果不希望这样,您可以在不更改默认设置的情况下对一个因素设置对比:

代码语言:javascript
复制
contrasts(WorkplaceCO$Area) <- contr.sum
contrasts(WorkplaceCO$Time) <- contr.sum

不管你怎么做,它都需要在aov()之前完成。

票数 2
EN
页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接:

https://stackoverflow.com/questions/29359192

复制
相关文章

相似问题

领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档