《条件异方差模型》PPT课件.ppt
《《条件异方差模型》PPT课件.ppt》由会员分享,可在线阅读,更多相关《《条件异方差模型》PPT课件.ppt(103页珍藏版)》请在三一办公上搜索。
1、1,第六章 条件异方差模型,EViews中的大多数统计工具都是用来建立随机变量的条件均值模型。本章讨论的重要工具具有与以往不同的目的建立变量的条件方差或变量波动性模型。我们想要建模并预测其变动性通常有如下几个原因:首先,我们可能要分析持有某项资产的风险;其次,预测置信区间可能是时变性的,所以可以通过建立残差方差模型得到更精确的区间;第三,如果误差的异方差是能适当控制的,我们就能得到更有效的估计。,2,6.1 自回归条件异方差模型 自回归条件异方差(Autoregressive Conditional Heteroscedasticity Model,ARCH)模型是特别用来建立条件方差模型并对
2、其进行预测的。ARCH模型是1982年由恩格尔(Engle,R.)提出,并由博勒斯莱文(Bollerslev,T.,1986)发展成为GARCH(Generalized ARCH)广义自回归条件异方差。这些模型被广泛的应用于经济学的各个领域。尤其在金融时间序列分析中。按照通常的想法,自相关的问题是时间序列数据所特有,而异方差性是横截面数据的特点。但在时间序列数据中,会不会出现异方差呢?会是怎样出现的?,3,恩格尔和克拉格(Kraft,D.,1983)在分析宏观数据时,发现这样一些现象:时间序列模型中的扰动方差稳定性比通常假设的要差。恩格尔的结论说明在分析通货膨胀模型时,大的及小的预测误差会大量
3、出现,表明存在一种异方差,其中预测误差的方差取决于后续扰动项的大小。,4,从事于股票价格、通货膨胀率、外汇汇率等金融时间序列预测的研究工作者,曾发现他们对这些变量的预测能力随时期的不同而有相当大的变化。预测的误差在某一时期里相对地小,而在某一时期里则相对地大,然后,在另一时期又是较小的。这种变异很可能由于金融市场的波动性易受谣言、政局变动、政府货币与财政政策变化等等的影响。从而说明预测误差的方差中有某种相关性。为了刻画这种相关性,恩格尔提出自回归条件异方差(ARCH)模型。ARCH的主要思想是时刻 t 的ut 的方差(=t2)依赖于时刻(t 1)的扰动项平方的大小,即依赖于 t2-1。,5,6
4、.1.1 ARCH模型 为了说得更具体,让我们回到k-变量回归模型:(6.1.1)如果 ut 的均值为零,对 yt 取基于(t-1)时刻的信息的期望,即Et-1(yt),有如下的关系:(6.1.2)由于 yt 的均值近似等于式()的估计值,所以式()也称为均值方程。,6,在这个模型中,变量 yt 的条件方差为()其中:var(yt Yt-1)表示基于(t-1)时刻的信息集合Yt-1=yt-1,yt-2,y1的 yt 的条件方差,假设在时刻(t 1)所有信息已知的条件下,扰动项 ut 的条件分布是:(6.1.7)也就是,ut 遵循以0为均值,(0+1u2t-1)为方差的正态分布。,7,由于(6.
5、1.7)中 ut 的方差依赖于前期的平方扰动项,我们称它为ARCH(1)过程:通常用极大似然估计得到参数0,1,2,k,0,1的有效估计。容易加以推广,ARCH(p)过程可以写为:(6.1.8)这时方差方程中的(p+1)个参数0,1,2,p也要和回归模型中的参数0,1,2,k一样,利用极大似然估计法进行估计。,8,如果扰动项方差中没有自相关,就会有 H0:这时 从而得到扰动项方差的同方差性情形。恩格尔曾表明,容易通过以下的回归去检验上述虚拟假设:其中,t 表示从原始回归模型()估计得到的OLS残差。,9,在 ARCH(p)过程中,由于 ut 是随机的,ut2 不可能为负,所以对于 ut 的所有
6、实现值,只有是正的,才是合理的。为使 ut2 协方差平稳,所以进一步要求相应的特征方程()的根全部位于单位圆外。如果 i(i=1,2,p)都非负,式()等价于 1+2+p 1。,10,6.1.2 ARCH的检验,下面介绍检验一个模型的残差是否含有ARCH效应的两种方法:ARCH LM检验和残差平方相关图检验。1.ARCH LM检验 Engle在1982年提出检验残差序列中是否存在ARCH效应的拉格朗日乘数检验(Lagrange multiplier test),即ARCH LM检验。自回归条件异方差性的这个特殊的设定,是由于人们发现在许多金融时间序列中,残差的大小与最近的残差值有关。ARCH本
7、身不能使标准的OLS估计无效,但是,忽略ARCH影响可能导致有效性降低。,11,ARCH LM检验统计量由一个辅助检验回归计算。为检验原假设:残差中直到q阶都没有ARCH,运行如下回归:式中 t 是残差。这是一个对常数和直到 q 阶的滞后平方残差所作的回归。这个检验回归有两个统计量:(1)F 统计量是对所有残差平方的滞后的联合显著性所作的一个省略变量检验;(2)TR2 统计量是Engles LM检验统计量,它是观测值个数 T 乘以回归检验的 R2;,12,普通回归方程的ARCH检验都是在残差检验下拉列表中进行的,需要注意的是,只有使用最小二乘法、二阶段最小二乘法和非线性最小二乘法估计的方程才有
8、此项检验。,Breusch-Pagan-GodfreyHarveyGlejserARCHWhiteCustom Test Wizard,图6.4 普通方程的ARCH检验列表,13,2.残差平方相关图 显示直到所定义的滞后阶数的残差平方t2的自相关系数和偏自相关系数,计算出相应滞后阶数的Ljung-Box统计量。残差平方相关图可以用来检查残差自回归条件异方差性(ARCH)。如果残差中不存在ARCH,在各阶滞后自相关和偏自相关系数应为0,且Q统计量应不显著。可适用于LS,TSLS,非线性LS方程。在图6.4中选择Residuals Tests/Correlogram Squared Residua
9、ls项,它是对方程进行残差平方相关图的检验。单击该命令,会弹出一个输入计算自相关和偏自相关系数的滞后阶数设定的对话框,默认的设定为36,单击OK按钮,得到检验结果。,14,例6.1 沪市股票价格指数波动的ARCH检验 为了检验股票价格指数的波动是否具有条件异方差性,本例选择了沪市股票的收盘价格指数的日数据作为样本序列,这是因为上海股票市场不仅开市早,市值高,对于各种冲击的反应较为敏感,因此,本例所分析的沪市股票价格波动具有一定代表性。在这个例子中,我们选择的样本序列sp是1996年1月1日至2006年12月31日的上海证券交易所每日股票价格收盘指数,为了减少舍入误差,在估计时,对sp进行自然对
10、数处理,即将序列ln(sp)作为因变量进行估计。,15,由于股票价格指数序列常常用一种特殊的单位根过程随机游动(Random Walk)模型描述,所以本例进行估计的基本形式为:(6.1.12)首先利用最小二乘法,估计了一个普通的回归方程,结果如下:(6.1.13)(2.35)(951)R2=0.997,16,可以看出,这个方程的统计量很显著,而且,拟合 的程度也很好。但是需要检验这个方程的误差项是否存在条件异方差性,。,17,图6.1 股票价格指数方程回归残差,观察上图,该回归方程的残差,我们可以注意到波动的“成群”现象:波动在一些较长的时间内非常小,在其他一些较长的时间内非常大,这说明残差序
11、列存在高阶ARCH效应。,18,因此,对式(6.1.26)进行条件异方差的ARCH LM检验,得到了在滞后阶数p=3时的ARCH LM检验结果如下。此处的P值为0,拒绝原假设,说明式()的残差序列存在ARCH效应。,可以计算式()的残差平方t2的自相关(AC)和偏自相关(PAC)系数,结果说明式()的残差序列存在ARCH效应。,19,例6.2 中国CPI模型的ARCH检验 本例建立CPI模型,因变量为中国的消费价格指数(上年同月=100)减去100,记为cpit;解释变量选择货币政策变量:狭义货币供应量M1的增长率,记为m1rt;3年期贷款利率,记为Rt,样本期间是1994年1月2007年12
12、月。由于是月度数据,利用X-12季节调整方法对 cpit 和 m1rt 进行了调整,结果如下:t=(19.5)(-5.17)(2.88)(-2.74)R2=0.99 对数似然值=-167.79 AIC=2.045 SC=2.12,20,这个方程的统计量很显著,拟合的程度也很好。但是观察该回归方程的残差图,也可以注意到波动的“成群”现象:波动在一些时期内较小,在其他一些时期内较大,这说明误差项可能具有条件异方差性。,21,从自相关系数和偏自相关系数可以看出:残差序列存在着一阶ARCH效应。再进行条件异方差的ARCH LM检验,得到了在滞后阶数p=1时的ARCH LM检验结果:,因此计算残差平方t
13、2的自相关(AC)和偏自相关(PAC)系数,结果如下:,22,从自相关系数和偏自相关系数可以看出:残差序列存在着一阶ARCH效应。因此利用ARCH(1)模型重新估计模型(),结果如下:均值方程:z=(12.53)(-1.53)(4.72)(-3.85)方差方程:z=(5.03)(3.214)R2=0.99 对数似然值=-151.13 AIC=1.87 SC=1.98 方差方程中的ARCH项的系数是统计显著的,并且对数似然值有所增加,同时AIC和SC值都变小了,这说明ARCH(1)模型能够更好的拟合数据。,23,再对这个方程进行条件异方差的ARCH LM检验,得到了残差序列在滞后阶数p=1时的统
14、计结果:此时的相伴概率为0.69,接受原假设,认为该残差序列不存在ARCH效应,说明利用ARCH(1)模型消除了式()的残差序列的条件异方差性。式()的残差平方相关图的检验结果为:自相关系数和偏自相关系数近似为0。这个结果也说明了残差序列不再存在ARCH效应。,24,6.1.3 GARCH模型 扰动项 ut 的方差常常依赖于很多时刻之前的变化量(特别是在金融领域,采用日数据或周数据的应用更是如此)。因此 必须估计很多参数,而这一点很难精确的做到。但是如果我们能够意识到方程(6.1.8)不过是 t2 的分布滞后模型,我们就能够用一个或两个 t2 的滞后值代替许多 ut2 的滞后值,这就是广义自回
15、归条件异方差模型(generalized autoregressive conditional heteroscedasticity model,简记为GARCH模型)。在GARCH模型中,要考虑两个不同的设定:一个是条件均值,另一个是条件方差。,25,在标准化的GARCH(1,1)模型中:均值方程:(6.1.17)方差方程:(6.1.18)其中:xt 是(k+1)1维外生变量向量,是(k+1)1维系数向量。(6.1.17)中给出的均值方程是一个带有扰动项的外生变量函数。由于t2是以前面信息为基础的一期向前预测方差,所以它被称作条件方差,式()也被称作条件方差方程。,26,(6.1.18)中给
16、出的条件方差方程是下面三项的函数:1常数项(均值):2用均值方程(6.1.11)的扰动项平方的滞后来度量从前期得到的波动性的信息:ut2-1(ARCH项)。3上一期的预测方差:t2-1(GARCH项)。GARCH(1,1)模型中的(1,1)是指阶数为1的GARCH项(括号中的第一项)和阶数为1的ARCH项(括号中的第二项)。一个普通的ARCH模型是GARCH模型的一个特例,GARCH(0,1),即在条件方差方程中不存在滞后预测方差t2-1的说明。,27,在EViews中ARCH模型是在扰动项是条件正态分布的假定下,通过极大似然函数方法估计的。例如,对于GARCH(1,1),t 时期的对数似然函
17、数为:(6.1.19)其中(6.1.20)这个说明通常可以在金融领域得到解释,因为代理商或贸易商可以通过建立长期均值的加权平均(常数),上期的预期方差(GARCH项)和在以前各期中观测到的关于变动性的信息(ARCH项)来预测本期的方差。如果上升或下降的资产收益出乎意料地大,那么贸易商将会增加对下期方差的预期。这个模型还包括了经常可以在财务收益数据中看到的变动组,在这些数据中,收益的巨大变化可能伴随着更进一步的巨大变化。,28,有两个可供选择的方差方程的描述可以帮助解释这个模型:1如果我们用条件方差的滞后递归地替代()式的右端,就可以将条件方差表示为滞后扰动项平方的加权平均:()我们看到GARC
18、H(1,1)方差说明与样本方差类似,但是,它包含了在更大滞后阶数上的,扰动项的加权条件方差。,29,2设 vt=ut2 t2。用其替代方差方程()中的方差并整理,得到关于扰动项平方的模型:()因此,扰动项平方服从一个异方差ARMA(1,1)过程。决定波动冲击持久性的自回归的根是 加 的和。在很多情况下,这个根非常接近1,所以冲击会逐渐减弱。,30,方差方程的回归因子 方程(6.1.18)可以扩展成包含外生的或前定回归因子 z 的方差方程:()注意到从这个模型中得到的预测方差不能保证是正的。可以引入到这样一些形式的回归算子,它们总是正的,从而将产生负的预测值的可能性降到最小。例如,我们可以要求:
19、,31,高阶GARCH(p,q)模型 高阶GARCH模型可以通过选择大于1的 p 或 q 得到估计,记作GARCH(q,p)。其方差表示为:()这里,q 是GARCH项的阶数,p是ARCH项的阶数,p0并且,(L)和(L)是滞后算子多项式。,32,为了使GARCH(q,p)模型的条件方差有明确的定义,相应的ARCH()模型()的所有系数都必须是正数。只要(L)和(L)没有相同的根并且(L)的根全部位于单位圆外,那么当且仅当0=0/(1-(L),(L)=(L)/(1-(L)的所有系数都非负时,这个正数限定条件才会满足。例如,对于GARCH(1,1)模型()这些条件要求所有的3个参数都是非负数。,
20、33,6.1.4 IGARCH模型 如果限定GARCH模型的方差方程中的参数和等于1,并且去掉常数项:()其中()这就是Engle和Bollerslev(1986)首先提出的单整GARCH模型(Intergrated GARCH Model,IGARCH)。,34,6.1.5 约束及回推 1约束 在估计一个GARCH模型时,有两种方式对GARCH模型的参数进行约束(restrictions)。一个选择是IGARCH方法,它将模型的方差方程中的所有参数之和限定为1。另一个就是方差目标(variance target)方法,它把方差方程()中的常数项设定为GARCH模型的参数和无条件方差的方程:(
21、)这里的 是残差的无条件方差。,35,2回推 在计算GARCH模型的回推初始方差时,首先用系数值来计算均值方程中的残差,然后计算初始值的指数平滑算子()其中:t 是来自均值方程的残差,是无条件方差的估计:()平滑参数 为0.1至1之间的数值。也可以使用无条件方差来初始化GARCH过程:(),36,6.1.6 GARCH模型的残差分布假设 在实践中我们注意到,许多时间序列,特别是金融时间序列的无条件分布往往具有比正态分布更宽的尾部。为了更精确地描述这些时间序列分布的尾部特征,还需要对误差项ut的分布进行假设。GARCH模型中的扰动项的分布,一般会有3个假设:正态(高斯)分布、学生t-分布和广义误
22、差分布(GED)。给定一个分布假设,GARCH模型常常使用极大似然估计法进行估计。下面分别介绍这3种分布,其中的 代表参数向量。1对于扰动项服从正态分布的GARCH(1,1)模型,它的对数似然函数为()这里的t2是ut的条件方差。,37,2如果扰动项服从学生t分布,GARCH(1,1)模型的对数似然函数的形式就是()这样,参数的估计就变成了在自由度k2的约束下使对数似然函数()最大化的问题。当k时,学生t-分布接近于正态分布。注 式()和()中的()代表 函数:若N是偶整数,则(N/2)=12 3(N/2)-1,有(2/2)=1;若N是奇整数,则,有。,38,3扰动项的分布为广义误差分布(GE
23、D)时,GARCH(1,1)模型的对数似然函数的形式为()这里的参数r 0。如果r=2,那么GED就是一个正态分布。,39,ARCH-M模型 金融理论表明具有较高可观测到风险的资产可以获得更高的平均收益,其原因在于人们一般认为金融资产的收益应当与其风险成正比,风险越大,预期的收益就越高。这种利用条件方差表示预期风险的模型被称为ARCH均值模型(ARCH-in-mean)或ARCH-M回归模型。在ARCH-M中我们把条件方差引进到均值方程中:()ARCH-M模型的另一种不同形式是将条件方差换成条件标准差:()或取对数(),40,ARCH-M模型通常用于关于资产的预期收益与预期风险紧密相关的金融领
24、域。预期风险的估计系数是风险收益交易的度量。例如,我们可以认为某股票指数,如上证的股票指数的收益率(returet)依赖于一个常数项及条件方差(风险):这种类型的模型(其中期望风险用条件方差表示)就称为GARCH-M模型。,41,在EViews中估计ARCH模型,估计GARCH和ARCH模型,首先选择Object/New Object/Equation,然后在Method的下拉菜单中选择ARCH,得到如下的对话框。,图6.5 ARCH模型定义对话框,42,与选择估计方法和样本一样,需要指定均值方程和方差方程。一、均值方程(Mean equation)在因变量编辑栏中输入均值方程形式,均值方程的
25、形式可以用回归列表形式列出因变量及解释变量。如果方程包含常数,可在列表中加入C。如果需要一个更复杂的均值方程,可以用公式的形式输入均值方程。,43,如果解释变量的表达式中含有ARCHM项,就需要点击对话框右上方对应的按钮。EViews中的ARCH-M的下拉框中,有4个选项:1.选项None表示方程中不含有ARCHM项;2.选项Std.Dev.表示在方程中加入条件标准差;3.选项Variance则表示在方程中含有条件方差 2。4.选项Log(Var),表示在均值方程中加入条件方差的对数ln(2)作为解释变量。,44,二、方差设定和分布设定(Variance and distribution sp
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 条件异方差模型 条件 方差 模型 PPT 课件
链接地址:https://www.31ppt.com/p-5584603.html