作者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/cn.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