二元配置分散分析(主効果と交互作用)

二元配置分散分析(two-way ANOVA)は、2つの要因(因子)によって分類されたデータに対して、それぞれの要因の主効果(main effect)および両者の交互作用(interaction effect)を同時に検定する手法である。これは一元配置分散分析の自然な拡張であり、複数の要因が応答変数に与える影響を分離して評価することを目的とする。

本手法の本質は、観測データの総変動を「要因Aによる変動」「要因Bによる変動」「交互作用による変動」「誤差による変動」に分解し、それぞれの寄与を偏差平方和として定量化する点にある。

基本設定

要因Aが $a$ 水準、要因Bが $b$ 水準を持ち、各組合せ $(i,j)$ において $n_{ij}$ 個の観測があるとする:

\[x_{ijk}, \quad i=1,\dots,a,\quad j=1,\dots,b,\quad k=1,\dots,n_{ij}\]

総標本数は

\[N = \sum_{i=1}^{a} \sum_{j=1}^{b} n_{ij}\]

である。以下の平均を定義する:

\[\bar{x}_{ij} = \frac{1}{n_{ij}} \sum_{k=1}^{n_{ij}} x_{ijk}\]\[\bar{x}_{i\cdot} = \frac{1}{\sum_j n_{ij}} \sum_{j=1}^{b} \sum_{k=1}^{n_{ij}} x_{ijk}\]\[\bar{x}_{\cdot j} = \frac{1}{\sum_i n_{ij}} \sum_{i=1}^{a} \sum_{k=1}^{n_{ij}} x_{ijk}\]\[\bar{x} = \frac{1}{N} \sum_{i,j,k} x_{ijk}\]

線形モデルによる定式化

二元配置分散分析は、以下の線形モデルで表される:

\[x_{ijk} = \mu + \alpha_i + \beta_j + (\alpha\beta)_{ij} + \varepsilon_{ijk}\]

ここで、

  • $\mu$:全体平均
  • $\alpha_i$:要因Aの主効果
  • $\beta_j$:要因Bの主効果
  • $(\alpha\beta)_{ij}$:交互作用効果
  • $\varepsilon_{ijk}$:誤差項

識別性のために以下の制約を課す:

\[\sum_i \alpha_i = 0,\quad \sum_j \beta_j = 0\]\[\sum_i (\alpha\beta)_{ij} = 0,\quad \sum_j (\alpha\beta)_{ij} = 0\]

誤差項は独立同分布で

\[\varepsilon_{ijk} \sim \mathcal{N}(0, \sigma^2)\]

と仮定する。

偏差平方和の分解

総平方和(SST)は

\[\mathrm{SST} = \sum_{i,j,k} (x_{ijk} - \bar{x})^2\]

で定義され、以下のように分解される:

\[\mathrm{SST}= \mathrm{SSA} + \mathrm{SSB} + \mathrm{SSAB} + \mathrm{SSE}\]

ここで、

要因Aの平方和:

\[\mathrm{SSA}= \sum_{i=1}^{a} n_{i\cdot} (\bar{x}_{i\cdot} - \bar{x})^2\]

要因Bの平方和:

\[\mathrm{SSB}= \sum_{j=1}^{b} n_{\cdot j} (\bar{x}_{\cdot j} - \bar{x})^2\]

交互作用平方和:

\[\mathrm{SSAB}= \sum_{i=1}^{a} \sum_{j=1}^{b} n_{ij}(\bar{x}_{ij} - \bar{x}_{i\cdot} - \bar{x}_{\cdot j} + \bar{x})^2\]

誤差平方和:

\[\mathrm{SSE}= \sum_{i,j,k} (x_{ijk} - \bar{x}_{ij})^2\]

である。

この分解は、加法モデルからの逸脱(交互作用)を明示的に取り出す点に特徴がある。

自由度と平均平方

  • SSA:$a - 1$
  • SSB:$b - 1$
  • SSAB:$(a - 1)(b - 1)$
  • SSE:$N - ab$(各セルに1以上の観測がある場合)

総自由度は

\[N - 1 = (a-1) + (b-1) + (a-1)(b-1) + (N-ab)\]

と分解される。

平均平方は

\[\mathrm{MSA} = \frac{\mathrm{SSA}}{a-1},\quad\mathrm{MSB} = \frac{\mathrm{SSB}}{b-1}\]\[\mathrm{MSAB} = \frac{\mathrm{SSAB}}{(a-1)(b-1)},\quad\mathrm{MSE} = \frac{\mathrm{SSE}}{N-ab}\]

である。

F検定

それぞれの効果について以下のF統計量を定義する:

\[F_A = \frac{\mathrm{MSA}}{\mathrm{MSE}},\quad F_B = \frac{\mathrm{MSB}}{\mathrm{MSE}},\quad F_{AB} = \frac{\mathrm{MSAB}}{\mathrm{MSE}}\]

帰無仮説はそれぞれ:

  • $H_0^A : \alpha_i = 0$(要因Aの主効果なし)
  • $H_0^B : \beta_j = 0$(要因Bの主効果なし)
  • $H_0^{AB} : (\alpha\beta)_{ij} = 0$(交互作用なし)

これらのもとで、各F統計量は対応する自由度を持つF分布に従う。

交互作用の解釈

交互作用とは、「ある要因の効果が、他の要因の水準によって変化する」ことを意味する。すなわち、加法モデル

\[x_{ijk} \approx \mu + \alpha_i + \beta_j\]

では説明できない構造を表す項である。

交互作用が有意である場合、主効果の単独解釈は一般に適切ではなく、要因の組合せごとの比較(単純主効果の解析)が必要となる。

幾何学的解釈

観測ベクトル $\mathbf{x} \in \mathbb{R}^N$ に対し、二元配置分散分析は以下の直交分解として理解できる:

\[\mathbf{x} = \mathbf{x}_{A} + \mathbf{x}_{B} + \mathbf{x}_{AB} + \mathbf{e}\]

これに対応して、ノルムの分解

\[\|\mathbf{x} - \bar{x}\mathbf{1}\|^2= \|\mathbf{x}_A\|^2 + \|\mathbf{x}_B\|^2 + \|\mathbf{x}_{AB}\|^2 + \|\mathbf{e}\|^2\]

が成立する。この直交性が平方和分解の基盤である。

実務上の注意

  • 設計が不均衡($n_{ij}$ が不等)な場合、平方和の定義に注意(Type I, II, III)
  • 交互作用が有意な場合、主効果の解釈は慎重に行う
  • 等分散性が破れる場合、ロバスト手法や変換を検討

まとめ

二元配置分散分析は、2つの要因の主効果および交互作用を同時に評価する手法であり、総変動を4つの成分に分解することでその構造を明らかにする。線形モデルおよび直交射影として統一的に理解できる点において、統計的推論の中核的枠組みをなすものである。

Mathematics is the language with which God has written the universe.





















数理統計学 機械学習