Re: [問題] R 迴圈 x[i,]+x[i-1,]的問題

看板Statistics作者 (........)時間17年前 (2008/05/09 22:59), 編輯推噓2(200)
留言2則, 1人參與, 最新討論串2/3 (看更多)
看不懂 t=1 時, x[t-1] 是多少... 我這邊先把 x[1] + x[0] 這項不列出來 故程式碼可以寫成下面這樣: x <- matrix(rnorm(100*3), ncol=3) t(apply(x, 1, function(x) x[-3] + x[-1])) 多善用 R 向量計算的功能, 可以簡化程式與增加執行效率... ※ 引述《bcs (= ="frailty..gggg XD)》之銘言: : 請問我想在對100*3的matrix做lag一期的加總 x[t]+x[t-1] : 但是當t=1時 x[0,] => col names。 : 若是用lag(x),dim仍是100*3。 : 下面的語法出現nothing to replace, 請問我該怎麼改語法才能 : 順利完成。 : 謝謝 : #j:col, t:row : x <- matrix(rnorm(100*3),100,3) : k <- matrix(0,100,3) : for(j in 1:3){ : for(t in 1:100){ : k[t,j] <- x[t,j] + x[t-1,j] : }} -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 211.74.96.39 ※ 編輯: KnightX 來自: 211.74.96.39 (05/09 23:01)

05/09 23:29, , 1F
謝啦,我要的是apply(x, 2, function(x) x[-3] + x[-1])
05/09 23:29, 1F

05/10 01:11, , 2F
apply(x,2,function(x) x[-100]+x[-1])
05/10 01:11, 2F
文章代碼(AID): #1896ROlY (Statistics)
文章代碼(AID): #1896ROlY (Statistics)