异方差时间序列模型PPT课件.ppt

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

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

2、Mandelbrot,1963、Fama,1965)。波动集群性表明股票报酬波动是时变的，表明是异方差。异方差虽然不会影响回归系数的最小二乘估计的无偏性，但是将影响到回归系数估计的标准差和置信区间。例如，汇率，股票价格常常用随机游走过程描述，第三张，PPT共四十九页，创作于2022年6月 其中ut为白噪声过程。1995-2000年日元兑美元汇率时间序列及差分序列见图1和图2。图1 日元兑美元汇率序列JPY(1995-2000)图2日元兑美元汇率差分序列(收益)D(JPY)第四张，PPT共四十九页，创作于2022年6月图3 收益绝对值序列(1995-2000)图4 D(JPY)的平方(1995-

3、2000)这种序列的特征是（1）过程的方差不仅随时间变化，而且有时变化得很激烈。（2）按时间观察，表现出“波动集群”（volatility clustering）特征，即方差在一定时段中比较小，而在另一时段中比较大。（3）从取值的分布看表现的则是“高峰厚尾”（leptokurtosis and fat-tail）特征，即均值附近与尾区的概率值比正态分布大，而其余区域的概率比正态分布小。图5给出高峰厚尾分布示意图。图6给出一个高峰厚尾分布实例。第五张，PPT共四十九页，创作于2022年6月第六张，PPT共四十九页，创作于2022年6月显然现期方差与前期的“波动”有关系。描述这类关系的模型称为自回

4、归条件异方差（自回归条件异方差（ARCH）模型）模型（Engle 1982年提出）。使用ARCH模型的理由是：（1）通过预测xt或ut的变化量评估股票的持有或交易对收益所带来的风险有多大，以及决策的代价有多大；（2）可以预测xt的置信区间，它是随时间变化的；（3）对条件异方差进行正确估计后可以使回归参数的估计量更具有有效性。第七张，PPT共四十九页，创作于2022年6月为了说得更具体，让我们回到k-变量回归模型：如果ut的均值为零，对yt取基于(t-1)时刻的信息的期望，即Et-1(yt)，有如下的关系：由于yt的均值近似等于式（1）的估计值，所以式（1）也称为均值方程。(1)(2)第八张，P

5、PT共四十九页，创作于2022年6月假设在时刻(t 1)所有信息已知的条件下，扰动项ut的条件分布是：也就是，ut遵循以0为均值，(0+1u2t-1)为方差的正态分布。由于(3)中ut的方差依赖于前期的平方扰动项，我们称它为ARCH(1)过程：(3)通常用极大似然估计得到参数0,1,2,k,0,1的有效估计。第九张，PPT共四十九页，创作于2022年6月第二节 ARCH模型一、一、ARCH模型的定义模型的定义若一个平稳随机变量xt可以表示为AR(p)形式，其随机误差项的方差可用误差项平方的q阶分布滞后模型描述，（1）（2）则称ut 服从q阶的ARCH过程，记作utARCH(q)。其中(1)式称

6、作均值方程，(2)式称作ARCH方程。第十张，PPT共四十九页，创作于2022年6月(1)和(2)式还应满足如下条件。对于(1)式，为保证平稳性，特征方程的根应在单位圆之外。xt 的条件期望是xt 的无条件期望（T 时）是对于(2)式，由于 的非负性，对i应有如下约束，当全部i=0,i=1,2,q时，条件方差 。因为方差是非负的，所以要求0 0。为保证 是一个平稳过程，(2)式的特征方程第十一张，PPT共四十九页，创作于2022年6月的根应在单位圆之外。对i,i=1,2,q的另一个约束是对（2）式求期望的：则无条件方差为：可见若保证 是一个平稳过程，应该有约束0 (1+2+q)1。因为Var(

7、xt)=Var(ut)=，所以上式可以用来预测xt的方差。当T时：第十二张，PPT共四十九页，创作于2022年6月 二、二、ARCH模型的极大似然估计模型的极大似然估计ARCH模型经常应用在回归模型中。其中=(0 1,k-1),xt=(1 x1,xk-1)（xt的分量也可以包括yt的滞后变量），utARCH(q)。为计算方便，假定已知yt,xt的T+q组观测值，从而保证估计参数所用的样本容量为T。utARCH(q)可以表示为：其中vtIID(0,1)，vt与xt相互独立，所以有 ，E(ut)=0。yt服从正态分布，概率密度函数为第十三张，PPT共四十九页，创作于2022年6月其中：用参数和=(

8、0 1 2 q)组成参数向量，对数似然函数是:求的极大似然估计量就是求 使 logL()在=处获得极大值。求log L()对 的偏导数，第十四张，PPT共四十九页，创作于2022年6月第十五张，PPT共四十九页，创作于2022年6月在上式为零条件下求到的 即是 的极大似然估计量。具有一致性。第十六张，PPT共四十九页，创作于2022年6月 三、三、ARCH模型检验模型检验在均值方程（回归模型或时间序列模型）的误差项中是否存在自回归条件异方差应该进行假设检验。检验ARCH可以使用F、LM、LR、W统计量。下面介绍F、LM检验。1、自回归条件异方差的LM检验。建立原假设H0：1=2=q=0 （不存

