如何纠正两个滞后期？答案

【问题标题】：How to right two lags periods?如何纠正两个滞后期？
【发布时间】：2020-07-01 19:10:20
【问题描述】：

我有一个关于 R 的问题；我想在我的模型中使用两个滞后期而不是一个（请检查我下面的代码），但我不知道如何用 R 编写它。有人可以帮忙吗？

以下是我的 R 代码的详细信息：

library(plm)
fixed = plm(sp ~lag(debt)+lag(I(debt^2))+outgp+gvex+vlimp+vlexp+bcour+infcpi, data=pdata, index=c("country", "year"), model="within")

滞后必须在可变债务上。

【问题讨论】：

标签： r statistics

【解决方案1】：

这应该会给 debt 变量带来 2 个滞后。

library(plm)
    
fixed = plm(sp ~lag(debt, k=1:2)+lag(I(debt^2))+outgp+gvex+vlimp+vlexp+bcour+infcpi, data=pdata, index=c("country", "year"), model="within")

例如：

data("Grunfeld", package = "plm")
lags2mod <- plm(inv ~ lag(value, k=1:2) + capital, data = Grunfeld, model = "within")

summary(lags2mod)

Oneway (individual) effect Within Model

Call:
plm(formula = inv ~ lag(value, k = 1:2) + capital, data = Grunfeld, 
    model = "within")

Balanced Panel: n = 10, T = 18, N = 180

Residuals:
      Min.    1st Qu.     Median    3rd Qu.       Max. 
-272.21434  -19.24168    0.42825   18.09930  260.85548 

Coefficients:
                      Estimate Std. Error t-value  Pr(>|t|)    
lag(value, k = 1:2)1  0.078234   0.015438  5.0677 1.059e-06 ***
lag(value, k = 1:2)2 -0.018754   0.016078 -1.1664    0.2451    
capital               0.352658   0.021003 16.7910 < 2.2e-16 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Total Sum of Squares:    2034500
Residual Sum of Squares: 617850
R-Squared:      0.69631
Adj. R-Squared: 0.67449
F-statistic: 127.633 on 3 and 167 DF, p-value: < 2.22e-16

【讨论】：

@Pradanna S，非常感谢您的帮助，再次感谢。
@EmilioDas，当然。请接受我的回答并点赞！！