如何在lm模型的summary（）中放入循環？

Question

我正在嘗試分析不同分布的影響。 我創建了兩種類型的分布：mu =累積[epsilon N（0,1）] + eta N（0,1）mu =累積[epsilon N（0,1）] + eta N（0,4）

我已經創建了線性回歸模型，並希望輸出摘要，因為我對R ^ 2感興趣。 由於每次創建循環時數據都是不同的，因此摘要要進行多次，並且打印結果是所有摘要的平均值。

我根本不知道如何創建該循環。

eps = rnorm(200,0,1)
eta200_1 = rnorm(200,0,1)
eta200_4 = rnorm(200,0,4)

y0 = cumsum(eps) + eta200_1
y1 = cumsum(eps) + eta200_4

#The model output that shall be repeated
modely0_y1 = lm(y0 ~ y1)
summary(modely0_y1)

Answer 1

像這樣

iterations <- 100
r_squared <- c()

for(i in 1:iterations) {

  eps = rnorm(200,0,1)
  eta200_1 = rnorm(200,0,1)
  eta200_4 = rnorm(200,0,4)

  y0 = cumsum(eps) + eta200_1
  y1 = cumsum(eps) + eta200_4

  modely0_y1 = lm(y0 ~ y1)
  r_squared <- c(r_squared, summary(modely0_y1)$r.squared)


}

r_squared
# [1] 0.8989347 0.7321244 0.5185411 0.7102153 0.4646995 etc

mean(r_squared)
[1] 0.6025012