分散分析

分散分析（ぶんさんぶんせき、英: analysis of variance、略称: ANOVA）は、観測データにおける変動を誤差変動と各要因およびそれらの交互作用による変動に分解することによって、要因および交互作用の効果を判定する、統計的仮説検定の一手法である。

分散分析の最も単純な形は，2つ以上の標本の母集団が等しいかどうかを判断するt-検定であり、分散分析ではこれを一般化したものである。

統計学者で遺伝学者のロナルド・フィッシャーによって1920年代から1930年代にかけて基本手法が確立された。そのため「フィッシャーの分散分析」「フィッシャーのANOVA法」とも呼ばれる。

基本的な手法として、まず、データの分散成分の平方和を分解し、誤差による変動から要因効果による変動を分離する。次に、平方和を自由度で割ることで平均平方を算出する。そして、要因効果（または、交互作用）によって説明される平均平方を分子、誤差によって説明される平均平方を分母とすることでF値を計算する（F検定）。各効果の有意性については有意水準を設けて判定する。

交互作用の性質を詳しく調べるには、単純主効果の検定や交互作用対比を行うとよい。また、3つ以上の水準を持つ要因の効果が有意であったとき、具体的にどの群とどの群の間に差があったかを知るためには、多重比較を行う必要がある。 したがって、分析の目的によっては、分散分析のみから結論が導かれるものではなく、これらの手法と組み合わせて用いることが肝要である。

分散分析には各種のモデルがあり、データの性質や要因計画の型、検証したい仮説に応じてそれらを使い分けることが適切な利用法である（一元配置分散分析・回帰分散分析・共分散分析など）。現在では、分散分析は一般線形モデル、構造方程式モデリングの一部として扱えることが判明しており、さらなる拡張も可能である（潜在変数に対する分散分析など）。

分散方法は、実験データの分析やモデルの開発によく用いられる。相関関係と比べると、数値変数でない変数を扱えるなどの利点がある。

基本概念と用語

分散分析は，実験データの解析に多用される統計的仮説検定の一形態である。 帰無仮説が正しいと仮定して、帰無仮説と標本から計算された検定結果が、偶然に発生した可能性が低いと判断される場合、”統計的に有意”と呼ばれる。 帰無仮説の発生確率が低い場合、確率(P値)が事前に指定されたしきい値(有意水準)よりも小さい場合、統計的に有意な結果によって帰無仮説は棄却される。

分散分析の典型的な使い方は、「すべてのグループが同じ母集団からのランダムなサンプルである」という帰無仮説である。 たとえば，患者の類似したサンプルに対する異なる治療法の効果を研究するとき，帰無仮説を「すべての治療法が同じ効果を持つ（何も起こらないという結果になる）」とする。 帰無仮説を棄却されると、「治療グループ間で観察された効果の違いが、無作為な偶然によるものである可能性が低い（＝効果の違いは治療に原因がある）」ことを意味する。

仮説検定では、第一種の過誤（偽陽性）の発生率を有意水準に制限する。実験者は，第二種の過誤（偽陰性）も制限したいと考える。 第二種の過誤の発生率は、サンプルサイズ、有意水準、効果量に大きく依存する。 （サンプルが小さいほど発生率は大きくなる。有意水準が高いと、差異を見落とす可能性も高くなる。効果量が小さいほど第2種エラーが発生しやすくなる。）

分散分析の用語は、多くが実験計画法に由来する。 実験者は、効果を決定するために、要因を調整し反応を測定する。 結果の妥当性を保証するために，ランダム化とブロッキングの組み合わせによって因子は実験ユニットに割り当てられる。 また実験対象によっては盲検化を行い、結果の公平性を確保する。 応答は、効果による偏りとランダムなエラーの複合で変動する。

分散分析は、いくつかのアイデアの合成であり、多くの目的に使用され，簡潔または正確に定義することは困難である。

バランスのとれたデータのための "古典的 "分散分析は、次の3つを同時に行う。

1.探索的データ分析として、分散分析は加法的データ分解を採用し、その二乗和は分解の各成分（または同等に、線形モデルの各項のセット）の分散を示す。 2.平均二乗の比較とF検定により、モデルの入れ子構造を検証することができる。 3.係数推定値と標準誤差を伴う線形モデルのフィッティングに密接に関連する。

