Me gustaría trazar regiones de confianza 2D (a 1-sigma, 2-sigma) para un modelo que he ajustado a los datos. He utilizado PyMC para generar 50k muestras posteriores 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 los contornos de isodensidad 3.) a partir de un punto de partida seleccionado (por ejemplo, la media) integrar hacia fuera perpendicularmente a los contornos de isodensidad hasta que la fracción deseada de puntos de muestra esté contenida en la región.
¿Existe una función conveniente en el mundo de numpy/scipy/pymc/pylab/etc que cree el gráfico de la región de confianza en 2D? Alternativamente, ¿dónde puedo encontrar un algoritmo codificado, o una herramienta independiente, que compute los contornos para su posterior trazado?
