1、语句说明:(1) DATA步中的命令 ” ,这是指令系统对变量 x进行1阶差分,差分后的序列值赋值给变量 difx。其中dif ()是差分函数,假如要差分的变量名为 x,常见的几种差分表示为:1阶差分:dif(x)2阶差分:(dif (x)k步差分:difk(2)我们在GPLOT过程中添加绘制了一个时序图“ difx*t ”,这是为了直观考察 1阶时序图显示差分后序列 difx没有明显的非平稳特征。(3)“ identify var=x ( 1); ”,使用该命令可以识别查分后序列的平稳性、纯随机性和适当的拟合模型阶数。其中 x( 1)表示识别变量x的1阶差分后序列。SAS支持多种形式的差分序
2、列识别:var=x (1),表示识别变量 x的1阶查分后序列 xt ;var=x (1,1 ),表示识别变量x的2阶查分后序列 2xt ;var=x ( k),表示识别变量 x的k步差分后序列 kxt ;var=x ( k, s),表示识别变量x的k步差分后,再进行 s步查分后序列 sA kxt。识别部分的输出结果显示 1阶查分后序列difx为平稳非白噪声序列,且具有显著的自相关系数不截尾、偏自相关系数 1截尾的性质。(4) “ estimate p=1; ”对1阶差分后序列 xt拟合AR( 1)模型。输出拟合结果显示 常数项不显著,添加或修改估计命令如下:estimate p=1 non i
3、t ;这是命令系统不要常数项拟合 AR(1)模型,拟合结果显示模型显著且参数显著。如图 5-51所示。Eat limtloaPtriHeterHUH41 .1AR i p 11. oao-4.4*4-Q, QSiAuto匚广广uiom Check ofRes IduI *图5-51序列difx模型拟合结果 xt=艾普龙t或等阶记为:xt=艾普龙t(5) “forecast lead=5 id=t ;,利用拟合模型对序列 xt作5期预测。建立数据集,绘制时序图data example5_2;input x;lagx=lag(x); cards;proc gplot data=example5_2
4、;plot x*t=1;symboll c=black i=join v=star;输出时序图如5-52所示。:图5-52 序列x时序图时序图显示,序列X有一个明显的随时间线性递增的趋势, 同时又有一定规律性的波动,所以不妨考虑使用误差自回归模型拟合该序列的发展。二、因变量关于时间的回归模型proc autoreg data =example5_2;model x=t/ dwprob;(1) “ proc autoreg data =example5_2;指令 SAS系统对临时数据集 example5_2 进行 回归程序分析。(2) “ model x=t/ dwprob; ”指令SAS系统以
5、变量t作为自变量,变量 x作为因变量,建立线性模型:xt a bt ui并给出残差序列DW检验统计量的分为点。本例中,序列x关于变量t的线性回归模型最小二乘估计输出结果如图 5-53所示。Ordi niry Leist Squares imates图5-53 序列x关于变量t的线性回归模型最小二乘估计结果本例输出结果显示,DW统计量的值等于,输出概率显示残差序列显著正相关。所以应该考虑对残差序列拟合自相关模型,修改 AUTORE程序如下:model x=t/ nlag =5 backstep method=ml;Model语句是指令系统对线性回归模型 的残差序列显示延迟 5阶的自相关图,并拟合
6、延迟5阶自相关模型,特别注意,SAS输出的自回归模型结构为:Ut 1Ut 1 5Ut 5 t即输出的自相关回归参数值与我们习惯定义的自回归参数值相差一个负号。由于自相关延迟阶数的确定是由我们尝试选择的,所以 nlag得阶数通常会指得大一些。这就导致残差自回归模型中可能有部分参数不显著,因而添加逐步回归选项 backstep,指令系统使用逐步回归的方法筛选出显著自相关因子, 并 使用极大似然的方法进行参数估计。输出如下四方面的结果:1因变量说明如图 5-54所示:The RUTOREG ProcedureDependent Uar i abIe x图5-54 因变量说明2普通最小二乘估计结果该部
7、分输出信息包括差平方和(SSE)自由度(DFE)均方误差(MSE)根号均方误差(Root MSE) SBC信息量。回归部分相关系数平方(regress rsquare)总的相关系数平方(totel rsquare),DW 统计量(durbin watson)及所有待参数的自由度、估计值、标准差、t值和统计量的P值。如图5-55所示Ord I nary Lo-ast Eat I hjSSE1IMSM EtHDFEMe$E冏KNI96EH8SBC37H.BEBBIBAIC3TI .353HT9HOQFG0B R-SqUdrQ.930*Total08300Jur bin-UtBIPr OH DU1f
8、i*hOTEi 户 rcou i* the Hue for i nia five 祁 d the p-v*me for* incgaTivb utocarrBIion.(jr IoDFE*t i nst-c乃口* CErrort UlUBApprox |f|Intereapt1-20.91ae6. I7B0i-30.00IGt3 .M9T70 .50313 .9T图5-55普通最小二乘估计结果3回归误差分析该部分共输出四方面的信息:残差序列自相关图、逐步回归消除的不显著项报告、初步均方误差(MSE)自回归参数估计值。本例该部分输出结果如图5-56所示Est i mat as of AiutoE
9、arre I at i ans:LScoMir ItncecorrelttlcnHl 9 9 7 G 5 H 3 I D1Z3M56I89HT3 .T1 .000040|333 -E0 .70047 |18,30 .H61871 II5H0,396651 OS .30 .E304B156 .90 .331 It!Bdckuard EIim i nat i on ofAutoreress i ue TermsEstI mates of fiutoregress i ve ParametersStandardLs9 Cocff i s i ent Error t Iuc1 -0 .70H007 0
10、112E91 -6.27图5-56自回归误差分析输出结果本例输出的残差序列自相关图显示残差序列有非常显著的 1阶正相关性。逐步回归消除报告显示除了 1阶的序列值显著自相关外, 延迟其他阶数的序列值均不具有显著的自相关性, 因此延迟2阶-5阶的自相关项被剔除。初步均方误差为ut 0.602573ut 1 t4、最终拟合模型该部分包括三方面的汇总信息:收敛状况、极大似然估计结果和回归系数估计。 本例该部分输出结果如图5-57所示。H&ximum Liike 1 i hood Est imates10037.807ME40MSE250.91$Roct MSE16 閑 125357.818944ilC
11、363.843805Re&ress R-Square0.8876Tciftl R-Square0J533IXirtin-Tatson1.6875NCTE:No intencept term is used. R-squar&sare redef i ned.St anda,rdftppraxVariable匚st i巾応et Vlue lt|0.29403-38c.ooaiARI-0.6S880.1124-6.12.uoaiAutoreressi ve parameters Assumed ;gien.St ardeirdEst i mat et Yalue III2.75360-23468J900(11图5-57最终拟合模型输出结果本例得到最终拟合模型为:xt 2.763& utut 0.6883ut 1 t为了得到直观的拟合效果,我们可以利用 OUTPU命令将拟合结果存入SAS数据集中,并对输出结果作图,相关命令如下:output out =out p=xp pm=trend;proc autoreg data =ex
copyright@ 2008-2022 冰豆网网站版权所有
经营许可证编号:鄂ICP备2022015515号-1