分析の例

例としてドッグショーに出品される犬について、その特徴に基づいて犬の体重を予測することを考える。 ドッグショーとは純血犬種の品評会であり、その犬種の基本的な犬種標準（スタンダード）に最も近いことを競うものである。

ドッグショーに出品される全部の犬の体重のヒストグラムは、かなり複雑なものになる。 犬を特徴に基づきグループ分けすることで、体重の分布を予測する。 グループ分けが上手くいったときには（a）各グループの分散が小さく（グループが比較的均質さを意味する）、（b）各グループの平均値が異なる（2つのグループの平均値が同じであればグループ分けが意味をなしていない）。

・若い犬と年老いた犬、短毛の犬と長毛の犬という2つの2値グループの積（相互作用）によって犬をグループ化することを考える。 ⇒　あまり予測として役に立たなさそう。このとき各グループ内の犬の体重分布は比較的大きな分散を持ち、平均値はグループ間で非常によく似ている。これらの特徴で犬をグループ化しても、犬の体重の変動を説明する効果的な方法にはならない。

・体重分布をペット犬と作業犬、運動量の少ない犬と多い犬でグループ化することを考える。 ⇒　ある程度うまくいきそう。最も重い犬達は、大きくて強い作業犬種である可能性が高く、ペットとして飼われている犬種は小さくて軽い傾向がある。 分散は1つ目のケースよりもかなり小さく、グループ間での平均値の差もよりはっきりする。しかし分布の重なりは依然大きく、予測の精度はまだ低い。

・体重を犬種でグループ化することを考える。 ⇒非常によく予測が当てはまりそう。 チワワはみんな軽く、セントバーナードは総じて重い。 セッターとイングリッシュ・ポインターで体重の違うと言っても、直感的な説明では説得力に欠くとき、分散分析を用いることで論理的に違うことを説明できる。

ソフトウェア

SASやSPSSといった主要な統計パッケージで、分散分析も実行可能である。R言語にも、分散分析に関わる関数がある。また、分散分析やそれに伴う多重比較に特化したソフトウェアもあり、多くはフリーソフトである。

js-STAR^[1]

田中敏（信州大学教授）作成による"STAR"をJavaScriptに移植したもの。3要因までの分散分析、単純主効果の検定および多重比較（LSD法、HSD法、Bonferroni法、Holm法）が一度にできる。また、その他にχ²検定や相関係数なども扱うことができる。ウェブ上でそのまま使うことができ、ダウンロードすることもできる。インターフェイスがシンプルで、使い方も分かりやすい。仕様の理論的背景は、田中敏・山際勇一郎による『ユーザーのための教育・心理統計と実験計画法－方法の理解から論文の書き方まで－』（教育出版、1992年、新訂版）に基づくと思われる^[誰?]。『実践データ心理解析－問題の発想・データ処理・論文の作成－』（新曜社、2006年、改訂版）には、js-STARの使用法、分散分析表の読み取り方、論文への記述の仕方などが詳しく解説されている。

ANOVA4 on the Web^[2]

桐木建始（広島女学院大学教授）が作成。4要因までの分散分析および多重比較（Ryan法）が一度に可能。ブラウザ上でそのまま動作し、インストール不要。また、『わかって楽しい心理統計法入門』（北大路書房、2007年）では、ANOVA4の使い方が解説されている。

分散分析プログラム（AIST-ANOVA）^[3]

独立行政法人産業技術総合研究所計測標準総合センター計測標準部門物性統計科応用統計研究室で開発された、Excelのアドイン。10要因までの分散分析を行うことができる。

Excel NAG 統計解析アドイン^[4]

Excelに分散分析関数を追加するアドイン。アカデミック版、無料の試供版もあり。

ezANOVA^[5]

単独で動作するソフトなので、ブラウザやExcelを要しない。

MAANOVA^[6]

MATLAB上で分散分析を行う。

Origin

Origin上で分散分析を行う。一元配置／二元配置、繰返しのある配置、Post-hoc検定などを含む。

その他、汎用言語であるC言語のプログラム^[7]、R言語で独自に書かれた関数^[8]などもある。それぞれのプログラミング言語が使える計算機環境と操作能力が必要である。

脚注

関連項目

• 推計統計学
• t検定
• 相関係数
• 回帰分析
• 一元配置分散分析

