Aus einem frequentistischen Ergebnis einen Bayesianischen Prior machen

13

Wie sollte man ein frequentistisches Ergebnis in einen bayesianischen Prior verwandeln?

Stellen Sie sich das folgende ziemlich allgemeine Szenario vor: In der Vergangenheit wurde ein Experiment durchgeführt und ein Ergebnis für einen Parameter ϕ gemessen. Die Analyse wurde mit einer frequentistischen Methodik durchgeführt. Ein Konfidenzintervall für ϕ ist in den Ergebnissen angegeben.

Ich führe jetzt ein neues Experiment durch, bei dem ich einige andere Parameter messen möchte, z. B. θ und ϕ . Mein Experiment unterscheidet sich von der vorherigen Studie - es wird nicht mit der gleichen Methodik durchgeführt. Ich würde gerne eine Bayes'sche Analyse durchführen und muss daher θ und ϕ mit Prioren versehenϕ .

Keine vorherigen Messungen von θ durchgeführt, daher setze ich eine uninformative (sagen wir ihre Uniform) voran.

Wie bereits erwähnt, gibt es ein vorheriges Ergebnis für ϕ , das als Konfidenzintervall angegeben ist. Um dieses Ergebnis in meiner aktuellen Analyse zu verwenden, müsste ich das vorherige Ergebnis des Frequentisten in einen informativen Prior für meine Analyse übersetzen.

Eine Möglichkeit , die in diesem Szenario aus nicht verfügbar ist , ist die vorherige Analyse zu wiederholen , dass auf die führte ϕ in einem Bayes - Mode - Messung. Wenn ich dies tun könnte, ϕ würde ein posterior aus dem vorherigen Experiment hat , dass ich dann als meine vorherigen verwenden würde, und es gäbe keine Frage sein.

Wie soll ich das frequentistische CI für meine Analyse in eine bayesianische Vorabverteilung übersetzen? Oder mit anderen Worten, wie könnte ich ihr häufigstes Ergebnis am ϕ in ein posteriores am ϕ übersetzen ϕ , dass ich dann als vor meiner Analyse verwenden würde?

Alle Erkenntnisse oder Referenzen, die diese Art von Problem behandeln, sind willkommen.

bill_e
quelle
Vorherige oder nachherige Verteilung?
Tim
der Übersichtlichkeit halber bearbeitet, besser?
bill_e
Können Sie eine Uniform von-unendlich bis + unendlich
mdewey
Ich bin nicht sicher, was dies mit der Metaanalyse zu tun hat. Können Sie
erklären
3
Sie suchen nach passenden Priors im Welch- und Peers-Stil. Schauen Sie sich diese Rezension an: projecteuclid.org/euclid.lnms/1215091929
Zen

Antworten:

3

Kurzfassung: Nehmen Sie einen Gaußschen Wert, der auf der vorherigen Schätzung zentriert ist, mit std. dev. gleich dem CI.

Lange Version: Let der wahre Wert des Parameters, und sei φ die Schätzung , dass Sie haben. Man nehme a priori eine einheitliche Zahl vor P ( ϕ ) = c t . Sie möchten die Verteilung der wissen , φ 0 gegeben , dass eine Schätzung φ bereits erzielt:ϕ0ϕ^P(ϕ)=ctϕ0ϕ^

Nundie einzige Abhängigkeit vonφ0istin dem AusdruckP (

P(ϕ0|ϕ^)=P(ϕ^|ϕ0)P(ϕ0)P(ϕ^)=P(ϕ^|ϕ0)ctP(ϕ^)
ϕ0, der Resteine Normierungskonstante. UnterAnnahmedie φ wird einen MaximumLikelihoodSchätzer (oder einen anderen konsistenten Schätzer), können wir die folgenden Tatsachen verwenden:P(ϕ^|ϕ0)ϕ^
  1. Wenn die Anzahl der Beobachtungen zunimmt, ist der MLE asymptotisch
  2. Es ist asymptotisch unbefangen (zentriert auf den wahren Wert ),ϕ0
  3. Sie schwankt um mit einer Varianz, die der inversen Fisher-Information der vorherigen Beobachtungen entspricht, und das hätte ich als CI (Quadrat) verwendet.ϕ0

Anders ausgedrückt: Der Bayesianische Posterior und die Verteilung eines konsistenten und effizienten Schätzers werden asymptotisch gleich.

Alex Monras
quelle
Ich sollte hinzufügen, dass diese Lösung für 68% CI ist, was 1 Sigma ist. Wenn Ihr Konfidenzintervall 95% beträgt, sind Sie bei zwei Sigmen, also sollten Sie das CI durch 2 teilen, wenn sie bei 99,7% liegen, dann sind es 3 Sigmen, also sollten Sie durch 3 teilen. De.wikipedia.org/wiki/ 68% E2% 80% 9395% E2% 80% 9399.7_Rule
Alex Monras
Ich sollte genau das kommentieren, was in Ihrem Kommentar steht :-) Vielleicht sollten Sie das Ihrer Antwort hinzufügen. Ich würde ...
Rolazaro Azeveires
1

ttσ2S2(np)/σ21/σ2

Jarle Tufto
quelle