作者warex14 (航)
看板R_Language
标题[问题] 类似上三角矩阵
时间Thu May 15 18:44:49 2014
文章分类提示:
程式谘询(我想用R 做某件事情，但是我不知道要怎麽用R 写出来)
[软体熟悉度]:
入门(写过其他程式，只是对语法不熟悉)
[问题叙述]:
对不起，当初问得太快了，以下补得详细一点。
目前我有一个大概这得像这样的矩阵(n x n)
a a a a a a ...a (n个)
b b b b b b ...b
c c c c c c ...c
...
...
...
n n n n n n ...n
我想要把它变成
0 0 0 0 0 0 ... a
0 0 0 0 0 0...b b
0 0 0 0 ... c c c
...
...
...
n n n n ... n n n
所以
我需要产生一个类似长得像这样的矩阵
A=
0 0 0 0 ... 1
0 0 0 0...1 1
0 0... 1 1 1
...
...
...
1 1 1 1...1 1
我想要利用A来跟原本的矩阵相乘
这样应该就可以得到我想要的矩阵
本来是想利用lower.tri()或是upper.tri()函式，但好像没办法得到像A这样的矩阵
目前还想不出什麽好方法。
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 140.112.25.121
※ 文章网址: http://webptt.com/cn.aspx?n=bbs/R_Language/M.1400150691.A.F31.html
1F：→ andrew43:所以你的得到什麽？ 05/15 19:30
2F：→ andrew43:我想呼吁大家问问题可不可以仔细点，不然对自己对别人都 05/15 19:32
3F：→ andrew43:没什麽好处。 05/15 19:32
5F：推 Wush978:其实你大概只差一个column index的rev而已 05/15 20:00
※ 编辑: warex14 (140.112.231.6), 05/15/2014 22:59:28
6F：推 koai:n <- 4; A <- (lower.tri(diag(n))+diag(n))[,n:1]; A 05/16 10:55
谢谢W大跟k大
※ 编辑: warex14 (140.114.136.190), 05/16/2014 14:15:52
7F：→ showfeb:1*(outer(1:n,n:1,'-')>=0) 05/16 22:23