差分

このページの2つのバージョン間の差分を表示します。

--- メモ:分散分析 [2017/09/07 09:56] – Wiki Editor
+++ メモ:分散分析 [2017/09/08 18:01] – [備忘録] Wiki Editor
@@ 行 11: / 行 11: @@
 ===== 備忘録 =====
-ネストした分散分析について
+==== ネストした分散分析について ====
 ネスト構造にあるデータの分散分析では、親要素のみで分散分析をした場合と、子要素もモデルに含めて分析した場合では、親要素の平方和は変わらない。子要素を投入すると、誤差の平方和が子要素で説明できる分だけ減少し、その分が子要素の平方和となる。
@@ 行 67: / 行 68: @@
 こうすることで、子は親にネストされているものとして平方和が計算される。具体的には、子要素の平方和は親要素との交互作用で計算される。そのため、親要素をまたがって子要素のコードが重なっていても、親要素が異なれば計算上は無関係となるので、期待通りの計算結果が得られる。
+==== 単純主効果が有意な場合の多重比較について ====
+要因の分散分析において交互作用が有意になった場合、単純主効果の検定を行う。
+ここでは竹内理・水本篤『外国語教育研究ハンドブック 【改訂版】 研究手法のより良い理解のために 』2014年8月の第7章を例に、筆者が提供してくれる[[http://mizumot.com/handbook/?page_id=246|データ]]（繰り返しなしの二元配置分散分析）を用いてRで同様の分析する方法をメモしておく。ポイントは単純主効果が有意な場合の多重比較の方法であるが、ついでなので7章の分析方法を全て再現しておく。
+SPSSでは一般線形モデルで分析を行うようである。分散分析表で平方和のTypeがIIIとなっているので、RでもType IIIの平方和を計算するためにcarパッケージを読込、Anova関数を利用する。
+  > library(car) # パッケージの読込
+  # Type III平方和を計算する際には、factorのダミーコードの方法を「contr.sum」に変更しないと計算結果がおかしくなる
+  > options(contrasts=c("contr.sum","contr.sum"))
+  # 筆者のサイトからダウンロードしてデータを読み込む
+  > DAT <- read.csv("データ")
+  # 数値でコーディングされているのでfactorとして認識させる
+  > DAT$教室 <- factor(DAT$教室);
+  > DAT$クラスサイズ <- factor(DAT$クラスサイズ);
+  # carパッケージのAnova関数を利用して分散分析
+  > Anova(lm(点数 ~ 教室 * クラスサイズ, data=DAT), type=3);
+  Anova Table (Type III tests)
+  Response: 点数
+                    Sum Sq  Df   F value    Pr(>F)
+  (Intercept)       664301   1 6575.8307 < 2.2e-16 ***
+  教室                  26   1    0.2618    0.6095
+  クラスサイズ       11238   2   55.6219 < 2.2e-16 ***
+  教室:クラスサイズ   4057   2   20.0823 1.421e-08 ***
+  Residuals          17578 174
+  # 偏イータ2乗の計算
+  # 対象の変数の平方和 / （対象の変数の平方和 + 誤差平方和）
+  # Anovaの結果を変数に投入
+  > res.anova <- Anova(lm(点数 ~ 教室 * クラスサイズ, data=DAT), type=3);
+  # res.anovaから平方和を取り出して計算
+  # 教室の偏イータ2乗
+  > res.anova$"Sum Sq"[2] / (res.anova$"Sum Sq"[2] + res.anova$"Sum Sq"[5])
+  [1] 0.001502481
+  # クラスサイズの偏イータ2乗
+  > res.anova$"Sum Sq"[3] / (res.anova$"Sum Sq"[3] + res.anova$"Sum Sq"[5])
+  [1] 0.3899955
+  # 交互作用の偏イータ2乗
+  > res.anova$"Sum Sq"[4] / (res.anova$"Sum Sq"[4] + res.anova$"Sum Sq"[5])
+  [1] 0.187541