Fisher-Informationen in einem hierarchischen Modell

20

Ausgehend von dem folgenden hierarchischen Modell sind und wobei ist eine Normalverteilung. Gibt es eine Möglichkeit, einen genauen Ausdruck für die Fisher-Information der Randverteilung von wenn . Das heißt, was ist die Fisher-Information von: Ich kann einen Ausdruck für die Randverteilung von gegebenem c erhalten , Aber es scheint sehr schwierig zu sein, zwischen wrt c und den Erwartungen zu unterscheiden. Vermisse ich etwas Offensichtliches? Jede Hilfe wäre dankbar.μ L a p l a c e ( 0 , c ) N ( , ) X c p ( x | c ) = p ( x | μ ) p ( μ | c ) d μ X c c

XN(μ,1),
μLaplace(0,c)
N(,)Xc
p(x|c)=p(x|μ)p(μ|c)dμ
Xcc
emakalisch
quelle
Ich habe es selbst ausprobiert, aber es ist jenseits meiner Fähigkeiten. Absolutwertfunktionen ruinieren alles! Sie sind im Grunde mit numerischen Methoden fest.
Wahrscheinlichkeitslogik
3
@probability Sie können einen Ausdruck für den Integranden erhalten, indem Sie das Integral einfach in die Bereiche μ0 und μ<0 . Es werden keine absoluten Werte benötigt. Aber das Ergebnis ist eine chaotische rationale Funktion von x , exp(x2) und Fehlerfunktionen, und es ist unwahrscheinlich, dass sie in geschlossener Form integriert werden kann.
whuber
1
@whuber - das habe ich mit "hoffnungslos" gemeint. Nicht dass das Integral unmöglich ist, aber die Fischerinformation ist unmöglich. Weil Sie den erwarteten Wert über X eines Verhältnisses von zwei dieser Arten von Integralen nehmen müssen
Wahrscheinlichkeitsrechnung
1
Eine Untergrenze für die Fisher-Information ist in diesem Fall 1/(1+2c2) . Ist es möglich, eine engere Obergrenze für die Fisher-Informationen zu erhalten als für die allgemeine 1+1/c2 ?
Emakalic
Während eine analytische Lösung eine Herausforderung hinsichtlich der menschlichen Traktabilität darstellt (außerhalb einer mathematischen Disziplin), gibt es eine Empfänglichkeit für eine ungefähre rechnerische Lösung? Man könnte eine stochastische Simulation durchführen und dann die Näherungen für die Anpassung betrachten.
EngrStudent - Wiedereinsetzung von Monica

Antworten:

2

Es gibt keinen geschlossenen analytischen Ausdruck für die Fisher-Informationen für das von Ihnen angegebene hierarchische Modell. In der Praxis können Fisher-Informationen nur für exponentielle Familienverteilungen analytisch berechnet werden. Für exponentielle Familien ist die Log-Wahrscheinlichkeit in den ausreichenden Statistiken linear, und die ausreichenden Statistiken haben bekannte Erwartungen. Bei anderen Distributionen vereinfacht sich die Log-Wahrscheinlichkeit auf diese Weise nicht. Weder die Laplace-Verteilung noch das hierarchische Modell sind exponentielle Familienverteilungen, daher ist eine analytische Lösung unmöglich.

Gordon Smyth
quelle
0

Die beiden Normalen und Laplace stammen aus der Exponentialfamilie. Wenn Sie die Verteilung in Exponentialform schreiben können, ist die Fisher-Informationsmatrix der zweite Gradient des Log-Normalisierers der Exponentialfamilie.

A. Yazdiha
quelle
Ich glaube nicht, dass der übliche Zwei-Parameter-Laplace in der Exponentialfamilie liegt. Wenn der Standortparameter bekannt ist, wird er in der Exponentialfamilie liegen, aber ich glaube, dass nicht in Exponentialfamilienform geschrieben werden kann. 12exp(|xμ|)
Glen_b