使用Stata中的数据进行R中的pgrangertest

我有一个不平衡的面板数据集，需要对其进行（面板）格兰杰因果关系测试。我已将数据文件从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？见所附图片。

我的问题是：这会导致错误吗？如果是这样，我该如何解决/翻转呢？如果没有，您认为问题出在哪里？