如何在 R 中处理泊松分布

Question

考虑泊松分布x1, x2, ... ~ pois(1)，其中 lambda=1。 我想编写一个函数，它接收一个数字作为输入（将其视为 a）并将我们作为输出 最小的 n（最小 n），这对于 sum(xi)>=a, i=1:n 是正确的. 我认为在这种情况下使用while 循环可能很好（但我不确定这是不是最好的方法）。也许可以使用其他循环来完成，例如for 循环。我不知道如何处理在 R 中包含泊松分布的这种情况？

Answer 1

while 循环足够简单，可以编写代码。它在 s 中为 sum 而在 s < a 中累积 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