作者windguide (夜很深)
看板Statistics
標題[問題] 用R分析2^4 design data
時間Wed Apr 28 12:08:39 2010
如果是跟統計軟體有關請重發文章
如果跟論文有關也煩請您重發文章
文章類別是為了幫助大家搜尋資料與解答，造成不便之處請見諒
其實我遇到的問題是這樣的
我用R分析2^4的Factorial Design，在做ANOVA分析的時候沒有問題
問題是出在做回歸分析的時候：
我記得回歸分析的係數就代表了那個因子效應大小的一半
可是我做出來的回歸分析係數跟我實際計算的因子效應並沒有關係
例如：
回歸分析： x_1 的係數 = A因子效應大小的一半 = 8.2526
實際用手算： A因子大小 = 3.0188
請問為什麼會出現這樣的問題，是我的Model Fit錯誤嗎?
麻煩各位幫我看一下，謝謝
以下是CODE
y =
c(7.037,6.376,14.707,15.219,11.635,12.089,17.273,17.815,10.403,10.151,
4.368,4.098,9.360,9.253,13.440,12.923,8.561,8.951,16.867,17.052,
13.876,13.658,19.824,19.639,11.846,12.337,6.125,5.904,11.190,10.935,
15.653,15.053)
trtA <- gl(2,2,32)
trtB <- gl(2,4,32)
trtC <- gl(2,8,32)
trtD <- gl(2,16,32)
out = aov(y~trtA*trtB*trtC*trtD)
summary(out)
out1 = lm(y~trtA*trtB*trtC*trtD)
summary(out1)
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.112.75.121
1F：推 kazuki:在two-level design factorial effects跟 regression 04/28 21:03
2F：→ kazuki:係數的估值會呈比例，也就是你說得1/2還是兩倍。 04/28 21:04
3F：→ kazuki:在但mixed-level experiments裡就不一定了。以上是我猜的^^ 04/28 21:05