后验分布MCMC的Fisher信息

Question

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

但是，有人建议我使用optim(..., Hessian=TRUE)但是我没有看到优化例程如何帮助。