Berechnung von 2D-Konfidenzregionen aus MCMC-Proben


8

Ich möchte 2D-Konfidenzbereiche (bei 1-Sigma, 2-Sigma) für ein Modell zeichnen, das ich an Daten angepasst habe. Ich habe PyMC verwendet, um 50.000 MCMC-Posterior-Samples für mein Modell mit 6 Parametern zu generieren.

Ich weiß, dass der Prozess zum Erstellen von Vertrauensbereichen ähnlich ist wie: 1.) Erstellen eines Histogramms der Proben im 2D-Raum 2.) Identifizieren von Isodichtekonturen 3.) von einem ausgewählten Startpunkt (z. B. dem Mittelwert) nach außen integrieren senkrecht zu Isodichtekonturen, bis der gewünschte Anteil der Probenpunkte in der Region enthalten ist.

Gibt es eine praktische Funktion in der Welt von numpy / scipy / pymc / pylab / etc, mit der das 2D-Konfidenzbereichsdiagramm erstellt wird? Wo finde ich alternativ einen codierten Algorithmus oder ein eigenständiges Tool, das die Konturen für das spätere Plotten berechnet?

Antworten:


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.