0

我希望提高下面代码的性能。我的想法是尝试通过以某种方式使 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
}
4

1 回答 1

0

我不确定如何在应得的情况下给予信用,但@DWin 是正确的......复制是解决方案。

于 2013-04-03T00:50:16.493 回答