r - "Non conformable arguments"pgmm 错误(plm 库)

Question

我尝试使用 plm 包中的 pgmm 进行 Arellano 和 Bond (1991) 估计，但没有成功。为了查看问题是否出在我的数据中，我改为使用 plm 库提供的数据，但在使用“summary”命令时遇到了相同的问题:

Error in t(y) %*% x : non-conformable arguments

但是可以获得模型的系数。

我自己的数据有T=3，N=290。据我了解，T=3 是最小值，但应该足够了。当使用 Arellano 和 Bond 数据时，当 T=4 时我得到同样的错误。

data("EmplUK", package = "plm")
library(sqldf)
UK<-sqldf("select * from EmplUK where year in ('1982','1981',
'1980','1979')")


z1 <- pgmm(log(emp) ~ lag(log(emp), 1) + log(wage) +
log(capital) + log(output) | lag(log(emp), 2),
       data = UK, effect = "twoways", model = "twosteps")
summary(z1)

我理解估计方法和 R 公式的方式是，左边项是因变量的差值，右边第一个项是滞后差值。后一项由 (t-2) 中因变量的水平来衡量

我已经验证我使用的子集是 T=4 的平衡面板。当我考虑更多的年份时，一切都会顺利。所以一定是面板的长度造成了麻烦。

任何帮助将不胜感激。

Answer 1

提出了类似的问题here 。建议该错误与 mtest 有关，这是由 pgmm 摘要方法执行的串行相关测试。单独运行该函数似乎证实了这一点

>mtest(z1, order = 2)
Error in t(y) %*% x : non-conformable arguments

T=3 足以估计模型，但这只能给您留下最后一个时期的估计值。二阶 mtest 要求残差至少包含 3 个周期，即模型的 T=5。