如何在R中使用boot包保存并复用回归模型的自助法系数？

当然有办法啦！你用boot包生成的results对象里，所有自助法迭代得到的系数都存在$t这个矩阵里，咱们一步步操作就能拿到wt的1000个系数：

第一步：提取wt的自助法系数 #

你的回归公式是mpg~wt+disp，对应系数顺序是截距（t1*）、wt（t2*）、disp（t3*）。results$t矩阵的第二列就是wt的1000个自助法系数，提取代码如下：

# 提取所有wt的自助法系数
wt_boot_coefs <- results$t[, 2]

第二步：查看与复用这些系数 #

提取完成后，你就可以自由使用这个系数向量了：

• 查看前10个系数：

head(wt_boot_coefs, 10)

• 分析系数的分布特征：

# 计算均值
mean(wt_boot_coefs)
# 计算标准差
sd(wt_boot_coefs)
# 绘制直方图观察分布
hist(wt_boot_coefs, main="Bootstrap Coefficients for wt", xlab="Coefficient Value")

• 计算95%百分位数置信区间：

quantile(wt_boot_coefs, c(0.025, 0.975))

小提示 #

如果不确定列对应的变量，可以用colnames(results$t)来确认：

colnames(results$t)
# 输出会显示："(Intercept)" "wt" "disp"

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