stata笔记教学提纲

1、1 .一般检验假设系数为0,t比较大则拒绝假设，认为系数不为0.假设系数为0,P比较小则拒绝假设，认为系数不为0.假设方程不显著，F比较大则拒绝假设，认为方程显著。2.小样本运用OLS进行估计的前提条件为：(1)线性假定。即解释变量与被解释变量之间为线性关系。这一前提可以通过将非线性转换为线性方程来解决。(2)严格外生性。即随机扰动项独立于所有解释变量：与解释变量之间所有时候都是正交关系，随机扰动项期望为0。(工具变量法解决)(3)不存在严格的多重共线性。一般在现实数据中不会出现，但是设置过多的虚拟变量时，可能会出现这种现象。Stata可以自动剔除。(4)扰动项为球型扰动项，即随即扰动项同方差

2、，无自相关性。3.大样本估计时，一般要求数据在30个以上就可以称为大样本了。大样本的前提是(1)线性假定(2)渐进独立的平稳过程(3)前定解释变量，即解释变量与同期的扰动项正交。4 4)E(XiX让)为非退化矩阵。5 5)gt为鞅差分序列，且其协方差矩阵为非退化矩阵。与小样本相比，其不需要严格的外生性和正太随机扰动项的要求。4.命令稳健标准差回归：regyx1x2x3,robust回归系数与OLS一样，但标准差存在差异。如果认为存在异方差，则使用稳健标准差。使用稳健标准差可以对大样本进行检验。只要样本容量足够大，在模型出现异方差的情况下，使用稳健标准差时参数估计、假设检验等均可正常进行，即可以

