※ブログ記事の商品・サービスリンクにはアフィリエイトリンクが含まれます。

ワイブル分布の期待値と分散【統計検定1級対策】

別記事で代表的な確率分布についてまとめていますが、文量が増えそうなのでワイブル分布の平均と分散の計算はここでまとめてみます。

そういえばワイブル分布については今までほとんど問題も解いたことなかったですし、馴染みがあまりありませんでした。

ワイブル分布の式は『現代数理統計学の基礎』に倣ってf(x)=abx^{b-1}exp(-ax^b)としておきます。

期待値の導出

定義通りに計算します。

\int_0^\infty xabx^{b-1}exp(-ax^b)dx\\=\int_0^\infty abx^bexp(-ax^b)dx

ここでax^b=tとして置換します。すると

\int_0^\infty(\frac{t}{a})^{\frac{1}{b}}exp(-t)dt\\=a^{-\frac{1}{b}}\Gamma(\frac{1}{b}+1)

となります。

分散の導出

これも二次のモーメントからゴリゴリ計算します。

\int_0^\infty x^2abx^{b-1}exp(-ax^b)dx\\=\int_0^\infty abx^bexp(-ax^b)dx\\=\int_0^\infty(\frac{t}{a})^{\frac{2}{b}})exp(-t)dt\\=a^{-\frac{2}{b}}\Gamma(\frac{2}{b}+1)

これをみるとワイブル分布のn次モーメントは

a^{-\frac{n}{b}}\Gamma(\frac{n}{b}+1)

という単純な形で推測できることが分かります。

よって分散は

V(X)=E[X^2]-(E[X])^2\\=a^{-\frac{2}{b}}\Gamma(\frac{2}{b}+1)-a^{-\frac{1}{b}}\Gamma^2(\frac{1}{b}+1)\}\\=a^{-\frac{2}{b}}\{\Gamma(\frac{2}{b}+1)-\Gamma^2(\frac{1}{b}+1)\}\}

となります。

コメントを残す

メールアドレスが公開されることはありません。 が付いている欄は必須項目です

CAPTCHA


日本語が含まれない投稿は無視されますのでご注意ください。(スパム対策)