tecktom.com

繰り返しのない二元配置

2026年10月17日9分で読めますAI生成教科書『資格試験』第15章 F 分布・等分散性の検定と分散分析
目次
  1. モデルと平方和
  2. ブロックとしての要因
  3. 確認テスト
  4. 木下の層

AI 生成・木下が未勉強

統計学の復習 第15章 §15.6 ★★ 各組合せで1回だけ測ると、交互作用と誤差を分けることができない。交互作用が無いと仮定し、残差を誤差とみなす。一方の要因をブロックと見れば、ブロックの違いを誤差から取り除く方法になる

モデルと平方和

各組合せ AiBjA_iB_j で1回だけデータを取る。r=1r=1 なので、定理 15.11 の SeS_e は 00 になり、自由度も ab(1−1)=0ab(1-1)=0 になる(確認テスト 15.5-3)。組合せの中のばらつきが観測されないので、誤差の分散 σ2\sigma^2 を直接推定する手段が無い。そこで交互作用が無いと仮定して、モデルを

xij=μ+ai+bj+εij,∑iai=0, ∑jbj=0,εij∼N(0,σ2) 独立(15.18)x_{ij}=\mu+a_i+b_j+\varepsilon_{ij},\qquad \sum_ia_i=0,\ \sum_jb_j=0,\qquad \varepsilon_{ij}\sim\mathrm N(0,\sigma^2)\ \text{独立}\tag{15.18}

とする。偏差の分解は式 (15.16) で r=1r=1 としたもので、

ST=SA+SB+Se,SA=b∑i(xˉi⋅−xˉ)2,SB=a∑j(xˉ⋅j−xˉ)2,Se=∑i,j(xij−xˉi⋅−xˉ⋅j+xˉ)2S_T=S_A+S_B+S_e,\qquad S_A=b\sum_i(\bar x_{i\cdot}-\bar x)^2,\quad S_B=a\sum_j(\bar x_{\cdot j}-\bar x)^2,\quad S_e=\sum_{i,j}(x_{ij}-\bar x_{i\cdot}-\bar x_{\cdot j}+\bar x)^2

である。交差項が消える証明は定理 15.11 と同じである。SeS_e は、繰り返しのある場合の交互作用の平方和と同じ式で、自由度は (a−1)(b−1)(a-1)(b-1) である。

モデル (15.18) のもとでは E[Ve]=σ2E[V_e]=\sigma^2 で、FA=VA/VeF_A=V_A/V_e を Fa−1,(a−1)(b−1)F_{a-1,(a-1)(b-1)} と、FB=VB/VeF_B=V_B/V_e を Fb−1,(a−1)(b−1)F_{b-1,(a-1)(b-1)} と比べる。しかし本当は交互作用 (ab)ij(ab)_{ij} があると、残差には誤差と交互作用が混ざり

E[Ve]=σ2+∑i,j(ab)ij2(a−1)(b−1)E[V_e]=\sigma^2+\frac{\sum_{i,j}(ab)_{ij}^2}{(a-1)(b-1)}

となる。分母が大きくなる分だけ FAF_A・FBF_B は小さくなり、検定は本当にある主効果を見逃しやすくなる。繰り返しが無いと交互作用を誤差と分けることができないので、交互作用が無いことを、データでなく実験の対象についての知識から判断して使う。

ブロックとしての要因

繰り返しのない二元配置がよく使われるのは、要因 BB が調べたい処理でなく、実験の日・材料のロット・被験者のように、ばらつきの原因として分かっているが、それ自体には関心が無いものの場合である。このような要因をブロックと呼ぶ。同じブロックの中ですべての処理を試せば、ブロックどうしの違いは SBS_B に分かれ、誤差 SeS_e に入らない。一元配置としてブロックを無視すると、ブロックの違いは群内変動に入り、VeV_e が大きくなって処理の違いが見えにくくなる(例題 15.7)。これは Fisher の3原則のうちの局所管理の考え方で、ブロックの中で処理を割り付ける順を無作為にした実験を乱塊法と呼ぶ。第16章と実験計画法の解説 第11章で扱う。

命題 15.12(2処理の繰り返しのない二元配置は対応のある tt 検定と同じである)

a=2a=2 の繰り返しのない二元配置で、差を dj=x1j−x2jd_j=x_{1j}-x_{2j}(j=1,…,bj=1,\dots,b)、その平均を dˉ\bar d、不偏分散を sd2s_d^2 とする。要因 AA の FAF_A は、第10章の対応のある差の tt 統計量 t=dˉ/(sd/b)t=\bar d/(s_d/\sqrt b) の2乗に等しい。

証明