3、很大程度上消除异方差带来的副作用对单个系数进行检验：testlnq=1线性检验：testnl_blnpl=_blnqA25.如果回归模型为非线性，不方便使用OLS则可以采取最大似然估计法(MLE)，或者非线性最小二乘法(NLS)6.违背经典假设，即存在异方差的情况。截面数据通常会出现异方差。因此检验异方差可以：(1)看残差图，但只是直观，可能并不准确。rvfplot(residual-versus-fittedplot)与拟合值的散点图rvpplotvarname(residual-versus-predictorplot)与解释变量的散点图扰动项的方差随观测值而变动，表示可能存在异方差。(2

4、)怀特检验：estatimtest,white(post-estimationinformationmatrixtest)P比较小，则拒绝同方差假设，表示存在异方差，不能用OLS反之则证明为同方差。(3)BP检验estathettest,iid(默认设置为使用拟合值yA)estathettest,rhsiid(使用方程右边的解释变量，而不是yA)estathettestvarlist,iid(使用某个指定的解释变量)P小，则拒绝原假设。如果存在异方差，则可以：(1)使用OLS+急健标准差robust(2)广义最小二乘法(GLS)(3)加权最小二乘法(WLS)predictel,res(预测残差

5、)ge2=e|A2辅助回归：glne2=log(e2)reglne2lnq,nocpredictlne2f计算辅助回归的拟合值ge2f=exp(lne2f)去掉对数即权重之倒数reglntclnqlnpllnpklnpfaw=1/e2fregyx1x2x3aw=1/var(aw表示analyticalweight,var表示随即扰动项的方差。)或者：predictu,residualspredictyf,xbgenlnu2=ln(uA2)genyf2=yfA2quietlyreglnu2yfyf2predictnlu2f=exp(xb()gensd=sqrt(u2f)vwlslntclnqln

6、pllnpflnpk,sd(sd)(4)可行广义最小二乘法(FGLSFGLS所做的过程和GLS一样， 只是GLS假设扰动项的方差已知， 若要用GLS,必须计算得到扰动项方差，而FGLS则是在未知方差的情况下求方差并最终通过将异方差转换为同方差后再运用OLS的结果。因此，GLS和FGLS在过程上是一致的。6.自相关时间序列中容易出现自相关，而截面数据也可能存在空间自相关。人为处理数据如移动平均等做法也可能导致自相关。检验自相关可以：(1)作图，但并不严格。定义滞后算子L.(只有时间序列数据和面板数据才能定义时间变量。)tssetyaear一阶差分：D.x=xt-xt-1D2.X=xt-xt-2L

7、D.表示一阶差分的滞后值画图：scatterelL.elacel(看自相关图)pacel(看偏相关图)(2)BG检验estatbgodfrey(默认p=1)estatbgodfrey,lags(p)estatbgodfrey,nomiss0(使用不添加0的BG检验)使用命令ac查看自相关图，或者设置较大的p值进行显著性检验，t期不显著了，则选择P=T-1统计检验P值小，则拒绝假设。(3)box-pierceQ检验/Ljung-BoxQregyx1x2x3predictel,residwntestqel(使用stata提供的默认滞后期)wntestqel,lags(p)(使用自己设定的滞后期)(

8、4)DW检验：现在已经不常用，因为其只能检验一阶自相关。estatdwatson自相关的处理方法：(1)使用OLS+异方差自相关稳健的标准差(HeteroskedasticityandAutocorrelationConsistentStandardError,HAC)neweyyx1x2x3,lag(p)(HAC标准差，必须制定滞后阶数p)滞后期数选择nA1/4(2)使用OLS噪类稳健的标准差(clusterrobuststandarderror)面板数据中经常使用聚类稳健的标准差。regyx1x2x3,cluster(state)(聚类稳健标准差，假设state”为聚类变量)(3)使用可行

9、广义最小二乘法(FGLSpraisyx1x2x3(使用默认的PW估计法)praisyx1x2x3,corc(使用CO估计法)(4)修改模型设定，可能自相关是由于遗漏了自相关的解释变量。7多重共线性在回归后，使用命令VIFestatvif经验表示，vif=1978)如希望将每个省设置为虚拟变量，则需要：tabulateprovince,generate(pr)回归简化为：regyx1x2x3pr2-pr3111.工具变量法这可以解决扰动项与自变量的相关问题， 设置的工具变量需要与扰动项无关而与内生解释变量相关。传统的工具变量法一般通过两阶段最小二乘法TSLS2SLS(twostageleasts

10、quare)。第一阶段，工具变量对内生解释变量回归；第二阶段，被解释变量对工具变量的拟合值进行回归。多个工具变量的线性组合仍然可以作为工具变量。命令：ivregress2slsdepvarvarlist1(varlist2=inslist)Depvar为被解释变量，varlist1为外生解释变量，varlist2为内生解释变量，instlist为工具变量。如：ivregress2slsyx1(x2=z1z2)ivregress2slsyx1(x2x3=z1z2z3z4),rfirst(r表示用异方差的标准差，first表示在结果中显示第一阶段的回归。)在面板数据中执行2SLSM以用：xtivr

11、egdepvarvarlistl（varlist_2=varlist_iv）（详见helpxtivreg）检验工具变量与解释变量的相关性：即检验工具变量是否为弱工具变量，命令：estatfirststage,allforcenonrobust（all表示显示每个内生变量的统计量，而非仅仅所有内生变量综合的统计量，forcenonrobust表示及时在进行工具变量法时用了稳健标准差，也仍然允许计算estatfirststage）解决弱工具变量的方法包括A.寻找更强的工具变量B.弱工具变量较多，则舍弃弱工具变量，C.用有限信息最大似然估计法（Limitedinformationmaximumlik

12、elihoodestimation,LIML）LIML与2SLSB进等价， 但在弱工具变量的情况下，LIML的小样本性质可能优于2SLS.命令为：ivregresslimldepvarvarlist1（varlist2=instlist）过度识别（即多余的工具变量的个数）命令为：estatoverid但并不能告诉哪些工具变量无效。使用工具变量的前提是存在内生解释变量（即解释变量与扰动项相关），这也需要检验。如果所有解释变量都是外生变量则用OLS比用工具变量法更有效，反之应该用工具变量法。豪斯曼检验就是假设所有解释变量都为外生变量。豪斯曼检验的stata命令：regyx1x2estimatess

13、toreols（存储OLS的结果）ivregress2slsyx1（x2=z1z2）（假设怀疑x2为内生变量）estimatesstoreiv（存储2SLS的结果）hausmanivols,constantsigmamore（根据存储的结果进行豪斯曼检验）但uguo存在异方差，则OLS并不是最有效的，传统额豪斯曼检验不适用于异方差的情形，解决方法是“自助法”；或者使用“杜宾-吴-豪斯曼DWH检验”也可以适用于存在异方差的情况。命令：estatendogenous在球型扰动项的假定下，2SLS是最有效的，但是如果扰动项存在异方差或者自相关，则广义矩估计（generalizedmethodofmo

14、ments,GMM）更有效。GMM与2SLS的关系就相当于GLSWOLS之间的关系。GMM过程：sccinstallivreg2（安装程序ivreg2）sccinstallranktest（安装另外一个在运行ivreg2时需要用到的辅助程序ranktest）usedataxtsetpanelvartimevar（设置面板变量和时间变量）ivreg2yx1（x2=z1z2）,gmm2s12.短面板：（固定效应和随机效应模型）N大T小的一般叫做短面板。面板数据的优点为：可以解决遗漏变量的问题，可以提供更多动态行为的信息，样本较多从而估计更准确。但面板数据一般不满足独立同分布的假设。处理面板数据的一

15、个方法是将面板数据当作横截面数据处理进行OLS回归，称为“混合回归”，但它忽略了同一个聚类存在的相关问题。虽然通常可以假定面板数据不同个体的扰动项的独立性，但是对于同一个体却存在固定效应和随机效应。固定效应模型通常可以采用组内估计法（FE）和一阶差分法（FD）。当T=2时，FD=FE;当T2,扰动项独立同分布时，FE优于FDo因此，实践用FE较多，而对动态面板，则用FD较多。随机效应的存在使得OLS估计是一致但无效的。可使用广义最小二乘法（FGLS来进行估计。也可以使用组间估计量（BE）用豪斯曼检验选择是选用固定效应模型还是随机效应模型，如果统计量大于零界值，则拒绝假设。如果假设成立，RE最有

16、效，但是不适用于异方差的情形。解决方法是自助法和辅助回归。非平衡面板经常会损失数据，导致破坏样本的随机性。xtsetpanelvatimevar（设定面板数据的面板个体变量和时间变量）encodecountry,gen（cntry）（为面板个体编号使得面板个体为整数）显示面板数据统计特征：xtdes（显示面板数据结构，是否为平衡面板）xtsum（显示组内、组间和整体的统计指标）xttabvarname（显示组内、组间与政体的分布频率）xtlinevarname（对每个个体分别显示该变量的时间序列图， 如果希望叠放，则选择overlay）regyx1x2x3,vce（clusterid）（混合回

17、归，VCE是以id为cluster的聚类稳健标准差，因为同一地区不同时间扰动项之间一般存在自相关。）VCE是考虑到了同一聚类之间的同方差现象。xtreg丫X1X2X3,fevce（clusterid）（使用固定效应模型回归）得到的回归中rho表示来自个体效应ui的比例。XtregYX1X2X3,fe中的F检验通过，则可以使用混合回归模型。若拒绝假设（p小），则FE更优，每个个体都有自己的固定效应。但此时的由于没有使用稳健回归，因此F检验并不有效，还需要进一步用LSDV法观察。LSDV去白stata命令为：xi:xtregyx1x2x3i.id,vce（clusterid）（xi为增添互动项in

18、teractionexpansion,i.id表示根据确定个体变量id生成的虚拟变量， 在这里是state）（P小，则说明个体虚拟变量很显著，因此，不满足无个体效应的假设，则不能用混合回归模型。）对于固定效应模型，也可以用一阶差分法FD的命令：xtserialyx1x2x3,output（我的stata不能识别xtserial咋回事？）一般认为FE比FD有效，故较少使用FD=也可以在固定效应中考虑时间效应，即双向固定效应：tabyeargen（year）定义年度虚拟变量xtregYX1X2X3year2-year7,fevce（clusterid） （明明是stata12,为什么说数据太多不能

19、计算？要把前面的计算保存等清空，那后面怎么对比？）testyear1year2year3。year7对年度虚拟变量的联合检验为检验存在时间效应和固定效益后，还可能存在随机效应，对随机效应进行检查：xtregyx1x2x3,revce（clusterid）（随机效应FGLS）xtregyx1x2x3,mle（随机效应MLE）检验个体随机效应存在一个LM检验，需要在随机效应估计之后进行：Xtregfatalspirconsunrateperincbeertax,revce（cluster,state）estimatesstoreRExttest0如果P很小，则拒绝假设，认为在“随即效应”与“混合回

20、归”之间，应该选择“随即效应”。如果数据质量不好，可以考虑使用组间估计量，但会损失较多信息量。regyx1x2x3,beestimatesstoreBE选用固定效应模型还是随机效应模型运用豪斯曼检验xtregyx1x2x3,feestimatesstoreFExtregyx1x2x3,reestimatesstoreREhausmanFERE,constantsigmamoreP小，则强烈拒绝使用随机效应模型，而应采用固定效应模型。但是如果聚类稳健标准差和普通标准差相差较大， 则传统的豪斯曼检验并不适用。这时需要进行辅助回归。目前在stata中需要手动进行。步骤如下：quietlyxtregf

21、atalbeertaxspirconsunrateperinck,rescalartheta=e(theata)(从回归中得到THETA)globalyandxforhausmanfatalbeertaxspirconsunrateperinck(使用yandxforhausman时，代表了所有使用的变量)foreachxofvarlist$yandxforhausmanbystate:egenmeanx=mean(x)至U底用什么表示？genmdx=x-mean父genredii-theata*meanxquietlyxtregredfatalredbeertaxredspirconsred

22、unrateredperinckmdbeertaxmdspirconsmdunratemdperinck,vce(clusterstate)testmdbeertaxmdspirconsmdunratemdperinck检验发现P小，则拒绝r=0的假设，应该选用固定效应模型。(疑问：那中间那个式子呢？)面板数据在回归之前要通过平稳性检验，其处理步骤简化为：面板数据单位根检验一协整检验一回归。一般为了方便，米用相同单位根检验LLC检验和不同单位根检验Fisher-ADF检验(非面板数据一般采用ADF检验)。两种检验均拒绝存在单位根的元假设，则认为序列平稳。13.长面板和动态面板(不能像短面板那样

23、假定独立同分布，而应该运用广义最小二乘法FGLS进行估计，解决组内和组间的自相关。)解释变量包含被解释变量的滞后项，则为动态面板，反之为静态面板。(1)仅解决组内自相关的FGLStabstate,gen(state)考虑个体效应，生成州虚拟变量。gent=year-1962考虑时间趋势，生成时间趋势变量reg丫X1X2X3State2-state10t,vce(clusterstate)用LSDV估计双向固定效应模型estimatesstoreOLS考虑组内自相关的情形：用命令stpcsestpcseYX1X2X3state2-state10t,corr(ar1)约束条件为自回归系数都相同的一

24、阶。考虑各组自回归系数不同的组内自相关情形stpcseYX1X2X3state2-state10t,corr(psar1)如果仅考虑不同个体扰动项存在异方差而忽略自相关，则可以用stpcseYX1X2X3state2-state10t,hetonly将以上各估计的系数及标准差列表便于比较：究竟应该采用OLS还是采用AR1,则需要检验自相关。(2)同时处理组内自相关和组间同期相关的FGLS的命令xtglsxtglsYX1X2X3state2-state10t,panels(iid/het/cor)corr(ar1/psar1)iid表示不同个体(组间)的扰动项独立且具有相同的方差，het表示不同

25、个体的扰动项独立但具有不同的方差，cor表示不同个体的扰动项同期相关且具有不同的方差。Corr表示组内，解释同前。如果加上igls则表示用的迭代FGLS而非两步FGLS.若执行xtgls或者xtpcse时没有使用个体虚拟变量，则为随机效应模型，若使用了则为固定效应模型。(3)长面板数据自相关和异方差检验组间异方差检验，运用似然比检验XtglsYX1X2X3,iglspanel(het)(允许异方差的迭代式FGLS估计)P=0,则强烈拒绝LR检验的“组间同方差”假设。另一种检验组间异方差的方法为xttest3,只能在xtreg,fe或者xtgls之后使用:Sscinstallxttest3(安装

26、命令xttest3)Xttest3(进行组间异方差检验)这个超级简单组内自相关检验通过下载命令xtserial来实现这个检验netinstallst0039(安装命令stserial)为什么下不了？xtserialyx1x2x3,output(进行组内组相关检验，若P小，则强烈决绝不存在一阶组内自相关的假设。)组间截面相关检验： 非官方xttest2。 其仅能在xtgls,ivreg2,和xtreg,fe后使用，且能用于场面版sscinstallxttest2(安装xttest2)xttest2(组间截面相关检验)当这一种方法无法检验时，可以采用另一种方法，xtcsd命令，长短面板Estima

27、tesstorehetroXtglsyx1x2x2,iglsEstimatesstorehomolocaldf=e(N_g)-1lrtestheterohomo,df(hetero)(同方差条件下的FGLS)(将同方差条件下的估计结果储存为(计算自由度，即约束条件个数，其中df)(制定自由度，进行似然比检验。homo)e(N_g)为个体个数)符号到底是怎样的？都能用，用在fe之后？sscinstallxtcsd(安装命令xtcsd)xtcsd,pesaranabsshow(pesaran的检验，服从正态分布)xtcsd,friedmanabsshow(friedman的检验，服从喀方分布)xt

28、csd,freesabsshow(frees的检验)对于长面板除了可以让不同个体拥有不同的截距项外还可以使其拥有不同的斜率，这成为“变系数模型”，如果其系数为常数，则可分别回归，但如果各个体扰动项相关，则应该把所有个体回归方程叠放，然后使用“似不相关回归(SUR),但由于参数估计较多，会损失自由度。如果考虑“部分变系数模型，则SUR不再适用，而应用LSDV法，引入虚拟变量，以及虚拟变量与可变细数解释变量的互动项(参见helpxi)如果将系数看作是随机的，则FGLS古计模型，即利用OLS残差估计协方差矩阵中的参数，然后再使用GLS命令为：Xtrcyx1x2x3,betas(betas表示显示对每

29、一组系数的估计),其附带参数稳定性检验虽然面板数据能在一定程度上解决遗漏变量问题， 但如果存在内生解释变量，还是需要用工具变量法，先解决遗漏问题，再使用2SLS对于动态面板，就算是组内估计量(FE)也是不一致的。差分GMM和水平GMM结合就成了系统GMM估计。差分GMM的stata命令为：xtabonddepvarindepvars,lags(p)maxldep(q)pre(varlist)endogenous(varlist)inst(varlist)twostepvce(robust)如：xtabondlwageoccsouthsmsaind,lags(2)maxldep(3)pre(wk

30、s,lag(1,2)endogenous(ms,lag(0,2)endogenous(union,lag(0,2)twostepvce(robust)差分GMM存在的前提是扰动项不存在自相关，对此要进行检验，需要扰动项的差分不存在二阶或更高阶自相关即可。命令为：estatabond,p小，则存在自相关，P大，则不存在自相关。更高阶自相关的检验为：estatabondartests(3)由于使用过多的工具变量，还需要进行过度识别检验先回归，但要去掉最后的vce(robust)estatsargan系统GMM的stata命令为：xtdpdsysdepvarindepvars,lags(p)maxl

31、dep(q)pre(varlist)endogenous(varlist)inst(varlist)twostepvce(robust)若要对比差分GMM和系统GMM的系数和标准差，则用命令：estimatestableDGMMSGMM,bse14.离散被解释变量（通常不适合用OLS而应该用二值模型，probitlogit和多值选择模型）如果F为标准正态的累积分布函数，则P模型为Probit模型，若F为逻辑分布的累积分布函数，则P模型为Logit模型。计算Logit模型通常比Probit模型更方便。二值模型的Stata命令为：probityx1x2x3logityx1x2x3probit和lo

32、git的分布函数不同，因此参数不能直接比较，需要分别计算二者的边际效应，再进行比较，STAT配行处理的命令为：mfx（计算在样本均值处的边际效应）mfx,at（X1=0）（计算在X1=0时，X2,X3取值样本均值处的边际效应）mfx,eyex（计算在样本均彳I处的弹性）predictyhat（计算发生概率的预测值，并记为yhat）,对于Logit模型，系数表示解释变量x增加一个单位将引起的“对数几率比”的边际变化。衡量二值模型的拟合优度采用“准R2”，判断拟合优度还可以通过计算预测准确的百分比，相应STA怖令为：estatclaslogti和probit模型虽然估计系数不同，但其估计系数没有可

33、比性，其m仅计算出的边际效应及准R2与正确预测比几乎一致。如果接受似然比检验（LR，则可用同方差probit模型估计。在存在异方差的情况下进行probit模型估计为：hetprobYX1X2X3,het（varlist）那自相关呢？多值选择模型：多值选择模型的stata命令为：mlogitYX1X2X3,base（#）（多值logit选择模型，base的用来指定参照组）mlogitYX1X2X3,rrrbase（#）（多值logit选择模型，回报relativeriskratio）mprobitYX1X2X3,base（#）（多值probit选择模型）排序数列模型的stata命令：ordere

34、dprobit/logitoprobitYX1X2X3ologitYX1X2X3oprobit预测：predictp2p3p4p5列出第一个观测值的预测结果：listp2p3p4p5in1/1,ologit预测：predictr2r3r4r5列出第一个观测值的预测结果：listr2r3r4r5in1/1,计数模型（1）有些被解释变量仅能取非负整数，如金牌数量，看病次数等，一般用泊松回归。Stata命令为poissonyx1x2x2,r（稳健标准差）estatgof（goodnessoffitnes拟合优度检验）泊松回归的局限是其期望和方差一定相等， 但有些被解释变量的方差明显大于期望，即存在过

35、度分散”,则可以考虑“负二项回归（negativebinomialregression）,使用MLE估计。Nbregyx1x2x2,r如果计数数据中含有大量的零值， 则可以使用“零膨胀泊松回归这可以Vuong统计量来检验， 如果Vuong统计量很大， 则选择零膨胀泊松回归或者零膨胀负二项泊松回归。零膨胀泊松回归的stata命令为：zipyx1x2x3,inflate(varlist)vuong(零膨胀泊松回归)zipnbyx1x2x3,inflate(varlist)vuong(零膨胀负二项回归)15.受限制的被解释变量(断尾回归、截取回归)断尾回归stata命令：truncregyx1x2x

36、3,ll(#)(左边断尾)truncregYX1X2X3,ul(#)(右边断尾)truncregYX1X2X3,ll(#)ul(#)(双边断尾)截取回归stata命令：tobit模型tobityx1x2x3,ll(#)tobityx1x2x3,ul(#)tobityx1x2x3,ll(#)ul(#)样本选择模型的stata命令：Heckmanyx1x2x3,select(z1z2)(默认使用MLE,选择方程的被解释变量为y),最下的似然比检验P很小，则认为样本选择模型适用。Heckmanyx1x2x3,select(z1z2)twostep(两步法,选择方程的被解释变量为y)Heckmanyx

37、1x2x3,select(w=z1z2)(默认使用MLE选择方程的被解释变量为w)16时间序列(平稳时间序列、非平稳时间序列。AR自回归模型、MA为移动平均模型结合起来为ARMA模型、自回归分布滞后模型ADL、向量自回归模型VAR向量移动平均过程VMA、格兰杰因果检验：条件是变量协整)(1)自相关和偏自相关corrgramsy,lags(#)第1至#阶ACF和PACFacy,lags(#)pacy,lags(#)(2)ARIMAarimay,ar(1/#),ma(1/#)或者arimay,arima(#p,#d,#q)#p表示#阶自回归，#q表示#阶移动平均，#d表示#阶差分达到平稳过程。检验

38、残差是否存在自相关：predicte1,rescorrgrame1,lags(#)(检查残差是否存在第1到第#阶自相关)（3） ADL和ARMAXARMAX的stata命令arimayx1x2x3,ar(#)ma(#)对于时间序列，首先要设置时间项，tsettimevar对于非平稳时间序列，首先要用差分使其平稳：gdx1=d.x1(当然差分后会缺失一个变量)如果自相关系数和偏自相关系数都存在断尾，则考虑ARMA模型：理由计算信息准则estatic（目的是什么？根据信息准则判断哪一个模型更优）还可以去掉不显著的阶，回归后再运用信息准则比较。（4） VARvarsocxyz,maxlag（#）时算

39、信息准则，最大?t后期默认值为4）,根据信息准则可以估计VAR系统的阶数varxyz（进行VAR估计，默认?t后期为2）varxyz,lags（73）（滞后期为一至三期）varxyz,lags（3）（滞后期为第三期）varxyz,dfk（如果样本容量过小，可以用dfk进行自由度调整）varxyz,small（显示小样本统计量）varxyz,exog（w1,w2）（引入外生变量w1,w2）varbasicxyz,irf（估atVAR模型，画脉冲响应图，未正交化）varwle（进彳TVAR估计后，对每个方程及所有方程的各阶系数联合显著性进行wald检验，wle表示waldlag-exclusion

40、statistics）varlmar（估计后，对残差是否自相关进行LM检验）varnorm（检验残差是否服从正态分布）varstable,graph（进行VAR估计后，检查VAR系统是否为平稳过程，如果平稳则所有特征值都在单位圆内。）vargranger（估计后，进行格兰杰因果检验）irfcreatefilename,set（filename）step（#）replace（建立脉冲文件，是脉冲文件filename”成为当前的脉冲文件，step表示考察几期的脉冲响应函数，默认为8,replace代表替代已有的文件名）irfgraphirf（画脉冲响应图，未正交化）irfgraphoirf（画正交

41、化的脉冲响应图）fcastcomputeprefix,step（#）（估atVAR后，计算被解释变量的未来#期的预测值，并把预测值赋予被解释变量加上前缀“prefix”之后的变量名，如可以加F-）fcastgraphvarlist,observed（执行命令fcastcompute”后， 将变量“varlist”所代表的预测值画图，其中“observed”表示与实际观测值比较）17非平稳的时间序列有确定性趋势存在的时候，则成为不平稳，去掉时间因素，则成为平稳序列的我们叫它为趋势平稳序列；存在结构变动则为非平稳序列；存在随机趋势则也为非平稳序列。称平稳时间序列为零阶单整，如果时间序列的一阶差分为

42、平稳，则称为一阶单整，也称为单位根过程，D阶差分为平稳，则成为“d阶单整”。单位根是用来检验是否平稳的。如果时间序列存在单位根，则为非平稳序列。如果序列非平稳，则可能带来A.自回归系数的估计值向左偏向于0;B.t检验失效；C两个相互独立的变量可能出现伪回归或伪相关。进行单位根检验的方法有：（1）DF检验，使用一阶自回归来检验单位根，要求扰动项为独立白噪声，故扰动项无自相关，若有自相关，则可以引入更高阶的滞后项来控制。（2） ADF检验，多阶。在进行ADF检验时，确定滞后阶数P的大小，采用最大滞后阶数Pmax=12.（T/100）1/4,（stata命令为:di12*（t/100）A（1/4）然

43、后使用由大到小的序贯t规则，看ADF检验中最后一阶回归系数是否显著，也可以使用信息准则。ADF检验的Stata命令为：dfullery（DF检验，不包含?t后差分项），若DF统计量较小，则拒绝单位根存在假设。dfullery,lags（p）（包含p阶滞后差分项）dfullery,noconstant（不带常数项）dfullery,trend（带时间趋势项）dfullery,regress（显示回归结果）e（3） phillips-perron单位根检验（PP检验）PP检验使用一阶自回归，但使用异方差自相关稳健的标准差对DF统计量进行修正。pperrony（默认设置带常数项，不带时间趋势项）pp

44、errony,hocbnstant（不带常数项）pperrony,trend（带时间趋势项）pperrony,regress（显示回归结果）（4） DF-GLS单位根检验，是去趋势后再使用ADF检验dfglsy（默认带时间趋势项，自动根据信息准则选择最佳滞后期数）dfglsy,notrend（不带时间趋势项）（5） KPSSII位根检验，将假设改为“时间序列为平稳”。sscinstallkpss（下载安装KPS纲序）kpssy（默认包括时间趋势）kpssynotrend（不带时间趋势项）如果进行单位根检验认为是非平稳的，则要判断为1（1）或者是1（2）。协整概念和stata操作：协整检验就是为

45、了考察几个变量之间的长期均衡关系检验协整秩的检验，即确定有多少个线性无关的协整向量，（检验中打星号的是线性无关的协整向量）。stata命令：vecranky1y2yn,lags（#）（默认设置包括常数项，但不包括时间趋势项）vecranky1y2yn,lags（#）trend（none）（不包括常数项，也不包括时间趋势项）vecranky1y2yn,lags（#）trend（trend）（包括常数项，也包括时间趋势项）vecranky1y2yn,max（显示最大特征值统计量及其临界值）做完协整秩检验，并确认h=1后，就可以对VECM模型（误差修正模型）进行最大似然估计。做VECM模型的目的是？vecy1y2yn,lags（#）（默认设置包括常数项，但不包括时间趋势项）vecy1y2yn,lags（#）trend（none） （不包括常数项，也不包括时间趋势项）vecy1y2yn,lags（#）trend（trend）（包括常数项，也包括时间趋势项）后面加rank（1）表示什么？回归后表格里面的数据的含义是什么？最后一列为协整方程代表的长期均衡关系以