Ich bin kürzlich auf die bivariate Poisson-Verteilung gestoßen, bin jedoch ein wenig verwirrt, wie sie abgeleitet werden kann.
Die Verteilung erfolgt durch:
P ( X = x , Y = y ) = e - ( θ 1 + θ 2 + θ 0 ) θ x 1x ! θ y 2y ! m i n ( x , y ) Σ i=0 ( xi ) ( yi ) ich! ( θ 0θ 1 θ 2 )i
Nach allem, was ich zusammenfassen kann, ist der Term θ 0
Vor diesem Hintergrund beruht meine Verwirrung auf dem Summationsterm - ich gehe davon aus, dass dieser Term die Korrelation zwischen X
Es scheint mir, dass der Summand eine Art Produkt binomialer kumulativer Verteilungsfunktionen darstellt, bei dem die Wahrscheinlichkeit des "Erfolges" durch (θ0θ1θ2)
Könnte jemand eine Hilfestellung geben, wie diese Verteilung abgeleitet werden kann? Auch wenn es in einer Antwort enthalten sein könnte, wie dieses Modell auf ein multivariates Szenario erweitert werden könnte (sagen wir drei oder mehr Zufallsvariablen), wäre das großartig!
(Schließlich habe ich bemerkt, dass es zuvor eine ähnliche Frage gab ( Verständnis der bivariaten Poisson-Verteilung ), aber die Ableitung wurde tatsächlich nicht untersucht.)
quelle
Antworten:
In einer Folienpräsentation definieren Karlis und Ntzoufras ein bivariates Poisson als die Verteilung von wobei die unabhängig Poisson Verteilungen haben. Denken Sie daran, dass eine solche Verteilung bedeutet(X,Y)=(X1+X0,X2+X0)(X,Y)=(X1+X0,X2+X0) XiXi θiθi
Pr(Xi=k)=e−θiθkik!
für k=0,1,2,….k=0,1,2,….
Das Ereignis (X,Y)=(x,y)(X,Y)=(x,y) ist die disjunkte Vereinigung der Ereignisse
(X0,X1,X2)=(i,x−i,y−i)
für alle , die alle drei Komponenten zu nicht negativen ganzen Zahlen machen, aus denen wir ableiten können . Da die unabhängig sind, multiplizieren sich ihre Wahrscheinlichkeitenii 0≤i≤min(x,y)0≤i≤min(x,y) XiXi
F(θ0,θ1,θ2)(x,y)=Pr((X,Y)=(x,y))=min(x,y)∑i=0Pr(X0=i)Pr(X1=x−i)Pr(X2=y−i).
Dies ist eine Formel; wir sind fertig. Um jedoch zu sehen, dass es der Formel in der Frage entspricht, verwenden Sie die Definition der Poisson-Verteilung, um diese Wahrscheinlichkeiten in Bezug auf die Parameter ; zu schreiben und (vorausgesetzt, dass keines von & thgr; 1, & thgr ; 2 Null ist) algebraisch zu überarbeiten so viel wie möglich aussehen wie das Produkt \ Pr (X_1 = x) \ Pr (X_2 = y) :θiθi θ1,θ2θ1,θ2 Pr(X1=x)Pr(X2=y)Pr(X1=x)Pr(X2=y)
F(θ0,θ1,θ2)(x,y)=min(x,y)∑i=0(e−θ0θi0i!)(e−θ1θx−i1(x−i)!)(e−θ2θy−i2(y−i)!)=e−(θ1+θ2)θx1x!θy2y!(e−θ0min(x,y)∑i=0θi0i!x!θ−i1(x−i)!y!θ−i2(y−i)!).
Wenn Sie es wirklich wollen, können Sie die Terme in der Summe mit den Binomialkoeffizienten (xi)=x!/((x−i)!i!)(xi)=x!/((x−i)!i!) Und \ binom erneut ausdrücken {y} {i}(yi)(yi) , nachgebend
F(θ0,θ1,θ2)(x,y)=e−(θ0+θ1+θ2)θx1x!θy2y!min(x,y)∑i=0i!(xi)(yi)(θ0θ1θ2)i,
genau wie in der frage.
Die Verallgemeinerung zu multivariaten Szenarien kann je nach erforderlicher Flexibilität auf verschiedene Arten erfolgen. Das einfachste würde die Verteilung von erwägen
(X1+X0,X2+X0,…,Xd+X0)
für unabhängige verteilte Poisson-Variablen X_0 . Für mehr Flexibilität könnten zusätzliche Variablen eingeführt werden. Zum Beispiel verwenden unabhängige Poisson Variablen und die multivariate Verteilung des berücksichtigen ,X0,X1,…,XdX0,X1,…,Xd ηiηi Y1,…,YdY1,…,Yd Xi+(Yi+Yi+1+⋯+Yd)Xi+(Yi+Yi+1+⋯+Yd) i=1,2,…,d.i=1,2,…,d.
quelle
Hier ist eine Möglichkeit, die bivariate Poissonverteilung abzuleiten.
Sei unabhängige Poisson-Zufallsvariablen mit den Parametern . Dann definieren wir . Die Variable , beide gemeinsam ein , verursacht das Paar korreliert. Dann müssen wir die Wahrscheinlichkeitsmassenfunktion berechnen:X0,X1,X2X0,X1,X2 θ0,θ1,θ2θ0,θ1,θ2 Y1=X0+X1,Y2=X0+X2X0Y1Y2(Y1,Y2)
P(Y1=y1,Y2=y2)=P(X0+X1=y1,X0+X2=y2)=min(y1,y2)∑x0=0P(X0=x0)P(X1=y1−x0)P(X2=y2−y0)=min(y1,y2)∑x0=0e−θ0θ0x0x0!e−θ1θ1y1−x0(y1−x0)!e−θ2θ2y2−x0(y2−x0)!=e−θ0−θ1−θ2θ1y1θ2y2min(y1,y2)∑x0=0(θ0θ1θ2)x0x0!(y1x0)(y2x0)
hoffe, das hilft!
quelle