r - 在 R 中，optim 不会返回函数的最大值，并且 uniroot 返回错误的值

Question

功能是：

f1 = function(x) {
  -1.3 * (x-0.1)^2+0.5 * (x-0.1)^5
}

我试图找到区间 [-1, 1] 中的最大值。优化函数返回正确的值：

optimize(f1, interval = c(-1, 1), maximum = T)

这给出了结果（这是正确的）：

$maximum
[1] 0.09999769

$objective
[1] -6.942984e-12

然而，我被要求（这个问题是家庭作业）使用 optim 和 uniroot。尝试使用优化：

optim(par = c(-1, 1), fn = f1)

导致此错误消息（无论我做什么都会发生这种情况）：

Error in optim(par = c(-1, 1), fn = f1) : 
  objective function in optim evaluates to length 2 not 1

虽然 uniroot 给出了一个明显错误的答案：

uniroot(f1, lower = -1, upper = 1, f.lower = -1,  f.upper =1)

结果（明显错误）：

$root
[1] 0.9998795

$f.root
[1] -0.7576706

那么这里的问题是什么？

标签： r

Answer 1

你可以试试这个

> optim(0, f1, lower = -1, upper = 1, method = "L-BFGS-B", control = list(fnscale = -1))
$par
[1] 0.1

$value
[1] -4.806839e-26

$counts
function gradient
       5        5

$convergence
[1] 0

$message
[1] "CONVERGENCE: REL_REDUCTION_OF_F <= FACTR*EPSMCH"