独立 でない確率変数 でも「和の期待値 は期待値 の和」は成り立ちます。同時確率 の行と列を合計するだけで証明でき、独立性 を仮定する必要がありません。一方、分散の加法には共分散項が残るため独立性 が要ります。この非対称の理由を式で確かめます。
stats-06 の「期待値の線形性 」節では、E [ X + Y ] = E [ X ] + E [ Y ] E[X+Y] = E[X]+E[Y] E [ X + Y ] = E [ X ] + E [ Y ] が X X X と Y Y Y の独立性 を要求しないという事実を述べ、証明を <details> の折りたたみで 4 行にまとめました。本補足ではその 4 行を展開し、「なぜ独立性 が要らないのか」「完全従属の極端例でも式が崩れないか」「分散の線形性 は独立性 を要求するのに期待値 だけが要らないのはなぜか」という 3 点を式で確かめます。
親記事: 期待値
なぜこの補足が必要か
E [ X + Y ] = E [ X ] + E [ Y ] E[X+Y] = E[X]+E[Y] E [ X + Y ] = E [ X ] + E [ Y ] に「独立 でなくても成り立つ」という注記が添えられていると、直感に反します。X X X と Y Y Y が絡み合っていれば、その関係が和の期待値 に影響しそうに見えるからです。独立 な場合は P ( X = x i , Y = y j ) = P ( X = x i ) ⋅ P ( Y = y j ) P(X=x_i, Y=y_j) = P(X=x_i) \cdot P(Y=y_j) P ( X = x i , Y = y j ) = P ( X = x i ) ⋅ P ( Y = y j ) という式が成り立ち、計算が単純になると知っている分、「独立性 がないと計算が崩れるのでは」という直感が働きます。
独立性 は不要です。その理由は式変形の構造にあります。証明のどこにも独立性 を使う余地がなく、同時確率 を周辺分布 に集約する操作だけで E [ X ] + E [ Y ] E[X]+E[Y] E [ X ] + E [ Y ] に到達します。以下の節でこの構造を式で追います。
X X X が取りうる値を x i x_i x i (i = 1 , 2 , … i = 1, 2, \ldots i = 1 , 2 , … )、Y Y Y が取りうる値を y j y_j y j (j = 1 , 2 , … j = 1, 2, \ldots j = 1 , 2 , … )とします。P ( X = x i , Y = y j ) P(X=x_i, Y=y_j) P ( X = x i , Y = y j ) は X = x i X = x_i X = x i かつ Y = y j Y = y_j Y = y j が同時に起きる確率で、「同時確率 」と呼びます。これは 2 次元の確率の表で、X X X の値を行、Y Y Y の値を列に並べると全マスの和が 1 になります。
まず X + Y X + Y X + Y の期待値 の定義式を書きます。確率変数 X + Y X + Y X + Y が取りうる値は x i + y j x_i + y_j x i + y j で、その確率は P ( X = x i , Y = y j ) P(X=x_i, Y=y_j) P ( X = x i , Y = y j ) ですから
となります。i i i と j j j の二重和は「すべての ( x i , y j ) (x_i, y_j) ( x i , y j ) の組み合わせ」について足し合わせることを意味します。
括弧の中 ( x i + y j ) (x_i + y_j) ( x i + y j ) を分配すると、和を 2 つに分けられます。
左の塊は x i x_i x i と P ( ⋅ ) P(\cdot) P ( ⋅ ) の積を全 ( i , j ) (i, j) ( i , j ) ペアで足し合わせています。右の塊は y j y_j y j と P ( ⋅ ) P(\cdot) P ( ⋅ ) の積を同じ範囲で足し合わせています。この分配は確率の大小や X X X と Y Y Y の関係によらず成立します。
左の塊に注目します。x i x_i x i は j j j に依存しないので ∑ j \sum_j ∑ j の外に出せます。
内側の和 ∑ j P ( X = x i , Y = y j ) \sum_j P(X=x_i, Y=y_j) ∑ j P ( X = x i , Y = y j ) は「X = x i X = x_i X = x i という条件のもとで、Y Y Y がどんな値を取ろうと全部足す」操作です。Y Y Y の値によらず足し合わせると Y Y Y の添字が消え、X = x i X = x_i X = x i の確率だけが残ります。これが周辺化 で、数学的には
となります。同時確率 の表で「X = x i X = x_i X = x i の行を横に合計する」だけの操作であり、X X X と Y Y Y がどんな関係を持っていようと、行の合計は X = x i X = x_i X = x i の確率に等しくなります。独立性 を使う余地はここにありません。
右の塊も対称的に処理できます。y j y_j y j は i i i に依存しないので
∑ i P ( X = x i , Y = y j ) = P ( Y = y j ) \sum_i P(X=x_i, Y=y_j) = P(Y=y_j) ∑ i P ( X = x i , Y = y j ) = P ( Y = y j ) は「Y = y j Y = y_j Y = y j の列を縦に合計する」周辺化です。これも独立性 と無関係です。
2 つの塊を合わせると
に到達します。証明の中で P ( X = x i , Y = y j ) = P ( X = x i ) P ( Y = y j ) P(X=x_i, Y=y_j) = P(X=x_i) P(Y=y_j) P ( X = x i , Y = y j ) = P ( X = x i ) P ( Y = y j ) (独立性 の定義式)を一度も使っていません。使ったのは「分配則」「j j j に依存しない量を ∑ j \sum_j ∑ j の外に出す」「行を合計すると周辺確率 になる」の 3 つだけです。
E [ X + Y ] = E [ X ] + E [ Y ] E[X+Y] = E[X]+E[Y] E [ X + Y ] = E [ X ] + E [ Y ] の証明で使うのは「周辺化」だけです。同時確率 P ( X = x i , Y = y j ) P(X=x_i, Y=y_j) P ( X = x i , Y = y j ) の行を横に合計すれば X X X の周辺確率 が出ます。列を縦に合計すれば Y Y Y の周辺確率 が出ます。関係構造(独立 か従属か)は和の計算に入り込みません。
完全従属の極端例で確かめる
「独立性 なしで成り立つ」を最も鋭く試す状況は、完全従属の場合です。X X X の値が決まれば Y Y Y が完全に決まる関係、つまり片方を知れば他方が 100% 予測できる状況で、線形性 が保たれるかを確認します。
例 1: Y = X Y = X Y = X (完全正従属)
サイコロを 1 個振り、出た目を X X X とします。Y = X Y = X Y = X と定義すると、X X X が 3 なら Y Y Y も 3 で、X X X と Y Y Y は完全に一致します。この場合の同時確率 は
Y \ X Y \backslash X Y \ X 1 2 3 4 5 6 1 1/6 0 0 0 0 0 2 0 1/6 0 0 0 0 3 0 0 1/6 0 0 0 4 0 0 0 1/6 0 0 5 0 0 0 0 1/6 0 6 0 0 0 0 0 1/6
確率が対角線にだけ並んでいます。行を横に合計すると各行が 1 / 6 1/6 1/6 (X X X の周辺確率 )、列を縦に合計すると各列が 1 / 6 1/6 1/6 (Y Y Y の周辺確率 )になります。
Y=X の同時確率分布(6×6グリッド) 1 2 3 4 5 6 X 1 2 3 4 5 6 Y
横軸が X X X (1〜6)、縦軸が Y Y Y (1〜6)です。青い円は確率 1 / 6 1/6 1/6 の組み合わせで、対角線上の 6 マスにのみ存在します。行を横に合計すると各行の和が 1 / 6 1/6 1/6 (X X X の周辺確率 )になり、列を縦に合計すると各列の和が 1 / 6 1/6 1/6 (Y Y Y の周辺確率 )になります。対角線以外の 30 マスは確率 0 です。
直接計算すると:X + Y = 2 X X + Y = 2X X + Y = 2 X なので E [ X + Y ] = E [ 2 X ] = 2 × 3.5 = 7 E[X+Y] = E[2X] = 2 \times 3.5 = 7 E [ X + Y ] = E [ 2 X ] = 2 × 3.5 = 7 です。
線形性 で計算すると:E [ X ] = 3.5 E[X] = 3.5 E [ X ] = 3.5 、E [ Y ] = E [ X ] = 3.5 E[Y] = E[X] = 3.5 E [ Y ] = E [ X ] = 3.5 なので E [ X ] + E [ Y ] = 7 E[X]+E[Y] = 7 E [ X ] + E [ Y ] = 7 となります。一致します。
例 2: Y = − X Y = -X Y = − X (完全負従属)
X X X を + 1 +1 + 1 と − 1 -1 − 1 を各確率 1 / 2 1/2 1/2 で取る確率変数 とし、Y = − X Y = -X Y = − X と定義します。X = + 1 X = +1 X = + 1 なら Y = − 1 Y = -1 Y = − 1 になり、X = − 1 X = -1 X = − 1 なら Y = + 1 Y = +1 Y = + 1 になります。E [ X ] = ( + 1 ) × 1 / 2 + ( − 1 ) × 1 / 2 = 0 E[X] = (+1) \times 1/2 + (-1) \times 1/2 = 0 E [ X ] = ( + 1 ) × 1/2 + ( − 1 ) × 1/2 = 0 、E [ Y ] = E [ − X ] = − E [ X ] = 0 E[Y] = E[-X] = -E[X] = 0 E [ Y ] = E [ − X ] = − E [ X ] = 0 です。
X + Y = X + ( − X ) = 0 X + Y = X + (-X) = 0 X + Y = X + ( − X ) = 0 で、X + Y X+Y X + Y は常に 0 になります。したがって E [ X + Y ] = E [ 0 ] = 0 E[X+Y] = E[0] = 0 E [ X + Y ] = E [ 0 ] = 0 です。
線形性 で計算すると:E [ X ] + E [ Y ] = 0 + 0 = 0 E[X]+E[Y] = 0 + 0 = 0 E [ X ] + E [ Y ] = 0 + 0 = 0 となります。こちらも一致します。
ここで一歩踏み込みます。X + Y = 0 X+Y = 0 X + Y = 0 は定数で、その分散は V a r [ X + Y ] = V a r [ 0 ] = 0 \mathrm{Var}[X+Y] = \mathrm{Var}[0] = 0 Var [ X + Y ] = Var [ 0 ] = 0 です。一方、V a r [ X ] = E [ X 2 ] − ( E [ X ] ) 2 = 1 − 0 = 1 \mathrm{Var}[X] = E[X^2] - (E[X])^2 = 1 - 0 = 1 Var [ X ] = E [ X 2 ] − ( E [ X ] ) 2 = 1 − 0 = 1 なので、単純に足すと V a r [ X ] + V a r [ Y ] = 1 + 1 = 2 \mathrm{Var}[X]+\mathrm{Var}[Y] = 1 + 1 = 2 Var [ X ] + Var [ Y ] = 1 + 1 = 2 になります。V a r [ X + Y ] = 0 \mathrm{Var}[X+Y] = 0 Var [ X + Y ] = 0 なのに V a r [ X ] + V a r [ Y ] = 2 \mathrm{Var}[X]+\mathrm{Var}[Y] = 2 Var [ X ] + Var [ Y ] = 2 という状況です。
期待値の線形性 が独立性 なしで成り立つからといって、分散の線形性 も独立性 なしで成り立つわけではありません。Y = − X Y = -X Y = − X の例で V a r [ X + Y ] = 0 \mathrm{Var}[X+Y] = 0 Var [ X + Y ] = 0 ですが V a r [ X ] + V a r [ Y ] = 2 \mathrm{Var}[X]+\mathrm{Var}[Y] = 2 Var [ X ] + Var [ Y ] = 2 で、2 つは一致しません。この非対称が次の節の主題です。
分散の定義から出発します。μ X = E [ X ] \mu_X = E[X] μ X = E [ X ] 、μ Y = E [ Y ] \mu_Y = E[Y] μ Y = E [ Y ] として
括弧の中を ( X − μ X ) + ( Y − μ Y ) (X - \mu_X) + (Y - \mu_Y) ( X − μ X ) + ( Y − μ Y ) と書き直し、二乗を展開すると
最初の 2 項は V a r [ X ] \mathrm{Var}[X] Var [ X ] と V a r [ Y ] \mathrm{Var}[Y] Var [ Y ] です。第 3 項に出てくる E [ ( X − μ X ) ( Y − μ Y ) ] E[(X-\mu_X)(Y-\mu_Y)] E [( X − μ X ) ( Y − μ Y )] を共分散 と呼び、C o v ( X , Y ) \mathrm{Cov}(X, Y) Cov ( X , Y ) と書きます。整理すると
となります。
共分散 C o v ( X , Y ) = E [ ( X − μ X ) ( Y − μ Y ) ] \mathrm{Cov}(X, Y) = E[(X-\mu_X)(Y-\mu_Y)] Cov ( X , Y ) = E [( X − μ X ) ( Y − μ Y )] は「X X X が平均より大きいときに Y Y Y も平均より大きい傾向の強さ」を数値化した量です。( X − μ X ) (X-\mu_X) ( X − μ X ) と ( Y − μ Y ) (Y-\mu_Y) ( Y − μ Y ) の積を平均しているので、両者が同じ方向に動けば正、逆方向に動けば負、関係がなければ 0 に近くなります。
X X X と Y Y Y が独立 なら、同時確率 が P ( X = x i , Y = y j ) = P ( X = x i ) P ( Y = y j ) P(X=x_i, Y=y_j) = P(X=x_i)P(Y=y_j) P ( X = x i , Y = y j ) = P ( X = x i ) P ( Y = y j ) に因数分解できるため、積の期待値 も因数分解できます。つまり E [ ( X − μ X ) ( Y − μ Y ) ] = E [ X − μ X ] ⋅ E [ Y − μ Y ] = 0 ⋅ 0 = 0 E[(X-\mu_X)(Y-\mu_Y)] = E[X-\mu_X] \cdot E[Y-\mu_Y] = 0 \cdot 0 = 0 E [( X − μ X ) ( Y − μ Y )] = E [ X − μ X ] ⋅ E [ Y − μ Y ] = 0 ⋅ 0 = 0 となります。独立 なら共分散は 0 ですから C o v ( X , Y ) \mathrm{Cov}(X, Y) Cov ( X , Y ) 項が消え、V a r [ X + Y ] = V a r [ X ] + V a r [ Y ] \mathrm{Var}[X+Y] = \mathrm{Var}[X]+\mathrm{Var}[Y] Var [ X + Y ] = Var [ X ] + Var [ Y ] が成り立ちます。従属の場合は共分散が 0 でなく残ります。
前節の Y = − X Y = -X Y = − X で確認します。μ − X = − μ X = 0 \mu_{-X} = -\mu_X = 0 μ − X = − μ X = 0 なので
V a r [ X ] = 1 \mathrm{Var}[X] = 1 Var [ X ] = 1 なので C o v ( X , − X ) = − 1 \mathrm{Cov}(X, -X) = -1 Cov ( X , − X ) = − 1 となります。これを分散の式に代入すると
前節で直接計算した V a r [ X + Y ] = 0 \mathrm{Var}[X+Y] = 0 Var [ X + Y ] = 0 と一致しました。
期待値 と分散が非対称になる根本的な理由は、演算の構造の違いにあります。期待値 E [ X + Y ] E[X+Y] E [ X + Y ] の計算で中心になるのは ( x i + y j ) (x_i + y_j) ( x i + y j ) という和 です。和を x i x_i x i と y j y_j y j に分けて別々に平均を取る操作は、X X X と Y Y Y の関係構造に干渉しません。一方、分散 V a r [ X + Y ] \mathrm{Var}[X+Y] Var [ X + Y ] の計算では ( X + Y − μ X − μ Y ) 2 (X+Y-\mu_X-\mu_Y)^2 ( X + Y − μ X − μ Y ) 2 を展開したときに ( X − μ X ) ( Y − μ Y ) (X-\mu_X)(Y-\mu_Y) ( X − μ X ) ( Y − μ Y ) という積 の項が出ます。積の項には X X X と Y Y Y が「一緒にどう変化するか」の情報が含まれるため、この項をゼロにするには独立性 の仮定が必要です。
期待値 は和だけで計算が完結します。和は分配・集約だけで整理でき、変数間の関係構造が入り込みません。分散は二乗(積)を含みます。二乗を展開したときに生まれる交差項が共分散であり、この項に X X X と Y Y Y の関係構造が凝縮されています。期待値 は独立性 を見ず、分散は独立性 を見ます。この非対称は演算の構造から必然的に生まれます。
結論として何が言えるか + 親記事に戻る
本補足で確認したことを整理します。
E [ X + Y ] = E [ X ] + E [ Y ] E[X+Y] = E[X]+E[Y] E [ X + Y ] = E [ X ] + E [ Y ] の証明は「周辺化」だけで完結し、独立性 を一度も使いません。同時確率 の行・列を合計するだけで、関係構造は計算に入りません。
Y = X Y = X Y = X (完全正従属)でも Y = − X Y = -X Y = − X (完全負従属)でも期待値の線形性 は成り立ちます。直接計算と線形性 の計算は常に一致します。
分散には共分散項 2 C o v ( X , Y ) 2\,\mathrm{Cov}(X, Y) 2 Cov ( X , Y ) が出ます。独立 なら共分散は 0 になり V a r [ X + Y ] = V a r [ X ] + V a r [ Y ] \mathrm{Var}[X+Y] = \mathrm{Var}[X]+\mathrm{Var}[Y] Var [ X + Y ] = Var [ X ] + Var [ Y ] が成り立ちますが、従属の場合は共分散が残るため分散の単純加法は成り立ちません。
「期待値 は独立性 不要」という性質の応用例として、ランダム順列の固定点期待値 があります。n n n 枚のカードをシャッフルしたとき、元の位置に留まるカード(固定点)の枚数の期待値 は 1 になります。18 世紀の数学者モンモール(1708 年)が問題として定式化したこの計算でも、各カードが固定点かどうかを示す指示変数 X k X_k X k (固定点なら 1、そうでなければ 0)を使い、E [ ∑ k X k ] = ∑ k E [ X k ] = ∑ k 1 / n = 1 E[\sum_k X_k] = \sum_k E[X_k] = \sum_k 1/n = 1 E [ ∑ k X k ] = ∑ k E [ X k ] = ∑ k 1/ n = 1 と線形性 で一行に導けます。n n n 枚のカードが互いに独立 でないにもかかわらず、期待値の線形性 が直接使えるためです。同じ仕組みは機械学習のアンサンブル学習にも現れます。複数モデルの予測が互いに独立 でなくても、それぞれの予測期待値 の和として全体の期待値 が計算できます(ただし分散の縮小効果は独立性 に依存します)。
Part 3 で扱う標準誤差では、独立 な確率変数 の和の分散が V a r [ X 1 + X 2 + ⋯ + X n ] = n V a r [ X 1 ] \mathrm{Var}[X_1+X_2+\cdots+X_n] = n\,\mathrm{Var}[X_1] Var [ X 1 + X 2 + ⋯ + X n ] = n Var [ X 1 ] となる事実を使います。これは分散の加法に独立性 を要求した上で成り立つ結果です。本補足で見た「分散には共分散項が残る」という制約が、標準誤差の導出で独立性 の仮定を必要とする直接の理由になります。
親記事: 期待値 (線形性 の節、本補足の展開元)
関連する補足: サンクトペテルブルクのパラドックス (期待値 が発散する極端例)
関連する本編: 標準誤差(Part 3 で扱う、分散の独立 加法が使われる先)