我想运行一系列回归模型，并将与特定协变量相关的输出放入数组。我已经能够使用plyr包中的laply来实现这一点，如下所示：

set.seed(1)
df <- data.frame(x1 = rnorm(100, 0, 1), 
                 x2 = rnorm(100, 0, 1),
                 y1 = rnorm(100, 0, 1), 
                 y2 = rnorm(100, 0, 1), 
                 t2 = sample(0:1, 100, replace = TRUE), 
                 t1 = sample(0:1, 100, replace = TRUE))

run_regressions <- function(model) {
  model <- glm(model, data = df)
  return(summary(model)$coefficients[2,])
}

models <- list("y1 ~ t1 + x1",
               "y2 ~ t2 + x1 + x2")

results <- laply(models, run_regressions)
results

字符串
现在我想扩展它，以便我可以将回归模型的不同元素分别传递给run_regressions函数。回归函数的定义如下

run_regressions2 <- function(outcome, treatment, covariates) {
  model <- glm(paste(outcome, " ~ ", treatment, " + ", covariates), data = df)
  return(summary(model)$coefficients[2,])
}

型
其中要运行的特定模型被指定为

models <- list(c("y1", "t1", "x1"),
               c("y2", "t2", "x1 + x2"))

型
有没有一些方法可以使用laply或其他函数来实现这一点？
我试着在plyr文档中搜索类似的例子，并在SO中搜索类似的问题，但没有找到任何结果。