私はRにはまったく慣れていませんが、最近、スーパーバイザーがSPSSで行った結果を複製するために、双方向の反復測定ANOVAを作成しようとしました。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
私のデータはバランスの取れた設計のものです。SPSS でも使用されているタイプ III の SS を使用しました。Sum は Df を 2 乗し、線形モデルは SPSS でも同じですが、異なるのは F 値と p 値だけです。したがって、Sum square の間違いではないはずです。
SPSS での結果は次のとおりです。
F Sig.
A 7.831 .000
B 2.681 .073
A:B 2.247 .014
私は少し迷っています。コントラストに関連する問題でしょうか？
ルーカス