外部リンク

• 私のための統計処理 分散分析
• Visual ANOVA 分散分析のイメージを説明するFlashプログラム
• 反復測定（度）分散分析/基礎と応用　計算法の解説、プログラムのダウンロード。多重比較や球面性の仮定・球面性検定にも詳しい。
• 心理生理学データの分散分析　反復測定分散分析の使用や論文中への書き方の解説。多重比較や球面性検定の説明も。

表 話 編 歴 統計学
標本調査	標本 母集団 無作為抽出 層化抽出法
要約統計量	連続確率分布 位置 平均 算術 幾何 調和 中央値 分位数 順序統計量 最頻値 階級値 分散 範囲 偏差 偏差値 標準偏差 標準誤差 変動係数 決定係数 相関係数 自己相関 共分散 自己共分散 分散共分散行列 百分率 統計的ばらつき モーメント 分散 歪度 尖度 カテゴリデータ 頻度 分割表
推計統計学	仮説検定 パラメトリック t検定 ウェルチのt検定 F検定 Z検定 二項検定 ジャック-ベラ検定 シャピロ–ウィルク検定 分散分析 共分散分析 ノンパラメトリック ウィルコクソンの符号順位検定 マン・ホイットニーのU検定 カイ二乗検定 イェイツのカイ二乗検定 累積カイ二乗検定 フィッシャーの正確確率検定 尤度比検定 G検定 アンダーソン–ダーリング検定 コルモゴロフ–スミルノフ検定 カイパー検定 マンテル検定 コクラン・マンテル・ヘンツェルの統計量 その他 帰無仮説 対立仮説 有意 棄却 区間推定 信頼区間 予測区間 モデル選択基準 AIC BIC WAIC MDL その他 偏り 偏りと分散 過剰適合 推定量 点推定 最尤推定 尤度関数 尤度方程式 最小距離推定 メタアナリシス ブートストラップ法
ベイズ統計学	確率 主観確率 ベイズ確率 事前確率 事後確率 最大事後確率 その他 ベイズ推定 ベイズ因子
相関	交絡 ピアソンの積率相関係数 順位相関（スピアマンの順位相関係数 ・ケンドールの順位相関係数 ）
モデル	一般線形モデル 一般化線形モデル 混合モデル 一般化線形混合モデル
回帰	線形 リッジ回帰 ラッソ回帰 エラスティックネット 非線形 k近傍法 決定木 ランダムフォレスト ニューラルネットワーク サポートベクターマシン 射影追跡回帰 時系列 自己回帰モデル 自己回帰移動平均モデル ARCHモデル 対移動平均比率法 トレンド定常 傾向推定 共和分 構造変化
分類	線形 線形判別分析 ロジスティック回帰 <! -- 名前に回帰とついていますが確率を回帰する分類手法です --> 単純ベイズ分類器 単純パーセプトロン 線形サポートベクターマシン 二次 二次判別分析 非線形 k近傍法 決定木 ランダムフォレスト ニューラルネットワーク サポートベクターマシン ベイジアンネットワーク 隠れマルコフモデル その他 二項分類 多クラス分類 第一種過誤と第二種過誤
教師なし学習	クラスタリング k平均法 （k-means++法 ） DBSCAN 密度推定（英語版） カーネル密度推定 （ カーネル ） その他 主成分分析 独立成分分析 自己組織化写像
統計図表	棒グラフ バイプロット（英語版） 箱ひげ図 管理図 フォレストプロット ヒストグラム 円グラフ Q-Qプロット ランチャート 散布図 幹葉表示 バイオリン図 ドットプロット ヒートマップ 階級区分図
生存分析	生存関数 カプラン＝マイヤー推定量 ログランク検定 故障率 比例ハザードモデル
歴史	統計学の創始者 確率論と統計学の歩み
応用	社会統計学 疫学 生物統計学 系統学 統計力学 計量経済学 機械学習 実験計画法
出版物	統計学に関する学術誌一覧 重要な出版物
全般	統計 頻度主義統計学 統計学および機械学習の評価指標
その他	方向統計学 S言語 R言語 統計検定 社会調査士 JDLA Deep Learning For GENERAL JDLA Deep Learning for ENGINEER 実用数学技能検定 品質管理検定
カテゴリ