xˉ1⋅−xˉ2⋅=dˉ\bar x_{1\cdot}-\bar x_{2\cdot}=\bar d で、xˉ\bar x はその中点なので xˉ1⋅−xˉ=dˉ/2\bar x_{1\cdot}-\bar x=\bar d/2、xˉ2⋅−xˉ=−dˉ/2\bar x_{2\cdot}-\bar x=-\bar d/2 である。したがって SA=b(dˉ2/4+dˉ2/4)=bdˉ2/2S_A=b(\bar d^2/4+\bar d^2/4)=b\bar d^2/2 で、ϕA=1\phi_A=1 である。残差は、xˉ⋅j=(x1j+x2j)/2\bar x_{\cdot j}=(x_{1j}+x_{2j})/2 から x1j−xˉ⋅j=dj/2x_{1j}-\bar x_{\cdot j}=d_j/2 なので、x1j−xˉ1⋅−xˉ⋅j+xˉ=dj/2−dˉ/2x_{1j}-\bar x_{1\cdot}-\bar x_{\cdot j}+\bar x=d_j/2-\bar d/2 である。i=2i=2 の残差はその符号を変えたものである。よって Se=2∑j(dj−dˉ)2/4=(b−1)sd2/2S_e=2\sum_j(d_j-\bar d)^2/4=(b-1)s_d^2/2、ϕe=b−1\phi_e=b-1、Ve=sd2/2V_e=s_d^2/2 で

FA=bdˉ2/2sd2/2=dˉ2sd2/b=t2F_A=\frac{b\bar d^2/2}{s_d^2/2}=\frac{\bar d^2}{s_d^2/b}=t^2

である。

第10章で、対応のあるデータは差をとって分析すると個体どうしの違いが消えることを見た。命題 15.12 は、その操作が「個体をブロックとした二元配置」と同じものであることを示している。

例題 15.7 ブロックを取り除くと処理の違いが見える ★★★

オリジナル3種類の教え方 A1,A2,A3A_1,A_2,A_3 を、4つのクラス B1,…,B4B_1,\dots,B_4 でそれぞれ1回ずつ試し、小テストの平均点を記録した。

B1B_1B2B_2B3B_3B4B_4
A1A_1810711
A2A_21013912
A3A_311131014
  1. 繰り返しのない二元配置として分散分析表を作り、有意水準 1%1\% で AA と BB を検定せよ(F2,6(0.01)=10.92F_{2,6}(0.01)=10.92、F3,6(0.01)=9.780F_{3,6}(0.01)=9.780)。
  2. クラスを無視して一元配置(各水準4個)として分析し、有意水準 5%5\% で AA を検定せよ(F2,9(0.05)=4.256F_{2,9}(0.05)=4.256)。(1) と結論が違う理由を述べよ。
解答

(1) 行の合計は 36, 44, 4836,\ 44,\ 48、列の合計は 29, 36, 26, 3729,\ 36,\ 26,\ 37、総合計は 128128 で、CT=1282/12=1365.33CT=128^2/12=1365.33、∑x2=1414\sum x^2=1414 から ST=48.67S_T=48.67 である。

SA=362+442+4824−CT=1384−1365.33=18.67,S_A=\frac{36^2+44^2+48^2}{4}-CT=1384-1365.33=18.67,

SB=292+362+262+3723−CT=1394−1365.33=28.67S_B=\frac{29^2+36^2+26^2+37^2}{3}-CT=1394-1365.33=28.67

で、Se=48.67−18.67−28.67=1.33S_e=48.67-18.67-28.67=1.33 である。

要因平方和自由度平均平方F0F_01%1\% 点
AA18.6718.67229.3339.33342.042.010.9210.92
BB28.6728.67339.5569.55643.043.09.7809.780
ee1.331.33660.2220.222
計48.6748.671111

AA も BB も有意水準 1%1\% で有意である。教え方によって平均点は異なる。

(2) 一元配置では SA=18.67S_A=18.67(同じ式)、Se=ST−SA=30.00S_e=S_T-S_A=30.00、ϕe=9\phi_e=9、Ve=3.333V_e=3.333 で、F0=9.333/3.333=2.80<4.256F_0=9.333/3.333=2.80<4.256 となり、5%5\% でも有意でない。(1) と (2) で SAS_A は同じである。違うのは誤差の大きさで、(2) ではクラスどうしの違い SB=28.67S_B=28.67 が群内変動に入って VeV_e が 1515 倍になった。どのクラスでも A1<A2≤A3A_1<A_2\le A_3 の順はほぼ揃っているが、クラスの平均の違い(B3B_3 は低く B4B_4 は高い)がそれより大きいので、ブロックを取り除かないと処理の違いが誤差に埋もれる。

確認テスト

オリジナル

  1. 15.6-1 AA が4水準・BB が5水準の繰り返しのない二元配置で、各要因・誤差・計の自由度を求めよ。
  2. 15.6-2 4人の被験者に2つの条件で作業させ、所要時間を測った。条件 A1A_1 は 5, 7, 6, 95,\ 7,\ 6,\ 9、A2A_2 は 7, 8, 9, 117,\ 8,\ 9,\ 11(同じ順の被験者)であった。繰り返しのない二元配置の FAF_A と、対応のある差の tt 統計量を計算し、命題 15.12 を確かめよ。
  3. 15.6-3 モデル (15.18) が正しいとき、E[Se]=(a−1)(b−1)σ2E[S_e]=(a-1)(b-1)\sigma^2 となることを、a=b=2a=b=2 の場合に SeS_e を εij\varepsilon_{ij} で表して確かめよ。
  4. 15.6-4 本当は交互作用があるのに繰り返しのない二元配置で分析すると、要因 AA の検定の第1種の誤りの確率は名目の α\alpha より大きくなるか小さくなるか。理由とともに答えよ。

木下の勉強ステータス

木下はまだこの記事を勉強していません。

AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com

← AI生成教科書『資格試験』の記事一覧へ