二項分布の確率漸化式と最頻値の導出

公開日: 更新日:

【2023年3月3週】 【B000】数理統計学 【B030】離散型の確率分布

この記事をシェアする
  • B!
サムネイル画像

本稿では、確率漸化式を用いることで二項分布の最頻値を導出しています。離散型確率分布は、確率関数を微分するという概念がないため、確率漸化式を用いて、確率関数を最大化する確率変数の値を求めます。

なお、閲覧にあたっては、以下の点にご注意ください。

  • スマートフォンやタブレット端末でご覧の際、数式が見切れている場合は、横にスクロールすることができます。

【定理】二項分布の確率漸化式と最頻値

【定理】
二項分布の確率漸化式と最頻値
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

関連記事

自己紹介

自分の写真

yama

大学時代に読書の面白さに気づいて以来、読書や勉強を通じて、興味をもったことや新しいことを学ぶことが生きる原動力。そんな人間が、その時々に学んだことを備忘録兼人生の軌跡として記録しているブログです。

このブログを検索

ブログ アーカイブ

QooQ