如何使用 R 创建简单的分支流程？

Question

我正在尝试使用 R 进行简单的分支过程来计算 10 代以上的狐狸数量

我正在使用 lambda = 1.5，并且我知道可以使用rpois(1, 1.5)计算初始生成中的狐狸数量。

但是，我坚持如何创建一个 function 来计算下一代的人口，并将总人口加在一起。 如果有人能给出一些起始点，那就太好了。

Answer 1

我相信您正在寻找如下 function ：

pois_bprocess <- function(generations, lambda, trace = TRUE){
  z_i <- rpois(1, lambda)
  i <- 1
  
  if(trace) print(paste(i, ':', z_i))
  
  while(i < generations){
    # draw z_i samples from a poisson distribution with lambda and sum them
    z_i <- sum(rpois(z_i, lambda)) 
    i <- i + 1
    if(trace) print(paste(i, ':', z_i))
  }
  
  return(z_i)
}


# Example usage : pois_bprocess(10, 1.5)
# Example usage : pois_bprocess(10, 1.5, trace = F)

请注意，由于假设的 lambda 为 1.5 的泊松分布，导致早期绘制零并卡住的机会相当高。

HTH