Was ist der erwartete Wert einer modifizierten Dirichlet-Verteilung? (Integrationsproblem)

14

Es ist einfach, eine Zufallsvariable mit Dirichlet-Verteilung unter Verwendung von Gamma-Variablen mit demselben Skalenparameter zu erzeugen. Wenn:

XiGamma(αi,β)

Dann:

(X1jXj,,XnjXj)Dirichlet(α1,,αn)

Problem Was passiert, wenn die Skalenparameter nicht gleich sind?

XiGamma(αi,βi)

Wie ist dann die Verteilung dieser Variablen?

(X1jXj,,XnjXj)?

Für mich würde es ausreichen, den erwarteten Wert dieser Verteilung zu kennen.
Ich benötige eine ungefähre geschlossene algebraische Formel, die von einem Computer sehr schnell ausgewertet werden kann.
Nehmen wir an, eine Annäherung mit einer Genauigkeit von 0,01 ist ausreichend.
Sie können davon ausgehen, dass:

αi,βiN

Hinweis Kurz gesagt besteht die Aufgabe darin, eine Annäherung an dieses Integral zu finden:

f(α,β)=R+nx1jxjjβjαjΓ(αj)xjαj1eβjxjdx1dxn

Łukasz Lew
quelle
1
@ Łukasz Können Sie noch etwas zu den Parametern , und sagen ? Es ist möglich, exakte Ausdrücke für und damit die Erwartungen der Verhältnisse , aber für bestimmte Kombinationen der Parameter könnte man mit weniger Aufwand die Normal- oder Sattelpunktnäherungen verwenden. Ich denke nicht, dass es eine universelle Annäherungsmethode geben wird, weshalb zusätzliche Einschränkungen wünschenswert wären. α i β i j X jnαiβijXj
Whuber
X1 und sind korreliert, daher müssen wir das Integral selbst approximieren. ist oft eine kleine Zahl wie 1 oder 2 und manchmal auch 10000. Ähnlich wie aber normalerweise 10-mal größer als . jXjαiβiαi
Łukasz Lew
Das Problem ist mit small . Wenn alle groß sind, lautet die gute Näherung des gesamten Integrals:αiαiα1/β1jαj/βj
Łukasz Lew
@ Łukasz Wenn Sie den Ausdruck der Erwartung bewerten müssen, warum benötigen Sie eine algebraische Formel? Ich überlege, einen numerischen Trick anzuwenden, um die Erwartung zu erhalten, aber ich brauche ein Feedback :)
deps_stats
Ich muss es in meinem Programm oft auswerten. Es muss sehr schnell sein, dh keine Schleifen und möglichst nicht zu viele Unterteilungen.
Łukasz Lew

Antworten:

2

Nur eine erste Bemerkung: Wenn Sie Rechengeschwindigkeit wollen, müssen Sie normalerweise auf Genauigkeit verzichten. "Mehr Genauigkeit" = "Mehr Zeit" im Allgemeinen. Sowieso ist hier eine Annäherung zweiter Ordnung, sollte die "grobe" Annäherung verbessern, die Sie in Ihrem Kommentar oben vorgeschlagen haben:

=α j

E(XjiXi)E[Xj]E[iXi]cov[iXi,Xj]E[iXi]2+E[Xj]E[iXi]3Var[iXi]
=αjiβjβiαi×[11(iβjβiαi)+1(iαiβi)2(iαiβi2)]

BEARBEITEN Eine Erklärung für die obige Erweiterung wurde angefordert. Die kurze Antwort ist Wikipedia . Die lange Antwort ist unten angegeben.

schreibe . Jetzt brauchen wir alle Ableitungen zweiter Ordnung vonf. Die Derivate erster Ordnung "stornieren", da sie alle die VielfachenX-E(X)undY-E(Y) enthalten,die beide Null sind, wenn sie Erwartungen annehmen.f(x,y)=xyfXE(X)YE(Y)

2f

2fx2=0
2f
2fxy=1y2
2fy2=2xy3

Und so ist die Taylor-Reihe bis zur zweiten Ordnung gegeben durch:

xyμxμy+12(1μy22(xμx)(yμy)+2μxμy3(yμy)2)

Die erwarteten Erträge nehmen:

E[xy]μxμy1μy2E[(xμx)(yμy)]+μxμy3E[(yμy)2]

Which is the answer I gave. (although I initially forgot the minus sign in the second term)

probabilityislogic
quelle
This looks like exactly what I need. Can you explain how you got this expansion? I tried in a lot of ways and was unable to do that ...
Łukasz Lew