R中的线性对比和Anova

我试图理解这个问题（我的老师正在度假），我希望得到一些帮助。R中的线性对比和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适合在哪里？多谢。

来源

2015-03-31 igiari21

设置对比度需要在适合模型之前完成。所以如果你在拨打aov()之前运行options(contrasts=c("contr.sum", "contr.poly"))，那么你会得到你想要的模型。请注意，options()呼叫会将未来呼叫的默认值更改为contrasts();如果你不想，你可以在一个因素集对比不改变这样的默认值：

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

你做到这一点无论怎样，它需要aov()之前完成。

来源

2015-03-31 02:02:03 drammock

输出的其余部分是否正常？我只是不确定我的对比符合我的方差分析，因为它们似乎并不相互联系。 – igiari21 2015-03-31 02:15:54

对比度决定如何计算ANOVA系数。当你用一个没有指定'contrasts'的因子调用'aov'时，它会在模型拟合之前自动完成（使用默认的对比度方法）。你可以阅读更多关于[这里]（http://talklab.psy.gla.ac.uk/tvw/catpred/）和[这里]（http://www.ats.ucla.edu/stat/r/库/ contrast_coding.htm）。 – drammock 2015-03-31 02:50:51

R中的线性对比和Anova

回答

相关问题