使用 lapply 和两个参数运行滞后回归

我正在运行多个单变量回归，就像在这个可重现的示例中一样：

require(dynlm)
data(USeconomic)
US<-USeconomic
vars<-colnames(US)[-2]
a<-lapply(colnames(US),function(x) dynlm(log(GNP)~get(x),data=US))

a包含 3 个单变量回归的列表。假设现在我想运行相同的回归，但有 3 个滞后：l<-c(0,1,4)其中 0 当然是我已经得到的情况。有没有办法使用向量l直接，比如

# this did not work for me, I obtain multivariate regressions including all lags at once
lapply(colnames(US),function(x) dynlm(log(GNP)~L(get(x),l),data=US),l=l)

在这不起作用之后，我尝试了另一种方法并添加到以下向量：

lagged_vars <- paste("L(",rep(vars,each=3),",",l,")",sep="")

to get:

[1] "L(log(M1),0)" "L(log(M1),1)" "L(log(M1),4)" "L(rs,0)"      "L(rs,1)"     
[6] "L(rs,4)"      "L(rl,0)"      "L(rl,1)"      "L(rl,4)"

不幸的是，我无法使用新的字符向量运行它， get() 没有帮助。我不明白为什么它适用于vars但不是lagged_vars它们都是字符向量。

请注意，L()语法来自 dynlm 包。 附带问题：如果我只打印 a，回归结果中的系数仍标记为 get(x) – 如何更改它？

i,j 循环可能是可能的解决方案，但我宁愿使用 lapply 或这个系列中的其他东西......

EDIT: as.formula不与一起工作L() from dynlm。我收到此错误消息：

merge.zoo(log(GNP), L(log(M1), 0), retclass = "list", all = FALSE) 中的错误： 找不到功能“L”

编辑： 找到了一个Achim Zeileis 的有趣帖子提到这个问题。

要构建 R 公式，您必须粘贴它all在一起，而不仅仅是它的预测端。所以你需要类似的东西：

formula <- as.formula(
    paste("log(GNP)~",
        paste("L(",rep(vars,each=3),",",l,")",sep=""),
        sep = ""
    )
)

然后运行

dynlm(formula, data = ...)