我正在尝试估计一个面板数据模型，该模型对旅游目的地的市场份额数据具有随机效应 Swamy 和 Arora 转换。请在https://github.com/Joseperles/Statistical-questions中找到附加的数据集
该模型根据以下代码进行估计：
random<-plm(ld_Share~ld_Gdp+ld_Gdppc+ld_Gkf+ld_Cpi+ld_Fdi+ld_Exrate,
data=Panel, index=c("Country", "obs"), model="random")
但是我收到此错误消息：
swar(object, data, effect) 中的错误：个体效应的估计方差为负
然而，当我使用 Gretl、EViews 或 Stata 进行比较时，我得到了这个模型的估计（基于 GLS 估计，它们的输出相同）。（请在同一网址找到此计量经济学软件包的输出）。
使用固定效果或合并 ols，我也可以使用 plm 和所有其他软件获得相同的估计。
有人知道可能是什么问题，以及是否有可能使用 plm R 包为随机模型（相同的估计过程）重现 EViews 或 Stata 的输出？