如何使循环更快

时间:2019-04-17 07:51:15

标签: r loops

此循环大约需要一个小时。有什么办法可以使速度更快?

l <- 2

while (l < 300001) {
  VV[1,l] = log(b+(l-1) * 0.001-0.0005) + k + 
    beta * (1-delta) * (1-p * sum(Pr[1,(1:l)])) * V[1, l] + 
    beta  *(1-delta) * p * sum((V[1,] * Pr[1,])[(l+1):300001])
  l = l + 1
}

1 个答案:

答案 0 :(得分:0)

如果在* nix(Linux,Unix,macOS等)上,则可以使用parallel处理:

library(parallel)

VV <- matrix(ncol = 300001)
mclapply(2:300001, function(l){
  VV[1,l] <<- log(b+(l-1) * 0.001-0.0005) + k + beta * (1-delta) * (1-p * sum(Pr[1,(1:l)])) * V[1,l]+beta  *(1-delta) * p * sum((V[1,] * Pr[1,])[(l+1):300001])
  }, mc.cores = 5)