如何在R中重复执行1000次OLS过程并逐步保存结果？

问题背景 #

你需要重复执行1000次自助抽样+OLS参数估计的流程并保存结果，已有基础代码但不清楚如何将完整逻辑嵌入重复执行的函数中，纠结循环或replicate的选择。

你的原始代码：

x1 = runif(4000, min = 0, max = 1)
x2 = runif(4000, min = 0, max = 1)
y <- 1*x1 - 2*x2 + rnorm(4000)
df <- data.frame(y, x1, x2)

# 需要重复1000次的部分:
set.seed(2)
df2 <- df[sample(1 : nrow(df), 4000, replace = T ), ]
x= as.matrix(df2[,-1])
y= as.matrix(df2[,1])
OLS <- solve(t(x)%*%x)%*%t(x)%*%y

# 尝试的代码
set.seed(2)
n = 1000
out <- replicate(n, {df2 <- df[sample(1 : nrow(df), 4000, replace = T ), ]})

最优实现方式：使用replicate函数 #

replicate是R中专门用于重复执行表达式并收集结果的函数，底层基于向量化操作，比手动循环更简洁高效，是这类场景的首选。只需把完整的抽样、矩阵转换、OLS计算逻辑放进replicate的表达式块中，让它返回每次的OLS系数即可。

完整代码：

# 1. 构建原始数据集
x1 = runif(4000, min = 0, max = 1)
x2 = runif(4000, min = 0, max = 1)
y <- 1*x1 - 2*x2 + rnorm(4000)
df <- data.frame(y, x1, x2)

# 2. 执行1000次重复自助抽样与OLS估计
set.seed(2) # 设置随机种子，保证结果可复现
n_reps <- 1000

ols_results <- replicate(n_reps, {
  # 自助抽样
  df2 <- df[sample(nrow(df), size = 4000, replace = TRUE), ]
  # 转换为矩阵
  x_mat <- as.matrix(df2[, -1])
  y_mat <- as.matrix(df2[, 1])
  # 计算OLS系数
  solve(t(x_mat) %*% x_mat) %*% t(x_mat) %*% y_mat
})

# 3. 将结果整理为易读的数据框（可选）
ols_df <- as.data.frame(t(ols_results))
colnames(ols_df) <- c("x1_coef", "x2_coef")

说明：

• replicate会自动将每次返回的2维OLS系数向量收集为一个矩阵，每一列对应一次重复的结果；转置后得到每行对应一次重复的数据框，方便后续统计分析（比如计算系数的均值、标准差）。
• 如果担心手动矩阵计算的数值稳定性，可以改用lm函数简化代码，结果一致且内部有更稳定的数值处理：

  ols_results <- replicate(n_reps, {
    df2 <- df[sample(nrow(df), 4000, replace = TRUE), ]
    coef(lm(y ~ x1 + x2 - 1, data = df2)) # -1表示不拟合截距项
  })
  

备选：预先分配容器的for循环 #

如果需要更精细的流程控制，可以使用for循环，但必须预先分配结果容器（如矩阵），避免动态扩容带来的性能损耗，此时性能与replicate相当。

代码示例：

set.seed(2)
n_reps <- 1000
# 预先分配存储结果的矩阵
ols_coefs <- matrix(nrow = n_reps, ncol = 2)

for(i in 1:n_reps){
  df2 <- df[sample(nrow(df), 4000, replace = TRUE), ]
  x_mat <- as.matrix(df2[, -1])
  y_mat <- as.matrix(df2[, 1])
  ols_coefs[i, ] <- solve(t(x_mat) %*% x_mat) %*% t(x_mat) %*% y_mat
}

ols_df <- as.data.frame(ols_coefs)
colnames(ols_df) <- c("x1_coef", "x2_coef")

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