同時確率分布の分散、共分散の導出がわかる(その１ 離散系の場合)
「同時確率分布の分散、共分散の導出がわからない」、と困っていませんか？
こういう疑問に答えます。
本記事のテーマ
- ①【共通】２段サンプリングの分散公式を導出するために知っておくべき内容
- ➁【何度も復習しよう！】離散型確率分布の場合
- ➂【何度も復習しよう！】連続型確率分布の場合(その２で解説)
サンプリングの分散公式への道ですが、徐々に難しくなっていきます。１つずつ理解してクリアーしましょう。
①【共通】２段サンプリングの分散公式を導出するために知っておくべき内容
2段サンプリングの分散の式
「２段サンプリングの分散」の式があります。
E(\(\bar{\bar{x}}\))=μ
V(\(\bar{\bar{x}}\))=\(\frac{M-m}{M-1}・\frac{σ_b^2}{m}\)+\(\frac{N-n}{N-1}・\frac{σ_w^2}{mn}\)
・\(m\)：1次サンプルの大きさ
・\(n\)：2次サンプルの大きさ
・\(σ_b^2\)：1次単位間の特性xの分散
・\(σ_w^2\)：1次単位内の特性xの分散
・M:1次単位の総数
・N:1次単位の大きさ
・\(\frac{M-m}{M-1},\frac{N-n}{N-1}\)：有限修正項
となりますよね。
でも、
と困ってしまいますよね。QCプラネッツも苦労しました。
そこで、
という思いで、解説していきます。
2段サンプリングの分散の式に必要な内容
まとめると、以下を理解しておく必要があります。
- 条件付き確率
- ２変数の確率分布関数（同時確率質量関数）
- 同時確率分布の分散、共分散の導出
- 条件付き確率の期待値・分散
- 全分散の公式の導出
- ２段サンプリングの分散の公式導出
残念ながら、「Yes」です。
だから、「２段サンプリングの分散」の式を暗記して代入する問題だけ出ます。
だから、サンプリングの分散はみんな苦手なのです!
2段サンプリングの分散の式の丁寧な導出はQCプラネッツだけ
「だから、教科書やサイトに、2段サンプリングの分散の式を導出する内容が書いているはず」と懇願しても、残念、ありませんでした。
では、１つ１つ解説します。
本記事のテーマ(再掲)
第3弾として「同時確率分布の分散、共分散の導出」を解説します。
- 条件付き確率
- ２変数の確率分布関数（同時確率質量関数）
- 同時確率分布の分散、共分散の導出
- 条件付き確率の期待値・分散
- 全分散の公式の導出
- ２段サンプリングの分散の公式導出
●2変数の同時確率質量関数については、関連記事で解説しています。ご確認ください。
 |
２変数の確率分布関数（同時確率質量関数）がわかる ２変数の確率分布関数（同時確率質量関数）が説明できますか？本記事では、2変数の確率分布関数の基礎をわかりやすく解説します。サンプリングの分散、全分散の公式導出に必須です。 |
➁【何度も復習しよう！】離散型確率分布の場合
●まず、わかりやすい「離散型」の場合で、数列∑を使った計算を解説します。
例題
●2次元の確率変数(X,Y)が、下表のような分布を持っている。
| X/Y | 1 | 2 | 3 | 計 |
| 1 | \(\frac{2}{8}\) | \(\frac{1}{8}\) | \(\frac{1}{8}\) | \(\frac{1}{2}\) |
| 2 | \(\frac{1}{8}\) | \(\frac{1}{8}\) | \(\frac{2}{8}\) | \(\frac{1}{2}\) |
| 計 | \(\frac{3}{8}\) | \(\frac{2}{8}\) | \(\frac{3}{8}\) | 1 |
(2)分散V[X],V[Y],共分散COV[X,Y]を求めよ。
期待値と分散のフルセットを計算してみましょう。
解法に必要な公式集
離散系の場合の期待値と分散の解法に慣れるために必要な公式集をまとめます。以下の式を使って、解いていきます。なお、連続系の場合は∑を∫に変えればOKです。
期待値の公式
●E[Y]=∑X・Pr(Y)
●E[X+Y]=∑(X+Y)・Pr(X+Y)
●E[XY]=∑XY・Pr(XY)
分散の公式
●V[Y]=E[\((Y-μ_Y)^2\)]
●COV[X,Y]=E[\((X-μ_X)(Y-μ_Y)\)]
●Cov[X,Y]= E[XY]- E[X]E[Y]
解法(期待値)
では、解いていきましょう。
E[X]の解法
表から、X=1の確率が1/2、X=2の確率が1/2ですから期待値は、
E[X]=1×1/2+2×1/2=3/2
簡単ですね！
E[Y]の解法
表から、Y=1の確率が3/8、Y=2の確率が2/8、Y=3の確率が3/8ですから期待値は、
E[Y]=1×3/8+2×2/8+3×3/8=2
簡単ですね！
E[X+Y]の解法
X+Yの場合について下表を追加しましょう。
| X/Y | 1 | 2 | 3 |
| 1 | X+Y=2 | X+Y=3 | X+Y=4 |
| \(\frac{2}{8}\) | \(\frac{1}{8}\) | \(\frac{1}{8}\) | |
| 2 | X+Y=3 | X+Y=4 | X+Y=5 |
| \(\frac{1}{8}\) | \(\frac{1}{8}\) | \(\frac{2}{8}\) |
表から、
X+Y=2の確率が2/8、
X+Y=3の確率が2/8、
X+Y=4の確率が2/8、
X+Y=5の確率が2/8
ですから期待値は、
E[X+Y]=2×2/8+3×2/8+4×2/8+5×2/8=3.5
表を追加すれば簡単ですね！
E[XY]の解法
同様にXYの場合について下表を追加しましょう。
| X/Y | 1 | 2 | 3 |
| 1 | XY=1 | XY=2 | XY=3 |
| – | \(\frac{2}{8}\) | \(\frac{1}{8}\) | \(\frac{1}{8}\) |
| 2 | XY=2 | XY=4 | XY=6 |
| – | \(\frac{1}{8}\) | \(\frac{1}{8}\) | \(\frac{2}{8}\) |
表から、
XY=1の確率が2/8、
XY=2の確率が2/8、
XY=3の確率が1/8、
XY=4の確率が1/8、
XY=6の確率が2/8
ですから期待値は、
E[XY]=1×2/8+2×2/8+3×1/8+4×1/8+6×2/8=25/8
表を追加すれば簡単ですね！
期待値をまとめると、
E[X]=3/2、E[Y]=2、E[X+Y]=3.5、E[XY]=25/8
となります。
また、
E[X+Y]= E[X]+ E[Y] は成り立ちますが、
E[XY]= E[X] E[Y] は成り立ちません。
X,Yは互いに独立ではないからですね。
解法(分散)
V[X]の解法
●ここで、分散V[X]の式をおさえましょう。
V[X]=E[\((X-μ_X)^2\)]
ですね。
次に、Xの平均\(μ_X\)を求めましょう。
平均\(μ_X\)はX=1,2の平均ですから3/2ですね。
表から、X=1の確率が1/2、X=2の確率が1/2ですから分散は、
V[X]=E[\((X-μ_X)^2\)]=E[\((X-1.5)^2\)]
=\((1-1.5)^2\)×1/2+\((2-1.5)^2\)×1/2
=1/4
ちょっと難しいですね。
V[Y]の解法
●ここで、分散V[Y]の式をおさえましょう。
V[Y]=E[\((Y-μ_Y)^2\)]
ですね。
次に、Yの平均\(μ_Y\)を求めましょう。
平均\(μ_Y\)はY=1,2,3の平均ですから2ですね。
表から、Y=1の確率が3/8、Y=2の確率が2/8、Y=3の確率が3/8ですから分散は、
V[Y]=E[\((Y-μ_Y)^2\)]=E[\((Y-2)^2\)]
=\((1-2)^2\)×3/8+\((2-2)^2\)×2/8+\((3-2)^2\)×3/8
=3/4
ちょっと難しいですね。
共分散COV[X,Y]の解法
●ここで、共分散COV[X,Y]の式をおさえましょう。
COV[X,Y]=E[\((X-μ_X)(Y-μ_Y)\)]
ですね。
共分散は、
COV[X,Y]=E[\((X-μ_X)(Y-μ_Y)\)]= E[\((X-1.5)(Y-2)\)]
=(1-1.5)(1-2)×2/8+(1-1.5)(2-2)×1/8+(1-1.5)(3-2)×1/8+
(2-1.5)(1-2)×1/8+(2-1.5)(2-2)×1/8+(2-1.5)(3-2)×2/8
=1/8
なお、共分散Cov[X,Y]はもう１つ公式があり、
Cov[X,Y]= E[XY]- E[X]E[Y]
1/8=25/8-3/2・3
が成り立ちます。
ちょっと難しいですが、解き方は１パターンなので、何度も復習しましょう。
分散をまとめると、
V[X]=1/4、V[Y]=3/4、Cov[X,Y]=1/8
となります。
離散系で使った公式一覧
➂【何度も復習しよう！】連続型確率分布の場合
その２の記事で解説します。
まとめ
同時確率分布の分散、共分散の導出をわかりやすく解説しました。
- ①【共通】２段サンプリングの分散公式を導出するために知っておくべき内容
- ➁【何度も復習しよう！】離散型確率分布の場合
- ➂【何度も復習しよう！】連続型確率分布の場合(その２で解説)
Warning: count(): Parameter must be an array or an object that implements Countable in /home/qcplanets/qcplanets.com/public_html/wp-content/themes/m_theme/sns.php on line 119