引导lm的标准错误

Question

我是R的新手，我可能需要一些帮助。

我已经获得了标准误差和y~z的lm回归的p值。现在我想做一个bootstrap并比较结果。到目前为止，我所能找到的是以下用于引导R平方值的代码：

for(i in 1:B) {
   indices <- sample(1:N, N, replace = TRUE)
   bs.data <- data[indices,]
   bs.OLS <- lm(y~z, data = bs.data)
   stor.se[i] <- summary(bs.OLS)$r.squared
}    

我无法弄清楚如何返回z系数的标准误差而不是R平方。用$r.squared替换循环最后一行的$sd不会有效。 $sd(z)和类似的事情也不会起作用。