I'm trying to estimate a panel data model with random effects Swamy and Arora 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?