求助！这个模拟1000次的程序问题出在哪？

1767

收藏 2012-08-04

请大家帮我看看下面的这段程序，最后要返回的是满足N>1时的k的平均值，我做了1000次，所以应返回这1000个k的均值，每次符合条件的k都赋予n[i]，但为什么结果最后竟然是500多？？请问问题出在哪里？

kehou3=function(B)
{
  k=0
  N=0
  n=array()
  for(i in 1:B)
{

   repeat
      {
      k=k+1
      u=runif(1)
      N=N+u
      if(N>1)
         {
            n[i]=k
            break
         }
      }

}
  mean(n)
}
kehou3(1000)

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

全部回复

meishanjia1900

2012-8-5 12:11:01

每次for循环时你确定不用在repeat语句之前加上K=0,N=0么？

而且，就算照你的程序运行，n也确实含有1000个元素呀？

不可能出现只含有500多个元素的情况！

你这500多究竟是什么意思？——照上下文理解就是“只含500多个元素”的意思。

但是我估计你是想说“mean(n)的值，即n的均值是500多”，对么？

你认为这样很不合理，当然，问题就出在for循环里，你应该在repeat前加上k = 0及N = 0的语句！