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

回答

相關問題