Helfen Sie mir, die Bayes'schen Prior- und Posteriorverteilungen zu verstehen

125

In einer Gruppe von Studenten sind 2 von 18 Linkshändern. Finden Sie die posteriore Verteilung der Linkshänder in der Bevölkerung unter der Annahme, dass sie zuvor nicht informativ waren. Fassen Sie die Ergebnisse zusammen. Laut Literatur sind 5-20% der Menschen Linkshänder. Berücksichtigen Sie diese Informationen in Ihrem vorherigen und berechnen Sie den neuen Seitenzahn.

Ich weiß, dass die Beta-Distribution hier verwendet werden sollte. Erstens mit α und β Werten als 1? Die Gleichung, die ich im Material für posterior gefunden habe, ist

π(r|Y)r(Y+1)×(1r)(NY+1)

Y=2 ,N=18

Warum ist das r in der Gleichung? ( r bezeichnet den Anteil der Linkshänder). Es ist unbekannt, wie kann es in dieser Gleichung sein? Mir scheint es lächerlich zu berechnen r gegeben Y und verwenden Sie diese r in der Gleichung geben r . Nun, mit der Probe r=2/18 Das Ergebnis war 0,0019 . Das f soll ich daraus ableiten?

Die Gleichung, die einen erwarteten Wert von R bei bekanntem Y und N funktionierte besser und ergab 0,15 was ungefähr richtig klingt. Die Gleichung ist wobei der Wert und zugewiesen ist . Welche Werte sollte ich für und , um vorherige Informationen zu berücksichtigen?E(r|X,N,α,β)=(α+X)/(α+β+N)1αβαβ

Einige Tipps wären sehr dankbar. Ein allgemeiner Vortrag über frühere und spätere Verteilungen würde auch nicht schaden (ich verstehe vage, was sie sind, aber nur vage) Fortgeschrittene Mathematik wird wahrscheinlich über meinen Kopf fliegen.

Bob
quelle
4
Haben Sie sich diese Frage und Antwort angesehen ?
David Robinson
7
Der Satz " Finden Sie die hintere Verteilung von Linkshändern " macht keinen Sinn. Zufallsvariablen haben Verteilungen, und "Linkshänder" ist vermutlich nicht das Ziel von " Finden Sie die hintere Verteilung des Anteils der Linkshänder ". Es ist wichtig, solche Details nicht zu beschönigen, sondern klar zu sagen, wovon Sie tatsächlich sprechen.
Glen_b
2
Beim Lesen Ihrer Frage scheint es mir, dass Ihr Problem nicht so sehr Bayes'sche Statistiken sind, als dass Sie lediglich die Wahrscheinlichkeitsverteilungen verstehen. Es ist immer so, dass das Argument einer Verteilungsfunktion (oder einer Wahrscheinlichkeitsfunktion, wie Sie sie haben) eine Funktion eines Unbekannten (der Zufallsvariablen) ist. Das ist ganz und gar der Punkt von ihnen.
Glen_b
Kommentare sind nicht für eine längere Diskussion gedacht. Diese Unterhaltung wurde in den Chat verschoben .
gung

Antworten:

234

Lassen Sie mich zunächst erklären, was ein konjugierter Prior ist. Ich werde dann die Bayesianischen Analysen anhand Ihres konkreten Beispiels erläutern. Die Bayes'schen Statistiken umfassen die folgenden Schritte:

  1. Definieren Sie die vorherige Verteilung , die Ihre subjektiven Überzeugungen über einen Parameter enthält (in Ihrem Beispiel ist der interessierende Parameter der Anteil der Linkshänder). Der Prior kann "nicht informativ" oder "informativ" sein (aber es gibt keinen Prior, der keine Informationen enthält, siehe die Diskussion hier ).
  2. Daten sammeln.
  3. Aktualisieren Sie Ihre vorherige Verteilung mit den Daten unter Verwendung des Bayes-Theorems, um eine posteriore Verteilung zu erhalten. Die hintere Verteilung ist eine Wahrscheinlichkeitsverteilung, die Ihre aktualisierten Überzeugungen über den Parameter darstellt, nachdem Sie die Daten gesehen haben.
  4. Analysieren Sie die posteriore Verteilung und fassen Sie sie zusammen (Mittelwert, Median, SD, Quantile, ...).

Die Grundlage aller Bayes'schen Statistiken ist der Bayes'sche Satz

posteriorprior×likelihood

