Bayesianischer Modellvergleich in der High School

8

Ich unterrichte Gymnasiasten in Physik und möchte, dass meine Schüler einen rudimentären Bayes'schen Modellvergleich für Daten aus ihren Experimenten durchführen. Ich habe einen Weg gefunden, wie sie das tun können (siehe unten), aber ich bin nicht sicher, ob es richtig ist. Ich würde mich über Feedback (insbesondere negatives Feedback!) Und Vorschläge zur besseren Verbesserung sehr freuen.

Ich möchte eine lineare Theorie mit den Parametern Steigung und Achsenabschnitt mit einer Nullhypothese einer Konstanten vergleichen, dh Steigung = 0. In beiden Fällen gehe ich von einem symmetrischen Gaußschen Rauschen aus.b aaba

Die Schüler können mithilfe von Excel die Schätzungen der maximalen Wahrscheinlichkeit für die Steigung und den Achsenabschnitt ( und ) sowie deren Fehler und . b deinedba^b^dadb

  1. Für den Prior auf der Steigung betrachte ich einen breiten Gaußschen Wert, der auf der Maximum-Likelihood-Schätzung ( ) zentriert ist und eine zehnfache Standardabweichung aufweist. Meine Argumentation ist, dass ich realistisch erwarte, dass sie die "richtigen" Linienparameter mindestens innerhalb einer Größenordnung finden, und in der Praxis werden sie diejenigen noch näher finden. Wenn ich also die "richtige" Steigung durch ihre MLE ersetze, werde ich die nicht ändern Zahlen zu viel.a^
  2. Für die Wahrscheinlichkeit der Evidenz bei einer bestimmten linearen Theorie betrachte ich die multivariate Standard-Gauß-Verteilung mit einer Standardabweichung ( ), die sich auf die Summe der quadratischen Residuen bezieht.σe
  3. Die Wahrscheinlichkeit des Beweises für die lineare Theorie im Allgemeinen, dh das Integral des obigen Prior und der Wahrscheinlichkeit, wird daher als der Prior und die Wahrscheinlichkeit am MLE-Punkt multipliziert mit dem Fehler in der Steigung geschätzt .da
  4. Die Wahrscheinlichkeit des Beweises bei gegebener Nullhypothese wird als ein weiterer multivariater Gaußscher angenommen, der nun die Gesamtstandardabweichung ( ) verwendet, basierend auf der Differenz zum Durchschnitt-Y.σT
  5. Dies ist der Teil, dessen ich mir am wenigsten sicher bin: Ich schätze den Bayes-Faktor als das Verhältnis der beiden oben genannten Wahrscheinlichkeiten (3 und 4 oben), wodurch ich die folgende Formel finden kann:

    B10=da(10|a^|2π)(σT/σe)Ne

Würde uns dies vernünftige Schätzungen für den Bayes-Faktor geben? Jedes Feedback ist willkommen.

Physik Lehrer
quelle
Ich habe Ihre Formel mit MathJax bearbeitet, um ein optimierteres Erscheinungsbild zu erzielen, das leichter zu lesen ist. Fühlen Sie sich frei, es zu bearbeiten, wenn ich es falsch übersetzt habe
Marquis de Carabas
Vielen Dank! Die letzten beiden Terme (das s-Verhältnis und die Quadratwurzel von e) sollten jedoch außerhalb des Bruchs oder im Zähler liegen.
Physiklehrer
1
Oh! Es ist genau wie bei LaTex! Ich habe die Formeln korrigiert; Danke noch einmal.
Physiklehrer

Antworten:

1

Lassen Sie mich zunächst sagen, dass ein vernünftiges Testen einer scharfen Hypothese wie eine sorgfältige vorherige Verteilung für a erfordert , da der Bayes-Faktor entscheidend von diesem Prior abhängt. Viele Bayesianer werden keine scharfe Hypothese prüfen, aber ich werde es tun.a=0a

Bevor ich fortfahre, muss ich Ihnen sagen, dass ich nicht wirklich verstehe, was Sie sagen, dass Sie tun, und daher kann ich Ihnen Ratschläge geben, nach denen Sie nicht suchen. Ich hoffe du kannst der Notation folgen.

Die Daten seien Beobachtungen: y = ( ( x 1 , y 1 ) , , ( x n , y n ) ) , wobei (gemäß dem allgemeineren Modell und einschließlich der Steigung) p ( y i | a , b , σ 2 ) = N ( y i | b + any=((x1,y1),,(xn,yn)) (Ich unterdrücke die unabhängige Variable x i aus der Liste der Konditionierungsargumente zur Vereinfachung der Notation.) Die Wahrscheinlichkeit ist gegeben durch p ( y | a , b , σ 2 ) = n i = 1 p ( y i | a , b , σ 2 ) . Gegeben ein Prior für ( a , b , σ 2 )

p(yi|a,b,σ2)=N(yi|b+axi,σ2).
xi
p(y|a,b,σ2)=i=1np(yi|a,b,σ2).
(a,b,σ2)ist die hintere Verteilung wobei die Wahrscheinlichkeit der Daten gemäß dem allgemeineren Modell p ( y ) ist
p(a,b,σ2|y)=p(y|a,b,σ2)p(a,b,σ2)p(y),
wo ichp(a,b,σ2)=p(a|b,σ2) verwendet habe
p(y)=p(y|a,b,σ2)p(a,b,σ)dσ2dbda=(p(y|a,b,σ2)p(b,σ2)dσ2db)p(a|b,σ2)da=p(y|a)p(a|b,σ2)da,
. Man beachte, dass p ( y | a ) die (marginale) Wahrscheinlichkeit für a ist und p ( a | b , σ 2 ) die bedingte Priorität für a ist . Wenn der Prior für a unabhängig von ( b , σ 2 ) ist , dann ist p ( a | b , σ 2 ) = p ( a )p(a,b,σ2)=p(a|b,σ2)p(b,σ2)p(y|a)ap(a|b,σ2)aa(b,σ2)p(a|b,σ2)=p(a). Ich werde annehmen, dass das wahr ist.

a

p(a|y)=p(y|a)p(a)p(y).
p(y|a)p(y)=p(a|y)p(a).
aa=0
p(y|a=0)p(y)=p(a=0|y)p(a=0).
a=0

a=0a=0aaa=0a=0

Sie können sich vorstellen, dass Sie diesem Problem nicht ausgesetzt sind, wenn Sie die von mir beschriebenen Schritte nicht befolgen, aber Sie würden sich irren. Die von mir vorgestellte Logik gilt unabhängig von dem von Ihnen angewendeten "Algorithmus".

p(b,σ2)1/σ2.
(b,σ2)ap(y|a)atyta

Ich hoffe, Sie finden etwas in dem, was ich gesagt habe, nützlich.

mef
quelle
p(a)=110|a^|2πe(aa^)22(10|a^|2p(a|y)=1σa|2πe(aa^)22σa2
Ich verstehe Ihren Prior nicht, da er die Daten über die Maximum-Likelihood-Schätzung einzubeziehen scheint.
Mef
a^hata
Ich verstehe die Gründe für Ihre Annahmen über den Prior nicht. Trotzdem lautet die Antwort auf die Frage in Ihrem ersten Kommentar "Ja". Ich denke, Sie werden feststellen, dass der Bayes-Faktor (BF) sehr empfindlich auf Ihre Wahl der vorherigen Varianz reagiert. Wenn Sie (zum Beispiel) 10 auf 20 ändern, werden Sie vermutlich eine große Änderung im BF erhalten. Und das ist der Punkt, den ich anstrebte.
Mef
a^