3

constrOptim.nl我正在使用来自alabama包的非线性优化。但是，我的问题更多地与传递参数有关（以及点-点-点（省略号/“...”），也许还有do.call）-所以我首先给出一个一般示例，然后再参考该constrOptim.nl函数。

假设，我有以下功能 - 我只能编辑second和third但不能first编辑.

first<-function (abc, second, third, ...){
second(abc,...)
third(abc,...)
}

second<- function(abc, ttt='nothing special'){
print(abc)
print(ttt)
}

third<- function(abc, zzz="default"){
print(abc)
print(zzz)
}

我想要的输出与我刚运行时得到的输出相同

second("test", ttt='something special')
third("test", zzz="non-default")

这是

"test"
"something special"
"test"
"non-default"

但是，下面的代码无法做到这一点。

first("test",second=second, third=third, ttt='something special',zzz="non-default")

如何更改调用或secondandthird函数以使其正常工作？

http://www.r-bloggers.com/r-three-dots-ellipsis/ 在这里我找到了一些可以帮助我的建议，do.call但目前我无法理解它应该如何工作。

我无法更改first函数，因为这是constrOptim.nl我的特定问题 - 它旨在能够将更多参数传递给不同的函数。但是，我可以更改secondandthird功能 - 因为它们是我想要最小化的限制和功能。显然我也可以改变函数的调用。

所以更具体地说，这是我的具体问题：我使用非线性限制执行最大似然估计：

minimize <- function(Param,VARresiduals){ 
   #Blahblah
   for (index in 1:nrow(VARreisduals)){
   #Likelihood Blahbla
   }
  return(LogL)
}
heq<-function(Param,W){   
  B<-Param[1:16]
  restriction[1]<-Lrestriction%*%(diag(4)%x%(solve(W))%*%as.vector(B))
  restriction[2:6]<-#BlablaMoreRestrictions
  return(restriction)
}

现在我打电话给constrOptim.nl...

constrOptim.nl(par=rnorm(20), fn=minimize,hin=NULL heq=heq,VARresiduals,W)

...但是得到相同的错误，正如我在调用first上面的函数时收到的那样 - 类似：“第二个错误（abc，...）：未使用的参数（zzz =“非默认”）”。

我该如何更改minimize和heq/或通话？：） 提前致谢

帖子被标记为重复后更新：相关帖子的答案更改了first我的示例中的函数 - 因为它在那里实现了一个 do.call，它调用了其他函数。但是，我无法更改示例中的第一个功能，因为我想保持constrOptim.nl各种不同功能的工作。还有其他方法吗？

4

1 回答 1

0

我想出的解决方案不是很优雅，但它确实有效。

second_2<- function(abc, extras){
  a<-extras[[1]]
  print(abc)
  print(a)
}

third_2<- function(abc, extras){
  a<-extras[[2]]
  print(abc)
  print(a)
}

extras<-list()
extras[[1]]<-'something special'
extras[[2]]<-"non-default"
extras


first("test",second=second_2, third=third_2, extras)

令人惊讶的是，以下代码也有效，但结果略有不同

 first("test",second=second, third=third, extras)

毕竟，现在设置默认值有点笨拙，但并非不可行。

于 2016-05-20T12:16:22.020 回答