短面板数据分析的基本程序

1、精选ppt 短面板数据分析的基本程序 方红生 浙江大学经济学院 2013年秋精选ppt参考书o计量经济学导论第四版（伍德里奇）中文版或英文版 o用用Stata学计量经济学学计量经济学 o高级计量经济学及stata应用（陈强）精选ppt内容安排第1讲 短面板数据分析 第2讲 长面板数据分析（PPT 第3讲 内生性与工具变量法第4讲 动态面板数据模型 第5讲 双重差分模型及其应用 第6讲 基于DID的权威文献做对了吗？ （学生报告与讨论） 第7讲 PSMDID 第8讲 如何识别核心变量的作用机制？ 精选ppt短面板数据o面板数据（panel data）是同时在时间和截面上取得的二维数据，也称时间序

2、列与截面混合数据（pooled time series and cross section data）。o是在一段时间内跟踪同一组个体的数据。既有横截面的维度（n个个体），又有时间维度（T个时期）。精选ppt StataStata中面板数据结构中面板数据结构companyyearinvestmvalue11951755.9483311952891.24924.9119531304.46241.7119541486.75593.621951588.22289.521952645.52159.4219536412031.321954459.32115.531951135.21819.4319521

3、57.32079.731953179.52371.631954189.62759.9精选pptregioncodeyearrgdpinflation北京12000北京12001北京12002北京12003北京12004北京12005北京12006北京12007北京12008天津22000天津22001精选ppto短面板：NT;反之为长面板。o平衡面板数据（balanced panel data）：如果每个个体在相同的时间内都有观测值记录。 For any i, there are T observations.o非平衡面板数据（unbalanced panel）:T may different

4、 over i.精选pptBenefits of panel data analysis精选pptittiititititituperinckunratespirconsbeertaxatal43210fuuse traffic.dta udes第一步:构造计量模型精选ppt面板数据模型p 非观测效应模型（unobserved effects model）固定效应模型(Fixed Effects Model, FE)随机效应模型(Random Effects Model, RE)o混合回归模型（Pooled Regression Model）精选ppt固定效应模型(Fixed Effects

5、Model, FE)则为固定效应模型与某个解释变量相关，如果为不可观测的个体效应其中iiitiitTtniuxuu,.,1;,.1yit精选ppt随机效应模型(Random Effects Model, RE)，则为随机效应模型与所有解释变量不相关如果为不可观测的个体效应其中iiitiitTtniuxuu,.,1;,.1yit精选ppt混合回归模型（Pooled Regression Model）Ttnixitit,.,1;,.1y, 0uiti为混合回归模型：即不存在个体效应，则如果精选ppt模型的估计o固定效应模型固定效应变换（Fixed Effects Transformation）（组

