Random effects estimation error with plm package d

Question

I'm trying to estimate a panel data model with random effects Swamy and Ahora transformation on market share data of tourism destinations. Please, find attached the dataset in https://github.com/Joseperles/Statistical-questions

The model is estimated according the following code:

random<-plm(ld_Share~ld_Gdp+ld_Gdppc+ld_Gkf+ld_Cpi+ld_Fdi+ld_Exrate,
 data=Panel, index=c("Country", "obs"), model="random")

However I get this error message:

Error in swar(object, data, effect) : the estimated variance of the individual effect is negative

However, when I use Gretl, Eviews or Stata to compare, I get an estimation for this model (the same output in all of them based on GLS estimation). (Please, find the output from this econometric packages at the same web adress).

With fixed effect or pooled ols I also get the same estimations with plm and all other softwares.

Somebody knows what could be the problem and if there is a possbility to reproduce this output of eviews or stata for the random model (the same estimation procedure) with the plm R-package?

Thanks in advance and sorry for any inconvenience.