9、在ARCH）H1：1,2,q 不全为零在原假设成立条件下，OLS估计量是一致的、有效的；在备择假设成立条件下，OLS估计量是一致的，但不是有效的。先介绍使用LM统计量检验H0。因为计算LM统计量的值，只需估计原假设成立条件下的方程。具体步骤是：第十七张，PPT共四十九页，创作于2022年6月 估计 ，求 ，计算 。估计辅助回归式 用第3步得到的可决系数R2构成统计量LM=T R2。其中T表示辅助回归式的样本容量。在原假设成立条件下有若LM ，接受H1。注意：注意：辅助回归式中要有常数项0。第十八张，PPT共四十九页，创作于2022年6月 2、自回归条件异方差的F检验。建立原假设 H0：1=2=

10、q=0（不存在ARCH）H1：1,2,q不全为零 估计yt=xt+ut，求 ，计算 。用 估计2个辅助回归式 构造F统计量，在原假设成立条件下有 注意：注意：辅助回归式中要有常数项0。若F F,(q,T-q-1)，接受H1。（约束模型，同方差）（约束模型，同方差）（非约束模型，存在（非约束模型，存在ARCH）第十九张，PPT共四十九页，创作于2022年6月3、自回归条件异方差的LR检验。建立原假设 H0：1=2=q=0（不存在ARCH）H1：1,2,q不全为零 估计yt=xt+ut，求 ，计算 。用 估计2个辅助回归式，并计算极大似然函数值LogLr和LogLu 用LogLr和LogLu构造L

11、R统计量，在原假设成立条件下有 （约束模型，同方差）（约束模型，同方差）（非约束模型，存在（非约束模型，存在ARCH）若 ，接受H0。若 ，接受H1。如果结论是应该建立ARCH模型，则进一步应该对ARCH模型的阶数q进行检验。对此可采用t检验。第二十张，PPT共四十九页，创作于2022年6月4、自回归条件异方差的Q检验。残差平方意味着方差，若存在自相关，说明存在自回归条件异方差。第二十一张，PPT共四十九页，创作于2022年6月四、四、ARCH模型检验（模型检验（EViews操作案例）操作案例）1995.1-2000.8日元兑美元汇率值（1427个）序列（JPY）见图。极小值为81.12日元，

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)建立时间序列模型。第二十二张，PPT共四十九页，创作于2022年6月通过相关图和偏自相关图分析，应该建立AR(3)或MA(3)模型。建立AR(3)模型如下：(2.0)(-3.3)R2=0.01，DW=1.91，Q(15)=8.6第二十三张，PPT共四十九页，创作于2022年6月 方法方法1：通过Q检验考察AR(3

13、)模型中是否存在自回归条件异方差。第二十四张，PPT共四十九页，创作于2022年6月 方法方法2：ARCH的LM检验。在均值方程估计窗口，选在均值方程估计窗口，选ARCH的的LM检验。用检验。用1阶检验式检验。阶检验式检验。(9.4)(9.9)R2=0.0643,T=1423第二十五张，PPT共四十九页，创作于2022年6月 方法方法3、4：自回归条件异方差的F检验和LR检验。用参差平方序列用参差平方序列1阶自回归检验式做参数约束的阶自回归检验式做参数约束的F检验和检验和LR检验检验第二十六张，PPT共四十九页，创作于2022年6月 例例 中国中国CPI模型的模型的ARCH检验检验 本例建立C

14、PI模型，因变量为中国的消费价格指数（上年同月=100）减去100，记为cpit；解释变量选择货币政策变量：狭义货币供应量M1的增长率，记为m1rt；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 第二十七张，PPT共四十九页，创作于2022年6月 这个方程的统计量很显著，拟合的程度也很好。但是观察该回归方程的残差图，也可以注意到波动的“成群”现象：波

