R的箱子模拟中的球

Question

我想模拟固定数量的球分布m的固定数目个二进制位n在R.到现在为止我一直在使用泊松逼近rpois() 对于n箱中的大量球来说，这是一个不错的近似值。

但是， rpois()仅允许您指示速率lambda ，即m/n 。 结果，正箱的数量通常小于球的数量。

有人会知道允许我随机将球分配到箱子中的功能或脚本吗？

最终，我试图通过自举来计算-log(empty bins/total bins)的置信区间。 这个问题可以说是“打破我的球”。

Answer 1

我想你想要多项分布。

这是一个快速的功能 - 我们在n个箱子中取m个球，然后给出x个结果，为每个x试验返回一个度量向量：

myfunc <- function(m,n,x){
  out <- rmultinom(x,m,rep(1,n))
  -log(colSums(out == 0)/n)
}

myfunc(10,40,10)
[1] 0.1923719 0.2548922 0.2231436 0.2548922 0.2876821 0.2876821 0.2231436 0.2231436 0.2231436 0.2548922

然后，您可以获得分位数/置信区间：

out = myfunc(10,40,1000)
quantile(out, c(0.05,0.95))
       5%       95% 
0.1923719 0.2876821