stata操作介绍之时间序列-(四)

四、平稳时间序列分析2021/6/271平滑分析自相关分析ARIMA模型2021/6/272平滑分析

大多数时间序列数据都具有上下起伏的波动，对于时间序列数据的识别十分困难。平滑分析可以把数据分为两个部分：一部分逐渐发生变化，便于分析处理；另一部分则含有突变的成份。时间平滑：用某时刻及其前后若干时刻的值进行加权平均，所得值作为该时刻的替代值以滤去小扰动的方法。2021/6/273平滑分析

该部分以1983年1月到7月Milford城镇的自来水消费量为例。文件：ch51.dta2021/6/274导入数据：.use"C:\Users\Administrator\Desktop\时间序列\数据\ch51.dta",clear2021/6/275平滑分析描述性统计：.describe2021/6/276平滑分析生成日期变量（一）：.gendate=mdy(month,day,year).listin1/52021/6/277平滑分析设置时间（二）：.tssetdate,format(%d).listin1/52021/6/278平滑分析趋势图：.graphtwolinewaterdate,ylabel(300(100)900)2021/6/279平滑分析生成移动平均值（1）：.genwater3=(water[_n-1]+water[_n]+water[_n+1])/32021/6/2710平滑分析生成移动平均值（2）：.tssmoothmawater5=water,window(212)注：tssmooth:表示移动平均值平滑(加权或不加权);window(212):表示使用该值的前两个值、该值与该值的后两个值进行平均计算；2021/6/2711平滑分析趋势图：.graphtwoline

water5date,clwidth(thick)||linewaterdate,clwidth(thin)clpattern(solid)2021/6/2712平滑分析波动幅度：.gench=water-water5.listin1/52021/6/2713平滑分析波动幅度：.graphtwolinechdate2021/6/2714平滑分析-滞后变量导入数据：.use"C:\Users\Administrator\Desktop\时间序列\数据\ch52.dta",clear描述性统计：.describe时间设置：.tssetyear,yearly2021/6/2715平滑分析-滞后变量时间平滑：.tssmoothmafylln=fylltemp,window(212)2021/6/2716平滑分析-滞后变量趋势图：.graphtwospikefylltempyear,base(1.67)yline(1.67)||linefyllnyear,clpattern(solid)2021/6/2717平滑分析-滞后变量生成n阶滞后变量的两种方法：.genwNAO_n=wNAO[_n-1].genwNAO_n=Ln.wNAO

注：第二种方法中的''Ln''表示Lag(n)；2021/6/2718平滑分析-滞后变量生成一阶滞后变量：.genwNAO_1=wNAO[_n-1]生成二阶滞后变量：.genwNAO_2=L2.wNAO2021/6/2719平滑分析-滞后变量2021/6/2720平滑分析-滞后变量三种滞后回归方法:.regfylltemp

wNAOwNAO_1wNAO_2ifyear>=1970&year<=1990.regfylltempwNAOL.wNAOL2.wNAOifyear>=1970&year<=1990.regfylltempL(0/2).wNAOifyear>=1970&year<=19902021/6/2721平滑分析-滞后变量

回归结果：2021/6/2722相关分析

自相关系数是对变量自身与其滞后变量之间相关关系的估计。偏相关系数是在消除其他变量影响的条件下，所计算的某两变量之间的相关系数。交叉相关是分析两个时间序列之间的关系。注：本部分继续使用ch52.dta数据。2021/6/2723自相关分析自相关表：.corrgramfylltemp,lags(9)说明:大多数P值都小于0.05，故认为fylltemp具有显著的自相关性；相关关系或偏相关关系越强，相应的线条越长；2021/6/2724自相关分析自相关图：.acfylltemp,lags(9)注：阴影部分是95%的置信区间；2021/6/2725自相关分析偏相关图：.pacfylltemp,lags(9)yline(0)ciopts(bstyle(outline))注：ciopts(bstyle(outline))表示将偏相关图中阴影部分改为矩阵区域。2021/6/2726自相关分析交叉相关图:.xcorrwNAOfylltemp

ifyear>=1970&year<=1990,lags(7)xlabel(-9(1)9,grid)说明：由图可知，lag为0时，交叉相关性最强（线条最长），且为负。2021/6/2727自相关分析交叉相关表:.xcorrwNAOfylltempifyear>=1970&year<=1990,lags(7)table2021/6/2728ARIMA模型

时间序列中的自相关集成移动平均模型（autoregressiveintegratedmovingaverage简称ARIMA），是指将非平稳时间序列转化为平稳时间序列，然后将因变量仅对它的滞后值以及随机误差项的现值和滞后值进行回归所建立的模型。ARIMA模型根据原序列是否平稳以及回归中所含部分的不同，包括移动平均过程（MA）、自回归过程（AR）、自回归移动平均过程（ARMA）以及ARIMA过程。2021/6/2729ARIMA模型ARIMA模型操作步骤：1、对变量进行检验，若变量具有稳定性，进行第二步，否则，不能使用ARIMA模型；2、做出变量的自相关图；3、根据变量的自相关图，选择合适的模型；4、根据选定的模型进行分析，并检验系数是否显著。若有的系数不显著，所选择的模型可能存在问题；若所有系数都显著，进行第五步；5、检验残差是否具有自相关性。若残差具有自相关性，则所选择的模型存在问题；若残差不具有自相关性，则所选择的模型是合适的。2021/6/2730ARIMA模型三种单位根检验方法：pperron检验：.pperronfylltemp,lag(3)说明：P值为0.0003，小于0.05，故拒绝原假设，说明该变量不满足稳定性检验；2021/6/2731ARIMA模型dfuller检验：.dfullerfylltemp,lag(3)说明：P值为0.089，大于0.05，故不能拒绝原假设，说明该变量满足稳定性检验；2021/6/2732ARIMA模型dfgls检验：.dfglsfylltemp,maxlag(3)notrend说明：P值为0.0，小于0.05，故拒绝原假设，说明该变量不满足稳定性检验；2021/6/2733

虽然pperron检验和dfgls检验拒绝了变量fylltemp具有稳定性的假设，但是dfuller检验不能拒绝原假设，还是可以认为该变量具有稳定性。2021/6/2734ARIMA模型自相关表：.corrgramfylltemp,lags(4)说明：该变量的自相关关系随着滞后期的增加而减少，偏自相关关系在一期自后滞后消失，故适合模型AR(1)来分析该变量。2021/6/2735ARIMA模型AR(1):.arimafylltemp,arima(1,0,0)nolog说明：由上图可知，常数项与一期滞后变量系数都是统计显著的，卡方检验也显著。2021/6/2736故上图可得到模型为：在得出结论之前，还需对残差进行检验，检验残差是否存在自