生物統計学 学習ノート概念・導出から実践的な解釈まで

CHAPTER 06 · TOPIC 04

二元配置ANOVAとMANOVA

二因子の主効果・交互作用と、複数目的変数の平均ベクトルを扱うMANOVAを学びます。

このページの内容
  1. 二元配置ANOVA:二因子を同時に調べる
  2. 二元配置ANOVAモデル
  3. 交互作用を先に見る理由
  4. 二元配置ANOVAの条件
  5. MANOVA:複数の目的変数を同時に調べる
  6. MANOVAの検定の作り方
  7. MANOVAの条件と限界
  8. 有意なMANOVAの追跡
  9. 三つの方法を区別する

二元配置ANOVAと多変量分散分析(MANOVA)は異なる方向の拡張です。二元配置ANOVAは一つの量的目的変数に対する二つのカテゴリー因子を扱い、MANOVAは関連する複数の量的目的変数を同時に扱います。

二元配置ANOVA:二因子を同時に調べる#

治療法を因子A、性別を因子B、血圧を目的変数とすれば、Aの主効果、Bの主効果、A×B交互作用を同じモデル内で検定できます。

効果問う内容
因子Aの主効果Bの水準で平均したとき、Aの水準間で平均が異なるか
因子Bの主効果Aの水準で平均したとき、Bの水準間で平均が異なるか
A×B交互作用Aの効果がBの水準によって変わるか

二元配置ANOVAモデル#

Yijk=μ+αi+βj+(αβ)ij+εijkY_{ijk}=\mu+\alpha_i+\beta_j+(\alpha\beta)_{ij}+\varepsilon_{ijk}
記号意味
μ全条件の総平均
αᵢ因子Aの第i水準の効果
βⱼ因子Bの第j水準の効果
(αβ)ᵢⱼ水準組合せの交互作用
εᵢⱼₖモデルで説明されない誤差
A主効果H0,A:α1==αa=0H_{0,A}:\alpha_1=\cdots=\alpha_a=0
B主効果H0,B:β1==βb=0H_{0,B}:\beta_1=\cdots=\beta_b=0
交互作用H0,AB:(αβ)ij=0for all i,jH_{0,AB}:(\alpha\beta)_{ij}=0\quad\text{for all }i,j

交互作用を先に見る理由#

交互作用が有意なら、一方の因子効果は他方の水準に依存します。たとえば薬剤が一方の性別でだけ有効なら、性別を平均した薬剤主効果だけではデータ構造を隠します。通常は交互作用を先に確認し、その後に主効果を解釈します。

有意な交互作用は、各性別内で治療法を比べるような単純主効果へ追跡します。これらの比較にも多重性調整が必要です。

二元配置ANOVAの条件#

  • 観測値が独立。反復測定なら反復測定法または混合モデルを使う
  • 各因子組合せ内の誤差が概ね正規
  • 各組合せで誤差分散が等しい
  • 交互作用を推定できる十分な観測が各セルにある

MANOVA:複数の目的変数を同時に調べる#

収縮期血圧、拡張期血圧、心拍数のように、同じ対象から関連する複数の量的結果を測る場合、MANOVAは各群の単一平均ではなく母平均ベクトルを比較します。

第j群の平均ベクトルμj=(μj1,μj2,,μjp)T\boldsymbol{\mu}_j=(\mu_{j1},\mu_{j2},\ldots,\mu_{jp})^{\mathsf T}
帰無仮説H0:μ1=μ2==μkH_0:\boldsymbol{\mu}_1=\boldsymbol{\mu}_2=\cdots=\boldsymbol{\mu}_k

目的変数ごとに別々のANOVAを行うと多重検定の第一種過誤が増え、結果の組合せとしての差も直接検定できません。MANOVAは目的変数間の共分散を利用します。

MANOVAの検定の作り方#

単変量ANOVAの群間・群内平方和を、平方和・交差積行列H(効果)とE(誤差)へ拡張し、その固有値から多変量統計量を作ります。

統計量判読
Pillai's trace仮定違反に比較的頑健。大きいほど群差が大きい
Wilks' lambda小さいほど平均ベクトルの群差が大きい
Hotelling–Lawley trace複数の判別方向の群間情報を累積
Roy's largest root最大差を持つ一つの判別方向に注目

MANOVAの条件と限界#

  • 観測値が独立
  • 各群の目的変数ベクトルが概ね多変量正規
  • 各群の共分散行列が等しい
  • 目的変数は関連するが、ほぼ重複するほど共線的ではない
  • 複数目的変数の共分散を推定できる十分な標本数がある

目的変数がほぼ無関係なら統合の利点が乏しく、相関が極端に高ければ同じ概念を重複測定している可能性があります。目的変数はデータを見た後ではなく研究課題から事前に選びます。

有意なMANOVAの追跡#

  1. どの因子または交互作用の多変量検定が有意か確認
  2. 研究課題に沿って各目的変数のANOVAや計画対比を行い多重性を調整
  3. 群平均ベクトル、信頼区間、効果量から方向と実質的意味を評価
  4. 必要なら判別関数や正準変量を調べ、原尺度の目的変数へ戻して解釈

三つの方法を区別する#

方法カテゴリー因子量的目的変数主な問い
一元配置ANOVA11一因子の群平均が異なるか
二元配置ANOVA21二つの主効果と交互作用があるか
MANOVA1以上2以上群の平均ベクトルが異なるか