15、动在一些时期内较小，在其他一些时期内较大，这说明误差项可能具有条件异方差性。第二十八张，PPT共四十九页，创作于2022年6月 从自相关系数和偏自相关系数可以看出：残差序列存在着一阶ARCH效应。再进行条件异方差的ARCH LM检验，得到了在滞后阶数p=1时的ARCH LM检验结果：因此计算残差平方t2的自相关（AC）和偏自相关（PAC）系数，结果如下：第二十九张，PPT共四十九页，创作于2022年6月 从自相关系数和偏自相关系数可以看出：残差序列存在着一阶从自相关系数和偏自相关系数可以看出：残差序列存在着一阶ARCH效效应。因此利用应。因此利用ARCH(1)模型重新估计模型模型重新估计模型，

16、结果如下：，结果如下：均值方程：均值方程：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)模型能够更好的拟合数模型能够更好的拟合数据。据。第三十张，PPT共四十九页，创作于2022年6月 再对这个方程进行条件异方差的再对这个方程进行条件异方差的ARCH L

17、M检验，得到了残差检验，得到了残差序列在滞后阶数序列在滞后阶数p=1时的统计结果：时的统计结果：此时的相伴概率为此时的相伴概率为0.69，接受原假设，认为该残差序列不存在，接受原假设，认为该残差序列不存在ARCH效应，说明利用效应，说明利用ARCH(1)模型消除了残差序列的条件异方差性。模型消除了残差序列的条件异方差性。残差平方相关图的检验结果为：残差平方相关图的检验结果为：自相关系数和偏自相关系数近似为自相关系数和偏自相关系数近似为0。这个结果也说明了残差序。这个结果也说明了残差序列不再存在列不再存在ARCH效应。效应。第三十一张，PPT共四十九页，创作于2022年6月第三节第三节 GARC

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

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

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

21、易商将会增加对下期方差的预期。这个模型还包括了经常可以在财务收益数据中看到的变动组，在这些数据中，收益的巨大变化可能伴随着更进一步的巨大变化。第三十五张，PPT共四十九页，创作于2022年6月 有两个可供选择的方差方程的描述可以帮助解释这个模型：1如果我们用条件方差的滞后递归地替代方差方程的右端，就可以将条件方差表示为滞后扰动项平方的加权平均：我们看到GARCH(1,1)方差说明与样本方差类似，但是，它包含了在更大滞后阶数上的，扰动项的加权条件方差。第三十六张，PPT共四十九页，创作于2022年6月 2设vt=ut2t2。用其替代方差方程中的方差并整理，得到关于扰动项平方的模型：因此，扰动项平

