我正在尝试估计一个面板数据模型，该模型具有随机效应Swamy和Ahora对旅游目的地市场份额数据的转换 . 请在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（对象，数据，效果）中的错误：单个效果的估计方差为负
但是，当我使用Gretl，Eviews或Stata进行比较时，我得到了对该模型的估计（基于GLS估计的所有模型中的相同输出） . （请在同一网址上找到此计量经济学软件包的输出） .
对于固定效果或合并的ols，我也可以使用plm和所有其他软件进行相同的估算 .
有人知道可能是什么问题，是否有可能使用plm R-package为随机模型（相同的估算程序）重现eviews或stata的输出？
提前致谢，不便之处，敬请原谅 .