R函数生成新的协方差矩阵

Question

我正在运行多个线性回归模型，我想遍历结果以随后生成可靠的标准误差。我的代码当前看起来像这样，但是我想运行多个模型，而不必复制用于为每个模型计算鲁棒标准误差的代码。

# load data
data(mtcars)

# run models
m1 <- lm("mpg ~ wt", data = mtcars)
m2 <- lm("mpg ~ wt + hp", data = mtcars)

# calculate robust standard errors
cov1 <- vcovHC(m1, type = "HC3")
robust_se1 <- sqrt(diag(cov1))
cov2 <- vcovHC(m2, type = "HC3")
robust_se2 <- sqrt(diag(cov2))

我该如何编写函数来处理此任务。我计划使用连续的整数（例如m1，m2，m3）为每个模型编号。到目前为止，我还无法适应相关的SO答案，以便使用像this one这样的循环来生成变量。

编辑：已更改为可执行代码。