r - 如何运行滚动窗口面板回归？

Question

我有一个公司年面板数据，其中包含 10000 家公司的 30 年数据：

 firm       year        y        x1         x2
 A          1997       150      130        202
 A          1998       200      250        230
 A          1999       350       22         50
 A          2000        30      150         30
 B          1995        50      350         30
 B          1996       350      270         59
 B          1997       789      360        382
 B          1998       360      385         75
 B          1999       362      567        123

我想为面板回归运行宽度 = 6 的滚动窗口：

plm(y~ lag(x1, 1)+ lag(log(x2), 1), data= dt, 
    model= "within", effect= "twoways",
    index= c("firm", "year"))

我认为对于 10000 家公司，滚动窗口应该从第 1 年到第 30 年 10000 次。我想获得每个公司年度观察的滞后（x1, 1）的每个系数。