循环并将结果收集到表中

Question

我试图让 RStudio 执行以下操作：

1. 运行代码 20 次“回归”
2. 收集每次回归“Alpha、Beta 和 Pr 值”的结果，并以表格格式提供给我。

以下是我正在使用的回归代码：

xaxis = rnorm(500,0,1)
eval = rnorm(500,0,1)
beta = 0.5
intercept = 1.5
yaxis = intercept + beta*xaxis + eval
Regression_yxe = lm(yaxis ~ xaxis)
summary(Regression_yxe)

谢谢，

Answer 1

首先将你的代码变成一个函数：

regress <- function () {
    xaxis <- rnorm(500,0,1)
    eval <- rnorm(500,0,1)
    beta <- 0.5
    intercept <- 1.5
    yaxis <- intercept + beta*xaxis + eval
    Regression_yxe <- lm(yaxis ~ xaxis)
    summary(Regression_yxe)$coefficients[, c(1, 4)]
}

然后运行该函数 20 次：

results <- replicate(20, regress(), simplify=FALSE)

然后合并结果：

do.call(rbind, results)
#              Estimate      Pr(>|t|)
# (Intercept) 1.4764338 2.957601e-124
# xaxis       0.4839075  7.532197e-27
# (Intercept) 1.4521559 5.397893e-117
# xaxis       0.4727487  3.001808e-20
#    . . . .
# (Intercept) 1.5099205 1.406824e-131
# xaxis       0.4963504  1.923430e-24
# (Intercept) 1.4693233 2.519495e-123
# xaxis       0.5008844  5.461081e-25