如何计算后验样本的最大后验密度( HPD )间隔?我有四个参数,从后验参数分布产生1000个样本。现在如何在R软件中计算HPD。no applicable method for 'HPDinterval' applied to an object of class "list"
其中"winbugsresult“是一个包含后验样本的列表
我想在PyMC3中实现捕获和捕获的基本模型(您捕获100只动物并标记它们,然后在它们混合后释放它们并重新捕获其中的100种动物,并注释标记了多少动物)。这是我的密码:import pymc3 as pmK = 100 #marked animals in first roundobs = 10 #observed marked in second round
with pm.Model() as my_mod