如何處理 R 中的泊松分布

Question

考慮 lambda=1 的泊松分布x1, x2, ... ~ pois(1) 。 我想寫一個 function 接收一個數字作為輸入（考慮它a），並給我們作為 output最小的 n （最小 n），這對於sum(xi)>=a, i=1:n是正確的。 我認為在這種情況下使用while循環可能很好（但我不確定這是最好的方法）。 也許可以使用其他循環（如for循環）來完成。 我不知道如何處理 R 中包含泊松分布的這種情況？

Answer 1

while循環很容易編寫代碼。 它在s < a時為sum 累加s中的rpois的值，計算n中的迭代次數。

minpois <- function(a){
  n <- 0L
  s <- 0L
  while(s < a) {
    n <- n + 1L
    s <- s + rpois(1L, lambda = 1)
  }
  n
}

set.seed(2020)
minpois(10)
#[1] 12