2つの独立したランダム変数、正規およびカイ二乗の積のpdf


17

XとYが独立している場合、2つの独立したランダム変数XとYの積のpdfは何ですか?Xは正規分布、Yはカイ二乗分布です。

Z = XY

場合正規分布を有する およびは、自由度が カイ二乗分布 ここで、は単位ステップ関数です。X

XN(μx,σx2)
fX(x)=1σx2πe12(xμxσx)2
Yk
Yχk2
fY(y)=y(k/2)1ey/22k/2Γ(k2)u(y)
u(y)

とが独立している場合、のpdfはどうなりますか?ZXY

解決策を見つけるための一つの方法は、場合Rohatgiのよく知られた結果(1976、141頁)を使用することである連続的なRVののジョイントPDFである及びのPDF、である fXY(x,y)XYZ

fZ(z)=1|y|fXY(zy,y)dy

なぜなら、とは独立している 積分を解く問題に直面する場所。誰でもこの問題で私を助けることができます。XYfXY(x,y)=fX(x)fY(y)

fZ(z)=1|y|fX(zy)fY(y)dy
fZ(z)=1σx2π12k/2Γ(k2)01|y|e12(zyμxσx)2y(k/2)1ey/2dy
01|y|e12(zyμxσx)2y(k/2)1ey/2dy

これを解決する代替方法はありますか?


2
その最後のステップは完全に正しく見えません。「」はを意味しているようにますが、さらに重要なことには、下限を変更することはできません。積分をで2つに分けを負の範囲の1つに設定し、2つを組み合わせます。これにより、統合が扱いやすくなる可能性があります。一般化された超幾何関数の線形結合が得られるようです。fXYfX00yy
whuber

はい、それは間違いでしたはなければなりません。fZY(zy)fX(zy)
ロビン

しかし、は上の関数であり、単位ステップ関数示されるため、下限を0に変更することは有効だとます。fY(y)(0,)u(y)
ロビン

私はもはやこの種の計算の訓練を受けていません...しかし、閉じた式で終わる可能性はありません。実用的なアプリケーションにこれが必要な場合は、「これを効率的に計算する方法」に焦点を当てるべきだと思います。
エルビス

4
この質問の動機はありますか?割ったNormal はスチューデントのですが、なぜ乗算または割ったNormalを考慮するのでしょうか?χtχ2
西安

回答:


1

積分の項を単純化する

T=e12((zyμxσx)2y)yk/22

多項式検索ようにp(y)

[p(y)e12((zyμxσx)2y)]=p(y)e12((zyμxσx)2y)+p(y)[12((zyμxσx)2y)]e12((zyμxσx)2y)=T

次のようなを見つけることになりますp(y)

p(y)+p(y)[12((zyμxσx)2y)]=yk/22

または

p(y)12p(y)(zμxσx2y2z2σx2y31)=yk/22

すべてのべき乗を別々に評価して行うことができますy


コメント後に編集する

上記のソリューションは、発散するため機能しません。

しかし、他の何人かはこのタイプの製品に取り組んできました。

Fourrier変換の使用:

Schoenecker、Steven、Tod Luginbuhl。「2つのガウス確率変数の積とガウスおよびガンマ確率変数の積の特性関数。」IEEE信号処理レター23.5(2016):644-647。 http://ieeexplore.ieee.org/document/7425177/#full-text-section

およびの積に対して、特性関数を取得しました。Z=XYXN(0,1)YΓ(α,β)

φZ=1βα|t|αexp(14β2t2)Dα(1β|t|)

ウィテカーの機能(http://people.math.sfu.ca/~cbm/aands/page_686.htmDα

Mellin変換の使用:

スプリンガーとトムソンは、ベータ、ガンマ、およびガウスの分布確率変数の積の評価をより一般的に説明しています。

スプリンガー、MD、WEトンプソン。「ベータ、ガンマ、ガウス確率変数の積の分布。」SIAM Journal on Applied Mathematics 18.4(1970):721-737。 http://epubs.siam.org/doi/10.1137/0118065

Mellin積分変換を使用します。メリンは、変換のメリン変換の積である及び(参照http://epubs.siam.org/doi/10.1137/0118065又はhttps://projecteuclid.org/euclid.aoms/1177730201を)。研究された製品の場合、この製品の逆変換は、計算方法を提供および証明するメイジャーG関数として表現できます。ZXY

同じ手法を使用できる可能性はありますが、ガウス分布変数とガンマ分布変数の積は分析しませんでした。すぐにこれを行おうとすると、H関数(https://en.wikipedia.org/wiki/Fox_H-function)を取得することが可能になるはずですただし、G-機能するか、他の単純化を行います。

M{fY(x)|s}=2s1Γ(12k+s1)/Γ(12k)

そして

M{fX(x)|s}=1π2(s1)/2σs1Γ(s/2)

あなたが得る

M{fZ(x)|s}=1π232(s1)σs1Γ(s/2)Γ(12k+s1)/Γ(12k)

の分布は次のとおりです。Z

fZ(y)=12πicic+iysM{fZ(x)|s}ds

少なくともH関数として項を削除するために変数を変更した後)232(s1)

残っているのは、この逆メリン変換をG関数として表現するパズルです。と両方出現はこれを複雑にします。ガウス分布変数のみの積の別のケースでは、変数代入することにより、をに変換できます。しかし、カイ二乗分布の条件のために、これはもう機能しません。たぶんこれが、このケースの解決策を誰も提供していない理由です。ss/2s/2sx=w2


1
...これは...?
wolfies


この分析がどのような進歩を表すかは不明です。あなたは解決策を手に入れましたか?
whuber

多項式(解を閉じるの係数を見つけるのは退屈ですが、簡単ですが、私が開いたままにしました。いくつかの例をすぐに入力します。p(y)k
セクストゥスエンピリカス
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.