2
0

Delete article

Deleted articles cannot be recovered.

Draft of this article would be also deleted.

Are you sure you want to delete this article?

More than 3 years have passed since last update.

3.4(難) 超幾何分布と二項分布の関係(方針でスターリングの公式を解説)

Last updated at Posted at 2021-07-24

方針

超幾何分布の密度関数を階乗表現を使って書き下すと,

\begin{align}
f_X(x)&= \frac{{}_MC_x\times{}_{N-M}C_{n-x}}{{}_{N}C_n}\\
&= \frac{M!}{x!(M-x)!}\times\frac{(N-M)!}{(n-x)!(N-M-n+x)!}\times\frac{n!(N-n)!}{N!}
\end{align}

となる.これを頑張って変形して,

f_X(x)={}_nC_xp^x(1-p)^{n-x},\ \ \ (\frac{M}{N}\rightarrow p)

に持っていく,という問題.実数の階乗は,実数が十分に大ならばスターリングの公式を用いた近似を行うことができる.

スターリングの公式


実数kが十分大きい時,

k!\approx\sqrt{2\pi}k^{k+\frac{1}{2}}\exp(-k)

答案

\begin{align}
f_X(x)&= \frac{{}_MC_x\times{}_{N-M}C_{n-x}}{{}_{N}C_n}\\
&= \frac{M!}{x!(M-x)!}\times\frac{(N-M)!}{(n-x)!(N-M-n+x)!}\times\frac{n!(N-n)!}{N!}\\
&= \frac{n!}{x!(n-x)!}\frac{M!(N-M)!(N-n)!}{(M-x)!(N-M-n+x)!N!}\\
&\approx {}_nC_x\frac{(Np)!(N-Np)!(N-n)!}{(Np-x)!(N-Np-n+x)!N!}\tag{1}\\
&\ \ \ \ \ \ (\because M\approx Np) 
\end{align}

今,N→∞より,(1)式の分数部分にスターリングの公式を適用すると,指数部分などが消去でき,

今,N→∞より,(1)式の分数部分にスターリングの公式を適用すると,

\frac{(Np)^{Np+\frac{1}{2}}(N-Np)^{N-Np+\frac{1}{2}}(N-n)^{N-n+\frac{1}{2}}}{(Np-x)^{Np-x+\frac{1}{2}}(N-Np-n+x)^{N-Np-n+x+\frac{1}{2}}N^{N+\frac{1}{2}}}\tag{2}

べき乗の数が分母分子で等しいので,分母分子をNで割ることは各要素の中身をNで割ることに等しい.よって,

\begin{align}
(2)&= \frac{(p)^{Np+\frac{1}{2}}(1-p)^{N-Np+\frac{1}{2}}(1-\frac{n}{N})^{N-n+\frac{1}{2}}}{(p-\frac{x}{N})^{Np-x+\frac{1}{2}}(1-p-\frac{n}{N}+\frac{x}{N})^{N-Np-n+x+\frac{1}{2}}1^{N+\frac{1}{2}}}\\
&\approx p^{Np+\frac{1}{2}-Np+x-\frac{1}{2}}(1-p)^{N-Np+\frac{1}{2}-N+Np+n-x-\frac{1}{2}}\ \ \ \ (\because N\rightarrow\infty)\\
&= p^x(1-p)^{n-x}\\
\therefore (1)&\approx {}_nC_xp^x(1-p)^{n-x}
\end{align}

参考文献

2
0
0

Register as a new user and use Qiita more conveniently

  1. You get articles that match your needs
  2. You can efficiently read back useful information
  3. You can use dark theme
What you can do with signing up
2
0

Delete article

Deleted articles cannot be recovered.

Draft of this article would be also deleted.

Are you sure you want to delete this article?