1

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

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)   
    }
4

0 回答 0