给定一些参数Θ的后p(Θ| D),可以定义以下内容:
的最高后验密度区域是集合Θ的最可能值,在总构成后部100质量(1-α)%的.
换句话说,对于给定的α,我们寻找满足以下条件的p*:
然后获得最高后部密度区域作为集合:
使用与上述相同的表示法,可信区域(或区间)定义为:
根据分布,可能有许多这样的间隔.中心可信区间定义为每个尾部有(1-α)/ 2质量的可信区间.
对于一般分布,从分布中给出样本,是否有任何内置函数可以在Python或PyMC中获得上述两个量?
对于常见的参数分布(例如Beta,Gaussian等),是否有任何内置函数或库可以使用SciPy或statsmodels进行计算?
python statistics scipy statsmodels pymc
pymc ×1
python ×1
scipy ×1
statistics ×1
statsmodels ×1