matplotlib - pymc3 - 如何将 HDI 添加到 traceplot？

Question

是否有内置方法可以将最高密度间隔添加到 pymc3 中的跟踪图？

我想直接在我的跟踪图中显示 HDI，最好是带有标签。基本上，我想制作一个后验分布图，类似于“做贝叶斯数据分析”一书中的那些：

我不想重新发明可能已经在 pymc3 中实现的东西。

Answer 1

不是内置方式，但您可以在此处查看代码。特别是 plot_post

更新

plot_posteriorPyMC3 中添加了一个用于绘制此图的函数 ( )。