У меня есть эта функция, чтобы максимизировать:

funcToOpt <- function (theta1,theta2,theta3)
{
  inner <- rowSums((wb + (theta1*X1+theta2*X2+theta3*X3) / Nt) * r)  
  innerU <- u(inner, gamma)                                     # u() is another function 

  return (sum(innerU) / T)
}

fbb <- function(x) funcToOpt(x[1],x[2],x[3])
max <- optim(c(0,0,0),fbb,control=list(fnscale=-1))

Куда:wb,X1,X2,X3,r являютсяTxN матрицы иNt вектор с длинойN

Я просто хочу это

(wb + (theta1*X1+theta2*X2+theta3*X3) / Nt)

не могу предположитьnegative ценности.

Есть кто-то, кто может мне помочь? Спасибо!