In Ihrem Fall ist die Wahrscheinlichkeit binomisch. Wenn sich die vorherige und die hintere Verteilung in derselben Familie befinden, werden die vorherige und die hintere Verteilung als konjugierte Verteilungen bezeichnet. Die Betaverteilung ist ein konjugierter Prior, da der Posterior ebenfalls eine Betaverteilung ist. Wir sagen, dass die Beta-Verteilung die konjugierte Familie für die Binomialwahrscheinlichkeit ist. Konjugatanalysen sind praktisch, treten jedoch in der Praxis nur selten auf. In den meisten Fällen muss die posteriore Distribution über MCMC numerisch gefunden werden (mit Stan, WinBUGS, OpenBUGS, JAGS, PyMC oder einem anderen Programm).

Wenn die Prioritätswahrscheinlichkeitsverteilung nicht zu 1 integriert wird, wird sie als unkorrekter Prioritätswert bezeichnet. Wenn sie zu 1 integriert wird, wird sie als korrekter Prioritätswert bezeichnet. In den meisten Fällen ist ein falscher Prior für Bayes-Analysen kein großes Problem. Die posteriore Verteilung muss jedoch korrekt sein, dh der posteriore Wert muss zu 1 integriert werden.

Diese Faustregeln ergeben sich direkt aus der Natur des Bayes'schen Analyseverfahrens:

  • Wenn der Prior nicht informativ ist, wird der Posterior sehr stark von den Daten bestimmt (der Posterior ist datengesteuert).
  • Wenn der Prior informativ ist, ist der Posterior eine Mischung aus dem Prior und den Daten
  • Je informativer die vorherigen Informationen sind, desto mehr Daten müssen Sie "ändern", sozusagen, weil der hintere Teil stark von den vorherigen Informationen bestimmt wird
  • Wenn Sie viele Daten haben, dominieren die Daten die posteriore Verteilung (sie überwältigen die vorherige).

Einen hervorragenden Überblick über einige mögliche "informative" und "nicht informative" Prioritäten für die Betaverteilung finden Sie in diesem Beitrag .

Angenommen, Ihre vorherige Betaversion ist wobei der Anteil der Linkshänder ist. Um die vorherigen Parameter und anzugeben , ist es hilfreich, den Mittelwert und die Varianz der Beta-Verteilung zu kennen (z. B. wenn Sie möchten, dass Ihre vorherige einen bestimmten Mittelwert und eine Varianz hat). Der Mittelwert ist . Wenn also , ist der Mittelwert . Die Varianz der Beta-Verteilung ist . Das Praktische ist nun, dass Sie an und denken könnenBeta(πLH|α,β)πLHαβπ¯LH=α/(α+β)α=β0.5αβ(α+β)2(α+β+1)αβwie zuvor beobachtete (Pseudo-) Daten, nämlich -Linkshänder und -Rechtshänder aus einer (Pseudo-) Stichprobe der Größe . Die Verteilung von ist gleichmäßig (alle Werte von sind gleich wahrscheinlich) und entspricht der Beobachtung von zwei Personen Davon ist einer Linkshänder und einer Rechtshänder.αβneq=α+βBeta(πLH|α=1,β=1)πLH

Die hintere Beta-Verteilung ist einfach wobei die Größe der Stichprobe und die Anzahl der Linkshänder in der Stichprobe ist. Das hintere Mittel von ist daher . Um die Parameter der posterioren Beta-Verteilung zu finden, fügen wir einfach Linkshänder zu und Rechtshänder zu . Die hintere Varianz istBeta(z+α,Nz+β)NzπLH(z+α)/(N+α+β)zαNzβ(z+α)(Nz+β)(N+α+β)2(N+α+β+1). Beachten Sie, dass ein sehr informativer Prior auch zu einer geringeren Varianz der posterioren Verteilung führt (die folgenden Grafiken veranschaulichen den Punkt gut).

In deinem Fall ist und und dein Prior ist die Uniform, die nicht aussagekräftig ist, also ist . Ihre hintere Verteilung ist daher . Der hintere Mittelwert ist . Hier ist eine Grafik, die den Prior, die Wahrscheinlichkeit der Daten und den Posterior zeigtz=2N=18α=β=1Beta(3,17)π¯LH=3/(3+17)=0.15

Der Prior, die Wahrscheinlichkeit der Daten und die hintere Verteilung mit einem einheitlichen Prior

Sie sehen, dass Ihre nachträgliche Verteilung vollständig von den Daten abhängt, da Ihre vorherige Verteilung nicht aussagekräftig ist. Dargestellt ist auch das höchste Dichteintervall (HDI) für die posteriore Verteilung. Stellen Sie sich vor, Sie legen Ihre hintere Verteilung in ein 2D-Becken und beginnen, Wasser einzufüllen, bis sich 95% der Verteilung über der Wasserlinie befinden. Die Punkte, an denen sich die Wasserlinie mit der posterioren Verteilung schneidet, bilden den 95% -HDI. Jeder Punkt innerhalb des HDI hat eine höhere Wahrscheinlichkeit als jeder Punkt außerhalb. Außerdem enthält der HDI immer den Peak der posterioren Verteilung (dh den Modus). Der HDI unterscheidet sich von einem gleichschwänzigen, zu 95% glaubwürdigen Intervall, bei dem 2,5% von jedem Schwanz des Seitenzahns ausgeschlossen sind (siehe hier ).

