對於R中具有不同邊界條件的多個變量，Optim（）函數

Question

有人可以幫助我爲具有不同邊界條件的兩個變量設置函數`optim（）'嗎？

x = rnorm(1:100) 
    y = rnorm(1:100)*50+2 
    #ndays is a vector of two values 

    test=function(ndays, x, y){ sum(x*ndays[1] - (x^2)*ndays[2] + y^2) } 

    optim(par = c(2,10),test, x=x, y=y, lower=c(1,10),upper=c(20,30), method="L-BFGS-B") 
    $par 
    [1] 20 30 

顯然，這種優化沒有工作的ndays [2]應該是10

如果optim()是不是適用於這裏，你可以建議一些其他的功能，良好的功能？

Answer 1

optim使函數test最小化，當ndays[2]增加時函數test減小。

所以得到30是正確的。此外ndays[1]將等於20，如果mean(x) < 0和等於1，如果mean(x) > 0