ひよこ統計

反復

統計検定2級

一元配置実験 を計画するとき、条件を何種類そろえるかと並んで決めなければならないのが、1つの条件につき何個の実験単位を試すかです。同じ水準を複数の実験単位に割り付けることを反復といい、その個数を繰り返し数と呼びます。フィッシャーの3原則 では3本柱の1つとして挙がりますが、ここでは繰り返し数が分散分析の計算のどこに現れるかを追います。

水準の数を a a 、各水準の繰り返し数を r r とします。データは全部で ar ar 個集まり、そのうち a a 個ぶんは水準ごとの平均を推定するために使われるので、残りが残差の自由度になります。ここで r=1 r = 1 、つまり各水準を1個ずつしか試さないと、この自由度は 0 になります。残差平均平方 Ve V_e は残差平方和 Se S_e をその自由度で割った量なので計算できません。水準間の差のほうも平方和を自由度で割って水準間平均平方 VA V_A とし、判定はこの2つの比 F=VA/Ve F = V_A / V_e で行うため、分母が無ければ差を比べる物差しそのものが手に入らないのです(比の作り方そのものは F比 で扱います)。反復が要る理由として、これが最もはっきりしています。

n=ar,n−a=a(r−1)n = ar, \qquad n - a = a(r-1)
a a :水準の数(比べる条件の種類)、r r :各水準の繰り返し数、n n :総データ数、n−a n - a :残差の自由度

同じ実験単位を測り直すことは反復に数えません。1枚の布の汚れを4回測れば数値は4つ得られますが、そこに現れるのは測定のぶれだけで、布どうしの違いは含まれていません。これを r=4 r = 4 として計算すると Ve V_e は本来より小さくなり、水準間の差が実際以上に大きな F F 値として出てしまいます。

繰り返しを増やす利点は、判定の側にも及びます。Ve V_e は自由度 a(r−1) a(r-1) で推定されるため、r r が大きいほど推定が安定します。加えて水準平均の 標準誤差 は Ve/r \sqrt{V_e / r} の形で表されるので、r r を4倍にすればこの値は半分になります。水準間の差が同じ大きさでも F F は大きくなり、帰無仮説が棄却されやすくなる、すなわち 検出力 が上がるということです。ただし総データ数は r r に比例して増えるため、どこまで反復するかは手間との兼ね合いで決めます。

洗剤3種を1枚ずつ試した場合と4枚ずつ試した場合

洗剤A・B・Cで、同じ汚れを付けた布を洗い、残った汚れの量(mg)を測ります。まず、洗剤ごとに布を1枚ずつ、r=1 r = 1 で試した結果です。

洗剤A洗剤B洗剤C
残った汚れ (mg)687

Aが最も少ないように見えます。しかしこの3つの数字からは、洗剤を変えたせいなのか、布や汚れの付き方が偶然そうなっただけなのかを切り分けられません。残差の自由度は 3−3=0 3 - 3 = 0 で、同じ洗剤の中にどれだけの幅があるかを示すデータが1つも無いからです。

次に、洗剤ごとに布を4枚ずつ、r=4 r = 4 に増やします。

1枚目2枚目3枚目4枚目平均
洗剤A56676
洗剤B78898
洗剤C67787

水準内の偏差はどの洗剤も −1, 0, 0, 1 -1,\ 0,\ 0,\ 1 なので、2乗和は 1+0+0+1=2 1+0+0+1 = 2 ずつ、合わせて残差平方和は 2+2+2=6 2+2+2 = 6 。自由度は 12−3=9 12 - 3 = 9 なので、残差平均平方は 6÷9=23 6 \div 9 = \frac{2}{3} です。総平均は 7 で、水準平均の偏差は −1, 1, 0 -1,\ 1,\ 0 だから水準間平方和は 4×(1+1+0)=8 4 \times (1+1+0) = 8 、自由度 2 で VA=4 V_A = 4 。F比は 4÷23=6 4 \div \frac{2}{3} = 6 となり、F0.05(2, 9)=4.26 F_{0.05}(2,\ 9) = 4.26 を超えるので、洗剤による差があると判断できます。

同じ12個の数値でも、1枚の布を12回測って集めたものだったら結論は変わります。測り直しの数値は 6.0, 6.1, 5.9, 6.0 のようにほとんど動かず、Se S_e は極端に小さくなります。分母が小さければF比は大きくなりますが、それは洗剤の差が大きいからではなく、洗剤の差を布どうしの幅ではなく測定器のぶれと比べてしまったからです。

入門問題

全 4 問
問 1

一元配置実験における反復の説明として最も適切なものはどれですか。

問 2

水準数が 3、各水準の繰り返し数が 4 の一元配置実験を行いました。残差の自由度はいくつですか。

問 3

各水準を1個ずつしか試さない(繰り返し数が 1 の)実験について、正しい記述はどれですか。

問 4

水準数を変えずに繰り返し数を増やしたときに起こることとして、最も適切なものはどれですか。