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