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

看板Statistics作者 ( )時間7年前 (2018/05/30 00:12), 7年前編輯推噓0(0推 0噓 8→)
留言8則, 2人參與, 7年前最新討論串1/1
[軟體程式類別]: 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
→
05/30 00:46, 7年前 , 1F
google Helmert 就有一堆解釋了 關係到coding的問題
05/30 00:46, 1F
→
05/30 00:47, 7年前 , 2F
也可以找contrast matrix, L-matrix 更詳細
05/30 00:47, 2F
→
05/30 01:37, 7年前 , 3F
重點不是helmert，而是在不平衡設計時需要正交對比
05/30 01:37, 3F
→
05/30 01:51, 7年前 , 4F
其它軟體没強調這個是因為在變方拆解時它自動做這個動作
05/30 01:51, 4F
→
05/30 01:53, 7年前 , 5F
在看板R_Language有過簡單討論，查標題type III anova
05/30 01:53, 5F
感謝回答 想順便請教，改變contrast方式會影響GLM的residuals嗎？ ※ 編輯: MTIS (59.126.233.13), 05/30/2018 21:15:32
→
05/31 00:07, 7年前 , 6F
應該是不會。你可以利用resid(lm.aov)取得殘差值
05/31 00:07, 6F
→
05/31 09:30, 7年前 , 7F
如果factor有3-level以上有可能 要看背後原始碼怎麼
05/31 09:30, 7F
→
05/31 09:30, 7年前 , 8F
弄coding
05/31 09:30, 8F
‣ 更多 MTIS 的文章
文章代碼(AID): #1R3NlKHC (Statistics)