[問題] 類似上三角矩陣

[warex14 (航)]

文章分類提示:
程式諮詢(我想用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://www.ptt.cc/bbs/R_Language/M.1400150691.A.F31.html
→ andrew43:所以你的得到什麼？
→ andrew43:我想呼籲大家問問題可不可以仔細點，不然對自己對別人都
→ andrew43:沒什麼好處。
推 Wush978:https://gist.github.com/wush978/1b67a9fcab232957b1e0
推 Wush978:其實你大概只差一個column index的rev而已
※ 編輯: warex14 (140.112.231.6), 05/15/2014 22:59:28
推 koai:n <- 4; A <- (lower.tri(diag(n))+diag(n))[,n:1]; A
謝謝W大跟k大
※ 編輯: warex14 (140.114.136.190), 05/16/2014 14:15:52
→ showfeb:1*(outer(1:n,n:1,'-')>=0)