Bei Ihrer zweiten Aufgabe werden Sie gebeten, die Information, dass 5-20% der Bevölkerung Linkshänder sind, zu berücksichtigen. Dafür gibt es verschiedene Möglichkeiten. Der einfachste Weg ist zu sagen, dass die vorherige Beta-Verteilung einen Mittelwert von was dem Mittelwert von und . Aber wie wählt man und der vorherigen Betaverteilung aus? Zunächst möchten Sie, dass Ihr Mittelwert der vorherigen Verteilung aus einer Pseudo-Stichprobe mit der entsprechenden Stichprobengröße beträgt . Allgemeiner gesagt, wenn Sie möchten, dass Ihr Vorgänger einen Mittelwert mit einer Pseudo-Stichprobengröße , das entsprechende0.1250.050.2αβ0.125neqmneqαund -Werte sind: und . Jetzt müssen Sie nur noch die Pseudo-Stichprobengröße auswählen, die bestimmt, wie sicher Sie über Ihre vorherigen Informationen sind. Angenommen, Sie sind sich Ihrer vorherigen Informationen sehr sicher und setzen . Die Parameter Ihrer vorherigen Verteilung sind daher und . Die hintere Verteilung ist mit einem Mittelwert von etwa was praktisch dem vorherigen Mittelwert vonβα=mneqβ=(1m)neqneqneq=1000α=0.1251000=125β=(10.125)1000=875Beta(127,891)0.1250.125. Die Vorinformation dominiert den posterioren Bereich (siehe folgende Grafik):

Der Prior, die Wahrscheinlichkeit der Daten und die posteriore Verteilung mit stark informativem Prior

Wenn Sie sich bezüglich der vorherigen Informationen weniger sicher sind, können Sie das Ihres Pseudo-Samples auf beispielsweise , was zu und für Ihre vorherige Beta-Verteilung führt. Die hintere Verteilung ist mit einem Mittelwert von etwa . Der hintere Mittelwert liegt nun in der Nähe des Mittelwerts Ihrer Daten ( ), da die Daten den vorherigen überschreiten. Hier ist die Grafik, die die Situation zeigt:neq10α=1.25β=8.75Beta(3.25,24.75)0.1160.111

Der Prior, die Wahrscheinlichkeit der Daten und die posteriore Verteilung mit Beta-Prior entsprechen einer Pseudo-Stichprobengröße von 3

Eine fortgeschrittenere Methode zum Einbeziehen der vorherigen Informationen wäre, zu sagen, dass das Quantil Ihrer vorherigen Beta-Verteilung etwa und das Quantil etwa sollte . Dies entspricht der Aussage, dass Sie zu 95% sicher sind, dass der Anteil der Linkshänder an der Bevölkerung zwischen 5% und 20% liegt. Die Funktion im R-Paket berechnet die entsprechenden und Werte einer solchen Quantilen entsprechenden beta-Verteilung. Der Code ist0.0250.050.9750.2beta.selectLearnBayesαβ

library(LearnBayes)

quantile1=list(p=.025, x=0.05)     # the 2.5% quantile should be 0.05
quantile2=list(p=.975, x=0.2)      # the 97.5% quantile should be 0.2
beta.select(quantile1, quantile2)

[1]  7.61 59.13

Es scheint, dass eine Beta-Distribution mit den Parametern und die gewünschten Eigenschaften hat. Der vorherige Mittelwert liegt bei was in der Nähe des Mittelwerts Ihrer Daten ( ) liegt. Wiederum beinhaltet diese vorherige Verteilung die Information einer Pseudo-Stichprobe mit einer äquivalenten Stichprobengröße von ungefähr . Die posteriore Verteilung ist mit einem Mittelwert von der mit dem Mittelwert der vorherigen Analyse unter Verwendung eines hochinformativen vor vergleichbar ist. Hier ist die entsprechende Grafik:α=7.61β=59.137.61/(7.61+59.13)0.1140.111neq7.61+59.1366.74Beta(9.61,75.13)0.113Beta(125,875)

Der Prior, die Wahrscheinlichkeit der Daten und die hintere Verteilung mit dem Prior haben 0,05 und 0,975 Quantile von 0,05 und 0,2

Siehe auch diese Referenz für einen kurzen, aber sehr guten Überblick über das Bayes'sche Denken und eine einfache Analyse. Eine längere Einführung für Konjugatanalysen, insbesondere für Binomialdaten, finden Sie hier . Eine allgemeine Einführung in das Bayes'sche Denken finden Sie hier . Weitere Folien zu Aspekten der baysischen Statistik finden Sie hier .

