R语言中runif()函数结合for循环与手动设置观测数时回归结果不一致问题

核心问题：随机数流的消耗差异 #

你遇到的问题确实和set.seed()的使用以及循环中随机数的生成顺序有关，具体原因拆解如下：

1. 单次运行的随机数来源：当你单独设置i=1000并运行时，set.seed(69420)初始化随机数生成器后，两次runif(1000)会直接使用种子后的前2000个随机数（每次调用生成1000个）。

2. 循环中的随机数消耗：你的循环代码里set.seed(69420)只执行了一次，之后每次循环都会调用两次runif(i)生成J和M。由于循环是从100逐步递增到1000，在第10次生成1000个样本之前，已经消耗了大量随机数：

   • 前9次循环总共生成的随机数数量是：2*(100+200+...+900) = 9000个
   • 第10次循环才会取随机数流中接下来的2000个值
     这就导致循环中第10次的J、M和单独运行时的完全不是同一批随机数，回归系数自然不一样。

另外你的循环代码还有个逻辑小问题：内层for(i in x)里的x并未定义，这会让循环行为偏离预期，进一步加剧了随机数的消耗差异。

修正方案 #

方案1：让循环第10次与单独运行结果一致 #

如果你希望两者结果匹配，可以通过跳过前面消耗的随机数来对齐随机数流的位置：

修正后的循环代码

# 初始化结果矩阵
Vr <- matrix(nrow = 10, ncol = 2)
sample_sizes <- seq(100, 1000, by = 100)

set.seed(69420)
for (k in 1:length(sample_sizes)) {
  n <- sample_sizes[k]
  J <- runif(n, min = 0, max = 2)
  M <- runif(n, min = 0, max = 2)
  L <- pmax(J, M)
  W <- pmin(J, M)
  
  # 计算回归系数
  vW <- mean((W - mean(W))^2)
  cWL <- mean((L - mean(L)) * (W - mean(W)))
  a <- cWL / vW
  b <- mean(L) - a * mean(W)
  
  Vr[k, ] <- c(a, b)
}

print(Vr)

匹配结果的单独运行代码

set.seed(69420)
# 跳过前9次循环消耗的9000个随机数
dummy <- runif(9000)
# 生成n=1000的样本
n <- 1000
J <- runif(n, min = 0, max = 2)
M <- runif(n, min = 0, max = 2)
L <- pmax(J, M)
W <- pmin(J, M)

vW <- mean((W - mean(W))^2)
cWL <- mean((L - mean(L)) * (W - mean(W)))
a <- cWL / vW
b <- mean(L) - a * mean(W)

cat("a =", a, "b =", b)

此时得到的(a,b)会和循环中第10次的结果完全一致。

方案2：独立种子的模拟方式（更推荐） #

如果你的目标是做独立模拟（每次模拟不受前一次影响），可以为每个样本量设置独立种子：

Vr <- matrix(nrow = 10, ncol = 2)
sample_sizes <- seq(100, 1000, by = 100)

for (k in 1:length(sample_sizes)) {
  n <- sample_sizes[k]
  set.seed(69420 + k)  # 每个模拟用不同的种子，保证独立性
  J <- runif(n, min = 0, max = 2)
  M <- runif(n, min = 0, max = 2)
  L <- pmax(J, M)
  W <- pmin(J, M)
  
  vW <- mean((W - mean(W))^2)
  cWL <- mean((L - mean(L)) * (W - mean(W)))
  a <- cWL / vW
  b <- mean(L) - a * mean(W)
  
  Vr[k, ] <- c(a, b)
}

print(Vr)

总结 #

• R的随机数生成器是状态机，每次调用runif()都会推进状态，后续随机数基于当前状态生成。
• 循环和单独运行的差异，本质是使用了随机数流中不同位置的数值。
• 通过对齐随机数流位置或设置独立种子，就能解决结果不一致的问题。

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