估计多元回归系数

Question

我有以下问题：

我正在使用此库“ mosaicData”中的数据“ Railtrail”。 我已经有了以下线性回归模型的系数： lm(volume ~ hightemp + cloudcover + weekday, data = RailTrail)，为总体计算。

现在，我需要估计该模型与样本的系数，并建立一个置信区间（95％）。 因此，我需要计算先前生成的数据样本的所有系数。我被要求使用“ for”循环，但我不知道如何计算LR模型。我还需要存储获得的系数。

我试图这样做

trial <- list()
set.seed(101)

for(i in 1:100){
  trial[[i]] <- RailTrail %>% 
    lm(volume ~ hightemp + cloudcover + weekday, data = RailTrail)
}

但出现以下错误：

xj [i]中的错误：下标类型为“语言”无效

谢谢

如果我的要求不清楚，请随时询问进一步的精度。

Francisco

Answer 1

您是指模型参数的置信区间吗？如果是这样，我希望此示例说明一种简洁的方法：

model <- lm(mpg ~ cyl + disp + gear, data = mtcars)
bind_cols(broom::tidy(model), broom::confint_tidy(model))