2

パラメータの合計を 1 に設定して、1 つの関数を最適化したい

varxyz<-function(param){
a<-param[1]
b<-param[2]
c<-param[3]
a^2*matcov[1,1]+b^2*matcov[2,2]+c^2*matcov[3,3]+2*a*b*matcov[1,2]+2*a*c*matcov[1,3]+2*b*c*matcov[2,3]   
}  


parammv <- optim(c(0.3,0.3,0.4),varxyz,method="L-BFGS-B",lower=c(0.1,0.1,0.1))

ここで、matcov は共分散行列です

3 つのパラメーターの合計を 1 に設定するにはどうすればよいですか? 回答ありがとうございます

4

3 に答える 3

2

二次最適化問題なので、使えますquadprog。

# Sample data
n <- 3
matcov <- var(matrix(rnorm(2*n^2),2*n,n))

library(quadprog)
solve.QP(matcov, rep(0,n), matrix(1,nc=1,nr=n), 1, meq=1)

# With lower- and upper-bound constraints:
solve.QP(
  matcov, rep(0,n), 
  cbind( 
    rep(1,n),  # Equality constraint
    diag(n),   # Lower bound constraint
    -diag(n)   # Upper bound constraint
  ),
  c(
    1,         # Equality constraint RHS
    rep(.1,n), # Lower bound
    rep(-1,n)  # - Upper bound
  ),
  meq = 1      # The first constraint is an equality
)
于 2013-04-16T12:46:04.497 に答える