關於optim（）函數取決於循環的R問題

Question

這是關於R上的optim（）函數的一個簡單問題。我試圖根據循環的變量找到函數的最小值。

f<-function(x,i){ 
    i*x 
} 

M<-matrix(0,10,1) 

    for (i in 1:5) { 
    M[i,1]<-optim(c(1),f,lower=-2,upper=2,method="L-BFGS-B") 
    } 
M 

[R回答我：

Erreur dans i * x : 'i' est manquant 

我真的不知道如何處理這一個....

Answer 1

你要通過i到optim電話：

optim(1,f,i=i,lower=-2,upper=2,method="L-BFGS-B") 

...和optim返回一個列表，所以你需要選擇哪一個元素nt在你想要存儲的列表中M。我想，你希望par：

M[i,1] <- optim(1,f,i=i,lower=-2,upper=2,method="L-BFGS-B")$par 

UPDATE：

賈斯汀在評論中指出，更多的R-樣的方式來做到這一點是：

L <- list() 
for (i in 1:5) { 
    L[[i]] <- optim(1,f,i=i,lower=-2,upper=2,method="L-BFGS-B") 
} 
V <- sapply(L, `[[`, "par")