確率質量関数を用いた二項分布の期待値・分散の導出

更新日

公式の確認

まずは、公式の確認をしましょう。

確率質量関数

P(X=k)=(nk)pk(1−p)n−kP(X=k)=\begin{pmatrix}n \\ k\end{pmatrix} p^{k}{(1-p)}^{n-k}

期待値

E(X)=npE(X)=np

分散

V(X)=np(1−p)V(X)=np(1-p)

期待値の導出

E(X)=∑k=0nkP(X=k)=∑k=0nk(nk)pk(1−p)n−k=∑k=0nk <em>nC</em>kpk(1−p)n−k=∑k=0nkn!(n−k)!k!pk(1−p)n−k=∑k=0nn(n−1)!(n−k)!(k−1)!ppk−1(1−p)n−k=np∑k=0n(n−1)!((n−1)−(k−1))!(k−1)!pk−1(1−p)n−k=np\begin{equation*}\begin{split}E(X)&=\sum_{k=0}^{n}kP(X=k)\\ &=\sum_{k=0}^{n}k\begin{pmatrix}n \\ k\end{pmatrix} p^{k}{(1-p)}^{n-k}\\ &=\sum_{k=0}^{n}k\ <em>n C</em>k p^{k}{(1-p)}^{n-k}\\ &=\sum_{k=0}^{n}k\frac{n!}{{(n-k)}!k!}p^{k}{(1-p)}^{n-k}\\ &=\sum_{k=0}^{n}\frac{n {(n-1)}!}{{(n-k)}!{(k-1)}!}p p^{k-1}{(1-p)}^{n-k}\\ &=np\sum_{k=0}^{n}\frac{{(n-1)}!}{{((n-1)-(k-1))}!{(k-1)}!}p^{k-1}{(1-p)}^{n-k}\\ &=np\end{split}\end{equation*}

【補足】

(n−1)!((n−1)−(k−1))!(k−1)!pk−1(1−p)n−k\begin{equation*}\begin{split}\frac{{(n-1)}!}{{((n-1)-(k-1))}!{(k-1)}!}p^{k-1}{(1-p)}^{n-k}\end{split}\end{equation*}

この部分のn-1をn'、k-1をk'とおくと、

