异方差时间序列模型.ppt

《异方差时间序列模型.ppt》由会员分享，可在线阅读，更多相关《异方差时间序列模型.ppt（49页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、关于异方差时间序列模型现在学习的是第1页，共49页Contents第一节第一节 问题的提出问题的提出第二节第二节 ARCH模型模型第三节第三节 GARCH模型模型第四节第四节 其他其他GARCH模型模型现在学习的是第2页，共49页第一节 问题的提出 在自回归移动平均模型中，我们主要讨论平稳时间序列的建模问题，由于针对平稳序列，实际上假定任一时点的随机误差项的期望值是相同的，一般为0，同时假定任一随机误差项平方的期望值就是随机误差的方差，即同方差。但是在金融市场上，金融资产报酬序列具有这样的特性，大的报酬紧连着大的报酬，小的报酬紧连着小的报酬，称为波动集群性(Mandelbrot,1963、Fa

2、ma,1965)。波动集群性表明股票报酬波动是时变的，表明是异方差。异方差虽然不会影响回归系数的最小二乘估计的无偏性，但是将影响到回归系数估计的标准差和置信区间。例如，汇率，股票价格常常用随机游走过程描述，现在学习的是第3页，共49页 其中ut为白噪声过程。1995-2000年日元兑美元汇率时间序列及差分序列见图1和图2。tttuxx1图1 日元兑美元汇率序列JPY(1995-2000)图2日元兑美元汇率差分序列(收益)D(JPY)现在学习的是第4页，共49页图3 收益绝对值序列(1995-2000)图4 D(JPY)的平方(1995-2000)这种序列的特征是（1）过程的方差不仅随时间变化，

3、而且有时变化得很激烈。（2）按时间观察，表现出“波动集群”（volatility clustering）特征，即方差在一定时段中比较小，而在另一时段中比较大。（3）从取值的分布看表现的则是“高峰厚尾”（leptokurtosis and fat-tail）特征，即均值附近与尾区的概率值比正态分布大，而其余区域的概率比正态分布小。图5给出高峰厚尾分布示意图。图6给出一个高峰厚尾分布实例。现在学习的是第5页，共49页现在学习的是第6页，共49页显然现期方差与前期的“波动”有关系。描述这类关系的模型称为自回归条件异方差（自回归条件异方差（ARCH）模型）模型（Engle 1982年提出）。使用ARC

