问题描述
我有一个不平衡的面板数据集，需要对其进行（面板）格兰杰因果关系测试。我已将数据文件从Stata（.dta）导入到R以执行功能pgrangertest（来自软件包plm）（有关此测试的更多信息，请参见以下链接：https://www.rdocumentation.org/packages/plm/versions/2.2-4/topics/pgrangertest ）。
但是，当我执行此测试时，会收到以下错误消息：
pgrangertest(GDPGrowth5Y ~ ratingUp,data = CoreDataSet)
Error in rval[i,j,drop = drop.,...] : subscript out of bounds
为了调试，我做了以下事情：
options(error=recover)
reach_full_in <- reachability(krack_full,'in')
不幸的是，我对R的经验并不丰富，因为直到现在我才需要在Stata工作。
我有1980年至2020年期间131个国家/地区的数据，其中包含12个变量。该.dta文件在Stata中可以正常运行。
我不确定，但我认为此错误可能是由于R中的数据集无法以与Stata相同的方式识别Panel数据引起的。如果我是正确的，那么R当前将年份（41）视为N，将国家数量（131）视为T？见所附图片。
我的问题是：这会导致错误吗？如果是这样，我该如何解决/翻转呢？如果没有，您认为问题出在哪里？
解决方法
暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！
如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@）