我一直在使用 R 包 optimx 进行统计模型估计。我的代码运行良好,我得到了所有 optimx 汇总结果(参数估计、函数值、收敛代码等)。但是,我还需要这些 optimx 结果的 Hessian 矩阵。我在 optimx 文档、R 站点和论坛中寻求帮助,但没有成功。有人能帮我吗?先感谢您!!!
这是一个示例(不是统计模型,但足以说明我的意思):
library(optimx)
fr <- function(x) {
## Rosenbrock Banana function
x1 <- x[1]
x2 <- x[2]
100 * (x2 - x1 * x1)^2 + (1 - x1)^2
}
grr <- function(x) {
## Gradient of ’fr’
x1 <- x[1]
x2 <- x[2]
c(-400 * x1 * (x2 - x1 * x1) - 2 * (1 - x1), 200 * (x2 - x1 * x1))
}
ans1<-optimx(c(-1.2,1), fr, grr)