R 3.5.3与3.6版本set.seed()向后兼容性及跨版本复现问题

我之前也遇到过类似的版本回退导致随机结果不匹配的问题，核心原因是R 3.6.0对sample()的默认采样算法做了关键变更：从旧版的"Rounding"（舍入法）改成了"Rejection"（拒绝采样法），而R 3.5.3只支持"Rounding"这一种采样方式。要在3.5.3里复现3.6生成的样本，关键是模拟出3.6版本的sample()行为，具体可以用这两种方法：

方法1：手动实现R 3.6的拒绝采样逻辑 #

如果你的样本生成主要依赖sample()或sample.int()，可以在R 3.5.3里自己写一个函数模拟3.6的拒绝采样逻辑。下面是一个简化版的实现，能匹配3.6中无放回、等概率的采样结果：

# 模拟R 3.6的sample.int（无放回、等概率）
sample_int_rejection <- function(n, size) {
  if (size < 0 || size > n) stop("size must be between 0 and n")
  if (size == 0) return(integer(0))
  
  # 核心逻辑参考R 3.6的官方源码
  result <- integer(size)
  used <- logical(n)
  count <- 0L
  
  while (count < size) {
    # 生成候选随机数
    candidate <- ceiling(runif(1) * n)
    if (!used[candidate]) {
      count <- count + 1L
      result[count] <- candidate
      used[candidate] <- TRUE
    }
  }
  result
}

# 包装成和sample()用法一致的函数
sample_rejection <- function(x, size, replace = FALSE, prob = NULL) {
  if (!replace && !is.null(prob)) {
    stop("带权重的无放回拒绝采样需要更复杂的实现，如需可参考R官方源码")
  }
  if (replace) {
    # 有放回采样在两个版本中行为一致，直接用原生函数即可
    return(sample(x, size, replace = TRUE, prob = prob))
  }
  if (is.null(prob)) {
    x[sample_int_rejection(length(x), size)]
  } else {
    stop("带权重的无放回拒绝采样暂未实现，建议参考R 3.6的src/main/random.c源码扩展")
  }
}

使用的时候，把原来调用sample()的地方换成sample_rejection()，同时必须保证种子和R 3.6中用的完全一样：

# 设置和R 3.6生成数据时相同的种子
set.seed(123)
# 比如在R 3.6中你用了sample(1:10, 5)，在3.5.3里就这么调用：
sample_rejection(1:10, 5)

方法2：移植R 3.6的采样函数代码 #

如果需要完整支持带权重的采样等复杂场景，最好直接从R 3.6的源码里提取sample.int()的拒绝采样逻辑，做成一个小的C扩展包在3.5.3里用。步骤如下：

• 下载R 3.6的源码，找到src/main/random.c里关于sample_int的实现
• 提取其中拒绝采样的核心代码，编写一个独立的C函数
• 在R 3.5.3中通过Rcpp或者原生的R扩展机制编译加载这个函数，用它替代原生的sample.int()

几个关键注意点 #

• 种子必须完全一致：不管用哪种方法，R 3.5.3里设置的种子必须和R 3.6生成数据时的种子一模一样，否则结果肯定对不上
• 其他随机函数（比如rnorm()、runif()）在两个版本里的底层生成器是一样的，只要种子相同，这些函数的结果会自动匹配
• 如果涉及带权重的无放回采样，3.5.3原生sample()和3.6的差异很大，一定要参考R官方源码实现对应的逻辑，才能保证结果一致

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