作者MTIS ( )
看板Statistics
標題[程式] R和SPSS做ANOVA結果不一致？
時間Wed May 30 00:12:01 2018
[軟體程式類別]:
R 3.5
SPSS 20
[程式問題]:
ANOVA
[軟體熟悉度]:
熟悉
[問題敘述]:
https://www.r-bloggers.com/ensuring-r-generates-
the-same-anova-f-values-as-spss/
我已經研究過了，問題似乎出在contrast，
可是這邊還是不太懂，究竟用Helmert的好處是什麼？
難道R採用的預設設定是比較嚴謹的嗎？
因為我用同樣資料做出來...R的結果沒有顯著...
有什麼書籍有專門談contrast的嗎？
這結果真的太震撼了....
以前統計系老師都沒說過這問題...〒△〒
[程式範例]:
require("car")
lm.aov<-lm(Y ~ X1*X2, data=data)
Anova(lm.aov, type=3)
如果不改options(contrasts = c("contr.helmert", "contr.poly"))
結果會和SPSS不同，而且R算的結果p value比較大。
--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 59.126.233.13
※ 文章網址: https://webptt.com/m.aspx?n=bbs/Statistics/M.1527610324.A.44C.html
1F：→ recorriendo: google Helmert 就有一堆解釋了 關係到coding的問題 05/30 00:46
2F：→ recorriendo: 也可以找contrast matrix, L-matrix 更詳細 05/30 00:47
3F：→ andrew43: 重點不是helmert，而是在不平衡設計時需要正交對比 05/30 01:37
4F：→ andrew43: 其它軟體没強調這個是因為在變方拆解時它自動做這個動作 05/30 01:51
5F：→ andrew43: 在看板R_Language有過簡單討論，查標題type III anova 05/30 01:53
感謝回答
想順便請教，改變contrast方式會影響GLM的residuals嗎？
※ 編輯: MTIS (59.126.233.13), 05/30/2018 21:15:32
6F：→ andrew43: 應該是不會。你可以利用resid(lm.aov)取得殘差值 05/31 00:07
7F：→ recorriendo: 如果factor有3-level以上有可能 要看背後原始碼怎麼 05/31 09:30
8F：→ recorriendo: 弄coding 05/31 09:30