python - 最高后牙密度区和中央可信区

给定一些参数Θ的后验p(Θ| D)，则可以将以下内容表示为:
最高后密度区域:
最高后密度区域是Θ的最可能值的集合，这些值总计构成后部质量的100(1-α)％。
换句话说，对于给定的α，我们寻找满足的p *:

然后获取最高后密度区域作为集合:

中部可信地区:
使用与上述相同的符号，和可信区域(或间隔)定义为:

根据分布，可能会有很多这样的间隔。 中央可信区间定义为可信区间，其中在每个尾的上有(1-α)/2 质量。
计算方式:

对于常规发行版，从该发行版给出的样本中，是否有任何内置函数可使用Python或 define 获得上述两个数量？

对于常见的参数分布(例如Beta，Gaussian等)，是否有任何内置文件或库可使用PyMC或SciPy进行计算？

最佳答案

要计算HPD，您可以利用pymc3，这是一个示例

import pymc3
from scipy.stats import norm
a = norm.rvs(size=10000)
pymc3.stats.hpd(a)

关于python - 最高后牙密度区和中央可信区，我们在Stack Overflow上找到一个类似的问题：https://stackoverflow.com/questions/22284502/