本記事は「これなら分かる!はじめての数理統計学」シリーズに含まれます。np(1-p)\end{align}ただし,$p=M/N$であり,$\tilde{C}$は二項係数を拡張した記号である。(\text{その他})\end{cases}\end{align}アタリが$M$個,ハズレが$N-M$個入っているくじ引きから$n$個を引くとき,アタリの個数$X$は超幾何分布に従います。超幾何分布は再生性を持たず,ロードマップ上では二項分布を非復元抽出として拡張した場合に相当します。確率質量関数は,高校数学の場合の数を利用して導かれます。アタリが$M$個,ハズレが$N-M$個入っているくじ引きから$n$個を引くとき,アタリの個数$X$は超幾何分布に従います。分母には,$N$個の中から$n$個選ぶ全事象が入ります。分子には,$M$個のアタリから$x$個選び,$N-M$個のはずれから$n-x$個選ぶ場合の数が入ります。n$のとき,超幾何分布の確率質量変数は以下のように表されます。\tilde{\C}_{x}$を導入しました。平均と分散を求める方法は2通りあります。以下で説明しますが,超幾何分布の平均と分散を定義に沿って愚直に計算すると少し厄介ですので,ぜひこちらの「有限母集団から抽出した標本平均の期待値と分散」を利用する方法をおさえておいて下さい。有限母集団における標本平均の期待値と分散で説明した通り,有限母集団から抽出した標本平均の期待値と分散は以下のように表されます。\frac{N-n}{N-1}\frac{\sigma^2}{n}\end{align}ただし,有限母集団の母平均を$\mu$,母分散を$\sigma^2$とおきました。この結果を利用しましょう。0\end{alignat}とおきます。特性値が$1$であることをアタリとしました。\overline{X}\end{align}が従う分布が超幾何分布になります。アタリが$M$個,ハズレが$N-M$個入っているくじ引きから$n$個を引くとき,アタリの個数は超幾何分布に従うからです。ここからは$E[X]$と$V[X]$を求めることが目標です。現段階では$\mu$と$\sigma^2$が分かっていないため,$E[X]$と$V[X]$を求めることができません。そこで,改めて$\mu$と$\sigma^2$の定義を確認しておきましょう。有限母集団の母平均と母分散の定義より,$\mu$と$\sigma^2$は以下のように定義されます。V[X_1]\end{align}ただし,有限母集団から$1$個の観測値を無作為抽出した確率変数を$X_1$とおきました。そこで,ここからは$E[X_1]$と$V[X_1]$を求めることが目標になりました。アタリが$M$個,ハズレが$N-M$個入っている有限母集団から無作為に$1$個の観測値$X_1$を抽出するとき,$X_1$が$1$である確率は$M/N$ですので,$X_1$はパラメータ$p=M/N$のベルヌーイ分布に従っていることになります。したがって,$E[X_1]$と$V[X_1]$は以下のように表されます。p(1-p)\end{align}さて,ここで有限母集団から抽出した標本平均の期待値と分散の定理を利用します。有限母集団から抽出された観測値の標本平均の期待値は母平均と変わらず,分散は無限母集団の場合の結果に対して有限修正項を付け加えた結果となります。以下のように表されます。\frac{N-n}{N-1}np(1-p)\end{align}「無作為に$1$個の観測値$X_1$を抽出するとき」で少し違和感を覚える方がいるかもしれません。自分も例外ではなく,「なぜ$1$回の抽出に注目するだけでよいのか」という部分で引っかかりました。実際に非復元抽出しているのは$n$個のはずなのに,最初の$1$回だけに着目してベルヌーイ試行とみなすのは「非復元抽出とは言えないのではないか」と思っていました。この疑問に対する答えは「非復元抽出という意味合いは有限修正項によって付け加えられる。そもそも,母平均$\mu$の定義が$E[X_1]$であり,有限母集団から無作為に抽出された$1$回の観測値に関する期待値を表している。母分散も同様に,$1$回の観測値に関する分散を表している。有限修正項は,$1$回の無作為抽出に対して,有限母集団からの非復元抽出という意味合いを付け加えるために利用される。$E[X_1]$の定義自体が,$1$回の無作為抽出に関する結果を表していたのです。$2$回以上の試行を考慮したいのであれば,確率変数を条件付ける必要があります。,$E[X_2|X_1]$は$1$回目の試行の結果を踏まえたうえでの$2$回目の試行の期待値を表しています。通常,離散分布の平均と分散を求める際には「確率母関数の性質」を利用します。超幾何分布の確率母関数は複雑であるため,平均と分散を求める際は定義を利用します。新しく導入した二項係数に関して重要な性質を利用します。{}_{N-M}\tilde{\C}_{n-x}\label{equation:拡張二項係数の性質}\end{align}この等式の意味は,左辺は$N$個のくじから$n$本を引くという場合の数で,右辺はそれをアタリの本数$x$による場合わけでカウントしています。全てのアタリの本数$M$から$x$本引く場合と,ハズレから残りを引く場合を考えているだけです。平均から計算していきます。ポイントは,分子の掛け算のうち右側の項は変わらないということです。np\end{alignat}ただし,$p=M/N$とおきました。続いて,分散を計算していきます。基本的には平均と同じ方針で計算することができるのですが,二項係数に「階乗」の形が複数含まれていることから,計算を簡単にするために$E[X^2]$ではなく$E[X(X-1)]$を考えることにします。確率母関数から分散を求めるときと同じアイディアです。E[X]-E[X]^2\end{align}$E[X]$はすでに求められていますので,以下では$E[X(X-1)]$を計算していきます。途中まで先ほどと同様の計算が現れます。ポイントも先ほどと同じで,分子の掛け算のうち右側の項は変わらないということです。\frac{n(n-1)M(M-1)}{N(N-1)}\end{align}したがって,分散は以下のように計算できます。再生性を示すためには,再生性を示したい分布に従う独立な二つの確率変数を考え,その和の確率母関数を計算したときに,パラメータが和の形になっていることを示します。超幾何分布の確率母関数は複雑な形をしており,和の確率母関数を考えてもパラメータが和の形にはならないため,超幾何分布は再生性を持ちません。ロードマップ確率分布のロードマップさて,ロードマップに戻りましょう。超幾何分布は,二項分布の非復元拡張として導入されました。以下の内容も参考になるでしょう。執筆にあたり参考にした文献は,以下でリストアップしております。※コメントに返信があったらメールで通知する。私はロボットではありません。