私はRを初めて使用するので、質問が少しばかげている場合はご容赦ください。値関数の反復のための単純な while ループを作成しようとしています。私の関数 (optim.routine) は、ソルバー ipoptr を使用しています。これが私のコードです:
d<-1
old1<-0
old2<-0
num.iter<-0
i.esp<-1e-05
i.T<-100
lb<-0
ub<-10
while (d>i.eps & num.iter<i.T){
new1 <- optim.routine(old1, old2, eval_f=eval_f, eval_grad_f=eval_grad_f, lb=lb, ub=ub, update=FALSE)
d<-dist(c(old1, new1), method="euclidean")
num.iter<-num.iter+1
old1<-new1
}
ここで、optim.routine は次の関数です。
optim.routine<-function(old1, old2, eval_f=obj, eval_grad_f=obj.deriv, lb=lb, ub=ub, update){
if (isTRUE(update)){
var2<-old2
var1<-old1
var1.deriv<-deriv(var1)
optimize <- ipoptr(x0 = old2, eval_f = eval_f, eval_grad_f = eval_grad_f, lb = lb,
ub = ub)
new1<- optimize$objective
new2<- optimize$solution
old2<-new2
old1<-new1
}else{
var2<-old2
var1<-old1
var1.deriv<-vf.deriv(var1)
optimize <- ipoptr(x0 = old2, eval_f = eval_f, eval_grad_f = eval_grad_f, lb = lb,
ub = ub)
new1<- optimize$objective
new2<- optimize$solution
old1<-new1
}
}
deriv は導関数を計算する関数です。
コードを実行しようとすると、次のエラーが発生します。
source('/mnt/ide0/home/myname/Documents/optim.R')
Error in get(as.character(FUN), mode = "function", envir = envir) :
object 'fn' of mode 'function' was not found
関数をデバッグすると:
Browse[2]> n
Error in isTRUE(update) : argument "update" is missing, with no default
while ループを使用せずに関数のみをソースすると、エラーは表示されません。正直なところ、私には手がかりがありません。どんな助けでも大歓迎です。ありがとう!
クラウディア