Je voudrais tracer des régions de confiance 2D (à 1-sigma, 2-sigma) pour un modèle que j'ai adapté aux données. J'ai utilisé PyMC pour générer des échantillons postérieurs de 50k MCMC pour mon modèle avec 6 paramètres.
Je sais que le processus de création de régions de confiance est quelque chose de similaire à: 1.) créer un histogramme des échantillons dans l'espace 2D 2.) identifier les contours d'iso-densité 3.) à partir d'un point de départ sélectionné (par exemple, la moyenne) intégrer vers l'extérieur perpendiculaire aux contours d'iso-densité jusqu'à ce que la fraction souhaitée de points d'échantillonnage soit contenue dans la région.
Existe-t-il une fonction pratique dans le monde numpy / scipy / pymc / pylab / etc qui créera le tracé de la région de confiance 2D? Sinon, où puis-je trouver un algorithme codé, ou un outil autonome, qui calculera les contours pour un tracé ultérieur?
la source