Как сделать квадратичное программирование в R?

Question

Я пытаюсь решить задачу квадратичного в R. У меня есть эти данные:

x1   x2  y 
0.3858 0.4687  1 
0.4871 0.6110  -1 
0.9218 0.4103  -1 
0.7382 0.8936  -1 
0.1763 0.0579  1 
0.4057 0.3529  1 
0.9355 0.8132  -1 
0.2146 0.0099  1 

И мне нужно, чтобы выяснить, Л1 к λ8. Я пробовал использовать подпрограмму: solve.QP(Dmat, dvec, Amat, bvec, meq=0, factorized=FALSE) из пакета quadprog. Я смущен, что ввести для dvec, bvec, Dmat в моем случае. Я попытался следующие, но это не дает правильный результат:

Dmat <- matrix(0,8,8) 
Amat<- matrix(c(0.3858,0.4871,0.9218,0.7382,0.1763,0.4057,0.9355,0.2146,0.4687,0.611,0.4103,0.8936,0.0579,0.3529,0.8132,0.0099),8,2) 
dvec<-c(1,-1,-1,-1,1,1,-1,1) 
solve.QP(Dmat,dvec,Amat) 

Результат должен быть λ1 = 65,5261, λ2 = 65,5261, λ3 к λ8 = 0. Любая помощь будет оценена.

Answer 1

добавить diag (Dmat) < - 1 в ваш код.