假设我有一个定义如下的函数
myFunction = function(input1, input2, input3) {
# do something with input and then return
}
现在我想在修复其他参数的同时,仅对第一个输入(myFunction
)最小化input1
。
在R中,似乎有一些预先包装的函数,如nlm
,optim
等。但是文档并没有真正解释如何解决上述问题。例如,optim
似乎只能在一个输入上最小化一个函数?我可能错了,但请纠正我,并告诉我这种最小化问题的建议方法。
非常感谢!
答案 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
的起始值(猜测从哪里开始优化)。第二个是函数,然后您可以传入固定参数(例如input2
和input3
)。
在这个例子中,最小值只是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
如果要根据input2
或input3
(不是函数的第一个参数)进行最小化,可以将其他参数作为命名参数:
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
。