STATA面板数据模型操作命令.doc

1、STATA 面板数据模型估计命令一览表一、静态面板数据的 STATA 处理命令固定效应模型itxyitiitititit随机效应模型ititit（一）数据处理输入数据tsset code year 该命令是将数据定义为“面板”形式xtdes 该命令是了解面板数据结构summarize sq cpi unem g se5 ln 各变量的描述性统计（统计分析）gen lag_y=L.y / 产生一个滞后一期的新变量gen F_y=F.y / 产生一个超前项的新变量gen D_y=D.y / 产生一个一阶差分的新变量gen D2_y=D2.y / 产生一个二阶差分的新变量（二）模型的筛选和检验1、检

2、验个体效应（混合效应还是固定效应） （原假设：使用 OLS 混合模型）xtreg sq cpi unem g se5 ln,fe对于固定效应模型而言，回归结果中最后一行汇报的 F 统计量便在于检验所有的个体效应整体上显著。在我们这个例子中发现 F 统计量的概率为0.0000，检验结果表明固定效应模型优于混合 OLS 模型。2、检验时间效应（混合效应还是随机效应） （检验方法：LM 统计量）（原假设：使用 OLS 混合模型）qui xtreg sq cpi unem g se5 ln,re (加上“qui”之后第一幅图将不会呈现)xttest0可以看出，LM 检验得到的 P 值为 0.0000，

3、表明随机效应非常显著。可见，随机效应模型也优于混合 OLS 模型。3、检验固定效应模型 or 随机效应模型 （检验方法：Hausman 检验）原假设：使用随机效应模型（个体效应与解释变量无关）通过上面分析，可以发现当模型加入了个体效应的时候，将显著优于截距项为常数假设条件下的混合 OLS 模型。但是无法明确区分 FE or RE 的优劣，这需要进行接下来的检验，如下：Step1：估计固定效应模型，存储估计结果Step2：估计随机效应模型，存储估计结果Step3：进行 Hausman 检验qui xtreg sq cpi unem g se5 ln,feest store fequi xtreg

4、 sq cpi unem g se5 ln,reest store rehausman fe (或者更优的是 hausman fe,sigmamore/ sigmaless)可以看出，hausman 检验的 P 值为 0.0000，拒绝了原假设，认为随机效应模型的基本假设得不到满足。此时，需要采用工具变量法和是使用固定效应模型。（三）静态面板数据模型估计1、固定效应模型估计xtreg sq cpi unem g se5 ln,fe (如下图所示)其中选项 fe 表明我们采用的是固定效应模型，表头部分的前两行呈现了模型的估计方法、界面变量的名称（id） 、以及估计中使用的样本数目和个体的数目。第

5、 3 行到第 5 行列示了模型的拟合优度、分为组内、组间和样本总体三个层面，通常情况下，关注的是组内（within） ，第 6 行和第 7 行分别列示了针对模型中所有非常数变量执行联合检验得到的 F 统计量和相应的 P 值，可以看出，参数整体上相当显著。需要注意的是，表中最后一行列示了检验固定效应是否显著的 F 统计量和相应的 P 值。显然，本例中固定效应非常显著。2、随机效应模型估计若假设本例的样本是从一个很大的母体中随机抽取的，且 与解释变量均i不相关，则我们可以将 视为随机干扰项的一部分。此时，设定随机效应模型i更为合适。xtreg sq cpi unem g se5 ln,re (如下

6、图所示)3、时间固定效应（以上分析主要针对的是个体效应）如果希望进一步在上述模型中加入时间效应，可以采用时间虚拟变量来实现。首先，我们需要定义一下 T-1 个时间虚拟变量。tab year ,gen(dumt) (tab 命令用于列示变量 year 的组类别，选项 gen（dumt）用于生产一个以 dumt 开头的年度虚拟变量)drop dumt1 （作用在于去掉第一个虚拟变量以避免完全共线性）若在固定效应模型中加入时间虚拟变量，则估计模型的命令为：xtreg sq cpi unem g se5 ln dumt*,fe（四）异方差和自相关检验1、异方差检验 （组间异方差）本节主要针对的是固定效

7、应模型进行处理（1）检验原假设：同方差 需要检验模型中是否存在组间异方差,需要使用 xttest3 命令。qui xtreg sq cpi unem g se5 ln,fexttest3显然，原假设被拒绝。此时，需要进一步以获得参数的 GLS 估计量，命令为xtgls：xtgls sq cpi unem g se5 ln,panels(heteroskedastic)其中，组间异方差通过 panels（）选项来设定。上述结果是采用两步获得，即，先采用 OLS 估计不考虑异方差的模型，进而利用其残差计算。 。 。 ，并最终得到 FGLS 估计量。2、序列相关检验对于 T 较大的面板而言， 往往无

8、法完全反映时序相关性，此时 便可i it能存在序列相关，在多数情况下被设定为 AR(1)过程。原假设：序列不存在相关性。（1） FE 模型的序列相关检验对于固定效应模型，可以采用 Wooldridge 检验法，命令为 xtserial:xtserial sq cpi unem g se5 ln可以发现，这里的 P=0.0000，我们可以在 1%的显著性水平下爱拒绝不存在序列相关的原假设。考虑到样本，该检验的最后一步是用 对 进行 OLS 回归，eit1-ti，因此，输入以下命令得到 。检验该值是否显著异于-0.5，因为在原85.0假设下（不相关） ，可见本例中不相等，拒绝原假设，说明存在序列.

9、-相关。mat list e(b)（2） RE 模型的序列相关检验对于 RE 模型，可以采用 xttest1 命令来执行检验：qui xtreg sq cpi unem g se5 ln dumt*,rexttest1这里汇报了 4 个统计量，分别用于检验 RE 模型中随机效应（单尾和双尾） 、序列相关以及二者的联合显著性，检验结果表明存在随机效应和序列相关，而且对随机效应和序列相关的联合检验也非常显著。（3） 稳健型估计上述结果表明，无论是 FE 还是 RE 模型，干扰项中都存在显著的序列相关。为此，我们进一步采用 xtregar 命令来估计模型，首先考虑固定效应模型：xtregar sq

10、cpi unem g se5 ln dumt*,fe lbi3、 “异方差 序列相关”稳健型标准误 虽然上述估计方法在估计方差-协方差矩阵时考虑了异方差和序列相关的影响，但都未将两者联立在一起考虑，要获得“异方差-序列相关”稳健型标准误，只需在 xtreg 命令中附加 vce(robust)或者 vce(cluster)选项即可。例如，对于 FE 模型，我们可以执行如下命令：xtreg sq cpi unem g se5 ln,fe vce(robust)与之前未经处理的估计结果相比，附加命令 vce(robust)选项时的结果，虽然系数的估计值未发生变化，但此时得到的标准误明显增大了，致使得

11、到的估计结果更加保守。对于面板数据模型而言，STATA 在计算所谓的“robust”标准误时，是以个体为单位调整标准误的。因此，我们得到的“robust”标准误其实是同时调整了异方差和序列相关后的标准误。换言之，上述结果与设定vce(cluster)选项的结果完全相同。4、截面相关检验 原假设：截面之间不存在着相关性（1）FE 模型检验对于 FE 模型，可以利用 xttest2 命令来检验截面相关性：qui xtreg sq cpi unem g se5 ln,fexttest2(该命令主要针对的是大 T 小 N 类型的面板数据，在本例中无法使用，故图标略去。)（2）RE 模型检验对于 RE 模型，可以利用 xtcsd 命令来检验截面相关性：qui xtreg sq cpi unem g se5 ln,rextcsd,pesaran (下面命令是另一个检验指标)