r - 从 R 的优化函数中获取更多细节

Question

我对 optim 函数不是很熟悉，我想从它的结果中获取这些信息:a) 需要多少次迭代才能获得结果？ b) 绘制部分解的序列，即每次迭代结束时得到的解。

到目前为止，我的代码如下所示:

  f1 <- function(x) {
  x1 <- x[1]
  x2 <- x[2]
  x1^2 + 3*x2^2
}

res <- optim(c(1,1), f1, method="CG")

如何改进它以获取更多信息？

提前致谢

Answer 1

您可以修改您的函数以将传递给它的值存储到全局列表中。

i <- 0  
vals <- list()
f1 <- function(x) {
  i <<- i+1
  vals[[i]] <<- x

  x1 <- x[1]
  x2 <- x[2]
  x1^2 + 3*x2^2  
}

res <- optim(c(1,1), f1, method="CG")

现在，如果您在运行该函数后检查 i 和 vals，您可以看到发生了什么。如果您想在 optim 运行时查看值，也可以将打印语句放入函数中。