1、但考虑yt的条件方差:其中,表示基于t-1时刻信息集合Yt-1的yt的条件方差,若残差项存在自回归结构,则yt的条件方差不固定。假设在前p期所有信息的条件下,残差项平方服从AR(p)模型: (*)其中为0均值、方差的白噪声序列。则残差项服从条件正态分布:残差项的条件方差:由两部分组成:(1)常数项;(2)ARCH项变动信息,前p期的残差平方和注:未知参数和利用极大似然估计法估计。方差非负性要求都非负。为了使协方差平稳,需进一步要求方程的根都位于单位圆外。若都非负,上式等价于. 若扰动项的条件方差不存在自相关,则有,此时,即残差的条件方差同方差性情形。3. GARCH(p,q)模型ARCH(p)
2、模型在实际应用中,为了得到较好的拟合效果,往往需要很大的阶数p,从而增加了待估参数个数、引发多重共线性、非限制估计违背非负性要求。1986年,Bollerslev将ARCH(p)模型推广为广义自回归条件异方差模型GARCH(p, q):残差的条件方差表示为由三项组成,(2)ARCH项;(3)GARCH项前q期预测方差.未知参数用极大似然法估计,通常残差的假设分布有正态分布、t分布、广义误差分布;该模型也要求非负;若要求是平稳过程,需要限制. 实际上,GARCH(p, q)模型是将残差平方用ARMA(q,p)模型描述。4. ARCH检验检验模型的残差是否具有ARCH效应有两种方法:(1). AR
3、CH LM检验拉格朗日乘数检验检验原假设H0:残差序列直到p阶都不存在ARCH效应;需进行如下回归:检验回归有两个统计量:F统计量检验回归系数是否显著为0.TR2统计量LM统计量,其中T为观察值个数,R2为回归拟合优度,该统计量渐近服从分布。(2). 残差平方相关图残差平方相关图显示残差平方序列,直到任意指定的滞后阶数的自相关函数(AC)和偏自相关函数(PAC),并计算相应滞后阶数的QLB统计量。若不存在ARCH效应,则任意滞后阶数的自相关函数(AC)和偏自相关函数(PAC)都近似为0. 5. GARCH-M模型一般风险越大,预期收益越大。在回归模型中加入一项“利用条件方差表示的预期风险”:称
4、为GARCH-M模型。另外,还有非对称冲击模型:TARCH、EGARCH、PARCH等(略)。二. SAS实现 -PROC AutoReg过程SAS中的AutoReg过程,是用于估计和预测误差项自相关或异方差的时间序列数据的线性回归模型。自回归误差模型被用来校正自相关系数和广义自回归条件异方差模型GARCH,并且其变体如广义的ARCH(GARCH)、方差无穷的GARCH(IGARCH)、指数的GARCH(EGARCH)和依均值的GARCH(GARCH-M)被用于异方差的建模和校正。自回归过程AutoReg可以拟合任意阶的自回归误差模型,并且可以拟合子集自回归模型。为了诊断自相关性,过程产生广义
5、Durbin-Watson(DW)统计量和其边缘概率。普通回归分析假定误差方差对于所有观察是相同的,但当误差方差不相同时,数据被称为异方差,此时普通最小二乘法估计不是有效的,同时也影响预测值置信区间的精确性。Autoreg过程能检验异方差,并且提供GARCH模型族来估计和校正数据易变性。对于带有自相关扰动和随时间变化的条件异方差模型,过程输出条件均值和条件方差的预测值。基本语法:Proc AutoReg data=数据集 ;model 因变量=独立回归变量列表 ;output out=数据集 说明:1. proc autoreg语句可选项 outest=数据集把估计参数输出到指定数据集; co
6、vout把估计参数的协方差阵输出到outest=指定数据集;该选项只有在指定了outest=选项后才有效.2. model语句 center通过减去均值,中心化因变量并且取消模型的均值参数; noint取消模型的均值参数; nlag数值 | (数值列表)指定自回归误差的阶或者自回归误差的时间间隔的子集。例如,nlag=3与nlag=(1 2 3)作用相同,但与nlag=(1 3)等不同; garch=(q=数值, p=数值,type=选择值, mean, noint, tr)指定广义条件异方差GARCH模型的类型。例如,定义GARCH(2,1)回归模型:model y=x1 x2 / garc
7、h=(q=2, p=1);注意:SAS系统的自回归参数符号q和p与前文所述公式中的符号p和q正好相反。定义GARCH-M(1,1)回归模型:model y=x1 x2 / garch=(q=2,p=1,mean);type=选择值指定GARCH模型的类型:默认为noineq表示无约束GARCH模型; nonneg表示非负约束GARCH模型;stn表示约束GARCH模型系数的和小于1; integ表示IGARCH模型; exp表示EGARCH模型;noint取消条件异方差模型中的均值参数;trGARCH模型的估计使用信赖区域方法,缺省值为对偶拟牛顿法; archtest要求用portmantea
8、 Q检验统计量和Engle的拉格朗日乘子LM检验是否存在条件异方差情况,即是否有ARCH效应; coef输出前几条观察的变换系数; corrb输出参数估计的估计相关系数; covb输出参数估计的估计协方差; dw=n输出直到n阶的DW统计量,默认n=1; dwprob输出DW统计量的p值,当误差自由度大于300时dwprob选项被忽略; ginv输出Yule-Walker解的自协方差的Toeplitz矩阵的逆; itprint输出每步迭代的目标函数和参数估计; lagdetp输出DW t统计量,它用于检验存在时滞因变量时残差的自相关性; lagdep=回归变量输出DW h统计量,它用于检验一阶
9、自相关性; partial输出偏自相关; backstep去掉非显著自回归参数,参数按最小显著性的次序被去掉; slstay=数值指定被backstep选项使用的显著水平,默认为0.05; converge=数值指定在迭代自回归参数估计时参数的变化量的最大绝对值小于此数值,那么认为收敛,默认为0.001; maxiter=数值指定允许迭代的最大次数,默认为50; method=ml | ols | yw | ityw指定估计的方法,分别为:最大似然估计、无条件最小二乘法、Yule-Walker估计、迭代Yule-Walker估计; nomiss使用没有缺失值的第一个连贯时间序列数据集,进行模型
10、拟合估计。否则,跳过数据集开始的任何缺失值,使用独立回归变量和因变量都不带缺失值的所有数据。请特别注意,为了保持时间序列中正确的时间间隔,必须要增加时间刻度值,这样就会产生因变量缺失值的观察。当因变量缺失时,过程可以产生预测值。如果缺失值很多,则应使用ML估计。3. output语句 out=数据集指定包含预测值和变换值的输出数据集; alphacli=数值设置时间序列预测值置信区间的显著水平,默认为0.05; alphaclm=数值设置模型结构部分预测值置信区间的显著水平,缺省值为0.05; cev=变量把条件误差方差写入到输出数据集的指定变量中,仅GARCH模型被估计时才使用; cpev变
11、量把条件预测误差方差写入到输出数据集的指定变量中,仅GARCH模型被估计时才使用; constant变量把被变换的均值写入到输出数据集的指定变量中; lcl变量把预测值的置信下限写入到输出数据集的指定变量中; ucl变量把预测值的置信上限写入到输出数据集的指定变量中; lclm变量把模型结构部分预测值的置信下限写入到输出数据集的指定变量中; uclm变量把模型结构部分预测值的置信上限写入到输出数据集的指定变量中; p变量把预测值写入到输出数据集的指定变量中; rm变量把来自模型结构部分预测的残差写入到输出数据集的指定变量中; transform变量把被变换的变量写入到输出数据集的指定变量中。三
12、. 例子对模拟方法生成的时间趋势加二阶自回归误差模型的时间序列数据,用PROC AUTOREG过程进行分析和建模,以便于比较和判断各种求解模型和运算结果的好坏。 模型:(一)按照模拟模型生成数据集代码:data randar; e1=0; e11=0; do t=-10 to 36; e=1.3*e1-0.5*e11+2*rannor(12346); x=10+0.5*t+e; e11=e1; e1=e; if t0 then output; end;run;proc print data=randar;运行结果:Obse1e11tex1-2.73816-4.991867.761820.0367411.036730.6635612.16364-0.8623311.13775-4.301658.19846-3.889579.11047-3.566359.93378-3.6890810.31099-2.5026311.997410-1.80957-2
copyright@ 2008-2022 冰豆网网站版权所有
经营许可证编号:鄂ICP备2022015515号-1