Question

我正在尝试复制我在文章中看到的回归：

Y_t=a(L) Y_(t-1)

L是滞后算子中的多项式，它表示

“滞后因变量[滞后：1-24]输入为四分之一移动平均分别滞后1、5、9、13、17和21”。

I saw，我可以执行以下操作：

lm(Y ~ I(lag(Y)) + I(lag(Y, 5)) + I(lag(Y, 9)) + I(lag(Y, 13)) + I(lag(Y, 17)) + I(lag(Y, 21)) + I(lag(Y, 24))

有更短的方法吗？

谢谢。

Answer 1

也许是这样的

lags <- c(1, 5, 9, 13, 17, 21)
xnam <- paste0("I(lag(Y, ", lags, "))")
(fmla <- as.formula(paste("y ~ ", paste(xnam, collapse= "+"))))