如何最小化R中一个输入参数的函数

时间:2013-06-09 19:58:56

标签: r minimize

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

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

现在我想在修复其他参数的同时,仅对第一个输入(myFunction)最小化input1

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

非常感谢!

1 个答案:

答案 0 :(得分:7)

为了在输入为矢量时最小化输出,您可以使用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