我想为我适合数据的模型绘制 2D 置信区域（在 1-sigma、2-sigma 处）。我使用 PyMC 为我的模型生成了 50k MCMC 后验样本，有 6 个参数。

我知道创建置信区域的过程类似于：1.) 在 2D 空间中创建样本的直方图 2.) 识别等密度轮廓 3.) 从选定的起点（例如，平均值）向外整合垂直于等密度轮廓，直到所需的样本点分数包含在该区域中。

在 numpy/scipy/pymc/pylab/etc 世界中是否有一个方便的函数可以创建 2D 置信区域图？或者，我在哪里可以找到一个编码算法或独立工具来计算轮廓以供以后绘图？