11

假设我有一个定义如下的函数

myFunction = function(input1, input2, input3) {
    # do something with input and then return
}

现在我只想最小化myFunction第一个输入,即input1,同时修复其他参数。

在 R 中,似乎有一些预先打包的函数,如nlm,optim等。但文档并没有真正解释如何解决上述问题。例如,似乎optim只能通过 1 个输入最小化一个函数?我可能是错的,但请纠正我,并告诉我做这种最小化问题的建议方法。

非常感谢!

4

1 回答 1

11

为了在输入是向量时最小化输出,您可以使用optim.

myFunction = function(input1, input2, input3) sum(abs(input1 - 3))
o = optim(7:10, myFunction, input2=5, input3=6)
o$par
# [1] 2.999989 2.999995 3.000000 3.000001

optim(在我的示例中)的第一个参数是(猜测从哪里开始优化)7:10的起始值。input1第二个是函数,然后你可以传入固定的参数(如input2input3)。

在这个例子中,最小值变成了一个只有 3s 的向量(因为它最小化了sum(abs(input1 - 3)))。


如果您只执行一维最小化,您将使用该optimize函数,并在传递函数后将其他输入传递给它:

myFunction = function(input1, input2, input3) {
    return(abs(input1 - 3) + input2 + input3)
}

o = optimize(myFunction, 2, 3, interval=c(-100, 100))
print(o$minimum)
# [1] 3.000003

如果您想基于input2input3(不是函数的第一个参数)最小化,您可以将其他参数作为命名参数:

myFunction = function(input1, input2, input3) {
    return(abs(input1 - 3) + abs(input2 - 12) + input3)
}

print(optimize(myFunction, input1=2, input3=3, interval=c(-100, 100)))
print(o$minimum)
[1] 11.99998

如果您想知道为什么这是可能的,那是因为optimize定义为:

function (f, interval, ..., lower = min(interval), upper = max(interval), 
          maximum = FALSE, tol = .Machine$double.eps^0.25) 

...允许使用任何数量的附加参数,然后将其传递给myFunction.

于 2013-06-09T20:03:04.803 回答