n′!(n′−k′)!k′!pk′(1−p)n′−k′\begin{equation*}\begin{split}\frac{n'!}{{(n'-k')}!k'!}p^{k'}{(1-p)}^{n'-k'}\end{split}\end{equation*}

となる。これはパラメータがn'とk'の二項分布の確率質量関数である。よって、

∑k=0n(n−1)!((n−1)−(k−1))!(k−1)!pk−1(1−p)n−k\begin{equation*}\begin{split}\sum_{k=0}^{n}\frac{{(n-1)}!}{{((n-1)-(k-1))}!{(k-1)}!}p^{k-1}{(1-p)}^{n-k}\end{split}\end{equation*}

は確率質量関数を確率変数がとりうる値において全て足しあわせた値であるため、1である。

(ある事象における全ての確率を足すと1になることと同義)

分散の導出

E(X2)=∑k=0nk2P(X=k)=∑k=0nk2(nk)pk(1−p)n−k=∑k=0nk2 <em>nC</em>kpk(1−p)n−k=∑k=0n(k(k−1)+k)n!(n−k)!k!pk(1−p)n−k=∑k=0nk(k−1)n!(n−k)!k!pk(1−p)n−k+∑k=0nkn!(n−k)!k!pk(1−p)n−k\begin{equation*}\begin{split}E(X^2)&=\sum_{k=0}^{n}k^2P(X=k)\\ &=\sum_{k=0}^{n}k^2\begin{pmatrix}n \\ k\end{pmatrix} p^{k}{(1-p)}^{n-k}\\ &=\sum_{k=0}^{n}k^2\ <em>n C</em>k p^{k}{(1-p)}^{n-k}\\ &=\sum_{k=0}^{n}(k(k-1)+k)\frac{n!}{{(n-k)}!k!}p^{k}{(1-p)}^{n-k}\\ &=\sum_{k=0}^{n}k(k-1)\frac{n!}{{(n-k)}!k!}p^{k}{(1-p)}^{n-k}+\sum_{k=0}^{n}k\frac{n!}{{(n-k)}!k!}p^{k}{(1-p)}^{n-k}\end{split}\end{equation*}

【補足】

∑k=0nkn!(n−k)!k!pk(1−p)n−k\begin{equation*}\begin{split}\\\sum_{k=0}^{n}k\frac{n!}{{(n-k)}!k!}p^{k}{(1-p)}^{n-k}\end{split}\end{equation*}

これは、先ほど導出した二項分布の期待値であるため、npとなる。

           =∑k=0nn!(n−k)!(k−2)!pk(1−p)n−k+np=∑k=0nn(n−1)(n−2)!(n−k)!(k−2)!p2pk−2(1−p)n−k+np=n(n−1)p2∑k=0n(n−2)!((n−2)−(k−2))!(k−2)!pk−2(1−p)n−k+np=n(n−1)p2+np\begin{equation*}\begin{split}\ \ \ \ \ \ \ \ \ \ \ &=\sum_{k=0}^{n}\frac{n!}{{(n-k)}!(k-2)!}p^{k}{(1-p)}^{n-k}+np\\ &=\sum_{k=0}^{n}\frac{n(n-1)(n-2)!}{{(n-k)}!(k-2)!}p^2p^{k-2}{(1-p)}^{n-k}+np\\ &=n(n-1)p^2\sum_{k=0}^{n}\frac{(n-2)!}{{((n-2)-(k-2))}!(k-2)!}p^{k-2}{(1-p)}^{n-k}+np\\ &=n(n-1)p^2+np\end{split}\end{equation*}

【補足】

(n−2)!((n−2)−(k−2))!(k−2)!pk−2(1−p)n−k\begin{equation*}\begin{split}\frac{{(n-2)}!}{{((n-2)-(k-2))}!{(k-2)}!}p^{k-2}{(1-p)}^{n-k}\end{split}\end{equation*}

この部分のn-2をn'、k-2をk'とおくと、

n′!(n′−k′)!k′!pk′(1−p)n′−k′\begin{equation*}\begin{split}\frac{n'!}{{(n'-k')}!k'!}p^{k'}{(1-p)}^{n'-k'}\end{split}\end{equation*}

となる。これはパラメータがn'とk'の二項分布の確率質量関数である。よって、

∑k=0n(n−2)!((n−2)−(k−2))!(k−2)!pk−2(1−p)n−k\begin{equation*}\begin{split}\sum_{k=0}^{n}\frac{{(n-2)}!}{{((n-2)-(k-2))}!{(k-2)}!}p^{k-2}{(1-p)}^{n-k}\end{split}\end{equation*}

は確率質量関数を確率変数がとりうる値において全て足しあわせた値であるため、1である。

(ある事象における全ての確率を足すと1になることと同義)

V(X)=E(X2)−(E(X))2=n(n−1)p2+np−(np)2=np(1−p)\begin{equation*}\begin{split}V(X)&=E(X^2)-{(E(X))}^2\\ &=n(n-1)p^2+np-{(np)}^2\\ &=np(1-p) \end{split}\end{equation*}

関連記事

二項分布のわかりやすいまとめ

積率母関数を用いた二項分布の期待値・分散の導出

【AVILEN 公式note】データサイエンティストの仕事とは

カテゴリ: 二項分布

データサイエンティストを募集しています AVILEN採用サイトへ

関連サービス

講座一覧ページ

記事一覧はこちら

無料で統計学を学ぶ

AVILEN 採用情報

データサイエンティストを募集しています AVILEN採用サイトへ