[程式] R和SPSS做ANOVA結果不一致？

[MTIS ( )]

[軟體程式類別]:
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://www.ptt.cc/bbs/Statistics/M.1527610324.A.44C.html
→ recorriendo: google Helmert 就有一堆解釋了 關係到coding的問題
→ recorriendo: 也可以找contrast matrix, L-matrix 更詳細
→ andrew43: 重點不是helmert，而是在不平衡設計時需要正交對比
→ andrew43: 其它軟體没強調這個是因為在變方拆解時它自動做這個動作
→ andrew43: 在看板R_Language有過簡單討論，查標題type III anova
感謝回答
想順便請教，改變contrast方式會影響GLM的residuals嗎？
※ 編輯: MTIS (59.126.233.13), 05/30/2018 21:15:32
→ andrew43: 應該是不會。你可以利用resid(lm.aov)取得殘差值
→ recorriendo: 如果factor有3-level以上有可能 要看背後原始碼怎麼
→ recorriendo: 弄coding