Der schnellste Weg, um das Bayes-Schätzerproblem zu lösen

7

Das folgende Problem stammt aus einer alten Doktorandenprüfung in unserer Abteilung. Meine eigene Lösung unten ist zeitaufwändig und möglicherweise falsch. Es hängt auch davon ab, dass eine weniger verbreitete Verteilung erkannt wird. Ich frage mich also: Gibt es einen schnelleren Weg, um dieses Problem zu lösen, und wie wird das gemacht? Ich vermute, dass es anwendbare Theoreme über Bayes-Schätzer gibt, die ich nicht kenne oder vergessen habe.

Problem

Gegeben θ>0, Lassen Y1,,Ynwerden iid aus der Verteilung der Dichte

f(yθ)=Iy>0cθexp{θy2/2},
wo c=(2π)1/2. Angenommen, die Verlustfunktion zum Schätzenf(0θ)=cθ ist gegeben durch L(δ,cθ):=(δcθ)2/(cθ)2 und der Prior auf θ ist die Gammaverteilung mit Parametern α>1,β>0.

Meine Lösung

Lassen γ=f(0θ). Dann

f(yγ)=γnexp{γ2iyi2/(2c2)}
und durch eine Dichtetransformation mit Jacobian 2γ wir bekommen

p(γ)γ2α1exp{βγ2/c2},

damit

p(γy)γn+2α1exp{γ2[β/c2+yy/(2c2)]}.
Jetzt minimieren EL(δ,γ) unter diesem hinteren wrt δ ist gleichbedeutend mit Minimieren E(δγ)2 (dh ergibt die gleiche Zielfunktion) unter dem Seitenzahn

p(γy)γn+2α3exp{γ2[β/c2+yy/(2c2)]}.

Wir kennen diese Einstellung δ=E(γy)löst ein solches Problem. Wir erkennen auch den neuesten Posterior als verallgemeinertes Gamma mit Parametern (gemäß der Wikipedia- Notation).d=n+2α2,p=2,1/a=β/c2+yy/(2c2). Der Mittelwert des verallgemeinerten Gammas istaΓ((d+1)/p)/Γ(d/p)und wir sind fertig.

ekvall
quelle
Entschuldigung, ich kann nicht anders, als die Frage zu lösen, aber aus Interesse, für welche Art von Doktorandenprogramm war das? War es in Mathematik oder Statistik? Ich frage nur als Student.
Chris C
@ ChrisC ist aus einem Statistikprogramm. Denken Sie nicht, dass man mit diesem Missbrauch der Notation in einem Mathematikprogramm davonkommt, aber ich bin mir nicht sicher :).
Ekvall
Ha, danke! Ich muss definitiv etwas auffrischen. Hoffe du findest eine Antwort!
Chris C

Antworten:

0

Schneller gefunden:

Wir wollen minimieren

r(δ): =E.(cθ- -δcθ)2=E.(c2θ- -2cθδ+δ2c2θ)
unter der posterioren Verteilung. Für das No-Data-Problem bekommen wir
r(δ)=1- -δ2c- -1E.θ- -1/.2+δ2c- -2E.θ- -1,

was bei minimiert wird

δ=cE.θ- -1/.2E.θ- -1.

Für eine gammaverteilte Zufallsvariable X.Geinmmein(ein,b)Wir haben für Parameterwerte, bei denen das Integral existiert:

E.X.- -c=beinΓ(ein)xein- -1- -ce- -bx=beinΓ(ein- -c)bein- -cΓ(ein).

Da der hintere auch Gamma ist, nämlich.

p(θy)θα+n/.2e- -θ(β+ichyich2/.2),

Daraus ergeben sich unmittelbar die erforderlichen Momente und geben δ(y).

ekvall
quelle