Eu gostaria de plotar regiões de confiança 2D (em 1-sigma, 2-sigma) para um modelo que eu ajustei aos dados. Eu usei o PyMC para gerar amostras posteriores de 50k MCMC para o meu modelo com 6 parâmetros.
Eu sei que o processo para criar regiões de confiança é algo semelhante a: 1.) criar um histograma das amostras no espaço 2D 2.) identificar contornos de iso-densidade 3.) a partir de um ponto de partida selecionado (por exemplo, a média) integrar para fora contornos perpendiculares a iso-densidade até que a fração desejada dos pontos de amostra esteja contida na região.
Existe uma função conveniente no mundo numpy / scipy / pymc / pylab / etc que criará o gráfico da região de confiança 2D? Como alternativa, onde posso encontrar um algoritmo codificado, ou ferramenta autônoma, que calculará os contornos para plotagem posterior?
