二変量ポアソン分布の導出


13

最近、2変量ポアソン分布に遭遇しましたが、その導出方法について少し混乱しています。

分布は次のとおりです。

P(X=x,Y=y)=e(θ1+θ2+θ0)θ1xx!θ2yy!i=0min(x,y)(xi)(yi)i!(θ0θ1θ2)i

私が収集できることから、θ0項はXYの間の相関の尺度Yです。したがって、XYが独立している場合、θ0=0あり、分布は2つの単変量ポアソン分布の積になります。

これを念頭に置いて、私の混乱は総和項に基づいています-この項はXYの間の相関を説明すると仮定していYます。

私には、被加数は「成功」の確率が\ left(\ frac {\ theta_ {0}} {\ theta_ {1} \ theta_ {2}で与えられる二項累積分布関数のある種の積を構成するように思われます} \ right)(θ0θ1θ2)および「失敗」の確率はi!^ {\ frac {1} {min(x、y)-i}}によって与えられます。i!1min(x,y)iなぜなら、(i!1min(x,y)i!)(min(x,y)i)=i!、しかしこれで大丈夫です。

誰かがこの分布をどのように導き出すことができるかについての支援を提供できますか?また、このモデルを多変量シナリオ(3つ以上のランダム変数など)に拡張する方法を回答に含めることができれば、それは素晴らしいことです!

(最後に、以前に投稿された同様の質問(2変量ポアソン分布を理解する)があったことに注意しましたが、その導出は実際には調査されませんでした。)


2
指数と最初の項はあってはならないの代わりに?e(θ1+θ2+θ0)eθ1+θ2+θ0
ジル

1
@Giles申し訳ありませんが、最初はあなたのコメントを読み間違えました-はい、あなたは正しいです。用語はます。それをキャッチしてくれてありがとう!e(θ1+θ2+θ0)
user9171 14

3
一般に、いくつかの従来の例外を除いて、単変量分布の多変量バージョンの「the」ではありません(「the」多変量正規分布など)。多機能拡張機能を取得する方法は、どの機能が最も重要かによって異なります。異なる著者は、一般的な単変量分布の異なる多変量バージョンを持っている場合があります。だから、一般的に、人は「何か言うかもしれない「ポアソンを多変量ポアソンを」、または「まあまあの二変量この1つはかなり自然なだけでなく、1です。。
Glen_b -Reinstateモニカ

2
(ctd)...例えば、一部の著者は、負の依存性を持つ多変量分布を探します。これは、この分布にはない機能です。
Glen_b -Reinstateモニカ

回答:


17

スライドプレゼンテーション、KarlisとNtzoufrasは分布として二変量ポアソンを定義する独立ポアソン有する分布。このような分布を持つことは、(X,Y)=(X1+X0,X2+X0)Xiθi

Pr(Xi=k)=eθiθikk!

以下のためk=0,1,2,.

イベントは、イベントの互いに素な結合です。(X,Y)=(x,y)

(X0,X1,X2)=(i,xi,yi)

3つのコンポーネントすべてを非負の整数にするすべてのについて、そこから推定できます。は独立しているため、確率は増加します。0 I X Y X Ii0imin(x,y)Xi

F(θ0,θ1,θ2)(x,y)=Pr((X,Y)=(x,y))=i=0min(x,y)Pr(X0=i)Pr(X1=xi)Pr(X2=yi).

これは式です。完了です。 しかし、それが問題の式と同等であることを確認するには、ポアソン分布の定義を使用して、これらの確率をパラメーターに関して記述し、(いずれもゼロではないと仮定して)代数的に再処理します可能な限り製品ように見えるようにするには:θ 1θ 2のPr X 1 = X のPr X 2 = Y θiθ1,θ2Pr(X1=x)Pr(X2=y)

F(θ0,θ1,θ2)(x,y)=i=0min(x,y)(eθ0θ0ii!)(eθ1θ1xi(xi)!)(eθ2θ2yi(yi)!)=e(θ1+θ2)θ1xx!θ2yy!(eθ0i=0min(x,y)θ0ii!x!θ1i(xi)!y!θ2i(yi)!).

あなたが本当にしたい場合-それはやや示唆的です-二項係数およびを使用して合計の項を再表現できます、降伏(xi)=x!/((xi)!i!)(yi)

F(θ0,θ1,θ2)(x,y)=e(θ0+θ1+θ2)θ1xx!θ2yy!i=0min(x,y)i!(xi)(yi)(θ0θ1θ2)i,

質問のとおりです。


多変量シナリオへの一般化は、必要な柔軟性に応じて、いくつかの方法で進めることができます。最も単純なのは、

(X1+X0,X2+X0,,Xd+X0)

独立ポアソン分布変量。柔軟性を高めるために、追加の変数を導入できます。例えば、ポアソン独立使用変数の多変量分布検討、X0,X1,,XdηiY1,,YdXi+(Yi+Yi+1++Yd)i=1,2,,d.


1
称賛!ところで、すべきではない第二最後のステップの前に大きな括弧内にも?eθ0eθ2
ジル

1
@Gilles誤字をキャッチしていただきありがとうございます-私はそれを修正しました。の初期指数は必要があり。括弧内のは正しいです。θ0+θ1θ1+θ2eθ0
whuber

@whuberどうもありがとう!それは完璧な答えです!
user9171 14

@whuber素晴らしい答え!イベントがイベントの互いに素な結合である理由はまだわかりません。これはのみ当てはまると思い。おそらく(コンポーネント単位)を意味していましたか?しかし、分布関数を特徴付けるのにそれで十分ですか?(X,Y)=(x,y)(X0,X1,X2)=(i,xi,yi)i=0(X,Y)(x,y)
vanguard2k

@ vanguard2kあなたのコメントがわかりません。これらのイベントはばらばらではないと断言していますか?(異なる値を持っているため、そうでなければなりません。)または、それらが網羅的ではないと断言していますか?(その場合、X0(X,Y)
のどの値

4

以下に、2変量ポアソン分布を導出する方法を示します。

ましょうパラメータを持つ独立したポアソンランダム変数である。次に、を定義します。と両方に共通の変数は、ペアを相関させます。次に、確率質量関数を計算する必要があります。X0,X1,X2θ0,θ1,θ2Y1=X0+X1,Y2=X0+X2X0Y1Y2(Y1,Y2)

P(Y1=y1,Y2=y2)=P(X0+X1=y1,X0+X2=y2)=x0=0min(y1,y2)P(X0=x0)P(X1=y1x0)P(X2=y2y0)=x0=0min(y1,y2)eθ0θ0x0x0!eθ1θ1y1x0(y1x0)!eθ2θ2y2x0(y2x0)!=eθ0θ1θ2θ1y1θ2y2x0=0min(y1,y2)(θ0θ1θ2)x0x0!(y1x0)(y2x0)

これが役に立てばです!

1
こんにちはKjetil-書式設定に関する問題を修正しました(ただし、できるだけ変更しないように、いくつかの入力ミスはそのまま残しました)。特に、最終結果が不正確になる重要な要因を途中で失った場合、なぜ以前の回答で派生物のレプリカを投稿するのか理解できません。あなたがしようとしている特定のポイントはありますか?TEX
whuber

1
whuber:私はあなたの答えが投稿される前に私の答えを書き始めました!そうでなければ、私はそれを書いていなかっただろう。
kjetil bハルヴォルセン14
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.