「統計学関連なんでもあり」の過去ログ--- 041

No.06566　Nested ANOVAについて　　【林】　2008/05/17(Sat) 16:37

Nested ANOVA.ppt -- Nested Designsを参考にRをやってみました。

X <- c(10,12,8,13,11,13,9,10,13,14,7,10,14,8,10,12,14,11,10,9,10,13,9,7,9,10,12,11,8,9,8,8,16,12,5,4)
Areas <- rep(rep(1:3, each=4), 3)
Sites <- rep(1:4, 9)
DF <- data.frame(F1=as.factor(Areas), F2=as.factor(Sites), X=X)
summary(aov (X ~ F1/F2, DF))

途中まで同じ結果になるのですが，
Df Sum Sq Mean Sq F value Pr(>F)
F1 2 4.50 2.25 0.5000 0.61271
F1:F2 9 128.25 14.25 3.1667 0.01156 *
Residuals 24 108.00 4.50

F1のF valueとPr(>F)が異なるのですが，
なぜでしょうか？

No.06568　Re: Nested ANOVAについて　　【青木繁伸】　2008/05/17(Sat) 19:02

> F1のF valueとPr(>F)が異なるのですが，なぜでしょうか？

わたしは，その分析法は知りませんが，引用されたパワーポイントファイルを見ると，スライド12 に F Area と F SiteA(B) の計算式が書いてあり，そのとおりすると，当たり前ですが，スライドに示したとおりになる。ということは，aov(X ~ F1/F2, DF) が間違えているというか，この検定を表現していないということでしょう？違うものの結果について「違うのはなぜか」といわれても？？？

No.06576　Re: Nested ANOVAについて　　【林】　2008/05/17(Sat) 20:36

早速コメントありがとうございました。やはりRへの指示が違っていたようです。
ヤマカンで，
summary(aov(X ~ F1+ Error(F1:F2), DF))
と入力したら同じ結果が出ました。

（Sitesについては出てこなくなりましたが，できれば一つの指示で済んでほしいものです）