abhängig von der Summe, wie ist die Verteilung der negativen Binome

9

Wenn ein negatives Binomial sind, wie ist dann die Verteilung von gegebenx1,x2,,xn(x1,x2,,xn)

x1+x2++xn=N ?

N ist fest.

Wenn x1,x2,,xn Poisson sind, ist (x_1, x_2, \ ldots, x_n) abhängig von der Summe (x1,x2,,xn)multinomial. Ich bin nicht sicher, ob es für negatives Binomial gilt, da es sich um eine Mischung aus Poisson handelt.

Falls Sie wissen möchten, ist dies kein Hausaufgabenproblem.

qkhhly
quelle
2
Angesichts des Zusammenhangs zwischen Gamma-Verteilungen und dem Dirichlet würde ich zunächst vermuten, dass es sich - zumindest bei angemessenen Einschränkungen der negativen Binome - in einigen Fällen als Dirichlet-Multinom herausstellen könnte.
Glen_b -State Monica
Wenn Sie die Begriffe in Ihrem Beitrag und in meinem Kommentar googeln, werden einige Treffer angezeigt, die darauf hindeuten, dass dies eine fruchtbare Linie sein könnte.
Glen_b -State Monica

Antworten:

7

Entschuldigung für die späte Antwort, aber das nervte mich auch und ich fand die Antwort. Die Verteilung ist in der Tat Dirichlet-Multinomial und das individuelle Neg. Binomialverteilungen müssen nicht einmal identisch sein, solange ihr Fano-Faktor (Verhältnis von Varianz zu Mittelwert) identisch ist.

Lange Antwort:

Wenn Sie NB wie folgt parametrieren:

p(X=x|λ,θ)=NB(x|λ,θ)=(θ1λ+x1x)(11+θ1)x(θ11+θ1)θ1λ

Dann ist und undE(X)=λVar(X)=λ(1+θ)

i:XiNB(λi,θ) impliziert

XiNB(λi,θ)

Nehmen Sie dann die Wahrscheinlichkeit unter Berücksichtigung der Summe:

NB(xi|λi,θ)NB(xi|λi,θ)=(11+θ1)xi(θ11+θ1)θ1λi(θ1λi+xi1xi)(11+θ1)xi(θ11+θ1)θ1λi(θ1λi+xi1xi)==Γ(xi+1)Γ(θ1λi)Γ(θ1λi+xi)Γ(θ1λi+xi)Γ(xi+1)Γ(θ1λi)=DM(x1,...,xn|θ1λ1,...,θ1λn)

wobei die Dirichlet-Multinomial-Wahrscheinlichkeit ist. Dies ergibt sich einfach aus der Tatsache, dass sich mit Ausnahme der Multinomialkoeffizienten viele Terme im Bruch auf der linken Seite aufheben und Sie nur die Gammafunktionsterme erhalten, die zufällig mit der DM-Wahrscheinlichkeit identisch sind.DM

Beachten Sie auch, dass die Parameter dieses Modells nicht als Zunahme von bei gleichzeitiger Abnahme von allen identifizierbar sind, was zu genau der gleichen Wahrscheinlichkeit führt.θλi

Die beste Referenz, die ich dafür habe, sind die Abschnitte 2 bis 3.1 von Guimarães & Lindrooth (2007): Kontrolle der Überdispersion in gruppierten bedingten Logit-Modellen: Eine rechnerisch einfache Anwendung der Dirichlet-multinomialen Regression - sie ist leider kostenpflichtig, aber ich konnte nicht Finden Sie eine nicht paywalled Referenz.

Martin Modrák
quelle