python - 后验分布 MCMC 的 Fisher 信息

Question

我使用的方法需要我计算后验分布的 Fisher 信息（关于所有超参数）。我目前拥有的是来自后验分布的蒙特卡洛样本。我觉得我可以分别通过二阶和一阶导数的样本均值和协方差来使用和近似这些，但我正在寻找一种更有效的方法。

但是，有人建议我使用，optim(..., Hessian=TRUE)但是我看不出优化例程有什么帮助。