r - R中的随机游走分布直方图

Question

我想研究分布（n = 10000 的 S 值大于 0 的直方图重复了 10000 次。但是，我没有得到正确的输出，我该如何获取直方图？这就是我所拥有的：

rwlength=function(nsims,n){
t=numeric(nsims)
  for( i in (1:nsims))
{
  t[i]=aboveaxis(n)
}
return(t)
}
hist(rwlength(10000,10000))

笔记：

aboveaxis = function(n) {
if (n<=0){ return(cat("n must be greater than 0"))}
else
step = c(1, -1)
S = c(0, cumsum(sample(step,n, prob=c(.5, .5), replace=T)))
above=sum(S[S > 0])
return(above)
}

aboveaxis 是返回大于 0 的 S 值之和的函数

谢谢！