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

統計

確率変数の和と比（和の分布・畳み込み・比の分布の変数変換）【統計検定1級対策】

最近は専門医試験のレポート期限が迫ってきたことや日本語で書かないといけない論文もあって、ちょっと更新ペースが落ちますが、ぼちぼちやっていきます。

確率変数の和と比について変数変換するとどのようになるか、といった問題も時々出ています。比についてはあまりまとめられているのを見たことがなかったので、まとめてみました。

目次：

確率変数の和・畳み込みを使う

Z=X+Yという確率変数がある時に、どうやってZの確率密度関数を求めるかを考えます。（　X、Yは互いに独立とします）

通常とられる方法として「畳み込み」という方法があります。

何をやるかと言えば

①ヤコビアンが1になるような2変数に変換する

②片方の変数を全範囲で積分して同時確率密度関数→目的の確率密度関数を出す

ということをやります。

まず①ですが、Z=X+Y, W=Yとすると

X=Z-W

Y=W

なので、ヤコビアンは1になります。

そうするとz,wの同時確率密度関数は

f_{Z,W}(z,w)=f_X(z-w)f_Y(w)

となります。

よって

\(f_Z(z)=\int f_X(z-w)f_Y(w)dw\)

が目的の確率密度関数といえます。

話としてはわかりやすいのですが、最後の積分で計算式が煩雑となるのが問題点です。

確率変数の和・モーメント母関数を使う

もう一つの方法がモーメント母関数を使う方法です。これは求めたい確率密度関数やX、Yの確率密度関数のモーメント母関数がわかっているときに有効な方法です。

M_Z(t)=E[e^{tZ}]\\=E[e^{t(X+Y)}]\\=E[e^{tX}]E[e^{tY}]

として求めることができます。

例として標準正規分布に従うX,Yの和を考えてみましょう。

E[e^{tX}]=E[e^{tY}]=e^{\frac{t^2}{2}}

なので

M_Z(t)=E[e^{tX}]E[e^{tY}]\\=e^{t^2}

となります。

よって、Z\sim N(0,2)であることがわかります。

正規分布における和の再生性と呼ばれる特性ですね。

確率変数の比

続いて、比について考えてみます。

Z=\frac{X}{Y}とした場合、Zの確率密度関数はどうなるでしょうか。

この手の問題で困ったら、分布関数に立ち戻って考えるとうまくいくことが多い気がします。

P(Z\lt z)=P(\frac{X}{Y}\lt z)

となります。ここでYの値が正か負かで場合分けをして

X\lt zY(Y\gt 0), X\gt zY(Y\lt 0)

そうすると、それぞれの値の範囲で同時確率密度関数f(x,y)を積分することでZの確率分布関数を求めることができます。

F_Z(z)=\int_{-\infty}^0dy\int_{zy}^\infty dxf(x,y)+\int_0^\infty dy\int_{-\infty}^{zy}dxf(x,y)

あとはこれをzで微分して

f_Z(z)=-\int_{-\infty}^0dyf(zy,y)y+\int_{0}^\infty dyf(zy,y)y\\=\int_{-\infty}^{\infty}|y|f(zy,y)dy

となります。

具体的な例として標準正規分布の比がコーシー分布となる、という問題があります。

現代数理統計学の基礎 4章 問12 – 脳内ライブラリアン

確率分布同士のいろんな関係性がわかると面白いですね。

参考文献：

Ratio distribution – Wikipedia

比の分布が説明された英語版Wikipediaです

比の分布関数

比の分布に関して書いてあるブログです、参考になりました

『現代数理統計学の基礎』

コメントを残す

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

CAPTCHA


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