我希望提高下面代码的性能。我的想法是尝试通过以某种方式使 st2 的每一列成为自己的矩阵来“矢量化”这个操作,这是数据框可以为我做的吗?
st<-matrix(nrow=257,ncol=10000)
st2<-matrix(nrow=257,ncol=10000)
st[1,]=1489
for(i in 1:10000)
{
st2[1,]<-st[(256+1),1]
for(j in 2:(256+1))
{
st2[j,]=st2[(j-1),] *exp(.5-(.2^2/2)) * rnorm(10000)
}
print(st2[256+1,]) #this for example, in reality I store these values for later
}