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

Question

我想研究分布（S值大于0的直方图，n = 10000重复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值之和

谢谢！