我想 (i) 计算和 (ii) 绘制 Distributions.jl 库中分布的中心可信区间和最高后验密度区间。理想情况下,可以编写自己的函数来计算 CI 和 HPD,然后使用 Plots.jl 绘制它们。但是,我发现实现非常棘手(免责声明:我是 Julia 的新手)。有什么关于库/gists/repo 的建议可以让计算和绘制它们更容易吗?
语境
using Plots, StatsPlots, LaTeXStrings
using Distributions
dist = Beta(10, 10)
plot(dist) # thanks to StatsPlots it nicely plots the distribution
# missing piece 1: compute CI and HPD
# missing piece 2: plot CI and HPD
预期的最终结果总结在下图或第页。BDA3的第 33 条。
目前找到的资源: