Andreas Jensen bentrevisor at gmail.com writes:
Hello.
I'm trying to solve a quadratic programming problem of the form min
||Hx - y||^2 s.t. x = 0 and x = t using solve.QP in the quadprog
package but I'm having problems with Dmat not being positive definite,
which is kinda okay since I expect it to be numerically semi-definite
in most cases. As far as I'm aware the problem arises because the
Goldfarb and Idnani method first solves the unconstrained problem
requiring a positive definite matrix. Are there any (fast) packages
that allows me to do QP with (large) semidefinite matrices?
Example:
t - 100
y - signalConvNoisy[1,]
D - crossprod(H,H)
d - crossprod(H,y)
A - cbind(rep(-1, nrow(H)), diag(ncol(H)))
b0 - c(t, rep(0, ncol(H)))
sol - solve.QP(Dmat=D, dvec = d, Amat = A, bvec = b0)$solution
Error in solve.QP(Dmat = D, dvec = d, Amat = A, bvec = b0) :
matrix D in quadratic function is not positive definite!
Thanks in advance,
Andreas Jensen
See the Optimization task view, there are several packages that promise to
handle quadratic programming problems. You may also try one of the nonlinear
optimization packages, they are not much slower nowadays.
Your code cannot be executed as signalConvNoisy is unknown.
Hans Werner
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.