A <- matrix(nrow = 4, ncol = 3,
+ data = c(3, 1, 2, 0, 2, 0, 0, 1, 1, 0, 2, 0))
B <- matrix(c(2, 1, 8, 3),4,1)

請問 使用Ｒ
1.上面例子Ａ是4X3矩陣  Ｂ是4X1矩陣
Ａx=B
求解 x（3X1矩陣 ）
而x內的三個元素數值都要大於0  且相加等於1 有什麼樣的package 和語法可以使用

如果看不懂問題
可以看以下網址中 2.3Replicating the mortality profile   裡面的2.8公式  我要計算死亡率的權重 此權重相加等於1
http://www.sciencedirect.com/science/article/pii/S0167668713001674?np=y

似乎可以用 或是
limsolve包 、Rdonlp2
                                                                                                                                                [size=8.000000pt]

[size=8.000000pt]

有查到limsolve包
                               
                       
               
在這個網址中的3.1. Equalities only 有說明到 可以用lsei 優化求解 min|Ａx-Ｂ|平方 算出x
https://cran.r-project.org/web/packages/limSolve/vignettes/limSolve.pdf
舉例
A <- matrix(nrow = 4, ncol = 3,
+ data = c(3, 1, 2, 0, 2, 0, 0, 1, 1, 0, 2, 0))
B <- c(2, 1, 8, 3)
lsei(A = A, B = B, fulloutput = TRUE, verbose = FALSE)

2.想問這個x矩陣的數值都是最優化最小的嗎?
可是矩陣中的數值有負的 可以怎樣改這個語法 或是其他的可用