積率母関数を用いた連続一様分布の期待値(平均)と分散の導出

更新日

公式の確認

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

積率母関数

etb−etat(b−a)\frac{{\mathrm{e}^{tb}}-{\mathrm{e}^{ta}}}{t(b-a)}

期待値

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

分散

V(X)=112(b−a)2V(X)=\frac{1}{12}(b-a)^2

連続一様分布の積率母関数の導出

mX(t)=E(etX)=∫−∞∞etxf(x)dx=∫−∞aetxf(x)dx+∫abetxf(x)dx+∫b∞etxf(x)dx=0+∫abetx1b−adx+0=1b−a[etxt]ab=etb−etat(b−a)\begin{equation*}\begin{split}m_X(t)&=E(\mathrm{e}^{tX})\\ &=\displaystyle \int_{ - \infty }^{ \infty }\mathrm{e}^{tx}f(x)dx\\ &=\displaystyle \int_{ - \infty }^{a} \mathrm{e}^{tx}f(x)dx+\displaystyle \int_{ a }^{ b } \mathrm{e}^{tx}f(x)dx+\displaystyle \int_{b}^{ \infty } \mathrm{e}^{tx}f(x)dx\\ &=0+\displaystyle \int_{ a }^{ b } \mathrm{e}^{tx}\frac{1}{b-a}dx+0\\ &=\frac{1}{b-a}\left[\frac{\mathrm{e}^tx}{t} \right]_{a}^{b}\\ &=\frac{{\mathrm{e}^{tb}}-{\mathrm{e}^{ta}}}{t(b-a)}\end{split}\end{equation*}

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

E(X)=dmX(t)dt∣t=0=(etb−eta)′t(b−a)−(etb−eta)(t(b−a))′(t(b−a))2∣t=0=(aetb−beta)t(b−a)−(etb−eta)(b−a)(t(b−a))2∣t=0\begin{equation*}\begin{split}E(X)&=\left.\frac{d{m_X}(t)}{dt}\right|_{t=0}\\ &=\left.\frac{{({\mathrm{e}^{tb}}-{\mathrm{e}^{ta}})}'t(b-a)-({\mathrm{e}^{tb}}-{\mathrm{e}^{ta}}){(t(b-a))}'}{{(t(b-a))}^2}\right|_{t=0}\\ &=\left.\frac{{({a\mathrm{e}^{tb}}-{b\mathrm{e}^{ta}})}t(b-a)-({\mathrm{e}^{tb}}-{\mathrm{e}^{ta}}){(b-a)}}{{(t(b-a))}^2}\right|_{t=0}\\ \end{split}\end{equation*}

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

ロピタルの定理

lim⁡x→af(x)g(x)=lim⁡x→af(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*}

これを前述の式に適応させると、

dmX(t)dt∣t=0=(aetb−beta)t(b−a)−(etb−eta)(b−a)(t(b−a))2∣t=0=lim⁡t→0(aetb−beta)t(b−a)−(etb−eta)(b−a)(t(b−a))2=lim⁡t→0f(t)g(t)(=lim⁡t→0f(t)′g(t)′=lim⁡t→0f(t)′′g(t)′′)\begin{equation*}\begin{split}\left.\frac{d{m_X}(t)}{dt}\right|_{t=0}&=\left.\frac{{({a\mathrm{e}^{tb}}-{b\mathrm{e}^{ta}})}t(b-a)-({\mathrm{e}^{tb}}-{\mathrm{e}^{ta}}){(b-a)}}{{(t(b-a))}^2}\right|_{t=0}\\ &=\lim_{t \to 0}\frac{{({a\mathrm{e}^{tb}}-{b\mathrm{e}^{ta}})}t(b-a)-({\mathrm{e}^{tb}}-{\mathrm{e}^{ta}}){(b-a)}}{{(t(b-a))}^2}\\ &=\lim_{t \to 0} \frac{f(t)}{g(t)}(=\lim_{t \to 0} \frac{f(t)'}{g(t)'}=\lim_{t \to 0} \frac{f(t)''}{g(t)''}) \end{split}\end{equation*}