COOLSerdash
quelle
1
Warum wählen wir hier die Beta-Distribution?
Metariat
1
@Metallica Der Hauptgrund ist, dass die Beta das Konjugat vor der Binomialverteilung ist. Dies bedeutet, dass, wenn wir eine Betaversion wie zuvor auswählen, die spätere Betaversion ebenfalls die Betaversion ist. Weitere Gründe sind, dass die Beta zwischen 0 und 1 liegt und sehr flexibel ist. Dazu gehört zum Beispiel die Uniform. Jede korrekte Verteilung mit Unterstützung in kann jedoch wie zuvor verwendet werden. Es ist nur so, dass der hintere Teil schwieriger zu berechnen ist. (0,1)
COOLSerdash
Haben Sie zufällig noch das Dokument für "Intro to Bayesian thinking"? Der Dropbox-Link ist nicht mehr verfügbar.
bs7280
@ bs7280 Ich habe die Links aktualisiert. Sie sollten jetzt wieder arbeiten.
COOLSerdash
1
@meduz Genau genommen gibt es keine wirklichen "uninformativen" Vorgänger. Ich möchte Sie auf die ausgezeichnete Antwort von Tim zu dieser Diskussion verweisen .
COOLSerdash
8

Eine Beta-Verteilung mit = 1 und = 1 entspricht einer Gleichverteilung. So ist es in der Tat uniformativ. Sie versuchen, Informationen über einen Parameter einer Verteilung zu finden (in diesem Fall den Prozentsatz der Linkshänder in einer Gruppe von Personen). Die Bayes-Formel lautet:αβ

P(r|Y1,...,n) =P(Y1,...,n|r)P(r)P(Y1,...,n|θ)P(r)

was Sie darauf hingewiesen haben, ist proportional zu:

α ( Y 1 , . . . , N | r ) * P ( r )P(r|Y1,...,n) (Y1,...,n|r)P(r)

Sie beginnen also im Grunde genommen mit Ihrer vorherigen Überzeugung, dass der Anteil der Linkshänder in der Gruppe (P (r), für die Sie einen einheitlichen Abstand verwenden), und berücksichtigen dann die Daten, die Sie sammeln, um Ihren Prior zu informieren (ein Binomial) in diesem Fall sind Sie entweder Rechts- oder Linkshänder, also ). Eine Binomialverteilung hat ein Beta-Konjugat vor, was bedeutet, dass die hintere VerteilungP ( r | Y 1 , . . . N ) αP(Y1,...,n|r)P(r|Y1,...n)gehört die Verteilung des Parameters nach Berücksichtigung der Daten zur selben Familie wie der vorherige. r hier ist am Ende nicht unbekannt. (Und ehrlich gesagt war es nicht vor dem Sammeln der Daten. Wir haben eine ziemlich gute Vorstellung vom Anteil der Linkshänder in der Gesellschaft.) Sie haben sowohl die vorherige Verteilung (Ihre Annahme von r) als auch die gesammelten Daten und füge die beiden zusammen. Der hintere Teil ist Ihre neue Annahme der Verteilung der Linkshänder nach Berücksichtigung der Daten. Sie nehmen also die Wahrscheinlichkeit der Daten und multiplizieren sie mit einer Uniform. Der erwartete Wert einer Beta-Distribution (das ist, was das Poster ist) ist . Also, als du angefangen hast, deine Annahme mit = 1 und αβ1αα+βαβ= 1 war, dass der Anteil der Linkshänder in der Welt . Jetzt haben Sie Daten gesammelt, bei denen 2 von 18 übrig sind. Sie haben einen Posterior berechnet. (noch eine Beta) Ihre Werte für und sind jetzt unterschiedlich und ändern Ihre Vorstellung vom Verhältnis von Linken zu Rechten. Wie hat es sich verändert? αβ12αβ

Eric Peterson
quelle
1

Im ersten Teil Ihrer Frage werden Sie aufgefordert, einen geeigneten Prior für "r" zu definieren. Mit den vorliegenden Binomialdaten wäre es ratsam, eine Betaverteilung zu wählen. Denn dann wird der hintere Teil eine Beta sein. Da es sich bei der gleichmäßigen Verteilung um einen speziellen Beta-Fall handelt, können Sie für "r" die gleichmäßige Verteilung vorziehen, sodass jeder mögliche Wert von "r" gleich wahrscheinlich ist.

Im zweiten Teil haben Sie die Informationen bezüglich der vorherigen Verteilung "r" bereitgestellt.

Mit dieser Antwort von @ COOLSerdash erhalten Sie die richtigen Anweisungen.

Vielen Dank, dass Sie diese Frage gestellt haben, und COOLSerdash für die richtige Antwort.

Nilupa Rupasinghe
quelle