重复测量ANOVA - SPSS与R的不同结果

时间:2015-03-31 09:28:27

标签: r spss anova

我正在尝试使用R运行重复测量ANOVA并将其与SPSS输出进行比较,结果差异很大!也许我在某处弄错了,但我无法弄明白 一些样本数据: id是主题。每个主题对三个项目(res_1,res_2和res_3)进行一次评级。我想比较一下项目的整体效果。

id<-c(1,2,3,4,5,6)
res_1<-c(1,1,1,2,2,1)
res_2<-c(4,5,2,4,4,3)
res_3<-c(4,5,6,3,6,6)
## wide format for spss
table<-as.data.frame(cbind(id, res_1, res_2, res_3))
## reshape to long format
library(reshape2)
table<-melt(table, id.vars="id")
colnames(table)<-c("id", "item", "rating")
aov.out = aov(rating ~ item+ Error(id/item), data=table) 
summary(aov.out)

这是我的SPSS代码(来自宽幅数据)

GLM item_1 item_2 item_3
/WSFACTOR=factor1 3 Polynomial 
/METHOD=SSTYPE(3)
/PRINT=DESCRIPTIVE 
/CRITERIA=ALPHA(.05)
/WSDESIGN=factor1.

我得到的结果 R:p值0.0526(误差:在)内 和SPSS:p值0.003(受试者效应内的测试)

有没有人有可能解释这种差异的建议? 如果我进行非参数Friedmann检验,我会在SPSS和R中得到相同的结果。 实际上,在查看我的数据时,摘要(aov.out)与SPSS&#34;对象内部的测试对比&#34; (但我学会了观察受试者内部效果的测试)。

谢谢!

1 个答案:

答案 0 :(得分:0)

那里有很多东西;我有点惊讶你的google for&#39; spss vs. R anova&#39;没有带您解释SPSS(III型)和R(I型)之间的平方和的差异以及对比度的差异。

这是我发现的前两个结果: http://myowelt.blogspot.ca/2008/05/obtaining-same-anova-results-in-r-as-in.htmlhttps://stats.stackexchange.com/questions/40958/testing-anova-hypothesis-with-contrasts-in-r-and-spss