GARCH模型与应用简介(文献)

GARCH模型与应用简介前言 (2)1.GARCH模型 (7)2.模型的参数估计 (16)3.模型检验 (27)4. 模型的应用 (32)5.实例 (42)6.某些新进展 (46)参考文献 (50)0. 前言(随机序列的条件均值与条件方差简介)考察严平稳随机序列{y t}, 且E|y t|<∞. 记其均值Ey t=μ, 协方差函数γk=E{(y t-μ)(y t+k-μ)}. 其条件期望(或条件均值): E(y t∣y t-1,y t-2,…)≡ϕ(y t-1,y t-2,…), (0.1) 依条件期望的性质有Eϕ(y t-1,y t-2,…)=E{E(y t∣y t-1,y t-2,…)}= Ey t =μ. (0.2) 记误差(或残差):e t≡ y t -ϕ(y t-1,y t-2,…). (0.3) 由(0.1)(0.2)式必有:Ee t=Ey t-Eϕ(y t-1,y t-2,…)=Ey t-Ey t=0, (0-均值性) (0.4) 及Ee t2=E[y t -ϕ(y t-1,y t-2,…)]2=E{(y t-μ)-[ϕ(y t-1,y t-2,…)-μ]}2 (中心化)=E(y t-μ)2+E[ϕ(y t-1,y t-2,…)-μ]2-2E(y t-μ)[ϕ(y t-1,y t-2,…)-μ]=γ0+Var{ϕ(y t-1,y t-2,…)}-2EE{(y t-μ)[ϕ(y t-1,y t-2,…)-μ]∣y t-1,y t-2,…}( 根据Ex=E{E[x∣y t-1,y t-2,…]} )=γ0+Var{ϕ(y t-1,y t-2,…)}-2E{[ϕ(y t-1,y t-2,…)-μ]E[(y t-μ)∣y t-1,y t-2,…]}( 再用E[x⨯ψ( y t-1,y t-2,…)∣y t-1,y t-2,…]=ψ( y t-1,y t-2,…) E[x∣y t-1,y t-2,…];并取x= (y t-μ), ψ( y t-1,y t-2,…)=[ϕ(y t-1,y t-2,…)-μ];由(0.1)(0.2)可得)=γ0+Var{ϕ(y t-1,y t-2,…)}-2E[ϕ(y t-1,y t-2,…)-μ]2=γ0-Var{ϕ(y t-1,y t-2,…)}. (0.5) 即有:γ0=Var(y t)=Var(ϕ(y t-1,y t-2,…))+Var(e t). (0.6) 此式表明, y t的方差(=γ0)可表示为: 回归函数的方差(Var(ϕ(y t-1,y t-2,…)), 与残差的方差(Var(e t))之和.下边讨论e t的条件均值与条件方差.为了符号简便, 以下记F t-1={y t-1,y t-2,…}.首先考虑e t的条件均值:E(e t∣F t-1)=E{y t-ϕ( y t-1,y t-2,…) ∣ F t-1}=E(y t∣ F t-1)- E{ϕ( y t-1,y t-2,…) ∣ F t-1}= ϕ( y t-1,y t-2,…)- ϕ( y t-1,y t-2,…)=0. (0.7) 再看条件方差:Var(e t∣F t-1)=E{[e t- E(e t∣F t-1)]2∣ F t-1}= E{e t2∣ F t-1} (用(0.7)式)≡S2(y t-1,y t-2,…). (0.8)此处S2(y t-1,y t-2,…)为条件方差函数. 注意, e t的条件均值是零, 条件方差是非负的函数S2(y t-1,y t-2,…), 它不一定是常数!依(0.3)式, 平稳随机序列{y t}总有如下表达式:y t = ϕ( y t-1,y t-2,…)+e t, (0.9)其中ϕ(y t-1,y t-2,…)被称为自回归函数, 不一定是线性的. {e t}可称为新息序列, 与线性模型的新息序列不同, 除非{y t}是正态序列. 顺便指出, 满足(0.4)式的{e t}为鞅差序列, 因为对它的求和是离散的鞅序列. 由于{y t}是严平稳随机序列, 且E|y t|<∞,上述推演是严格的, 从而{e t}是严平稳的鞅差序列. 当{y t}有遍历性时, 它也是遍历的. 此处所涉及的抽象概念可不必深究.现在将e t标准化, 即令εt≡ e t/S(y t-1,y t-2,…).则有,E(εt∣F t-1)=E[e t/S(y t-1,y t-2,…)∣F t-1]={1/S(y t-1,y t-2,…)}E[e t∣F t-1]=0. (依(0.7)式) (0.10) 以及E(εt2∣F t-1)=E[e t2/S2(y t-1,y t-2,…)∣F t-1]={1/S2(y t-1,y t-2,…)}E[e t2∣F t-1] (用(0.8))={S2(y t-1,y t-2,…)}/{S2(y t-1,y t-2,…)}=1. (a.s.) (0.11)由此可见, {εt}也是平稳鞅差序列, 与{e t}相比, {εt}的条件方差为常数1. 于是(0.9)式可写为:y t=ϕ( y t-1,y t-2,…) + S(y t-1,y t-2,…)εt, (0.12)此式可称为条件异方差自回归模型, 所谓条件异方差就是指: 条件方差S2(y t-1,y t-2,…)不为常数. 请注意, 条件异方差自回归模型与下文中的自回归条件异方差模型是不同的概念!* 还有一点很重要, 如果(0.9)模型具有可逆性, 那么, Var(e t∣F t-1)=Var(e t∣y t-1,y t-2,…)=Var(e t∣e t-1,e t-2,…)≡h(e t-1,e t-2,…). (0.13)因此, 模型(0.12)式又xie 可成:y t=ϕ( y t-1,y t-2,…) + h1/2(e t-1,e t-2,…)εt. (0.14)请注意, 模型(0.12)(0.14)式是普遍适用(或称万用)的模型!但是, 为便于研究建模理论, 在(0.12)式中还附加假定:εt与{y t-1,y t-2,…}相互独立!此假定是实质性的, 人为的.它对{y t}的概率分布有实质性的限制.还须指出: 若在(0.9)式中直接假定e t与{y t-1,y t-2,…}独立, 此假定除了上述的人为性含义外, 还增多了如下假定:Var(e t2∣y t-1,y t-2,…) =Var(e t2)=常数. (0.15)这里用了条件期望的一条性质, 即当X与Y独立时,E(X∣Y)= EX.大家要问, 为什么加这些人为的假定呢?让我们回顾一下这些假定演变的历程吧.在文献中(0.9)式e t先后被假定为:“i.i.d. 且N(0,σ2)”，（1943--）“i.i.d. 且0-均值-方差有穷”，（1960--）“鞅差序列，且条件方差S2(...)=常数”，（1970--）“e t=S(y t-1, y t-2, … )εt，但{εt}为i.i.d. N(0,σ2)序列,而且S(y t-1, y t-2, … )为有限参模型”,（1982--）“e t=S(y t-1, y t-2, … )εt , 但{εt}为i.i.d.序列而且S(y t-1, y t-2, … )为有限参模型”。

（2000--）究其根源, 主要是受时间序列统计理论知识的限制.以上专门讨论了{e t}的定义, 性质, 和人为限制的历程. 但是, 这里也顺便提一下自回归函数ϕ(y t-1,y t-2,…)的发展史, 大致如下(不细论):线性→非线性参数→半参数→非参数。

在以上的讨论中, 使用记号ϕ(y t-1,y t-2,…), 是为了突出普适性. 在文献中和实际应用中, 所考虑的ϕ(y t-1,y t-2,…)的形式很简单. 半个多世纪来, 虽说有了很大的改进, 但是, 与最一般的ϕ(y t-1,y t-2,…)还有很大差距.类似的讨论也适用S(y t-1,y t-2,…). 也是为了突出普适性, 才引入了记号S(y t-1,y t-2,…)和模型(0.12)(0.14). 在文献中和实际应用中, 直到近二十来年才考虑了不为常数的S(y t-1,y t-2,…)的简单情况---ARCH模型. 近几年来, 也在向着半参数, 非参数方面发展. 但是, 与最一般的S(y t-1,y t-2,…)也还相差甚远.1. ARCH与GARCH模型1.1. 概述在条件异方差模型问世以前, 时间序列分析主要讨论自回归结构, 或者说, 主要讨论ϕ(y t-1,y t-2,…)的有关内容. 当条件异方差模型问世后, 在时间序列分析中, 特别是建模分析中, 就包含了两个内容,一个与ϕ(y t-1,y t-2,…)有关; 另一个与S(y t-1,y t-2,…)有关. 如何统计分析它们, 是摆在我们面前的主要问题. 对此问题, 通常作法是: 分两步完成, 先按平稳序列建模方法, 对ϕ(y t-1,y t-2,…)建立适当的模型, 比如AR模型; 由此获得弥合的残差序列, 把它当做新息序列{e t}的样本值, 再对它进行条件异方差建模分析. 分两步完成有方便之处, 其一, 做第一步时, 由于{e t}是鞅差序列, 其建模有理论根据. 其二, 在介绍条件异方差建模时, 可以只讨论ϕ(y t-1,y t-2,…)=0的情况. 这并无损失, 还便于理解条件异方差概念. 其实, 还有一言, 在金融统计中, 专门考虑条件异方差建模问题, 也有一定的实际背景.综上所说, 我们将专门讨论如下的鞅差平稳序列, 即, E(y t∣y t-1,y t-2,…)≡ϕ(y t-1,y t-2,…)=0. (1.1)Var(y t∣y t-1,y t-2,…)≡S2(y t-1,y t-2,…)>0. (1.2) 换句话说, 考虑如下的(0.9)模型y t=e t, (1.3) 它的标准化的模型(0.12)为y t=S(y t-1,y t-2,…)εt. (1.4) 请注意, 这一模型几乎含盖了所有的条件异方差模型. 我们不可能泛泛地讨论它. 再请回看对鞅差序列{e t}的限制的历程, 以下我们要讲的恰好是:“e t=S(y t-1, y t-2, … )εt，但{εt}为i.i.d. N(0,σ2)序列, 而且S(y t-1, y t-2, … )为有限参模型’’, (1982--).再新的内容, 我们也将提到. 至此, 大家完全明白我们将要讨论什么样的序列.为说明该序列的某些特征, 先看一看序列{e t}的自协方差函数序列:γe(k)=Ee t+k e t= E[E(e t+k e t∣e t+k-1,e t+k-2,…)]= E[e t E(e t+k∣e t+k-1,e t+k-2,…)]= E{e t⨯0}=0, k≥1.可见, 平稳鞅差序列也是白噪声. 根据自协方差序列做平稳序列的建模和谱分析时, 除了判断ϕ(y t-1,y t-2,…)=0外, 几乎无话可说. 换句话说, 相关性分析和谱分析不能对(1.4)式的序列作出更深刻的分析. 为了进一步获得它的深入的结构特征, 必须引入新的概念和新的方法.1.2.ARCH(p)模型.(ARCH--- Autoregressive Conditional Heteroscedasticity)在金融界, 大量的数据序列呈现不可预报性, 相当于前面的(0.9)或(0.12)式中的ϕ(y t-1,y t-2, …)=0, 于是有兴趣研究(1.4)模型. Engle(1982)首先提出并使用了如下的有限参数模型:y t=S(y t-1, y t-2, … )εt≡ h t1/2 εt, (1.5)h t=α0+α1y t-12+α2y t-22+…+αp y t-p2, (1.6)α0>0, αi≥0, i=1,2,…,p.其中{εt}为i.i.d.的序列, εt~N(0, 1), 且εt与{y t-1, y t-2, …}独立, 为了简化记号, 记h t=S2(y t-1, y t-2, … ).此模型被称为自回归条件异方差模型, 简记ARCH(p),其中p表示模型的阶数.很明显, 此模型只是普遍适用的(1.4)式模型的子类, 因为, 在ARCH模型中对模型(1.4)添加了很多的人为限制.为了增进对ARCH模型的了解, 我们将作几点明, 以代替严格的推理论述.其一, 限定{εt}为i.i.d.序列! 这是很强的限制, 这是由于现有理论的基楚所限.其二, 限定条件方差有(1.6)式的简单形式, 即h t=S2(y t-1, y t-2, … )=α0+α1y t-12+α2y t-22+…+αp y t-p2,是为了统计分析方便.其三, 限定εt服从正态分布, 是为了求极大似然估计方便. 限制εt~N(0, 1), 而不用εt~N(0, σ2), 是因为{εt}满足标准化的模型(0.11)式.其四, 限制α0>0, αi≥0, i=1,2,…,p, 是为了保证条件方差函数h t=S2(y t-1, y t-2, … )>0. 限制α0>0, 而不是α0≥0, 这是为了保证模型(1.5)(1.6)有平稳解, 否则, 当α0=0时它没有平稳解! 这可从以下简单例子看出. 考查如下ARCH(1) 模型:h t=α1 y t-12,将它代入(1.5)式得y t=h t1/2 εt=(α1 y t-12)1/2 εt,将它两边平方得y t2=α1y t-12εt2,将它两边取对数得log(y t2)=log(α1)+log(y t-12)+log(εt2), (1.7) 记x t=log(y t2), c=log(α1), ηt=log(εt2)(仍为i.i.d.序列), 上式为x t = c+ x t-1+ ηt,这不是熟知的一元AR(1)模型吗? 而且不满足平稳性条件! 所以, 没有平稳解. 从而模型(1.5)也没有平稳解.其五, 为使ARCH模型有平稳解, 对系数αi(i=1,2,…,p)还要加限制. 较早的限制(也是较强)是α1+α2+…+αp<1. (1.8) 在此条件下, 不仅有平稳解, 还有有穷二阶矩. 后来, 也有人放宽条件, 只保证有平稳解, 不保证有有穷二阶矩. 所有这些结果的推理, 都要用到非线性时间序列分析的新成果.其六, Engle(1982)首次先提出ARCH模型时, 使用了如下叙述:y t⎥y t-1,y t-2,…,y1 ~ N(0, h t), (1.5)’h t=α0+α1y t-12+α2y t-22+…+αp y t-p2,α0>0, αi≥0, i=1,2,…,p.易见, (1.5)’式与(1.5)式是等价的.其七, ARCH模型有不同的变形形式. 仿(1.7)式的做法, 即将(1.5)式两边平方, 再将(1.6)式代入其中可得y t2=h tεt2=(α0+α1y t-12+α2y t-22+…+αp y t-p2)εt2=(α0+α1y t-12+α2y t-22+…+αp y t-p2)(1+εt2-1)=α0+α1y t-12+α2y t-22+…+αp y t-p2+(εt2-1)(α0+α1y t-12+α2y t-22+…+αp y t-p2)=α0+α1y t-12+α2y t-22+…+αp y t-p2+ h t(εt2-1)=α0+α1y t-12+α2y t-22+…+αp y t-p2+ w t , (1.9) 对序列{y t2}而言, 此式很像线性AR(p)模型, 其中w t=h t(εt2-1)是一个平稳的鞅差序列, 因为E{w t|y t-1,y t-2, …}==E{h t(εt2-1)|y t-1,y t-2, …}= E{h tεt2|y t-1,y t-2, …}-E{h t|y t-1,y t-2, …}= h t E{εt2|y t-1,y t-2, …}-E{h t|y t-1,y t-2, …} (依(1.6))= h t - h t =0. (1.10) 用(1.9)式和线性AR(p)模型的求解方法, 可得{y t2}的平稳解. 但是, 从原理上说, 得到了{y t2}的解, 还不能说就得到了原序列{y t}的解. 好在当我们只关心y t的条件方差时, 有了{y t2}的解也足够用了. (1.9)式的变形方式是严格的, 可放心地使用它. 所谓使用它, 就是将原数据平方后得到y12, y22, …, y T2, 对它们建立AR(p)模型, 便得到参数α0,α1,…,αp的一种估计.如果对y t2=h tεt2两边取对数可得log(y t2)=log(h t)+log(εt2)=log(α0+α1y t-12+α2y t-22+…+αp y t-p2)+log(εt2)记x(t)=log(y t2), c=Elog(εt2), ηt=log(εt2)-c, 于是上式可写成x(t)=c+log(α0+α1e x(t-1)+α2e x(t-2)+…+αp e x(t-p))+ ηt.于是又得到ARCH模型的另一种变形. 此式是关于序列{x(t)}的非线性自回归模型, 注意, 上式中的序列{ηt}是i.i.d.的. 此外, ARCH模型还有别的表示方法, 不再一一介绍了.其八, 根据数据y1,y2,…,y T, 要作自回归条件异方差模型的统计分析, 包含两项内容, 首先是用假设检验方法, 判别这些数据是否有条件异方差条件性, 即, S(y t-1, y t-2, … )=常数? 如果是否定回答, 第二项内容就是对ARCH模型未知参数的估计. 在第2节中, 我们将介绍参数的估计方法, 在第3节中, 介绍检验方法.1.3. GARCH(Generalized ARCH) 模型:在Engle(1982)提出ARCH模型后, 受到应用者的关注, 特别是金融界. 稍后几年, 也被时间序列分析理论研究所重视. 从前面对新息序列{e t}限制条件的放宽过程可见, 提出ARCH模型, 无疑是对时间序列分析理论和应用研究有开拓性的意义. 在对ARCH模型的理论研究和应用中, 人们自然会发问: 在(1.6)式中, y t的条件方差S2(y t-1, y t-2, … )≡ h t=α0+α1y t-12+α2y t-22+…+αp y t-p2,只依赖于p个历史值, 能否考虑依赖全部历史值的情况? Bollerslev(1986)给出了回答, 他提出了如下的更广的模型, 即GARCH模型:y t=S(y t-1, y t-2, … )εt≡ h t1/2 εt, (1.11)h t=α0+α1y t-12+α2y t-22+…+αp y t-p2+β1h t-1+…+βq h t-q, (1.12) α0>0, αi≥0, i=1,2,…,p; βj≥0, j=1,2,…,q. (1.13) 其中{εt}为i.i.d.的N(0,1)分布, 且εt与{ y t-1, y t-2, …}独立.对此GARCH模型作如下说明:其一, 利用(1.12)式反复迭代可得知, h t= S2(y t-1, y t-2, … )确实依赖序列的全部历史值, 但是, h t仅依赖有限个参数.其二, 在1997年诺贝尔经济学奖, 被两位研究期权定价理论的Black-Scholes方程的学者获得. 从理论上人们发现, Black-Scholes方程的解是连续时间变化的随机过程, 对它进行等间隔离散化采样, 所得到的序列, 恰好满足GARCH模型. 于是, GARCH模型更被认可, 而且, 金融界特别偏爱GARCH模型.其三, 如前所述, (1.13)式的条件α0>0, 仍不能放宽为α0≥0. 而且, (1.13)式中的条件αi≥0, i=1,2,…,p, 还应附加一个限制: α1+α2+…+ αp>0, 否则如果全部αi=0 (i=1,2,…,p)将导致(1.12)式的h t为常数(仍用迭代法可证明). 这一点未在文献中指出, 一个潜在原因是: 应用者默认p ≥1, 且αp>0.其四, 与对ARCH模型的说明中的其五很类似, 为使GARCH模型有平稳解, 对系数αi(i=1,2,…,p)和βj≥0, j=1,2,…,q. 还要加限制. 较早的限制(也是较强)是α1+…+αp+β1+…+ βq <1. (1.14) 在此条件下, 不仅有平稳解, 还有有穷二阶矩. 其余的叙述与ARCH情况相同, 从略.其五, 统计问题. 与对ARCH模型的说明中的其七很类似. 但是, 它比ARCH模型要复杂些, 具体如下: y t2= h tεt2= h t(εt2-1+1)= h t+ h t(εt2-1)=α0+α1y t-12+α2y t-22+…+αp y t-p2+β1h t-1+…+βq h t-q+w t=α0+α1y t-12+α2y t-22+…+αp y t-p2+β1h t-1(εt-12-εt-12+1)+…+βq h t-q(εt-12-εt-q2+1)+w t =α0+α1y t-12+α2y t-22+…+αp y t-p2+β1h t-1εt-12+…+βq h t-qεt-q2-β1h t-1(εt-12-1)-…-βq h t-q(εt-q2-1)+w t=α0+α1y t-12+α2y t-22+…+αp y t-p2+β1y t-12+…+βq y t-q2-β1w t-1-…-βq w t-q+w t=α0+(α1+β1)y t-12+(α2+β2)y t-22+…+(αm+βm)y t-p2-β1w t-1-…-βq w t-q+w t, (1.15) 其中m=max{p,q}, 而且, 当k>p时αk=0; 当k>q时βk=0, w t=h t(εt2–1). 如前所述{w t}是平稳鞅差序列, 所以, 以上表达式说明, {h t}是由{w t}驱动的平稳ARMA序列. 以上模型不仅表达了GARCH模型的结构特性, 而且, 依此可借助于平稳ARMA序列建模方法, 得到GARCH模型参数的一种简单的估计方法. 关于GARCH模型的参数估计和检验方法, 分别在第2节和第3节中介绍.2. GARCH模型的参数估计2.1. 概述在实际应用中, 人们拥有序列观测值y1,y2,…,y n, 如果要为它们建立GARCH模型, 将面对着下列问题: 为什么要建立GARCH模型? 用多少阶数的模型? 怎样获得模型的参数值? 回答了这些问题, 就解决了为GARCH模型建模的问题. 前两个问题将在下一节中讨论, 这一节只讨论模型的参数估计问题, 换言之, 讨论在模型阶数已知时, 如何根据观测值y1,y2,…,y n, 估计出GARCH(或者ARCH) 模型的参数. 在统计学中有多种方法可以用来解决这一问题, 这里只介绍两种估计方法. 一种是比较简单的方法, 另一种是熟知的极大似然估计方法. 前一种估计可能不如后者精细, 但是它可作为用迭代法求取后者时的初始值. 另外, 对ARCH和GARCH模型而言, 它们的参数估计方法的难易程度有明显差异, 所以, 我们将分别予以介绍.2.2. ARCH模型的参数估计2.2.1. 最小二乘法估计最小二乘法是非常熟悉的方法，此方法是基于最小二乘原理。