我对 R 很陌生,但我最近尝试进行双向重复测量方差分析,以复制我的主管在 SPSS 上所做的结果。我已经挣扎了好几天,阅读了几十篇文章来了解 R 中发生了什么,但我仍然没有得到相同的结果。
> mod <- lm(Y~A*B)
> Anova(mod, type="III")
Anova Table (Type III tests)
Response: Y
Sum Sq Df F value Pr(>F)
(Intercept) 0.000 1 0.0000 1.00000
A 2.403 5 8.6516 4.991e-08 ***
B 0.403 2 3.6251 0.02702 *
A:B 1.220 10 2.1962 0.01615 *
Residuals 51.987 936
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
我的数据来自平衡设计,我使用了 III 型 SS,因为它也是 SPSS 中使用的。Sum 平方 Df,线性模型在 SPSS 中是相同的,唯一不同的是 F 和 p 值。因此,它不应该是和平方错误。
SPSS中的结果是:
F Sig.
A 7.831 .000
B 2.681 .073
A:B 2.247 .014
我有点迷路了。这会是与对比度有关的问题吗?
卢卡斯