[問題] 類似上三角矩陣

作者: 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