[問題] 類似上三角矩陣

看板R_Language作者 (航)時間10年前 (2014/05/15 18:44), 10年前編輯推噓3(304)
留言7則, 4人參與, 最新討論串1/1
文章分類提示: 程式諮詢(我想用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

05/15 19:30, , 1F
所以你的得到什麼?
05/15 19:30, 1F

05/15 19:32, , 2F
我想呼籲大家問問題可不可以仔細點,不然對自己對別人都
05/15 19:32, 2F

05/15 19:32, , 3F
沒什麼好處。
05/15 19:32, 3F

05/15 20:00, , 5F
其實你大概只差一個column index的rev而已
05/15 20:00, 5F
※ 編輯: warex14 (140.112.231.6), 05/15/2014 22:59:28

05/16 10:55, , 6F
n <- 4; A <- (lower.tri(diag(n))+diag(n))[,n:1]; A
05/16 10:55, 6F
謝謝W大跟k大 ※ 編輯: warex14 (140.114.136.190), 05/16/2014 14:15:52

05/16 22:23, , 7F
1*(outer(1:n,n:1,'-')>=0)
05/16 22:23, 7F
文章代碼(AID): #1JT9gZyn (R_Language)
文章代碼(AID): #1JT9gZyn (R_Language)