これより、lim⁡t→0f(t)′′g(t)′′\begin{equation*}\begin{split}\lim_{t \to 0}\frac{f(t)''}{g(t)''}\end{split}\end{equation*}を求めることにより、平均を導出できる。以下でその過程について記述する。

f(t)=(betb−aeta)t(b−a)−(etb−eta)(b−a)f′(t)=(betb−aeta)(b−a)+(b2etb−a2eta)t(b−a)−(betb−aeta)(b−a)=(b2etb−a2eta)t(b−a)f′′(t)=(b2etb−a2eta)(b−a)−(b3etb+a3eta)t(b−a) g(t)=(t(b−a))2g′(t)=2t(b−a)2g′′(t)=2(b−a)2f(t)={({b\mathrm{e}^{tb}}-{a\mathrm{e}^{ta}})}t(b-a)-({\mathrm{e}^{tb}}-{\mathrm{e}^{ta}}){(b-a)}\\ \begin{equation*}\begin{split}f'(t)&={({b\mathrm{e}^{tb}}-{a\mathrm{e}^{ta}})}(b-a)+({b^{2}\mathrm{e}^{tb}}-{a^{2}\mathrm{e}^{ta}})t(b-a)-({b\mathrm{e}^{tb}}-{a\mathrm{e}^{ta}}){(b-a)}\\ &=({b^{2}\mathrm{e}^{tb}}-{a^{2}\mathrm{e}^{ta}})t(b-a)\\\end{split}\end{equation*}\\ f''(t)=({b^{2}\mathrm{e}^{tb}}-{a^{2}\mathrm{e}^{ta}})(b-a)-({b^{3}\mathrm{e}^{tb}}+{a^{3}\mathrm{e}^{ta}})t(b-a)\\\ g(t)={(t(b-a))}^2\\ g'(t)=2t(b-a)^2\\ g''(t)=2(b-a)^2\\\\\\
lim⁡t→0f(t)′′g(t)′′=lim⁡t→0(b2etb−a2eta)(b−a)−(b3etb+a3eta)t(b−a)2(b−a)2=b2−a22(b−a)=a+b2(=dmX(t)dt∣t=0=E(X))\begin{equation*}\begin{split}\lim_{t \to 0}\frac{f(t)''}{g(t)''}&=\lim_{t \to 0} \frac{({b^{2}\mathrm{e}^{tb}}-{a^{2}\mathrm{e}^{ta}})(b-a)-({b^{3}\mathrm{e}^{tb}}+{a^{3}\mathrm{e}^{ta}})t(b-a)}{2(b-a)^2}\\ &=\frac{b^{2}-a^{2}}{2(b-a)}\\ &=\frac{a+b}{2}(=\left.\frac{d{m_X}(t)}{dt}\right|_{t=0}=E(X))\end{split}\end{equation*}

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

E(X2)=d2mX(t)dt2∣t=0=((aetb−beta)t(b−a)−(etb−eta)(b−a)(t(b−a))2)′∣t=0\begin{equation*}\begin{split}E(X^2)&=\left.\frac{d^2{m_X}(t)}{dt^2}\right|_{t=0}\\ &=\left.(\frac{{({a\mathrm{e}^{tb}}-{b\mathrm{e}^{ta}})}t(b-a)-({\mathrm{e}^{tb}}-{\mathrm{e}^{ta}}){(b-a)}}{{(t(b-a))}^2})'\right|_{t=0} \end{split}\end{equation*}

※これを前述のロピタルの定理を用いて導くことができるが、煩雑であるため省略する。連続一様分布の分散は、通常、確率密度関数から求める手法が採用される。

その値は、

V(X)=112(b−a)2\begin{equation*}\begin{split}V(X)=\frac{1}{12}{(b-a)}^2\end{split}\end{equation*}

となる。

関連記事

確率密度関数を用いた連続一様分布の期待値(平均)と分散の導出

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

カテゴリ: 連続一様分布

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

関連サービス

講座一覧ページ

記事一覧はこちら

無料で統計学を学ぶ

AVILEN 採用情報

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