如何在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。 我嘗試使用以下例程：來自包quadprog solve.QP(Dmat, dvec, Amat, bvec, meq=0, factorized=FALSE) 。 在我的情況下，我對於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添加到您的代碼中。