• 统计学经济
  • 用plm做固定变量模型,R平方值不对,哪里出错了吗?

在用一个模型复制的练习。一个推测对照的50个国家的贸易自由化对其能源消耗量的影响的模型。
E = a0 + a1i + a2t + a3*LIBit + nit
其中E是能源消耗增长值,a0是常项,a1i是国家固定效应,a2t是时间固定效应,LIBit是贸易自由化的哑变量。其中原模型作者用了0代表贸易自由化前,1为贸易自由化中和贸易自由化后。时间t对照贸易自由化分为三个不同时间阶段。

我的数据库大概是这样的
countries t Eit LIB
country1 1 ... 0
country1 2 ... 1
country1 3 ... 1
country2 1 ... 0


我先做了一个"国家"固定效应的模型,这是我的代码

Regression <- plm(Eit ~ LIB, data = MyData,
index = c("countries","t"),
model = "within",
effect = "individual")
summary(Regression)

R平方结果出来是0.00几,原作者得出的结果是0.56。而且其他系数都不一样
请问哪里出错了呢