python - 事後分布MCMCのフィッシャー情報

Question

私が使用している方法では、事後分布のフィッシャー情報を計算する必要があります（すべてのハイパーパラメーターに関して）。現在私が持っているのは、事後分布からのモンテカルロサンプルです。これらを使用して、それぞれ2次および1次導関数のサンプル平均と共分散で近似できると思いますが、より効率的な方法を探しています。

ただし、使用することが提案されましたがoptim(..., Hessian=TRUE)、最適化ルーチンがどのように役立つかわかりません。