stata经验与疑惑

question:

陈强老师的教材：用stata16做短面板随机效应模型，试了不同数据，似乎当变量数小于时期数，回归结果就会显示sigma_u=0，rho=0，而当变量数大于时期数，却不会

answer1:

古扎拉蒂（2013）：如果长面板（T较大且N）较小，两种方法的估计值没什么差别。此时，固定效应模型更加可取。
在短面板（T小且N大）中，两种方法的估计值可能有显著差别。如果确信样本中横截面 单元不是从一个较大的总体中随机抽取的，那么固定效应模型是可取的。如果不是这种情 况，则随机效应模型是可取的。如果随机效应模型的假设成立，那么随机效应模型的估计 量比固定效应模型的估计量更加有效。

answer2:

默认情况下，stata采用GLS估计RE模型，此时 simga_e 采用固定效应模型的残差计算而得，而 sigma_u 则是利用 simga_e 的估计值，以及组间效应模型的估计系数计算而得。参见 [xt] p.464。计算过程中，由于解释变量中包含了 L.y，致使 sigma_u=sigma_e=0，并最终导致 OLS 和 RE 的估计结果完全相同，此时 RE 的估计结果是错的。

应该采用 MLE 来估计 RE 模型：
xtreg y L.y x1 x2, mle

此外，由于这个模型是动态面板模型，最好采用 GMM 估计，stata命令为 xtabond 或 xtdpdsys。