使用apply R的递归循环

Question

我有两个循环，我想使用一些Apply函数进行转换，以减少计算时间。 第一个似乎很容易做到，第二个问题就出现了，因为S在每次迭代中都用自己的值进行更新。

S0   = 100
a    = 0.00016
b    = 0.0126
sim    = 10000 
drifts = 1000
Si = rep(0,sim)

for(i in (1:sim))
{
 S =  S0
 for (j in (1:drifts))
  {
   z = rnorm(1, mean = 0, sd = 1)
   S = S * exp(a + b*z)
  }
  Si[i] =S 
}

有人可以帮忙吗？

Answer 1

calc_s <- function(S, i = 1) {
  S <- S * exp(a + b * rnorm(1, mean = 0, sd = 1))
  return (if (i < drifts) calc_s(S, i + 1) else S)}
S2 <- sapply(1:sim, function(x) {
  calc_s(S0)
})

虽然不快

Answer 2

下面的代码怎么样？ 我已经用产品替换了j循环，我认为数学是正确的。

sapply(1:sim,function(x) S0*prod(exp(a + b*rnorm(drifts))))

它也明显更快：

> system.time(for(i in (1:sim))
+ {
+   S =  S0
+   for (j in (1:drifts))
+   {
+     z = rnorm(1, mean = 0, sd = 1)
+     S = S * exp(a + b*z)
+   }
+   Si[i] =S 
+ }
+ )
   user  system elapsed 
  23.29    0.02   23.34 
> 
> system.time(Si<-sapply(1:sim,function(x) S0*prod(exp(a + b*rnorm(drifts)))))
   user  system elapsed 
   1.76    0.00    1.76