时间戳

首页 > 解决方案 > 求取两个输入变量的函数的局部最大值，每个变量的区间不同

r - 求取两个输入变量的函数的局部最大值，每个变量的区间不同

问题描述

您如何找到以下函数的局部最大值，两个输入具有不同的间隔？

f <- function(x, y) {
  y/50*(100*x)^0.9 + (50-y)/y*(80*(10-x))^0.8
}
# interval for y = (0, 50)
# interval for x = (0, 10)

我研究了使用 optim 函数，但我不知道如何为两个输入变量设置“par”、间隔和其他参数。

标签： roptimizationminimize

解决方案

optim期望一个带有vector参数（+可选参数）的函数，即要使问题中给出的函数与一个函数一起工作，optim必须稍微改变函数或使用某种包装函数，如下所示：

# Function of two scalar inputs
f_xy <- function(x, y) {
  y / 50 * (100 * x)^0.9 + (50 - y) / y * (80 * (10 - x))^0.8
}

# Wrapper or helper function with vector argument
f <- function(x) {
  y <- x[2]
  x <- x[1]

  f_xy(x, y)
}

# Default optim with starting value c(x = 5, y = 20)

optim(c(5, 20), fn = f)

注意：默认情况下optim执行最小化。

最大化：设置control = list(fnscale = -1)（感谢@alistaire 在评论中指出这一点。）

optim(c(5, 20), fn = f, control = list(fnscale = -1))

但是，对于问题中给出的函数，optim输出显示算法尚未收敛（$convergence != 0）（感谢@alistaire 在评论中指出这一点。）：

optim(c(5, 20), fn = f, control = list(fnscale = -1))
#> $par
#> [1] 1.591824e+00 3.861200e-34
#> 
#> $value
#> [1] 2.368542e+37
#> 
#> $counts
#> function gradient 
#>      501       NA 
#> 
#> $convergence
#> [1] 1
#> 
#> $message
#> NULL

在给定的情况下$convergence，1意味着达到了最大迭代次数。可以通过例如设置来增加最大迭代次数来解决这个问题control = list(..., maxit = 500))。但是，这并不能解决问题，因为optim仍然无法收敛。

推荐阅读