pymc3 - 如何将HDI添加到traceplot？

Question

是否有内置方法可以在pymc3中为traceplot添加最高密度区间？

我想直接在我的traceplot上显示HDI，理想情况下是标签。 基本上，我想制作一个后验分布图，类似于“做贝叶斯数据分析”一书中的那些：

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

Answer 1

不是内置方式，但您可以检查代码here。特别plot_post

<强>更新

将此图添加到PyMC3中的函数（plot_posterior）。