Me gustaría trazar regiones de confianza 2D (en 1-sigma, 2-sigma) para un modelo que he ajustado a los datos. He usado PyMC para generar 50k muestras posteriores de MCMC para mi modelo con 6 parámetros.
Sé que el proceso para crear regiones de confianza es algo similar a: 1.) crear un histograma de las muestras en el espacio 2D 2.) identificar contornos de iso-densidad 3.) desde un punto de inicio seleccionado (p. Ej., La media) integrar hacia afuera perpendicular a los contornos de iso-densidad hasta que la fracción deseada de puntos de muestra esté contenida en la región.
¿Existe una función conveniente en el mundo numpy / scipy / pymc / pylab / etc que creará el diagrama de región de confianza 2D? Alternativamente, ¿dónde puedo encontrar un algoritmo codificado o una herramienta independiente que calcule los contornos para su posterior trazado?
fuente
