[问题] 类似上三角矩阵

楼主: warex14 (航)   2014-05-15 18:44:49
文章分类提示:
程式咨询(我想用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这样的矩阵
目前还想不出什么好方法。
作者: andrew43 (讨厌有好心推文后删文者)   2014-05-15 19:30:00
所以你的得到什么?我想呼吁大家问问题可不可以仔细点,不然对自己对别人都没什么好处。
作者: Wush978 (拒看低质媒体)   2014-05-15 19:59:00
https://gist.github.com/wush978/1b67a9fcab232957b1e0其实你大概只差一个column index的rev而已
作者: koai (可笑的自以为...)   2014-05-16 10:55:00
n <- 4; A <- (lower.tri(diag(n))+diag(n))[,n:1]; A
作者: showfeb   2014-05-16 22:23:00
1*(outer(1:n,n:1,'-')>=0)

Links booklink

Contact Us: admin [ a t ] ucptt.com