R: censored cumsum (censored random walk)

【问题标题】：R: censored cumsum (censored random walk)R: censored cumsum (censored random walk)
【发布时间】：2016-02-28 06:19:45
【问题描述】：

我希望能够审查随机游走数据的生成，以便游走永远不会低于目标值（通常为 0）。下面的代码完成了我想要的，除了我宁愿有一个类似于 cumsum 的函数，我可以用它来快速处理数百万行这样的值，例如cumsum(x,min=0)：

x <- rnorm(1000)
y <- rep(0,length(x))
for(i in 2:length(x)) y[i] <- max(x[i]+y[i-1], 0)
plot(y, type='l')

【问题讨论】：

见this relevant QA

标签： r cumsum

【解决方案1】：

为什么不制作自己的函数（从您的代码中复制，但添加了 y[1] 的初始化以使其类似于 cumsum 行为）：

cumsum0<-function(x,min=0){
  y<-rep(0,length(x))
  y[1]<-max(x[1],0)
  for (i in 2:length(x)) y[i] <- max(x[i]+y[i-1], min)
  return(y)
}

x<-rnorm(1000)
plot(cumsum0(x),type="l")

【讨论】：