R语言rpois函数参数为数组时的行为差异咨询

核心结论 #

直接调用rpois(R, lambda_vec)（lambda为向量），本质是为每个lambda元素生成1个泊松随机数，最终输出长度与lambda向量一致的结果；而你分别用每个标量lambda调用rpois(R, lambda_i)，是为单个lambda生成R个泊松随机数。两者结果的关联完全由R的随机数生成器的状态顺序决定。

详细解释 #

1. rpois()的参数规则
   rpois(n, lambda)的隐含逻辑：

   • 当lambda是向量时，如果n等于lambda的长度，函数会为每个lambda元素生成恰好1个泊松数；如果n不等于lambda长度，会循环lambda来匹配n的长度。
   • 你案例中R=5，lambda_vec长度也是5，所以直接调用会生成5个泊松数，每个数对应lambda_vec中位置相同的lambda值。

2. 随机数生成的顺序性
   R的随机数生成器是状态机，每次调用随机函数都会消耗若干随机数（数量取决于具体算法，比如泊松生成的Knuth算法每个数可能需要多个均匀随机数）：

   • 当你分5次调用rpois(5, lambda_i)，每次调用都会一次性消耗生成5个泊松数所需的随机数，得到5个结果（对应你给出的5行数据）。
   • 当你直接调用rpois(5, lambda_vec)，函数会依次为每个lambda元素执行「生成1个泊松数」的逻辑，消耗的随机数序列和你**重新设置种子后依次调用5次rpois(1, lambda_i)**完全一致，因此结果相同。
   • 你分调用时每行的第一个结果，和单独调用rpois(1, lambda_i)的结果不同，是因为批量生成多个泊松数时，R可能使用了优化的批量算法，导致随机数的消耗路径和单个生成不同。

验证示例 #

重新设置种子后，依次调用单个生成：

set.seed(123)
R <- 5
X <- rexp(R, 2)
lambda_vec <- exp(X/4)
# 依次生成单个泊松数
sapply(lambda_vec, function(l) rpois(1, l))
# 输出结果：1 2 0 1 2，和直接调用rpois(R, lambda_vec)完全一致

内容的提问来源于stack exchange，提问作者Teodoro Marinucci