R泊松模拟

Question

我需要编写一个函数来计算在 n 次试验中直到时间 t 的到达次数。参数应该是 lambda（我可以假设它在 0.1 和 1 之间）、时间（我可以假设小于或等于 1）和要采样的计数。

我之前写过一个函数，它采用长度为 n 的向量，其中第一个 n-1 个元素作为事件间时间，第 n 个元素作为时间t，统计t之前发生的事件数。

inp <- readline(prompt="Input vector with each element seperated by a space")
inp <- strsplit(inp," ")
inp <- as.integer(as.vector(inp[[1]]))
t <- tail(inp, n=1)
c.e <- function(x) {
  inp = x
  stopped = NA
  for (i in seq_along(inp)) {
    runsum <- sum(inp[1:i])
    cat("The sum of the", i, "first elements is", runsum, "\n")

    if (runsum > tail(inp, 1)) {
      stopped = i - 1
      break()
    }
  }

  stopped
}
cat(c.e(inp), "events occur in", t, "time units")

（例如输入 1 2 3 4 7 将输出 3 个事件在 6 个时间单位内发生）

我认为我需要使用并可能编辑此功能，以便让它完成我需要它做的事情，但我真的不知道该怎么做。任何帮助将不胜感激:)

Answer 1

如果您想要诸如打印之类的东西以及您所拥有的良好输入/输出，您可以将其编辑到您的旧功能中，但至于您实际要求的，似乎您所需要的只是counts<-function(lambda, t, n) sum(rpois(n, lambda) < t).

rpois 生成试验结果，参数n（试验次数）和lambda，然后我们将它们与t（时间）和sum我们的结果进行比较。