推 bcs:謝啦，我要的是apply(x, 2, function(x) x[-3] + x[-1]) 05/09 23:29
推 bcs:apply(x,2,function(x) x[-100]+x[-1]) 05/10 01:11
看不懂 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)