ひよこ統計

3群以上の平均値の差(分散分析)

統計検定2級

3種類の接着剤A・B・Cで、接合部がはがれるまでの荷重に差があるかを確かめたいとします。手持ちの道具が2群を比べる 母平均の差の検定 だけなら、A対B・A対C・B対Cと3回に分けて検定することになります。ところがこの総当たりでは、1回あたりの有意水準を5%に取っても、3回のどれか1つで誤って帰無仮説を棄却してしまう確率は5%を大きく超えます。3回が互いに独立だと粗く見積もれば 1−0.953≈0.14 1 - 0.95^3 \approx 0.14 ですが、同じ標本を使い回す検定どうしは独立ではないので、これは目安であって正確な値ではありません。比べる回数は水準の数 a a に対して a(a−1)2 \frac{a(a-1)}{2} 回と急に増えるため、第1種の過誤と第2種の過誤 のうち第1種の過誤だけが積み上がっていきます。この積み上がりを避け、3つ以上の群の母平均を1回の検定でまとめて調べる方法が分散分析です。

平均の差を知りたいのに分散を分析する、という名前のねじれが引っかかるかもしれません。ここでは2種類の散らばりの大小関係を見ています。群と群の間の散らばりを、同じ群の中にある散らばりと比べ、前者が不釣り合いに大きければ、すべての群の母平均が等しいという帰無仮説を棄却します。一元配置実験 で分けた2つの平方和は、この2種類の散らばりを数値にしたものです。散らばりどうしの比を検定統計量へ組み立てる手順は F比 で扱います。

仮説の置き方も2群のときとは形が変わります。帰無仮説は H0:μ1=μ2=⋯=μa H_0: \mu_1 = \mu_2 = \dots = \mu_a 、つまりすべての水準で母平均が等しいという1つの主張にまとめ、対立仮説はその否定として、少なくとも1組の水準間で母平均が異なる、と置きます(帰無仮説と対立仮説)。この対立仮説を、すべての水準で母平均が異なる、と読み替えてはいけません。AとBだけがずれていてCはAと等しい、という状態でも帰無仮説は成り立っていないのですから、対立仮説の側がそれを含んでいなければ否定として不足します。

この判断は、群の中の散らばりを共通の物差しとして使うことで成り立っています。そのため、データの取られ方について次の3つを前提に置きます。前提が崩れても計算そのものは進んでしまうぶん、判断を誤ったことに気づきにくくなります。

帰無仮説が棄却されたとき言える範囲は、どこかに差がある、というところまでです。AとBなのか、BとCなのか、どの組に差があるかは分散分析だけでは分かりません。組まで特定したいときは、比較の回数を織り込んだ別の手続きへ進みます。逆に水準を2つ(a=2 a = 2 )まで減らすと、分散分析の結論は2群の両側の t t 検定と完全に一致し、F=t2 F = t^2 が成り立ちます(片側の t t 検定とは一致しません)。F分布 の t(ν)2∼F(1,ν) t(\nu)^2 \sim F(1, \nu) という関係がその根拠で、分散分析は t t 検定を置き換える別物ではなく、群の数を3つ以上へ広げた一般化だと捉えられます。

接着剤3種の比較を、総当たりの検定と分散分析で見比べる

接着剤A・B・Cのそれぞれで試験片を3枚ずつ作り、接合部がはがれるまでの荷重(kg)を測りました。

接着剤1枚目2枚目3枚目水準平均
接着剤A1353
接着剤B4686
接着剤C79119

総当たりで比べるなら A対B・A対C・B対C の3回の検定が要ります。分散分析はこの9個の数値をそのまま使い、H0:μ1=μ2=μ3 H_0: \mu_1 = \mu_2 = \mu_3 を1回で検定します。一元配置実験 と同じ手順で平方和を分けると SA=54 S_A = 54 ・Se=24 S_e = 24 が得られ、そこから作った F F は帰無仮説を棄却します(比の作り方と判定は F比)。検定を3回に分けなくても、同じ9個の数値から結論が出せたわけです。

ここで棄却されたのは、3つの母平均がすべて等しいという主張です。少なくとも1組の接着剤群の母平均に差があると言えるところまでは進みましたが、AとBの間に差があるのか、差があるのはBとCの間だけなのかは、この結果からは決まりません。また、この分析だけでは接着剤の違いが差の原因だとは言えず、因果解釈には無作為割付などの実験設計が必要です。

同じデータからAとCの2水準だけを取り出すと、水準数は2になります。各群の偏差平方和はどちらも 8 なのでプールした分散は sp2=(8+8)÷4=4 s_p^2 = (8 + 8) \div 4 = 4 、その平方根は sp=2 s_p = 2 です。t t 検定の統計量は t=(3−9)÷(21/3+1/3)≈−3.674 t = (3 - 9) \div (2\sqrt{1/3 + 1/3}) \approx -3.674 になります。一方この2水準の分散分析では SA=3(3−6)2+3(9−6)2=54 S_A = 3(3-6)^2 + 3(9-6)^2 = 54 (自由度1)と Se=8+8=16 S_e = 8 + 8 = 16 (自由度4)から VA=54 V_A = 54 ・Ve=4 V_e = 4 ・F=54÷4=13.5 F = 54 \div 4 = 13.5 が得られ、(−3.674)2=13.5 (-3.674)^2 = 13.5 と一致します。

入門問題

全 4 問
問 1

3つの群の母平均に差があるかを調べるために、2群の母平均の差の検定を有意水準5%で3回繰り返しました。この進め方の問題点として最も適切なものはどれですか。

問 2

一元配置分散分析の対立仮説として最も適切なものはどれですか。

問 3

接着剤A・B・Cを比べる一元配置分散分析で帰無仮説が棄却されました。この結果から言えることとして最も適切なものはどれですか。

問 4

2つの群の母平均の差を等分散を仮定した t t 検定で調べたところ、t=−2.5 t = -2.5 となりました。同じデータを一元配置分散分析にかけたときの F F 値として最も適切なものはどれですか。