Mittelwert und Varianz einer Null-Inflations-Poisson-Verteilung

11

Kann jemand zeigen, wie der erwartete Wert und die Varianz des Null-aufgeblasenen Poisson mit Wahrscheinlichkeitsmassenfunktion funktionieren

f((y)={π+((1- -π)e- -λ,wenn y=0((1- -π)λye- -λy!,wenn y=1,2 ....

wobei die Wahrscheinlichkeit ist, dass die Beobachtung durch einen Binomialprozess Null ist und der Mittelwert des Poisson ist, abgeleitet wird?πλ

Das Ergebnis ist der erwartete Wert und die Varianz ist .μ=((1- -π)λμ+π1- -πμ2

ADD: Ich suche einen Prozess. Können Sie beispielsweise eine Momenterzeugungsfunktion verwenden? Letztendlich würde ich gerne sehen, wie das geht, um das null aufgeblasene Gamma und andere besser zu verstehen.

B_Miner
quelle
1
Anscheinend kennen Sie ein Modell dafür, wie eine solche Wahrscheinlichkeitsverteilung entstehen würde. Können Sie das nutzen, um Ihnen zu helfen?
Kardinal

Antworten:

21

Methode 0 : Der faule Statistiker.

Man beachte , daß für haben wir f ( y ) = ( 1 - π ) p y wobei p y ist die Wahrscheinlichkeit , dass ein Poisson Zufallsvariable Wert annimmt y . Da der Ausdruck, der y = 0 entspricht, den erwarteten Wert nicht beeinflusst, sagt uns unsere Kenntnis des Poisson und der Linearität der Erwartung sofort, dass μ = ( 1 - π ) λ und E Y 2 = ( 1 -y0f((y)=((1- -π)pypyyy=0

μ=((1- -π)λ
EY2=(1π)(λ2+λ).

Eine kleine Algebra und die Identität ergeben das Ergebnis.Var(Y)=EY2μ2

Methode 1 : Ein probabilistisches Argument.

Es ist oft hilfreich, ein einfaches Wahrscheinlichkeitsmodell für die Entstehung einer Verteilung zu haben. Sei und Y P o i ( λ ) unabhängige Zufallsvariablen. Definiere X = Z Y.ZBer(1π)Y.P.Öich((λ) Dann ist leicht zu erkennen, dass X die gewünschte Verteilung f hat . Um dies zu überprüfen, ist zu beachten, dass P ( X = 0 ) = P ( Z = 0 ) + P ( Z = 1 , Y = 0 ) = π + ( 1 - π ) e - λ durch Unabhängigkeit. In ähnlicher Weise ist P ( X = k ) = P ( Z =

X.=Z.Y..
X.fP.((X.=0)=P.((Z.=0)+P.((Z.=1,Y.=0)=π+((1- -π)e- -λ für k 0 .P.((X.=k)=P.((Z.=1,Y.=k)k0

Daraus ist der Rest einfach, da durch die Unabhängigkeit von und Y , μ = E X = E Z Y = ( E Z ) ( E Y ) = ( 1 - π ) λZ.Y. Und V a r ( X ) = E X 2 - μ 2 = ( E Z ) ( E Y 2 ) - μ 2 = ( 1 - π ) ( λ 2 + λ ) - μ 2 = μ + π

μ=EX=EZY=(EZ)(EY)=(1π)λ,
Var(X)=EX2μ2=(EZ)(EY2)μ2=(1π)(λ2+λ)μ2=μ+π1πμ2.

Methode 2 : Direkte Berechnung.

λ

μ=k=1((1- -π)ke- -λλkk!=((1- -π)λe- -λj=0λjj!=((1- -π)λ.

E.X.2=((1- -π)k=1k2e- -λλkk!=((1- -π)λe- -λj=0((j+1)λjj!=((1- -π)((λ2+λ),

Nachtrag : Hier werden einige Tricks beschrieben, die in den obigen Berechnungen verwendet wurden.

k=0λkk!=eλ

k=0kλkk!=k=1kλkk!=k=1λk((k- -1)!=k=1λλk- -1((k- -1)!=λj=0λjj!=λeλ,
j=k- -1

E.X.((n)=E.X.((X.- -1)((X.- -2)((X.- -n+1)

eλE.X.((n)=k=nk((k- -1)((k- -n+1)λkk!=k=nλnλk- -n((k- -n)!=λnj=0λjj!=λneλ,
E.X.((n)=λnn0k<nk((k- -1)((k- -n+1)=0 da genau ein Term im Produkt Null ist.
Kardinal
quelle
Kardinal, das ist fantastisch. Würde es Ihnen etwas ausmachen, ein kurzes Detail zum Herausziehen des zu geben?λ? Meine Summe ist <sehr> rostig. Vielen Dank!
B_Miner
Nochmals vielen Dank dafür. Dies mag eine einfache Frage sein, aber was passiert mit dem oberen Teil des PDF (wenn y = 0)π+((1- -π)e- -λ Warum wird es nicht in die Berechnung für einbezogen? μ?
B_Miner
1
Erinnern Sie sich an die Definition des erwarteten Werts für eine diskrete Zufallsvariable: μ=E.Y.=y=0yP.((Y.=y). So füry=0ist der Begriff im erwarteten Wert 0((π+((1- -π)e- -λ)=0.
Kardinal