Analysis of Variance
（書きかけです）
分散分析とは、統計量がF分布に従うことを前提としたパラメトリック検定*1の一種です。３つ以上の群（水準）に対して、効果による分散（群間変動）と誤差による分散（群内変動）との比（F値）を F検定を用いて確認することで、各水準の母平均に誤差以上の効果があるか否かを判断する・・というものです。
分散分析を行うには、対象となるデータの分布が1) 正規分布に従っていること（正規性）、2) 分散の等質性があること、3) 観測値の独立性（無作為標本）があることが前提となります。
分散分析の帰無仮説と対立仮説を、3群の分散分析の場合を例に確認します。
以下、A・B・C の３社の製品の耐久性を比較する・・という例で、分散分析の用語について概説します。
データの値を変化させる原因を「要因」といいます。この場合「会社」です。要因の数は複数設定されることもあります。例えば、A・B・Cの３社の製品を、それぞれ、温暖な地域で使用した場合と、寒冷な地域で使用した場合とで比較する場合、要因は２つ。ひとつは「製造会社」もうひとつは「気候」です。この場合、3ｘ2 で６種類の評価実験が必要になります。
要因を構成する条件を「水準（群）」といいます。上の例で「会社」要因で言えば、A・B・Cの３つの「水準（群）」が存在します。
分散分析には大きく一元配置と二元配置（多元配置）の２種類があります。
要因と被験者の計悪を組み合わせて、「１要因被験者間計画」、「２要因被験者内計画」といった実験計画が行われます。
特定の要因単独で有意に差が認められるときは，主効果（または単純主効果）があるといいます。例えば「気候条件に関わらずA社の製品の評価が高い」や「製造会社に関わらず、温暖な地域の方が評価が高い」など
要因を組み合わせた場合の複合的な効果がある場合は、交互作用がある・・といいます。例えば「温暖な地域ではA社の製品の評価が高く、寒冷地ではC社の製品の評価が高い」など・・
多重比較（Multiple Comparison Procedure）とは水準間の平均値の差を比較する検定手法で、これを用いると、3水準以上の比較において、どの水準間に差があるかを調べることができます。
一般に、分散分析で有意と判断された場合、その下位検定（事後検定）として多重比較が行われますが、分散分析では有意なのに、多重比較では群間の有意差が見いだせない場合や、あるいはその逆の場合あり、分散分析と多重比較は独立に行うものと考えた方が良いとされています。
多重比較には様々な種類があって、分散分析を事前に行うことを前提とした Scheffe法や、分散分析とは独立に群間の有意差検定が可能な Bonferroni法、Tukey-Kramer法などがあります。
分散分析の結果は、一般に「分散分析表」の形で提示します。一元配置と二元配置では、異なるので、それぞれ事例を紹介します。
はじめに、分散分析表のイメージを紹介して、その後、その表ができるプロセスを具体的なデータを通して紹介します。
| 要因 | 平方和 S | 自由度 df | 平均平方 V | F値 |
| 群間 | 32 | 2 | 16.00 | 13.11 |
| 群内 | 11 | 9 | 1.22 | |
| 全体 | 43 | 11 |
この表は、以下のデータにもとづいて作成されたものです。
| メーカー | A社 | B社 | C社 |
| 3 | 7 | 5 | |
| 4 | 8 | 4 | |
| 4 | 9 | 6 | |
| 3 | 6 | 7 | |
| 群平均 | 3.5 | 7.5 | 5.5 |
| 全平均 | 5.5 | ||
被験者12名を４名づつ３つのグループ（群）に分けて、それぞれに A・B・C の３つのメーカーの製品の音質を評価してもらった結果です（一人につき一商品を、直感的に10段階で評価してもらうという単純なものです）。
表の左上から順に、どのような計算で値が決まるのかを概説します。
以下の関係があることをおさえると、理解しやすいかもしれません。
（ 各値 - 全平均 ）=（ 各値 - 群平均 ）+（ 群平均 - 全平均 ）
上記の３つの値は、この関係式の３つの項について、その平方（２乗）を、すべてのデータの数（12件）合計したものになっています（以下参考）。
（ 3 - 5.5 ）＝（ 3 - 3.5 ）+（ 3.5 - 5.5 ） ・・・１件目のデータ
（ 4 - 5.5 ）＝（ 4 - 3.5 ）+（ 3.5 - 5.5 ） ・・・２件目のデータ
：
（ 8 - 5.5 ）＝（ 8 - 7.5 ）+（ 7.5 - 5.5 ） ・・・6件目のデータ
：
（ 7 - 5.5 ）＝（ 7 - 5.5 ）+（ 5.5 - 5.5 ） ・・・12件目のデータ
↑ ↑ ↑
平方和(43) ＝ 平方和(11) + 平方和(32)
3 -1 = 2
（ 4 - 1 ）+ （ 4 - 1 ）+ （ 4 - 1 ）= 9
12 - 1 = 11
32 ÷ 2 = 16.00
11 ÷ 9 = 1.22
16.00 ÷ 1.22 = 13.11
という流れで表が埋まっていきます。
実際には、Python 等の統計ライブラリの利用で、一括出力されますが、なんとなく数字の意味がわかっていると、結果の考察で迷うことがありません。
F値は、群内変動よりも群間変動のほうが大きい場合にその値が大きくなります。F値大きな値をとる場合は、偶然による誤差よりも、群の要因による誤差のほうが大きい、すなわち群間に違いがある・・と判定できることになります。
F値は、F分布と呼ばれる確率分布に従うのですが、自由度によって分布の形状が変わります。p=0.05 の F値を列挙した F分布表では、分子の自由度（群間の自由度）が df1、分母の自由度（群内の自由度）が df2 として、縦横２軸のテーブルになっているので、その２つの値を追って、交差する部分を見つけます。
有意水準 α = 0.05 として、df1 = 2 と df2 = 9 が交差する部分の値を見ると 4.26、よってこの例では、棄却域は F ≧ 4.26 となります。Fの実現値は、13.11 なので、帰無仮説は棄却され「３つの群の平均値には、いずれかに有意な差がある」と言うことができます。
実際には Python 等の統計ライブラリの利用で、F値、P値が一括出力されるので、F分布表を目で追う必要はありません。出力された P値が、有意水準よりも小さな値であれば、有意と判定できます。
論文などではこの結果を、以下のように書きます。
\( F ( 2,9 ) = 13.11 ( p < .05 ) \)
\( F_{2,9} = 13.11 ( p < .05 ) \)
（書きかけです）