我对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")

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

提前感谢

回答：

你可以修改你的函数，将传入的值存储到一个全局列表中。

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运行时查看值，你也可以在函数中添加一个打印语句。