有没有办法在R中使用nlminb时返回Hessian矩阵？

Question

我想用R中的nlminb来优化几个参数，并且还要计算每个参数估计的标准误差以获得t值。但是，nlminb不会像其他优化命令那样返回Hessian矩阵。反正有吗？提前感谢您的帮助。

Answer 1

hessian在numDeriv包中可以用数值计算hessian。