4、H模型的理由是：（1）通过预测xt或ut的变化量评估股票的持有或交易对收益所带来的风险有多大，以及决策的代价有多大；（2）可以预测xt的置信区间，它是随时间变化的；（3）对条件异方差进行正确估计后可以使回归参数的估计量更具有有效性。现在学习的是第7页，共49页为了说得更具体，让我们回到k-变量回归模型：如果ut的均值为零，对yt取基于(t-1)时刻的信息的期望，即Et-1(yt)，有如下的关系：由于yt的均值近似等于式（1）的估计值，所以式（1）也称为均值方程。ttkkttuxxy110ktkttttxxxy221101)(E(1)(2)现在学习的是第8页，共49页假设在时刻(t 1)所有信息

5、已知的条件下，扰动项ut的条件分布是：也就是，ut遵循以0为均值，(0+1u2t-1)为方差的正态分布。由于(3)中ut的方差依赖于前期的平方扰动项，我们称它为ARCH(1)过程：)(,02110ttuNu(3)21102)var(tttuu通常用极大似然估计得到参数0,1,2,k,0,1的有效估计。现在学习的是第9页，共49页第二节 ARCH模型一、一、ARCH模型的定义模型的定义若一个平稳随机变量xt可以表示为AR(p)形式，其随机误差项的方差可用误差项平方的q阶分布滞后模型描述，222221102222110.)(.qtqtttttptptttuuuuEuxxxx（1）（2）则称ut 服

6、从q阶的ARCH过程，记作utARCH(q)。其中(1)式称作均值方程，(2)式称作ARCH方程。现在学习的是第10页，共49页(1)和(2)式还应满足如下条件。对于(1)式，为保证平稳性，特征方程0.1221ppLLL的根应在单位圆之外。xt 的条件期望是ptpttptttxxxxxxE.),.,(221101xt 的无条件期望（T 时）是ptxE101)(对于(2)式，由于 的非负性，对i应有如下约束，2tuqii,.,2,1,0,00当全部i=0,i=1,2,q时，条件方差 。因为方差是非负的，所以要求0 0。为保证 是一个平稳过程，(2)式的特征方程02t2t现在学习的是第11页，共4

7、9页0.1221qqLLL的根应在单位圆之外。对i,i=1,2,q的另一个约束是对（2）式求期望的：则无条件方差为：可见若保证 是一个平稳过程，应该有约束0 (1+2+q)1。因为Var(xt)=Var(ut)=，所以上式可以用来预测xt的方差。1.021q22222110222221102.)(.)()(qtqttqtqtttuEuEuE当T时：2222102.qt01211qii2t2t现在学习的是第12页，共49页 二、二、ARCH模型的极大似然估计模型的极大似然估计ARCH模型经常应用在回归模型中。其中=(0 1,k-1),xt=(1 x1,xk-1)（xt的分量也可以包括yt的滞后变

8、量），utARCH(q)。为计算方便，假定已知yt,xt的T+q组观测值，从而保证估计参数所用的样本容量为T。utARCH(q)可以表示为：ttuxyttttvhu 其中vtIID(0,1)，vt与xt相互独立，所以有 ，E(ut)=0。yt服从正态分布，概率密度函数为22222110.qtqtttuuuhttthuE)(22现在学习的是第13页，共49页)2)(exp(21),(2ttttitthyhxyf x其中：2222221110)(.)()(qtqtqtttttxyxyxyh用参数和=(0 1 2 q)组成参数向量，)(对数似然函数是:TttttTttTtthyhTyfrL1211)

9、(21)log(21)2log(2)(log)(log x求的极大似然估计量就是求 使 logL()在=处获得极大值。求log L()对 的偏导数，现在学习的是第14页，共49页TttttttttTtttttttttTttttttttTttyhhhhuhhuyhhhhhyyhhrL1222122212221)(121)(1121)()(121log21)(log xxxx22111221121201201200020000000010)(qttjtjtqjjjtjtqjjqttqjjtjqjjtjqjjtjtuuuxuxuuuuuh 现在学习的是第15页，共49页022112)(1221q1j

10、j222ttTttqttj-tj-ttttttuxhuuhhuxy0 x-ttuxu在上式为零条件下求到的 即是 的极大似然估计量。具有一致性。现在学习的是第16页，共49页 三、三、ARCH模型检验模型检验在均值方程（回归模型或时间序列模型）的误差项中是否存在自回归条件异方差应该进行假设检验。检验ARCH可以使用F、LM、LR、W统计量。下面介绍F、LM检验。1、自回归条件异方差的LM检验。建立原假设H0：1=2=q=0 （不存在ARCH）H1：1,2,q 不全为零在原假设成立条件下，OLS估计量是一致的、有效的；在备择假设成立条件下，OLS估计量是一致的，但不是有效的。先介绍使用LM统计量

11、检验H0。因为计算LM统计量的值，只需估计原假设成立条件下的方程。具体步骤是：现在学习的是第17页，共49页 估计 ，求 ，计算 。估计辅助回归式ttuxyttu 2tutqtqtttvuuuu222221102.用第3步得到的可决系数R2构成统计量LM=T R2。其中T表示辅助回归式的样本容量。在原假设成立条件下有)(22qTRLM若LM ，接受H1。)(2q)(2q注意：注意：辅助回归式中要有常数项0。现在学习的是第18页，共49页 2、自回归条件异方差的F检验。建立原假设 H0：1=2=q=0（不存在ARCH）H1：1,2,q不全为零 估计yt=xt+ut，求 ，计算 。用 估计2个辅助

12、回归式 构造F统计量，在原假设成立条件下有 注意：注意：辅助回归式中要有常数项0。tu 2tu2tuttttttvuuuvu.2222110202)1,()1/(/)(qTquurFqTSSEqSSESSEF 若F F,(q,T-q-1)，接受H1。（约束模型，同方差）（约束模型，同方差）（非约束模型，存在（非约束模型，存在ARCH）现在学习的是第19页，共49页3、自回归条件异方差的LR检验。建立原假设 H0：1=2=q=0（不存在ARCH）H1：1,2,q不全为零 估计yt=xt+ut，求 ，计算 。用 估计2个辅助回归式，并计算极大似然函数值LogLr和LogLu 用LogLr和LogL

13、u构造LR统计量，在原假设成立条件下有 tu 2tu2tuttttttvuuuvu.2222110202)()log(log22mLLLRur（约束模型，同方差）（约束模型，同方差）（非约束模型，存在（非约束模型，存在ARCH）若 ，接受H0。若 ，接受H1。)(2mLR)(2mLR如果结论是应该建立ARCH模型，则进一步应该对ARCH模型的阶数q进行检验。对此可采用t检验。现在学习的是第20页，共49页4、自回归条件异方差的Q检验。残差平方意味着方差，若存在自相关，说明存在自回归条件异方差。现在学习的是第21页，共49页四、四、ARCH模型检验（模型检验（EViews操作案例）操作案例）19

14、95.1-2000.8日元兑美元汇率值（1427个）序列（JPY）见图。极小值为81.12日元，极大值为147.14日元。其均值为112.93日元，标准差是13.3日元。1995年4月曾一度达到81.12日元兑1美元。1998年8月达到147.14日元兑1美元。JPY显然是一个非平稳序列。JPY的差分序列D(JPY)表示收益，见图9.2。因为D(JPY)是平稳序列，用D(JPY)建立时间序列模型。现在学习的是第22页，共49页通过相关图和偏自相关图分析，应该建立AR(3)或MA(3)模型。建立AR(3)模型如下：ttttuDJPYDJPYDJPY0859.00541.032(2.0)(-3.3

15、)R2=0.01，DW=1.91，Q(15)=8.6现在学习的是第23页，共49页 方法方法1：通过Q检验考察AR(3)模型中是否存在自回归条件异方差。现在学习的是第24页，共49页 方法方法2：ARCH的LM检验。在均值方程估计窗口，选在均值方程估计窗口，选ARCH的的LM检验。用检验。用1阶检验式检验。阶检验式检验。2122535.0685.0ttuu(9.4)(9.9)R2=0.0643,T=14234.910643.0142159.97)21423/(355.95191/)355.95191.10173()/(/)(2TRKTSSEmSSESSEFuur现在学习的是第25页，共49页

16、方法方法3、4：自回归条件异方差的F检验和LR检验。用参差平方序列用参差平方序列1阶自回归检验式做参数约束的阶自回归检验式做参数约束的F检验和检验和LR检验检验514.91)399.3371(656.3418(2),(log),(log28.3)21422,1(59.97)21423/(355.95191/)355.95191.10173()/(/)(2205.0LLLRFKTSSEmSSESSEFuur现在学习的是第26页，共49页 本例建立CPI模型，因变量为中国的消费价格指数（上年同月=100）减去100，记为cpit；解释变量选择货币政策变量：狭义货币供应量M1的增长率，记为m1rt；

17、3年期贷款利率，记为Rt，样本期间是1994年1月2007年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 ttttttuRrmcpicpicpi06.0168.236.035.12121现在学习的是第27页，共49页 这个方程的统计量很显著，拟合的程度也很好。但是观察该回归方程的残差图，也可以注意到波动的“成群”现象：波动在一些时期内较小，在其他一些时期内较大，这说明误差项可能具有条件异方差性。现在学习的是第

18、28页，共49页 从自相关系数和偏自相关系数可以看出：残差序列存在着一阶ARCH效应。再进行条件异方差的ARCH LM检验，得到了在滞后阶数p=1时的ARCH LM检验结果：因此计算残差平方t2的自相关（AC）和偏自相关（PAC）系数，结果如下：现在学习的是第29页，共49页 从自相关系数和偏自相关系数可以看出：残差序列存在着一阶从自相关系数和偏自相关系数可以看出：残差序列存在着一阶ARCH效应。因此利用效应。因此利用ARCH(1)模型重新估计模型模型重新估计模型，结果如下：，结果如下：均值方程：均值方程：z=(12.53)(-1.53)(4.72)(-3.85)方差方程：方差方程：z=(5.

19、03)(3.214)R2=0.99 对数似然值对数似然值=-151.13 AIC=1.87 SC=1.98 方差方程中的方差方程中的ARCH项的系数是统计显著的，并且对数似然值有所增项的系数是统计显著的，并且对数似然值有所增加，同时加，同时AIC和和SC值都变小了，这说明值都变小了，这说明ARCH(1)模型能够更好的拟合数据模型能够更好的拟合数据。ttttttuRrmcpicpicpi062.01098.313.0088.12121212648.0186.0ttu现在学习的是第30页，共49页 再对这个方程进行条件异方差的再对这个方程进行条件异方差的ARCH LM检验，得到了残差序列检验，得到

20、了残差序列在滞后阶数在滞后阶数p=1时的统计结果：时的统计结果：此时的相伴概率为此时的相伴概率为0.69，接受原假设，认为该残差序列不存在，接受原假设，认为该残差序列不存在ARCH效应，说明利用效应，说明利用ARCH(1)模型消除了残差序列的条件异方差模型消除了残差序列的条件异方差性。残差平方相关图的检验结果为：性。残差平方相关图的检验结果为：自相关系数和偏自相关系数近似为自相关系数和偏自相关系数近似为0。这个结果也说明了残差。这个结果也说明了残差序列不再存在序列不再存在ARCH效应。效应。现在学习的是第31页，共49页第三节第三节 GARCH模型模型一、一、GARCH模型定义模型定义 扰动项

21、ut的方差常常依赖于很多时刻之前的变化量（特别是在金融领域，采用日数据或周数据的应用更是如此）。因此 必须估计很多参数，而这一点很难精确的做到。但是如果我们能够意识到方差方程不过是t2的分布滞后模型，我们就能够用一个或两个t2的滞后值代替许多ut2的滞后值，这就是广义自回归条件异方差模型(Generalized AutoRegressive Conditional Heteroscedasticity model，简记为GARCH模型)。在GARCH模型中，要考虑两个不同的设定：一个是条件均值，另一个是条件方差。222221102ptptttuuu现在学习的是第32页，共49页 在标准化的GA

22、RCH(1,1)模型中：均值方程：方差方程：其中：xt是(k+1)1维外生变量向量，是(k+1)1维系数向量。均值方程是一个带有扰动项的外生变量函数。由于t2是以前面信息为基础的一期向前预测方差，所以它被称作条件方差，也被称作条件方差方程。tttuyx21212tttu现在学习的是第33页，共49页条件方差方程是下面三项的函数：1常数项（均值）：2用均值方程的扰动项平方的滞后来度量从前期得到的波动性的信息：ut-12（ARCH项）。3上一期的预测方差：t2-1（GARCH项）。GARCH(1,1)模型中的(1,1)是指阶数为1的GARCH项（括号中的第一项）和阶数为1的ARCH项（括号中的第二

23、项）。一个普通的ARCH模型是GARCH模型的一个特例，GARCH(0,1)，即在条件方差方程中不存在滞后预测方差t2-1的说明。现在学习的是第34页，共49页 在EViews中ARCH模型是在扰动项是条件正态分布的假定下，通过极大似然函数方法估计的。例如，对于GARCH(1,1)，t 时期的对数似然函数为：其中 这个说明通常可以在金融领域得到解释，因为代理商或贸易商可以通过建立长期均值的加权平均（常数），上期的预期方差（GARCH项）和在以前各期中观测到的关于变动性的信息（ARCH项）来预测本期的方差。如果上升或下降的资产收益出乎意料地大，那么贸易商将会增加对下期方差的预期。这个模型还包括了

24、经常可以在财务收益数据中看到的变动组，在这些数据中，收益的巨大变化可能伴随着更进一步的巨大变化。222/)(21ln21)2ln(21tttttylx2121212112)(ttttttuyx现在学习的是第35页，共49页 有两个可供选择的方差方程的描述可以帮助解释这个模型：1如果我们用条件方差的滞后递归地替代方差方程的右端，就可以将条件方差表示为滞后扰动项平方的加权平均：我们看到GARCH(1,1)方差说明与样本方差类似，但是，它包含了在更大滞后阶数上的，扰动项的加权条件方差。.12112jtjjtu现在学习的是第36页，共49页 2设vt=ut2t2。用其替代方差方程中的方差并整理，得到关

25、于扰动项平方的模型：因此，扰动项平方服从一个异方差ARMA(1,1)过程。决定波动冲击持久性的自回归的根是加的和。在很多情况下，这个根非常接近1，所以冲击会逐渐减弱。.1212ttttvvuu现在学习的是第37页，共49页 高阶高阶GARCH(p,q)模型模型 高阶GARCH模型可以通过选择大于1的p或q得到估计，记作GARCH(p,q)。其方差表示为：这里，p是GARCH项的阶数，q是ARCH项的阶数，p0并且,(L)和(L)是滞后算子多项式。22012122)()(ttpiitiqjjtjtLuLu现在学习的是第38页，共49页 为了使GARCH(q,p)模型的条件方差有明确的定义，相应的

26、ARCH()模型的所有系数都必须是正数。只要(L)和(L)没有相同的根并且(L)的根全部位于单位圆外，那么当且仅当0=0/(1-(L)，(L)=(L)/(1-(L)的所有系数都非负时，这个正数限定条件才会满足。例如，对于GARCH(1,1)模型这些条件要求所有的3个参数都是非负数。202)(ttuL21212tttu现在学习的是第39页，共49页 GARCH模型的检验模型的检验当原假设H0是ARCH(0)时，显然备择假设H1有两个。一个是一个是ARCH(r)，一个是，一个是GARCH(r,0)。若原假设 H0是ARCH(1)，则备择假设H1可以是ARCH(1+r)，也可以是GARCH(r,1)

27、。同理若原假设H0是ARCH(q)，备择假设H1可以有两个。ARCH(q+r)和GARCH(r,q)。LM统计量无法区别这两个备择假设。但仍然可以做LM检验。在实际应用中，备择假设既可以是ARCH，也可以是GARCH。对于q值很大的ARCH模型，建议使用GARCH模型。在实际应用中，GARCH(1,1)和GARCH(2,1)一般足以满足对自回归条件异方差的描述。现在学习的是第40页，共49页第四节第四节 其他其他GARCH模型模型一、一、IGARCH模型模型 对于ARCH(p)模型和GARCH(p,q)模型，在实际应用中，条件 有时不能得到满足。下面以GARCH(1,1)模型为例进行讨论。10

28、10111piiqiiqii（保证可以转换成无限阶的ARCH过程）（保证GARCH过程平稳）现在学习的是第41页，共49页用 分别表示对1,1的估计。有时会出现 21121102ttu、1甚至，。例如Engle-Chowdury（1992）对IBM收益率序列估计时，得如下结果，1212102953.0053.000056.0tttttuuy现在学习的是第42页，共49页 在计算在计算GARCH模型的回推初始方差时，首先用系数值来计模型的回推初始方差时，首先用系数值来计算均值方程中的残差，然后计算初始值的指数平滑算子算均值方程中的残差，然后计算初始值的指数平滑算子 （6.1.30）其中：是来自均

29、值方程的残差，是无条件方差的估计：其中：是来自均值方程的残差，是无条件方差的估计：（6.1.31）平滑参数平滑参数为为0.1至至1之间的数值。也可以使用无条件方差来初始之间的数值。也可以使用无条件方差来初始化化GARCH过程：过程：（6.1.32）)()1(02122020TjjTjTTuuTttuT122122020 u现在学习的是第43页，共49页 在实践中我们注意到，许多时间序列，特别是金融时间序列的无在实践中我们注意到，许多时间序列，特别是金融时间序列的无条件分布往往具有比正态分布更宽的尾部。为了更精确地描述这些时条件分布往往具有比正态分布更宽的尾部。为了更精确地描述这些时间序列分布的

30、尾部特征，还需要对误差项间序列分布的尾部特征，还需要对误差项ut的分布进行假设。的分布进行假设。GARCH模型中的扰动项的分布，一般会有模型中的扰动项的分布，一般会有3个假设：正态（高斯）个假设：正态（高斯）分布、学生分布、学生t-分布和广义误差分布（分布和广义误差分布（GED）。给定一个分布假设，）。给定一个分布假设，GARCH模型常常使用极大似然估计法进行估计。下面分别介绍这模型常常使用极大似然估计法进行估计。下面分别介绍这3种分布，其中的种分布，其中的 代表参数向量。代表参数向量。1对于扰动项服从正态分布的对于扰动项服从正态分布的GARCH(1,1)模型，它的对模型，它的对数似然函数为数

31、似然函数为 （6.1.33）这里的这里的 t2是是ut的条件方差。的条件方差。TttttTttyTL12212)(21ln21)2ln(2)(lnx现在学习的是第44页，共49页 2如果扰动项服从学生如果扰动项服从学生t分布，分布，GARCH(1,1)模型的对数似然函模型的对数似然函数的形式就是数的形式就是（6.1.34）这样，参数的估计就变成了在自由度这样，参数的估计就变成了在自由度k2的约束下使对数似然函数（的约束下使对数似然函数（6.1.34）最大化的问题。当）最大化的问题。当k时，学生时，学生t-分布接近于正态分布。分布接近于正态分布。注注 式（式（6.1.34）和（）和（6.1.35

32、）中的）中的()代表代表 函数：函数：若若N是偶整数，则是偶整数，则 (N/2)=1 2 3(N/2)-1，有，有(2/2)=1；若若N是奇整数，则是奇整数，则 ，有有 。1)2(252321)2(NN)21(TttttTttkykkkkTL1221222)2()(1ln2)1(ln212)1()2()2(ln2)(lnx 现在学习的是第45页，共49页 3扰动项的分布为广义误差分布（扰动项的分布为广义误差分布（GED）时，）时，GARCH(1,1)模模型的对数似然函数的形式为型的对数似然函数的形式为 （6.1.35）这里的参数这里的参数r 0。如果。如果r=2，那么，那么GED就是一个正态分

33、布。就是一个正态分布。TtrtttTttryrrrrTL12221223)1()(3(ln21)2)(3()1(ln2)(ln x现在学习的是第46页，共49页 金融理论表明具有较高可观测到风险的资产可以获得更高的平均金融理论表明具有较高可观测到风险的资产可以获得更高的平均收益，其原因在于人们一般认为金融资产的收益应当与其风险成正比收益，其原因在于人们一般认为金融资产的收益应当与其风险成正比，风险越大，预期的收益就越高。这种利用条件方差表示预期风险的，风险越大，预期的收益就越高。这种利用条件方差表示预期风险的模型被称为模型被称为ARCH均值模型均值模型(ARCH-in-mean)或或ARCH-

34、M回归模型。回归模型。在在ARCH-M中我们把条件方差引进到均值方程中中我们把条件方差引进到均值方程中:（6.1.38）ARCH-M模型的另一种不同形式是将条件方差换成条件标准差模型的另一种不同形式是将条件方差换成条件标准差：（6.1.41）或取对数或取对数 （6.1.42）ttttuy2xttttuyxttttuy)ln(2x现在学习的是第47页，共49页 ARCH-M模型通常用于关于资产的预期收益与预期风险紧模型通常用于关于资产的预期收益与预期风险紧密相关的金融领域。预期风险的估计系数是风险收益交易的度量密相关的金融领域。预期风险的估计系数是风险收益交易的度量。例如，我们可以认为某股票指数，如上证的股票指数的收益率。例如，我们可以认为某股票指数，如上证的股票指数的收益率（returet）依赖于一个常数项及条件方差依赖于一个常数项及条件方差(风险风险)：这种类型的模型（其中期望风险用条件方差表示）就称为这种类型的模型（其中期望风险用条件方差表示）就称为GARCH-M模型。模型。tttureture)ln(22112112tttu现在学习的是第48页，共49页感谢大家观看感谢大家观看现在学习的是第49页，共49页