0

我想研究分布(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 值之和的函数

谢谢!

4

0 回答 0