我使用了固定效果和第一次差异估计。为了确定哪种效率更高,Wooldridge提出了一种特定的测试,该测试通过以下功能包含在plm软件包中:
pwfdtest(Y~X1+X2+..., data=...)
但是,运行此操作会导致出现一条错误消息,提示我:
> pwfdtest(DepVar~ExplVar1+ExplVar2, data = data)
Error in `$<-.data.frame`(`*tmp*`, "FDres", value = c(-1.18517291896221, :
replacement has 521293 rows, data has 621829
In addition: Warning message:
In id - c(NA, id[1:(N - 1)]) :
longer object length is not a multiple of shorter object length
我试图查找是否有人在发布之前遇到此错误,但是我找不到答案。
有时候,我遇到了一些要求最低工作实例的人,但是我的工作根本没有。但是,plm软件包中的示例确实起作用。
请注意,这是我的首次研究,也是我第一次使用R。所以请多多包涵。
最良好的祝愿 亚历克斯
编辑:
我读到traceback()函数可能有点有用。但是,大多数情况下,它只是吐出各种各样的内容,我什至都无法到达顶部(?), 这些数字的最后几行是:
-1.65868856541809, 2.89084861854684, -1.68650260853188, 0.655681663187397,
-0.677329685017227, 0.993684102310348, 1.33441048058398, -2.0526651614649,
-1.64392358708552, 2.58673448155514, 0.952616064091869, -0.909754051474562,
0.815593306056627, -0.0542364686765445, 0.0184515528912868))
2: pwfdtest.panelmodel(fd1)
1: pwfdtest(fd1)
EDIT2:
我的第一个猜测是NA可能令人不快,因此我只将面板缩小为因变量和一个解释变量。事先,我检查是否有NA,而没有。还有一个错误消息:
Error in `$<-.data.frame`(`*tmp*`, FDres, value = c(-1.18517291896221, :
replacement has 521293 rows, data has 621829
In addition: Warning message:
In id - c(NA, id[1:(N - 1)]) :
longer object length is not a multiple of shorter object length
EDIT3:
我想我可能已经发现了问题:面板不平衡。我想这是有道理的……但是,传统意义上似乎没有解决方案,它根本行不通。
所以,如果有人对我所做的事情感兴趣:
我进一步将小组人数减少到300人以下。我将个人命名为1-300,并且 drumroll 起作用了。但是,将某些个人名称更改为例如555或556后,它给了我与以前相同的错误。
我对这些事情不是很精通,但是我没有学到的猜测是,该测试根本无法在不平衡的面板上进行。