本稿では、確率漸化式を用いることで二項分布の最頻値を導出しています。離散型確率分布は、確率関数を微分するという概念がないため、確率漸化式を用いて、確率関数を最大化する確率変数の値を求めます。
なお、閲覧にあたっては、以下の点にご注意ください。
- スマートフォンやタブレット端末でご覧の際、数式が見切れている場合は、横にスクロールすることができます。
【定理】二項分布の確率漸化式と最頻値
【定理】
二項分布の確率漸化式と最頻値
Recurrence Relation and Mode of Binomial Distribution
確率変数 $X$ が二項分布
\begin{align}
X \sim B \left(n,p\right)
\end{align}
に従うとき、
$P \left(X\right)$ と $P \left(X+1\right)$ の比について、
\begin{align}
\frac{P \left(X+1\right)}{P \left(X\right)}=\frac{ \left(n-x\right)p}{ \left(x+1\right) \left(1-p\right)}
\end{align}
が成り立ち、
最頻値は、
(i)$ \left(n+1\right)p$ が整数のとき
\begin{align}
X= \left(n+1\right)p-1,X= \left(n+1\right)p
\end{align}
(ii)$ \left(n+1\right)p$ が整数でないとき
\begin{align}
X= \left\lceil \left(n+1\right)p\right\rceil
\end{align}
ただし、$ \left\lceil x\right\rceil$ は、「$x$ を超えない最大の整数」を表す。
で与えられる。
証明
確率関数の定義式を用いて、比を計算すると、
\begin{align}
\frac{P \left(X+1\right)}{P \left(X\right)}&=\frac{{}_{n}C_{x+1}p^{x+1} \left(1-p\right)^{n-x-1}}{{}_{n}C_xp^x \left(1-p\right)^{n-x}}\\
&=\frac{\displaystyle\frac{n!}{ \left(x+1\right)! \left(n-x-1\right)!}p^{x+1} \left(1-p\right)^{n-x-1}}{\displaystyle\frac{n!}{x! \left(n-x\right)!}p^x \left(1-p\right)^{n-x}}\\
&=\frac{ \left(n-x\right)p}{ \left(x+1\right) \left(1-p\right)}
\end{align}
この漸化式は、2つの比が
1よりも小さいときは、$P \left(X\right)$ の方が大きい
1のときは、$P \left(X\right)$ と $P \left(X+1\right)$ が等しい
1よりも大きいときは、$P \left(X+1\right)$ の方が大きい
すなわち、
\begin{gather}
\frac{P \left(X+1\right)}{P \left(X\right)} \lt 1\Leftrightarrow P \left(X+1\right) \lt P \left(X\right)\\
1=\frac{P \left(X+1\right)}{P \left(X\right)}\Leftrightarrow P \left(X\right)=P \left(X+1\right)\\
1 \lt \frac{P \left(X+1\right)}{P \left(X\right)}\Leftrightarrow P \left(X\right) \lt P \left(X+1\right)
\end{gather}
ということを意味している。
比が1のときの $X$ の値を求めると、
\begin{gather}
\frac{ \left(n-x\right)p}{ \left(x+1\right) \left(1-p\right)}=1\\
\left(x+1\right) \left(1-p\right)= \left(n-x\right)p\\
x-xp+1-p=np-xp\\
x=np+p-1\\
x= \left(n+1\right)p-1
\end{gather}
したがって、$P \left(X\right)$ は、
$x \lt \left(n+1\right)p-1$ で単調増加
$ \left(n+1\right)p-1 \lt x$ で単調減少
すなわち
\begin{gather}
x \lt \left(n+1\right)p-1\Leftrightarrow P \left(X\right) \lt P \left(X+1\right)\\
\left(n+1\right)p-1 \lt x\Leftrightarrow P \left(X+1\right) \lt P \left(X\right)
\end{gather}
最頻値を $m$ とすると、
\begin{align}
P \left(0\right) \lt P \left(1\right) \lt \cdots \lt P \left(m-1\right) \lt P \left(m\right) \geq P \left(m+1\right) \gt \cdots \gt P \left(n\right)
\end{align}
となる。
したがって、最頻値の定義より、
(i)$ \left(n+1\right)p$ が整数のとき
\begin{align}
X= \left(n+1\right)p-1 \quad X+1= \left(n+1\right)p
\end{align}
において、確率が最大になる。
(ii)$ \left(n+1\right)p$ が整数でないとき、
\begin{align}
X+1= \left\lceil \left(n+1\right)p\right\rceil
\end{align}
すなわち、
$ \left(n+1\right)p$ を超えない最大の整数
のとき、確率が最大になる。
$\blacksquare$
参考文献
- 小寺 平治 著. 数理統計:明解演習. 共立出版, 1986, p.62
- 野田 一雄, 宮岡 悦良 著. 入門・演習数理統計. 共立出版, 1990, p.06-107 練習問題 ex.3.1.11
0 件のコメント:
コメントを投稿