22、方服从一个异方差ARMA(1,1)过程。决定波动冲击持久性的自回归的根是加的和。在很多情况下，这个根非常接近1，所以冲击会逐渐减弱。第三十七张，PPT共四十九页，创作于2022年6月 高阶高阶GARCH(p,q)模型模型 高阶GARCH模型可以通过选择大于1的p或q得到估计，记作GARCH(p,q)。其方差表示为：这里，p是GARCH项的阶数，q是ARCH项的阶数，p0并且,(L)和(L)是滞后算子多项式。第三十八张，PPT共四十九页，创作于2022年6月 为了使GARCH(q,p)模型的条件方差有明确的定义，相应的ARCH()模型的所有系数都必须是正数。只要(L)和(L)没有相同的根并且(L

23、)的根全部位于单位圆外，那么当且仅当0=0/(1-(L)，(L)=(L)/(1-(L)的所有系数都非负时，这个正数限定条件才会满足。例如，对于GARCH(1,1)模型这些条件要求所有的3个参数都是非负数。第三十九张，PPT共四十九页，创作于2022年6月 GARCH模型的检验模型的检验当原假设H0是ARCH(0)时，显然备择假设H1有两个。一个一个是是ARCH(r)，一个是，一个是GARCH(r,0)。若原假设 H0是ARCH(1)，则备择假设H1可以是ARCH(1+r)，也可以是GARCH(r,1)。同理若原假设H0是ARCH(q)，备择假设H1可以有两个。ARCH(q+r)和GARCH(r

24、,q)。LM统计量无法区别这两个备择假设。但仍然可以做LM检验。在实际应用中，备择假设既可以是ARCH，也可以是GARCH。对于q值很大的ARCH模型，建议使用GARCH模型。在实际应用中，GARCH(1,1)和GARCH(2,1)一般足以满足对自回归条件异方差的描述。第四十张，PPT共四十九页，创作于2022年6月第四节第四节 其他其他GARCH模型模型一、一、IGARCH模型模型 对于ARCH(p)模型和GARCH(p,q)模型，在实际应用中，条件 有时不能得到满足。下面以GARCH(1,1)模型为例进行讨论。（保证可以转换成无限阶的ARCH过程）（保证GARCH过程平稳）第四十一张，PP

25、T共四十九页，创作于2022年6月用 分别表示对1,1的估计。有时会出现 甚至，。例如Engle-Chowdury（1992）对IBM收益率序列估计时，得如下结果，第四十二张，PPT共四十九页，创作于2022年6月 2 2回推回推回推回推 在计算在计算GARCH模型的回推初始方差时，首先用系数值来计模型的回推初始方差时，首先用系数值来计算均值方程中的残差，然后计算初始值的指数平滑算子算均值方程中的残差，然后计算初始值的指数平滑算子 （6.1.30）其中：是来自均值方程的残差，是无条件方差的估计：其中：是来自均值方程的残差，是无条件方差的估计：（6.1.31）平滑参数平滑参数为为0.1至至1之间

26、的数值。也可以使用无条件方差来初始化之间的数值。也可以使用无条件方差来初始化GARCH过程：过程：（6.1.32）第四十三张，PPT共四十九页，创作于2022年6月6.1.6 GARCH6.1.6 GARCH模型的残差分布假设模型的残差分布假设模型的残差分布假设模型的残差分布假设 在实践中我们注意到，许多时间序列，特别是金融时间序列在实践中我们注意到，许多时间序列，特别是金融时间序列的无条件分布往往具有比正态分布更宽的尾部。为了更精确地描的无条件分布往往具有比正态分布更宽的尾部。为了更精确地描述这些时间序列分布的尾部特征，还需要对误差项述这些时间序列分布的尾部特征，还需要对误差项ut的分布进行

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

28、差。第四十四张，PPT共四十九页，创作于2022年6月 2如果扰动项服从学生如果扰动项服从学生t分布，分布，GARCH(1,1)模型的对数似然函数的模型的对数似然函数的形式就是形式就是（6.1.34）这样，参数的估计就变成了在自由度这样，参数的估计就变成了在自由度k2的约束下使对数似然函数的约束下使对数似然函数（6.1.34）最大化的问题。当）最大化的问题。当k时，学生时，学生t-分布接近于正态分布。分布接近于正态分布。注注 式（式（6.1.34）和（）和（6.1.35）中的）中的()代表代表 函数：函数：若若N是偶整数，则是偶整数，则 (N/2)=1 2 3(N/2)-1，有，有(2/2)=

29、1；若若N是奇整数，则是奇整数，则 ，有有 。第四十五张，PPT共四十九页，创作于2022年6月 3扰动项的分布为广义误差分布（扰动项的分布为广义误差分布（GED）时，）时，GARCH(1,1)模模型的对数似然函数的形式为型的对数似然函数的形式为 （6.1.35）这里的参数这里的参数r 0。如果。如果r=2，那么，那么GED就是一个正态分布。就是一个正态分布。第四十六张，PPT共四十九页，创作于2022年6月6.1.76.1.7 ARCH-MARCH-M模型模型模型模型 金金融融理理论论表表明明具具有有较较高高可可观观测测到到风风险险的的资资产产可可以以获获得得更更高高的的平平均均收收益益，其

30、其原原因因在在于于人人们们一一般般认认为为金金融融资资产产的的收收益益应应当当与与其其风风险险成成正正比比，风风险险越越大大，预预期期的的收收益益就就越越高高。这这种种利利用用条条件件方方差差表表示示预预期期风风险险的的模模型型被被称称为为ARCH均均值值模模型型(ARCH-in-mean)或或ARCH-M回回归归模模型型。在在ARCH-M中我们把条件方差引进到均值方程中中我们把条件方差引进到均值方程中:（6.1.38）ARCH-M模型的另一种不同形式是将条件方差换成条件标准差：模型的另一种不同形式是将条件方差换成条件标准差：（6.1.41）或取对数或取对数 （6.1.42）第四十七张，PPT

31、共四十九页，创作于2022年6月 ARCH-M模模型型通通常常用用于于关关于于资资产产的的预预期期收收益益与与预预期期风风险险紧紧密密相相关关的的金金融融领领域域。预预期期风风险险的的估估计计系系数数是是风风险险收收益益交交易易的的度度量量。例例如如，我我们们可可以以认认为为某某股股票票指指数数，如如上上证证的的股股票票指指数数的的收收益益率率（returet）依赖于一个常数项及条件方差依赖于一个常数项及条件方差(风险风险)：这种类型的模型（其中期望风险用条件方差表示）就称为这种类型的模型（其中期望风险用条件方差表示）就称为GARCH-M模模型。型。第四十八张，PPT共四十九页，创作于2022年6月感谢大家观看第四十九张，PPT共四十九页，创作于2022年6月