r - R での制約の使用

Question

ここの単純なもの、私はかつてこれを知っていましたが、何年にもわたって失われてきました.

Rでコーディングしやすい単純な方程式：

f(x,y) = 2x^2 + 4y^2 + 6x - 8y + 15

と の制約がx > 1ありy > -1ます。

Rで制約を適切に記述する方法を一生覚えられず、持っている本は役に立たない

どんな助けにも乾杯

最小値と最大値を求める

Answer 1

引数の単一ベクトルを取る関数を定義します。

myfun <- function(xy) { 
  x <- xy[1] 
  y <- xy[2] 

  2*x^2 + 4*y^2 + 6*x - 8*y + 15
}

に開始値を指定しoptim、 と の下限を指定しxますy。

starting_values <- c(0, 0)
optim(starting_values, myfun, lower=c(1, -1), method='L-BFGS-B')

optim出力：

$par
[1] 1 1

$value
[1] 19

$counts
function gradient 
       2        2 

$convergence
[1] 0

$message
[1] "CONVERGENCE: NORM OF PROJECTED GRADIENT <= PGTOL"