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