我正在尝试估算具有随机效应的面板数据模型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的输出？
提前致谢，不便之处。