limsolve用法 或其他package

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

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矩陣的數值都是最優化最小的嗎?
可是矩陣中的數值有負的 可以怎樣改這個語法 或是其他的可用








                       
               

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：package Solve LIMS Pack Age package matrix