如何对多个因变量进行回归并比较这些模型？

Question

下面，我设法使用多个因变量运行回归。我如何将这些估计的结果分配给循环内的模型对象，然后基于 AIC/BIC 度量比较这些模型，例如使用

rcompanion::compareLM()

？

depvar <- c("qsec", "hp")
indepvar <- paste(c("mpg", "wt", "am", "disp"), collapse = " + ")

for(i in 1:length(depvar)){
  est_model <- lm(paste(depvar[i],"~", indepvar), data=mtcars)
  res <- summary(est_model)
  print(list(depvar[i], res))
}