矩阵与每行中的另一个矩阵

时间:2013-04-02 21:50:23

标签: r

我希望改善下面代码的性能。我的想法是试图通过某种方式使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
}

1 个答案:

答案 0 :(得分:0)

我不确定如何在信用到期时给予信用,但@DWin是对的钱...复制是解决方案。