Wie leitet man die Poisson-Verteilung aus der Gamma-Verteilung ab?

8

Sei T1,T2, eine Folge exponentieller Zufallsvariablen mit dem Parameter λ . Die Summe Sn=T1+T2++Tn ist eine Gammaverteilung. Soweit ich weiß, wird die Poisson-Verteilung durch Nt wie folgt definiert:

Nt=max{k:Skt}

Wie zeige ich formal, dass eine Poisson-Zufallsvariable ist?Nt

Anregungen geschätzt. Ich habe versucht, eine Reihe von Beweisen auszuarbeiten, komme aber nicht zur endgültigen Gleichung.

Verweise

http://en.wikipedia.org/wiki/Exponential_distribution

http://en.wikipedia.org/wiki/Gamma_distribution

http://en.wikipedia.org/wiki/Poisson_distribution

user862
quelle
3
@ user862, die Beweise, die ich aus der Hand kenne, sind nicht besonders direkt. Durrett hat eine Ableitung in seinem Wahrscheinlichkeitsbuch, die ziemlich sauber ist. Es dauert 3-4 Seiten, denke ich; Was, wenn Sie eines seiner Bücher gelesen haben, nach seinen Maßstäben ein langer Beweis ist. Resnick geht in seinem Text über stochastische Prozesse etwas abstrakter vor. Durch den Bau und den Einsatz größerer Hämmer kann er jedoch allgemeinere Ergebnisse erzielen. Ross hat zweifellos eine Behandlung in seinem Buch über stochastische Prozesse, aber ich bin damit nicht so vertraut.
Kardinal
Fand den Beweis in Durretts Buch. Es wird wirklich klar erklärt. Danke für die Hinweise.
user862

Antworten:

6

Ich bin sicher, dass Durretts Beweis nett ist. Eine einfache Lösung für die gestellte Frage lautet wie folgt.

Fürn1

P(Nt=n)=0tP(Sn+1>tSn=s)P(Snds)=0tP(Tn+1>ts)P(Snds)=0teλ(ts)λnsn1eλs(n1)!ds=eλtλn(n1)!0tsn1ds=eλt(λt)nn!

Für gilt .n=0P(Nt=0)=P(T1>t)=eλt

Dies beweist nicht, dass ein Poisson-Prozess ist, der schwieriger ist, aber es zeigt, dass die Randverteilung von Poisson mit dem Mittelwert .(Nt)t0Ntλt

NRH
quelle
3
(+1) Das Ansprechen auf bedingte Dichten ist hier nicht unbedingt erforderlich. Beachten Sie, dass und wir müssen nur die Verbindungsdichte über den Bereich . Da und unabhängig sind, ist dies ein einfaches Unterfangen. { ( s n , t n + 1 ) : 0 s nt , t n +P(Nt=n)=P(Snt,Sn+1>t)=P(Snt,Sn+Tn+1>t)S n T n + 1{(sn,tn+1):0snt,tn+1>tsn}R2SnTn+1
Kardinal
@cardinal - und wie ist die Antwort von @ NRH nicht einfach? In der Tat würde ich sagen, es ist einfacher, weil nur 1 Integration erforderlich ist.
Wahrscheinlichkeitslogik
@probabilityislogic: Meine "einfache" Referenz war nur eine Bemerkung über die verbleibende Berechnung, die in meinem Kommentar nicht gezeigt wurde. Es war in keiner relativen Hinsicht in Bezug auf die Antwort von @ NRH gemeint.
Kardinal
2
@probabilityislogic: In den ersten beiden Zeilen des @ NRH-Beweises verbirgt sich eine Menge (zusätzlicher) Theorie. Der Punkt meines Kommentars war, dass wir das gleiche Ergebnis erzielen können, wenn wir nur die Barebones gemeinsamer Wahrscheinlichkeitsverteilungen verwenden, nämlich nur das Produktmaß. Dies ist meines Erachtens eine grundlegend einfachere Grundlage für die Berechnung als die Einführung einer bedingten Erwartung und der Begründung, die erforderlich ist, um von Zeile eins zu zwei der Antwort von @ NRH zu gelangen. Ich meine nicht, dass ich als Kritik in irgendeiner Weise nur eine alternative Methode anbieten wollte.
Kardinal