Was ist der Unterschied zwischen Metropolis Hastings, Gibbs, Importance und Rejection Sampling?


36

Ich habe versucht, MCMC-Methoden zu erlernen und bin auf Stichproben von Metropolis Hastings, Gibbs, Wichtigkeit und Ablehnung gestoßen. Während einige dieser Unterschiede offensichtlich sind, dh wie Gibbs ein Sonderfall von Metropolis Hastings ist, wenn wir die vollständigen Bedingungen haben, sind die anderen weniger offensichtlich, wenn wir MH in einem Gibbs-Sampler usw. verwenden möchten einfache Möglichkeit, den Großteil der Unterschiede zwischen diesen beiden zu sehen? Vielen Dank!


3
Iain Murray spricht dies in seinem Vortrag , zumindest in Bezug auf MCMC, gut an.
gwr

2
Ich stimme Xi'an zu, dass dies eine sehr breite Frage ist. Sie fragen effektiv nach einer Menge Informationen zu vier verschiedenen Dingen, von denen eine Diskussion (oder ein Kontrast zwischen zwei Dingen) zu einer etwas längeren Antwort führen würde. Wir könnten in der Lage sein, die Frage zu fokussieren, indem wir feststellen, dass alle vier Monte-Carlo-Methoden, wichtige Stichproben und Zurückweisungsstichproben keine MCMC sind (das heißt nicht, dass sie nicht in MCMC verwendet werden können).
Glen_b

Antworten:


47

Wie bereits in unserem Buch ausführlich mit George Casella, Monte Carlo statistischen Methoden werden diese Verfahren zur Herstellung von Proben aus einer bestimmten Verteilung verwendet wird , mit der Dichte sagen wir, entweder eine Vorstellung über diese Verteilung zu bekommen, oder eine Integration oder Optimierungsproblem im Zusammenhang zu lösen mit f . Zum Beispiel den Wert zu finden X h ( x ) f ( x ) d xff oder die Art der Verteilung von h ( X ), wenn X f ( x ) oder ein Quantil dieser Verteilung ist.

Xh(x)f(x)dxh(X)R
h(X)Xf(x)

Um die Monte-Carlo- und Markov-Ketten-Monte-Carlo-Methoden zu vergleichen, die Sie anhand relevanter Kriterien erwähnen, müssen Sie den Hintergrund des Problems und die Ziele des Simulationsexperiments festlegen, da die Vor- und Nachteile von Fall zu Fall unterschiedlich sind.

Hier einige allgemeine Bemerkungen , die die Komplexität des Problems mit Sicherheit nicht abdecken :

  1. fu1,u2,xfffxXf
  2. (xt)tfff
    f(x)Zf~(x,z)dz
    (xt)t(xt)txttfft
  3. g(x)
    f(x)/g(x).
    gfggf

I=Xh(x)f(x)dx,
I^=Xh(x)f(x)dx
f

f

Ich habe mich nur gefragt, was h(x)konkret h(x)f(x)dxin einem Bayes'schen Analyseszenario gemeint ist . Wir versuchen, den posterior zu bekommen, unter Berücksichtigung des Prior und der Daten. Es scheint jedoch, dass wir mit all diesen Stichprobenverfahren tatsächlich versuchen, uns anzunähern f(x). Kann man also sagen, dass dies f(x)bereits der Posterior ist, den wir suchen, und dass h(x)es sich nur um eine willkürliche Funktion handelt, die wir auch mit dem Posterior zusammenstellen könnten f(x)? Oder habe ich es nicht richtig verstanden. Vielen Dank.
xji

Xh(x)f(x)dx
fh
Durch die Nutzung unserer Website bestätigen Sie, dass Sie unsere Cookie-Richtlinie und Datenschutzrichtlinie gelesen und verstanden haben.
Licensed under cc by-sa 3.0 with attribution required.