6、内变换）(Within Transformation)LSDV (Least Square Dummy Variable 精选ppt( (式式1)1)给定第给定第i i 个个体，将个个体，将( (式式1)1)两边对时间取平均可得，两边对时间取平均可得，(式式2)固定效应变换o itiititXY1iiiiXY1精选ppt (式式1)(式式2)得：得：可以用可以用OLS方法估计方法估计，称为，称为“固定效应估计量固定效应估计量”（FixedEffectsEstimator），记为），记为FE由于由于主要使用了每个个体的组内离差信息，故主要使用了每个个体的组内离差信息，故也称为也称为“组内估计量组

7、内估计量”（withinestimator）。）。FE令令 ，则，则 精选pptStata 命令oxtreg ,feoxi: xtreg i.year, fe精选pptLSDV(Least Square Dummy Variable )o基本思想：将不可观测的个体效应基本思想：将不可观测的个体效应ai看做待估看做待估计的参数，计的参数，ai就是第就是第i个个体的截距。估计个个体的截距。估计n个截个截距的方法就是引入距的方法就是引入n1个虚拟变量（如果省略常个虚拟变量（如果省略常数项，则引入数项，则引入n个虚拟变量）个虚拟变量）。精选ppto例如：共有例如：共有7个州，方程可以写成：个州，方程可

8、以写成：01112233445566ititiYXDDDDDDu7个州的回归线斜率相同，但截距不同。个州的回归线斜率相同，但截距不同。第第1个州的截距是：个州的截距是：第第2个州的截距是：个州的截距是：第第3个州的截距是：个州的截距是：第第4个州的截距是：个州的截距是：10201302403精选pptStata 命令oxi: reg i.codeoxi: reg i.code i.year精选ppt随机效应模型估计oGLS The usual pooled OLS can give consistent estimators ,but as its standard errors ignore

9、 the positive serial correlation in the composite error term, they will be incorrect.itiitstisit),/(),(Corr222精选pptoSolution: GLS transformation to eliminate the serial correlation)()(.)()1 (y1110itiitikitkkiitixxxxy:These estimators can be based on the pooled OLS or fixed effects residuals.精选pptoRa

10、ndom Effects Estimator: The feasible GLS estimator that usesin place of精选pptRE,FEandPLS（相对重要时）不重要时相对，FEREOLSEii, 1)u(R001精选pptStata 命令oxtreg ,reoxi:xtreg i.year ,re精选ppt精选ppt进一步的解释o heteroscedasticity consistent or “White” standard errors are obtained by choosing option vce(robust) which is availabl

11、e for most estimation commands.精选pptoStatas estimation commands with option robust also contain a cluster() option and it is this option which allows the computation of so-called Rogers or clustered standard errors.精选ppt But oWhile all these techniques of estimating the covariance matrix are robust

12、to certain violations of the regression model assumptions, they do not consider cross-sectional correlation. However, due to social norms and psychological behavior patterns, spatial dependence can be a problematic feature of any microeconometric panel dataset even if the cross-sectional units(e.g.

13、individuals or firms) have been randomly selected.精选ppt精选ppt精选ppt精选ppt*引入了时间虚拟变量导致exper消失精选ppt第2步:描述性统计o 变量解释与变量的描述性统计usetraffic.dta xtset state year sum fatalbeertaxspirconsunrateperincko关键变量与被解释变量的散点图并画出回归直线 twoway (scatter fatalbeertax) (lfit fatalbeertax)精选pptoPLS or FE tab year, gen(year)1.xtre

14、gfatalbeertaxspirconsunrateperinckyear2-year7,fe 这里误差项可能存在自相关、异方差和截面相这里误差项可能存在自相关、异方差和截面相关问题，所以关问题，所以F F检验显示的结果可能不可靠检验显示的结果可能不可靠, ,所所以严格的话，首先要检验是否存在截面相关问以严格的话，首先要检验是否存在截面相关问题，命令如下：题，命令如下：xtcsd,pesxtcsd,frixtcsd,fre 第3步:模型选择精选pptTestingforCross-sectionalDependenceoxtcsd 短面板oxttest2 长面板oxtcsd is a pos

15、testimation command valid for use after running an FE or RE model.oxtcsd can also perform Pesarans CD test for unbalanced panels.精选pptPLS or FEo在使用命令“xtreg,fe ”时，如果不加选项cluster(state),则输出结果还包含一个F检验，其原假设为“H0:all ui=0”,即混合回归是可以接受的。精选ppt2.xi:xtsccfatalbeertaxspirconsunrateperinckyear2-year7i.state对州虚拟变量

16、做对州虚拟变量做F F检验检验 如果不存在截面相关，则如果不存在截面相关，则xi:regfatalbeertaxspirconsunrateperinckyear2-year7i.state,cluster(state)对州虚拟变量做对州虚拟变量做F检验检验精选pptoPLS or RElxtregfatalbeertaxspirconsunrateperinckyear2-year7,relxttest0/xttest1（AR(1)）精选pptPLS or RE精选pptoFE or REuHausman test1xtregfatalbeertaxspirconsunrateperinck

17、year2-year7,feeststoreFExtregfatalbeertaxspirconsunrateperinckyear2-year7,rehausmanFE,sigmamore精选ppt Hausman检验：检验：o基本思想：如果基本思想：如果 , Fe 和和 Re 都都是一致的，但是一致的，但Re更有效。更有效。 如果如果 , Fe 仍然一致，但仍然一致，但Re是有偏是有偏的。的。 因此：因此： 如果原假设成立，则如果原假设成立，则FE与与RE估计量将共同收敛于真估计量将共同收敛于真实的参数值，反之，两者的差距过大，则倾向于拒实的参数值，反之，两者的差距过大，则倾向于拒绝原假设

18、，选择绝原假设，选择FE( ,)0iitCovX( ,)0iitCovXFE or RE精选ppt精选ppt解决办法：构造一个辅助回归精选ppt继续o基于随机效应估计的自相关检验 xtserial fatalbeertaxspirconsunrateperinckyear2-year7精选pptHausman test2oquietlyxtregfatalbeertaxspirconsunrateperinckyear2year3year4year5year6year7,reoscalartheta=e(theta)oglobalyandxforhausmanfatalbeertaxspirc

19、onsunrateperinckyear2year3year4year5year6year7osortstateoforeachxofvarlist$yandxforhausmanbystate:egenmeanx=mean(x)genmdx=x-meanxgenredx=x-theta*meanxquietlyregredfatalredbeertaxredspirconsredunrateredperinckredyear2redyear3redyear4redyear5redyear6redyear7mdbeertaxmdspirconsmdunratemdperinckmdyear2m

20、dyear3mdyear4mdyear5mdyear6mdyear7,vce(clusterstate)otestmdbeertaxmdspirconsmdunratemdperinckmdyear2mdyear3mdyear4mdyear5mdyear6mdyear7精选pptHausman test 3精选ppto基于随机效应估计的截面相关检验 xtregfatalbeertaxspirconsunrateperinckyear2-year7,rextcsd,pesabsxtcsd,friabsxtcsd,freabs 精选pptHausman test3pquietlyxtsccredf

21、atalredbeertaxredspirconsredunrateredperinckredyear2redyear3redyear4redyear5redyear6redyear7mdbeertaxmdspirconsmdunratemdperinckmdyear2mdyear3mdyear4mdyear5mdyear6mdyear7otestmdbeertaxmdspirconsmdunratemdperinckmdyear2mdyear3mdyear4mdyear5mdyear6mdyear7精选ppt第4步：报告计量结果o假设Hausman test 选择FE，则xtregfatalbeertaxspirconsunrateperinckyear2-year7,fextcsd,pesabsxtcsd,friabsxtcsd,freabs精选ppto如果存在截面相关，则最终报告由如下命令估计的结果： xtscc fatalbeertaxspirconsunrateperinckyear2-year7,fe 其中标准误是其中标准