11

在尝试实现一些优化算法时，我注意到 R 参考类中有一些奇怪的行为。在特定的初始化方法中似乎有一些幕后的解析魔法，这使得使用匿名函数变得困难。这是一个说明困难的示例：我定义了一个要优化的函数（f_opt），一个在其上运行 optim 的函数，以及一个将这两个作为方法的引用类。奇怪的行为在代码中会更清楚

f_opt <- function(x) (t(x)%*%x)

do_optim_opt <- function(x) optim(x,f)
do_optim2_opt <- function(x)
  {
   f(x) #Pointless extra evaluation
   optim(x,f)
  }

optClass <- setRefClass("optClass",methods=list(do_optim=do_optim_opt,
                                 do_optim2=do_optim2_opt,
                                 f=f_opt))
oc <- optClass$new()
oc$do_optim(rep(0,2)) #Doesn't work: Error in function (par)  : object 'f' not found
oc$do_optim2(rep(0,2)) #Works. 
oc$do_optim(rep(0,2)) #Parsing magic has presumably happened, and now this works too.

只是我，还是其他人看起来也像一个错误？

4

1 回答 1

6

R-devel 中的这篇文章似乎是相关的，有解决方法

do_optim_opt <- function(x, f) optim(x, .self$f)

似乎值得在 R-devel 上发帖。

于 2011-09-07T13:33:03.250 回答