为了最小化一个关于 3 个参数的函数，我使用了 optim 函数和“L-BFGS-B”方法。 这是错误消息：
优化错误(c(3, 0.01, 0.75), fn = f, lower = c(0.6, 0.0001, 0.1), : L-BFGS-B 需要 'fn' 的有限值
当三个参数中的至少一个到达“边界”（即下限值或上限值）时，我已经检查了函数的值，但它从未给出无限值。
如何知道优化函数中哪些值给出了无限值？
【问题讨论】：
标签： r optimization
一种快速而肮脏的方法是包装你的函数并打印出传递给它的参数。这是执行此操作的示例会话。
f <- function(par, data){
if(0.1 < par & par < .9){
return(Inf)
}else{
sqrt(abs(par))
}
}
并使用 optim 调用它...
> optim(2, f, method = "L-BFGS-B", lower = -1)
Error in optim(2, f, method = "L-BFGS-B", lower = -1) :
L-BFGS-B needs finite values of 'fn'
现在包装代码...
wrap_f <- function(par, ...){
cat(par, "\n")
f(par, ...)
}
现在我们可以看到叫什么了...
> optim(2, wrap_f, method = "L-BFGS-B", lower = -1)
2
2.001
1.999
1.646447
1.647447
1.645447
1.256777
1.257777
1.255777
-0.3018999
-0.3008999
-0.3028999
0.7441084
Error in optim(2, wrap_f, method = "L-BFGS-B", lower = -1) :
L-BFGS-B needs finite values of 'fn'
所以在这个例子中，它最终试图评估 .7441084，它给出了函数定义方式的错误。
【讨论】：