• R语言

  • plm, fixed effect 的R square

当我用plm做fixed effect，得到的结果里的R square很低，应该是不包括a_i的贡献的R square.

比如：

firmfixed=plm(leverage~tangibility+markettobook+logsale+profitbility,data=panel,model="within",effect="individual");


我得到的结果的R square 比OLS的低，因为只包含了四个regressor的贡献

但是我在stata里运行同样的数据和模型，R square就特别高，因为stata默认的R square包括了a_i的贡献。

如果我想要stata里的那种R square，我该怎么办？
</p>

Stata代码是这样的么？

xtreg leverage tangibility markettobook logsale profitbility, fe


Stata会输出三个R^2: within, between and overall.

感觉plm()输出的R^2和Stata输出的within R^2很接近。
</p>

从某种意义上来说，xtreg输出的R2也不准确

应该用

xi:reg leverage tangibility markettobook logsale profitbility i.id

回复 第2楼 的 tryshy：没错，plm给的应该就是within的，我要是想得到overall的有办法吗

回复 第1楼 的 ahbuhorson：请问你会做空间面板数据模型吗？

回复 第4楼 的 ahbuhorson：Stata各种R^2的计算方法参见（http://www.stata.com/statalist/archive/2006-08/msg00400.html）。感觉plm()输出的within R^2是对的，fixed effects就应该用within R^2。

回复 第5楼 的 zhangbing4502431：在学习中啊。Lemmon(2008)的论文就是我的那个模型，比较within 和 overall R sqaure, within很小，overall很大，表明那四个变量根本没有解释力，能有解释力的是firm fixed,所以他说back to the begining。但是plm根本不给overall的R sqaure

回复 第6楼 的 tryshy：within表示的是模型中变量的解释力，但是我还想要firm fixed effect的解释力。

回复 第5楼 的 zhangbing4502431：看错了，空间面板数据没做过

回复 第5楼 的 zhangbing4502431：看错了，空间面板数据没做过