積率母関数を用いた離散一様分布の期待値・分散の導出

更新日

公式の確認

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

積率母関数

MX(t)=etN1etN1etM_X(t)=\frac{\mathrm{e}^t}{N}\frac{1-\mathrm{e}^{tN}}{1-\mathrm{e}^t}

期待値

E(X)=N+12E(X)=\frac{N+1}{2}

分散

V(X)=N2112V(X)=\frac{N^2-1}{12}

離散一様分布の積率母関数(モーメント母関数)の導出

MX(t)=E(etk)=k=1NetkP(X=k)=k=1Netk1N=1Nk=1Netk=1Net(1etN)1et=etN1etN1et\begin{equation*}\begin{split}M_X(t)&=E(\mathrm{e}^{tk})\\ &=\sum_{k=1}^{N}\mathrm{e}^{tk}P(X=k)\\ &=\sum_{k=1}^{N}\mathrm{e}^{tk}\frac{1}{N}\\ &=\frac{1}{N}\sum_{k=1}^{N}\mathrm{e}^{tk}\\ &=\frac{1}{N}\frac{\mathrm{e}^{t}(1-\mathrm{e}^{tN})}{1-\mathrm{e}^{t}}\\ &=\frac{\mathrm{e}^{t}}{N}\frac{1-\mathrm{e}^{tN}}{1-\mathrm{e}^{t}} \end{split}\end{equation*}

積率母関数を用いた期待値の導出

E(X)=dMx(t)dtt=0=(etN1etN1et)t=0=etN1etN1et+etN(NetN)(1et)(1etN)(et)(1et)2t=0=etN(1etN)(1et)+(NetN)(1et)(1etN)(et)(1et)2t=0\begin{equation*}\begin{split}E(X)&=\left.\frac{d{M_x}(t)}{dt}\right|_{t=0}\\ &=\left.(\frac{\mathrm{e}^{t}}{N}\frac{1-\mathrm{e}^{tN}}{1-\mathrm{e}^{t}})'\right|_{t=0}\\ &=\left.\frac{\mathrm{e}^{t}}{N}\frac{1-\mathrm{e}^{tN}}{1-\mathrm{e}^{t}}+\frac{\mathrm{e}^{t}}{N}\frac{(-N\mathrm{e}^{tN})(1-\mathrm{e}^{t})-(1-\mathrm{e}^{tN})(-\mathrm{e}^t)}{(1-\mathrm{e}^{t})^2}\right|_{t=0}\\ &=\left.\frac{\mathrm{e}^{t}}{N}\frac{(1-\mathrm{e}^{tN})(1-\mathrm{e}^{t})+(-N\mathrm{e}^{tN})(1-\mathrm{e}^{t})-(1-\mathrm{e}^{tN})(-\mathrm{e}^t)}{(1-\mathrm{e}^{t})^2}\right|_{t=0}\\ \end{split}\end{equation*}

【補足】

上式にt=0を代入すると、分母が0になる。そのため上式が定義できない。そこでロピタルの定理を用いる。

ロピタルの定理とは

limxaf(x)g(x)=limxaf(x)g(x)\begin{equation*}\begin{split}\lim_{x \to a} \frac{f(x)}{g(x)}=\lim_{x \to a} \frac{f(x)'}{g(x)'}\end{split}\end{equation*}

である。これを前述の式に適応させる。

=limt0etN(1etN)(1et)+(NetN)(1et)(1etN)(et)(1et)2=limt0(et((1etN)(1et)+(NetN)(1et)(1etN)(et)))(N(1et)2)=N2+N2N=N+12\begin{equation*}\begin{split}&=\lim_{t \to 0}\frac{\mathrm{e}^{t}}{N}\frac{(1-\mathrm{e}^{tN})(1-\mathrm{e}^{t})+(-N\mathrm{e}^{tN})(1-\mathrm{e}^{t})-(1-\mathrm{e}^{tN})(-\mathrm{e}^t)}{(1-\mathrm{e}^{t})^2}\\ &=\lim_{t \to 0}\frac{(\mathrm{e}^{t}((1-\mathrm{e}^{tN})(1-\mathrm{e}^{t})+(-N\mathrm{e}^{tN})(1-\mathrm{e}^{t})-(1-\mathrm{e}^{tN})(-\mathrm{e}^t)))''}{(N(1-\mathrm{e}^{t})^2)''}\\ &=\frac{N^2+N}{2N}\\ &=\frac{N+1}{2}\end{split}\end{equation*}

積率母関数を用いた分散の導出

E(X2)=d2mX(t)dt2t=0\begin{equation*}\begin{split}E(X^2)&=\left.\frac{d^2{m_X}(t)}{dt^2}\right|_{t=0}\end{split}\end{equation*}

【補足】

これを前述のロピタルの定理を用いて導くことができるが、計算が煩雑であるため省略する。離散一様分布の分散は、確率質量関数から求める手法が採用される。その値は、

V(X)=N2112\begin{equation*}\begin{split}V(X)=\frac{N^{2}-1}{12}\end{split}\end{equation*}となる。

関連記事

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

カテゴリ: 離散一様分布