一种计算R中矩阵滚动标准差的快速方法

时间:2013-04-09 13:50:09

标签: r

我想找到一种快速计算矩阵滚动标准偏差的方法。为了准确地解释我的意思,我在下面嘲笑了一些代码以准确显示我在寻找什么。如你所知,这将是非常缓慢的

aaa <- matrix(1:1000,10,100)
bbb <- aaa
bbb[] <- NA
tt.n <- 5

for (i in (tt.n:dim(aaa)[2]))
    {
    bbb[,i] <- apply(aaa[,(i-tt.n+1):i],1,sd)   
    }

0 个答案:

没有答案