首页 GARCH模型与应用简介

GARCH模型与应用简介

举报
开通vip

GARCH模型与应用简介PAGEPAGE7GARCH模型与应用简介(2006,5)前言……………………………………………..2GARCH模型………………………………………….7模型的参数估计………………………………………16模型检验………………………………………………274.模型的应用……………………………………………32实例……………………………….……………………42某些新进展……………………….…………………...46参考文献……………………………………………….500.前言(随机序列的条件均值与条件方差简介)考察严平稳...

GARCH模型与应用简介
PAGEPAGE7GARCH模型与应用简介(2006,5)前言……………………………………………..2GARCH模型………………………………………….7模型的参数估计………………………………………16模型检验………………………………………………274.模型的应用……………………………………………32实例……………………………….……………………42某些新进展……………………….…………………...46参考文献……………………………………………….500.前言(随机序列的条件均值与条件方差简介)考察严平稳随机序列{yt},且Eyt<.记其均值Eyt=,协方差函数k=E{(yt-)(yt+k-)}.其条件期望(或条件均值):E(ytyt-1,yt-2,…)(yt-1,yt-2,…),(0.1)依条件期望的性质有E(yt-1,yt-2,…)=E{E(ytyt-1,yt-2,…)}=Eyt=.(0.2)记误差(或残差):etyt-(yt-1,yt-2,…).(0.3)由(0.1)(0.2)式必有:Eet=Eyt-E(yt-1,yt-2,…)=Eyt-Eyt=0,(0-均值性)(0.4)及Eet2=E[yt-(yt-1,yt-2,…)]2=E{(yt-)-[(yt-1,yt-2,…)-]}2(中心化)=E(yt-)2+E[(yt-1,yt-2,…)-]2-2E(yt-)[(yt-1,yt-2,…)-]=0+Var{(yt-1,yt-2,…)}-2EE{(yt-)[(yt-1,yt-2,…)-]yt-1,yt-2,…}(根据Ex=E{E[xyt-1,yt-2,…]})=0+Var{(yt-1,yt-2,…)}-2E{[(yt-1,yt-2,…)-]E[(yt-)yt-1,yt-2,…]}(再用E[x(yt-1,yt-2,…)yt-1,yt-2,…]=(yt-1,yt-2,…)E[xyt-1,yt-2,…];并取x=(yt-),(yt-1,yt-2,…)=[(yt-1,yt-2,…)-];由(0.1)(0.2)可得)=0+Var{(yt-1,yt-2,…)}-2E[(yt-1,yt-2,…)-]2=0-Var{(yt-1,yt-2,…)}.(0.5)即有:0=Var(yt)=Var((yt-1,yt-2,…))+Var(et).(0.6)此式 关于同志近三年现实表现材料材料类招标技术评分表图表与交易pdf视力表打印pdf用图表说话 pdf 明,yt的方差(=0)可表示为:回归函数的方差(Var((yt-1,yt-2,…)),与残差的方差(Var(et))之和.下边讨论et的条件均值与条件方差.为了符号简便,以下记Ft-1={yt-1,yt-2,…}.首先考虑et的条件均值:E(etFt-1)=E{yt-(yt-1,yt-2,…)Ft-1}=E(ytFt-1)-E{(yt-1,yt-2,…)Ft-1}=(yt-1,yt-2,…)-(yt-1,yt-2,…)=0.(0.7)再看条件方差:Var(etFt-1)=E{[et-E(etFt-1)]2Ft-1}=E{et2Ft-1}(用(0.7)式)S2(yt-1,yt-2,…).(0.8)此处S2(yt-1,yt-2,…)为条件方差函数.注意,et的条件均值是零,条件方差是非负的函数S2(yt-1,yt-2,…),它不一定是常数!依(0.3)式,平稳随机序列{yt}总有如下表达式:yt=(yt-1,yt-2,…)+et,(0.9)其中(yt-1,yt-2,…)被称为自回归函数,不一定是线性的.{et}可称为新息序列,与线性模型的新息序列不同,除非{yt}是正态序列.顺便指出,满足(0.4)式的{et}为鞅差序列,因为对它的求和是离散的鞅序列.由于{yt}是严平稳随机序列,且Eyt<,上述推演是严格的,从而{et}是严平稳的鞅差序列.当{yt}有遍历性时,它也是遍历的.此处所涉及的抽象概念可不必深究.现在将et标准化,即令tet/S(yt-1,yt-2,…).则有,E(tFt-1)=E[et/S(yt-1,yt-2,…)Ft-1]={1/S(yt-1,yt-2,…)}E[etFt-1]=0.(依(0.7)式)(0.10)以及E(t2Ft-1)=E[et2/S2(yt-1,yt-2,…)Ft-1]={1/S2(yt-1,yt-2,…)}E[et2Ft-1](用(0.8))={S2(yt-1,yt-2,…)}/{S2(yt-1,yt-2,…)}=1.(a.s.)(0.11)由此可见,{t}也是平稳鞅差序列,与{et}相比,{t}的条件方差为常数1.于是(0.9)式可写为:yt=(yt-1,yt-2,…)+S(yt-1,yt-2,…)t,(0.12)此式可称为条件异方差自回归模型,所谓条件异方差就是指:条件方差S2(yt-1,yt-2,…)不为常数.请注意,条件异方差自回归模型与下文中的自回归条件异方差模型是不同的概念!*还有一点很重要,如果(0.9)模型具有可逆性,那么,Var(etFt-1)=Var(etyt-1,yt-2,…)=Var(etet-1,et-2,…)h(et-1,et-2,…).(0.13)因此,模型(0.12)式又可些成yt=(yt-1,yt-2,…)+h1/2(et-1,et-2,…)t.(0.14)请注意,模型(0.12)(0.14)式是普遍适用(或称万用)的模型!但是,为便于研究建模理论,在(0.12)式中还附加假定:t与{yt-1,yt-2,…}相互独立!此假定是实质性的,人为的.它对{yt}的概率分布有实质性的限制.还须指出:若在(0.9)式中直接假定et与{yt-1,yt-2,…}独立,此假定除了上述的人为性含义外,还增多了如下假定:Var(et2yt-1,yt-2,…)=Var(et2)=常数.(0.15)这里用了条件期望的一条性质,即当X与Y独立时,E(XY)=EX.大家要问,为什么加这些人为的假定呢?让我们回顾一下这些假定演变的历程吧.在文献中(0.9)式et先后被假定为:“i.i.d.且N(0,σ2)”,(1943--)“i.i.d.且0-均值-方差有穷”,(1960--)“鞅差序列,且条件方差S2(...)=常数”,(1970--)“et=S(yt-1,yt-2,…)t,但{t}为i.i.d.N(0,σ2)序列,而且S(yt-1,yt-2,…)为有限参模型”,(1982--)“et=S(yt-1,yt-2,…)t,但{t}为i.i.d.序列而且S(yt-1,yt-2,…)为有限参模型”。(2000--)究其根源,主要是受时间序列统计理论知识的限制.以上专门讨论了{et}的定义,性质,和人为限制的历程.但是,这里也顺便提一下自回归函数(yt-1,yt-2,…)的发展史,大致如下(不细论):线性→非线性参数→半参数→非参数。在以上的讨论中,使用记号(yt-1,yt-2,…),是为了突出普适性.在文献中和实际应用中,所考虑的(yt-1,yt-2,…)的形式很简单.半个多世纪来,虽说有了很大的改进,但是,与最一般的(yt-1,yt-2,…)还有很大差距.类似的讨论也适用S(yt-1,yt-2,…).也是为了突出普适性,才引入了记号S(yt-1,yt-2,…)和模型(0.12)(0.14).在文献中和实际应用中,直到近二十来年才考虑了不为常数的S(yt-1,yt-2,…)的简单情况---ARCH模型.近几年来,也在向着半参数,非参数方面发展.但是,与最一般的S(yt-1,yt-2,…)也还相差甚远.1.ARCH与GARCH模型1.1.概述在条件异方差模型问世以前,时间序列 分析 定性数据统计分析pdf销售业绩分析模板建筑结构震害分析销售进度分析表京东商城竞争战略分析 主要讨论自回归结构,或者说,主要讨论(yt-1,yt-2,…)的有关内容.当条件异方差模型问世后,在时间序列分析中,特别是建模分析中,就包含了两个内容,一个与(yt-1,yt-2,…)有关;另一个与S(yt-1,yt-2,…)有关.如何统计分析它们,是摆在我们面前的主要问 快递公司问题件快递公司问题件货款处理关于圆的周长面积重点题型关于解方程组的题及答案关于南海问题 .对此问题,通常作法是:分两步完成,先按平稳序列建模方法,对(yt-1,yt-2,…)建立适当的模型,比如AR模型;由此获得弥合的残差序列,把它当做新息序列{et}的样本值,再对它进行条件异方差建模分析.分两步完成有方便之处,其一,做第一步时,由于{et}是鞅差序列,其建模有理论根据.其二,在介绍条件异方差建模时,可以只讨论(yt-1,yt-2,…)=0的情况.这并无损失,还便于理解条件异方差概念.其实,还有一言,在金融统计中,专门考虑条件异方差建模问题,也有一定的实际背景.综上所说,我们将专门讨论如下的鞅差平稳序列,即,E(ytyt-1,yt-2,…)(yt-1,yt-2,…)=0.(1.1)Var(ytyt-1,yt-2,…)S2(yt-1,yt-2,…)>0.(1.2)换句话说,考虑如下的(0.9)模型yt=et,(1.3)它的标准化的模型(0.12)为yt=S(yt-1,yt-2,…)t.(1.4)请注意,这一模型几乎含盖了所有的条件异方差模型.我们不可能泛泛地讨论它.再请回看对鞅差序列{et}的限制的历程,以下我们要讲的恰好是:“et=S(yt-1,yt-2,…)t,但{t}为i.i.d.N(0,σ2)序列,而且S(yt-1,yt-2,…)为有限参模型’’,(1982--).再新的内容,我们也将提到.至此,大家完全明白我们将要讨论什么样的序列.为说明该序列的某些特征,先看一看序列{et}的自协方差函数序列:e(k)=Eet+ket=E[E(et+ketet+k-1,et+k-2,…)]=E[etE(et+ket+k-1,et+k-2,…)]=E{et0}=0,k1.可见,平稳鞅差序列也是白噪声.根据自协方差序列做平稳序列的建模和谱分析时,除了判断(yt-1,yt-2,…)=0外,几乎无话可说.换句话说,相关性分析和谱分析不能对(1.4)式的序列作出更深刻的分析.为了进一步获得它的深入的结构特征,必须引入新的概念和新的方法.ARCH(p)模型.(ARCH---AutoregressiveConditionalHeteroscedasticity)在金融界,大量的数据序列呈现不可预报性,相当于前面的(0.9)或(0.12)式中的(yt-1,yt-2,…)=0,于是有兴趣研究(1.4)模型.Engle(1982)首先提出并使用了如下的有限参数模型:yt=S(yt-1,yt-2,…)tht1/2t,(1.5)ht=0+1yt-12+2yt-22+…+pyt-p2,(1.6)0>0,i0,i=1,2,…,p.其中{t}为i.i.d.的序列,tN(0,1),且t与{yt-1,yt-2,…}独立,为了简化记号,记ht=S2(yt-1,yt-2,…).此模型被称为自回归条件异方差模型,简记ARCH(p),其中p表示模型的阶数.很明显,此模型只是普遍适用的(1.4)式模型的子类,因为,在ARCH模型中对模型(1.4)添加了很多的人为限制.为了增进对ARCH模型的了解,我们将作几点明,以代替严格的推理论述.其一,限定{t}为i.i.d.序列!这是很强的限制,这是由于现有理论的基楚所限.其二,限定条件方差有(1.6)式的简单形式,即ht=S2(yt-1,yt-2,…)=0+1yt-12+2yt-22+…+pyt-p2,是为了统计分析方便.其三,限定t服从正态分布,是为了求极大似然估计方便.限制tN(0,1),而不用tN(0,2),是因为{t}满足标准化的模型(0.11)式.其四,限制0>0,i0,i=1,2,…,p,是为了保证条件方差函数ht=S2(yt-1,yt-2,…)>0.限制0>0,而不是00,这是为了保证模型(1.5)(1.6)有平稳解,否则,当0=0时它没有平稳解!这可从以下简单例子看出.考查如下ARCH(1)模型:ht=1yt-12,将它代入(1.5)式得yt=ht1/2t=(1yt-12)1/2t,将它两边平方得yt2=1yt-12t2,将它两边取对数得log(yt2)=log(1)+log(yt-12)+log(t2),(1.7)记xt=log(yt2),c=log(1),t=log(t2)(仍为i.i.d.序列),上式为xt=c+xt-1+t,这不是熟知的一元AR(1)模型吗?而且不满足平稳性条件!所以,没有平稳解.从而模型(1.5)也没有平稳解.其五,为使ARCH模型有平稳解,对系数i(i=1,2,…,p)还要加限制.较早的限制(也是较强)是1+2+…+p<1.(1.8)在此条件下,不仅有平稳解,还有有穷二阶矩.后来,也有人放宽条件,只保证有平稳解,不保证有有穷二阶矩.所有这些结果的推理,都要用到非线性时间序列分析的新成果.其六,Engle(1982)首次先提出ARCH模型时,使用了如下叙述:ytyt-1,yt-2,…,y1N(0,ht),(1.5)’ht=0+1yt-12+2yt-22+…+pyt-p2,0>0,i0,i=1,2,…,p.易见,(1.5)’式与(1.5)式是等价的.其七,ARCH模型有不同的变形形式.仿(1.7)式的做法,即将(1.5)式两边平方,再将(1.6)式代入其中可得yt2=htt2=(0+1yt-12+2yt-22+…+pyt-p2)t2=(0+1yt-12+2yt-22+…+pyt-p2)(1+t2-1)=0+1yt-12+2yt-22+…+pyt-p2+(t2-1)(0+1yt-12+2yt-22+…+pyt-p2)=0+1yt-12+2yt-22+…+pyt-p2+ht(t2-1)=0+1yt-12+2yt-22+…+pyt-p2+wt,(1.9)对序列{yt2}而言,此式很像线性AR(p)模型,其中wt=ht(t2-1)是一个平稳的鞅差序列,因为E{wt|yt-1,yt-2,…}==E{ht(t2-1)|yt-1,yt-2,…}=E{htt2|yt-1,yt-2,…}-E{ht|yt-1,yt-2,…}=htE{t2|yt-1,yt-2,…}-E{ht|yt-1,yt-2,…}(依(1.6))=ht-ht=0.(1.10)用(1.9)式和线性AR(p)模型的求解方法,可得{yt2}的平稳解.但是,从原理上说,得到了{yt2}的解,还不能说就得到了原序列{yt}的解.好在当我们只关心yt的条件方差时,有了{yt2}的解也足够用了.(1.9)式的变形方式是严格的,可放心地使用它.所谓使用它,就是将原数据平方后得到y12,y22,…,yT2,对它们建立AR(p)模型,便得到参数0,1,…,p的一种估计.如果对yt2=htt2两边取对数可得log(yt2)=log(ht)+log(t2)=log(0+1yt-12+2yt-22+…+pyt-p2)+log(t2)记x(t)=log(yt2),c=Elog(t2),t=log(t2)-c,于是上式可写成x(t)=c+log(0+1ex(t-1)+2ex(t-2)+…+pex(t-p))+t.于是又得到ARCH模型的另一种变形.此式是关于序列{x(t)}的非线性自回归模型,注意,上式中的序列{t}是i.i.d.的.此外,ARCH模型还有别的表示方法,不再一一介绍了.其八,根据数据y1,y2,…,yT,要作自回归条件异方差模型的统计分析,包含两项内容,首先是用假设检验方法,判别这些数据是否有条件异方差条件性,即,S(yt-1,yt-2,…)=常数?如果是否定回答,第二项内容就是对ARCH模型未知参数的估计.在第2节中,我们将介绍参数的估计方法,在第3节中,介绍检验方法.1.3.GARCH(GeneralizedARCH)模型:在Engle(1982)提出ARCH模型后,受到应用者的关注,特别是金融界.稍后几年,也被时间序列分析理论研究所重视.从前面对新息序列{et}限制条件的放宽过程可见,提出ARCH模型,无疑是对时间序列分析理论和应用研究有开拓性的意义.在对ARCH模型的理论研究和应用中,人们自然会发问:在(1.6)式中,yt的条件方差S2(yt-1,yt-2,…)ht=0+1yt-12+2yt-22+…+pyt-p2,只依赖于p个历史值,能否考虑依赖全部历史值的情况?Bollerslev(1986)给出了回答,他提出了如下的更广的模型,即GARCH模型:yt=S(yt-1,yt-2,…)tht1/2t,(1.11)ht=0+1yt-12+2yt-22+…+pyt-p2+1ht-1+…+qht-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与{yt-1,yt-2,…}独立.对此GARCH模型作如下说明:其一,利用(1.12)式反复迭代可得知,ht=S2(yt-1,yt-2,…)确实依赖序列的全部历史值,但是,ht仅依赖有限个参数.其二,在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)式的ht为常数(仍用迭代法可证明).这一点未在文献中指出,一个潜在原因是:应用者默认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模型要复杂些,具体如下:yt2=htt2=ht(t2-1+1)=ht+ht(t2-1)=0+1yt-12+2yt-22+…+pyt-p2+1ht-1+…+qht-q+wt=0+1yt-12+2yt-22+…+pyt-p2+1ht-1(t-12-t-12+1)+…+qht-q(t-12-t-q2+1)+wt=0+1yt-12+2yt-22+…+pyt-p2+1ht-1t-12+…+qht-qt-q2-1ht-1(t-12-1)-…-qht-q(t-q2-1)+wt=0+1yt-12+2yt-22+…+pyt-p2+1yt-12+…+qyt-q2-1wt-1-…-qwt-q+wt=0+(1+1)yt-12+(2+2)yt-22+…+(m+m)yt-p2-1wt-1-…-qwt-q+wt,(1.15)其中m=max{p,q},而且,当k>p时k=0;当k>q时k=0,wt=ht(t2–1).如前所述{wt}是平稳鞅差序列,所以,以上表达式说明,{ht}是由{wt}驱动的平稳ARMA序列.以上模型不仅表达了GARCH模型的结构特性,而且,依此可借助于平稳ARMA序列建模方法,得到GARCH模型参数的一种简单的估计方法.关于GARCH模型的参数估计和检验方法,分别在第2节和第3节中介绍.2.GARCH模型的参数估计2.1.概述在实际应用中,人们拥有序列观测值y1,y2,…,yn,如果要为它们建立GARCH模型,将面对着下列问题:为什么要建立GARCH模型?用多少阶数的模型?怎样获得模型的参数值?回答了这些问题,就解决了为GARCH模型建模的问题.前两个问题将在下一节中讨论,这一节只讨论模型的参数估计问题,换言之,讨论在模型阶数已知时,如何根据观测值y1,y2,…,yn,估计出GARCH(或者ARCH)模型的参数.在统计学中有多种方法可以用来解决这一问题,这里只介绍两种估计方法.一种是比较简单的方法,另一种是熟知的极大似然估计方法.前一种估计可能不如后者精细,但是它可作为用迭代法求取后者时的初始值.另外,对ARCH和GARCH模型而言,它们的参数估计方法的难易程度有明显差异,所以,我们将分别予以介绍.2.2.ARCH模型的参数估计2.2.1.最小二乘法估计最小二乘法是非常熟悉的方法,此方法是基于最小二乘原理。我们先指出在此可以使用此原理的依据,为此不妨以ARCH(1)模型为例说明之。依(1.9)式知,满足ARCH(1)模型的序列{yt}必满足以下模型yt2=0+1yt-12+wt,(2.1)其中{wt}是鞅差序列,而且wt=ht(t2-1),于是有E{wt|yt-12}=E{ht(t2-1)|yt-12}=htE{(t2-1)|yt-12}=htE{t2|yt-12}-ht=htE{t2}-ht=ht-ht=0.(a.s.)(2.2)利用此式可得知,E{yt2-a0-a1yt-12}2=E{0+1yt-12+wt-a0-a1yt-12}2=E{(0-a0)+(1-a1)yt-12+wt}2=E{(0-a0)+(1-a1)yt-12}2+E{wt}2+2E{[(0-a0)+(1-a1)yt-12]wt}=E{(0-a0)+(1-a1)yt-12}2+E{wt}2+2EE{[(0-a0)+(1-a1)yt-12]wt|yt-12}=E{(0-a0)+(1-a1)yt-12}2+E{wt}2+2EE{[(0-a0)+(1-a1)yt-12]E{wt|yt-12}=E{(0-a0)+(1-a1)yt-12}2+E{wt}2(by(2.2))=E{(0-a0)+(1-a1)yt-12}2+E{ht(t2-1)}2=E{(0-a0)+(1-a1)yt-12}2+Eht2E(t2-1)2Eht2E(t2-1)2=c.(依平稳性)易见,上式中的=号成立,当且仅当(0-a0)=(1-a1)=0.此事实表明,min{E(yt2-a0-a1yt-12)2:a0,a1}=E{yt2-0-1yt-12}2.(2.3)此式表明,用所有可能的系数拟合(2.1)模型时,只有以其真系数拟合,才使拟合参差的方差最小!在实际应用时,我们没有(2.1)式中的确切的概率分布,但是,我们有序列{yt2}的观测数据y1,y2,…,yn,根据统计学的基楚性原理---大数定律,(2.3)式的最小化特征,用样本平均代替之,随着样本个数的增加将近似成立。换言之,求解以下最小化问题之解,即min{(n-1)-1t=2n(yt2-a0-a1yt-12)2:a0,a1}=(n-1)-1t=2n(yt2-a0*-a1*yt-12)2,显然,此问题等价于如下的最小化问题min{t=2n(yt2-a0-a1yt-12)2:a0,a1}=t=2n(yt2-a0*-a1*yt-12)2.(2.4)以其解(a0*,a1*)作为真参数(0,1)的估计,称它们为最小二乘估计。这就是使用最小二乘原理的依据。以上论述不难推广到一般的ARCH(p)模型,除了符号的繁琐外,并无本质差异。这里只强调一点:对ARCH(1)使用最小二乘原理时,残差项wt与yt-1相互独立且Ewt=0是常见的条件,至少也要满足条件E{wt|yt-12}=0(a.s.)。这一点对一般情况也适用。现在介绍ARCH(p)模型参数最小二乘估计方法。首先重新写出(1.9)式yt2=0+1yt-12+2yt-22+…+pyt-p2+wt,t=p+1,p+2,…,n.(2.5)在此特别强调足标t的取值范围,只是为了模型中的yt-p都落在我们的数据序列中。依前所述,未知参数=(0,1,2,…,p)的最小二乘估计*,就是如下的最小值问题的解,即min{t=2n(yt2-a0-a1yt-12-a2yt-22-…-apyt-p2)2:a0,a1,…,ap}=t=2n(yt2-a0*-a1*yt-12-…-ap*yt-p2)2,(2.6)取最小二乘估计*=(a0*,a1*,…,ap*)。欲给出(a0*,a1*,…,ap*)的表达方式,既可用分析方法,又可用代数方法。现在使用后一方法,为此将(2.5)式改写成Y=X+W,(2.7)其中Y=(yp+12,yp+22,…,yn2),W=(wp+12,wp+22,…,wn2),X=.当以a=(a0,a1,a2,…,ap)为自由参数向量时,于是有t=p+1n(yt2-a0-a1yt-12-a2yt-22-…-apyt-p2)2=||Y-Xa||2=(Y-Xa)t(Y-Xa)=YtY-YtXa-atXtY+atXtXa=(XtXa-XtY)t(XtX)-1(XtXa-XtY)+YtY–(XtY)t(XtX)-1(XtY)³YtY–(XtY)t(XtX)-1(XtY),其中用到了以下的矩阵性质(XtXa-XtY)t(XtX)-1(XtXa-XtY)³0.由前一式可知,(2.6)式的最小值解必满足XtXa-XtY=0.现在求解XtXa-XtY=0,即XtXa=XtY,其解为a*=(XtX)-1XtY.(2.8)注意,上式右边的矩阵X和向量Y,都是由已知数据量组成的,计算(XtX)-1和(XtX)-1XtY,有许多软件可供使用.当然,也可以自行编程序计算之.自回归模型(1.9)的系数的最小二乘估计,被(2.8)式明显的表达出,而且便于计算.这一优越性是自回归模型所特有的,因此,自回归模型在时间序列分析中问世最早.类似地,Engel(1982)最先引入的条件异方差模型,又是自回归型的条件异方差模型---ARCH模型,也是基于这一便于使用的优点.稍后几年才由Bollerslev(1986)提出更一般的GARCH模型.在时间序列分析中,自回归模型系数的最小二乘估计,有很多优良性质,这已经被研究得很完美了.但是,将它用于ARCH模型系数估计,这些优良性质不一定具有了.在此,我们仅指它的优缺点.其优点是:易理解,易计算;缺点是:欠精细,缺少某些优良性质.欠精细是相对极大似然估计而言的,详见后文.缺少某些优良性质,是指在使用最小二乘估计方法时,还需要条件E(yt2)2<,才具有相合性,然而此条件对ARCH模型而言,太强了.此外,使用最小二乘估计方法时,不能保证估计a*=(XtX)-1XtY的每个分量都是非负的,尽管其真值都是非负的.当然,对其它估计也存在同样的问题.在此顺便指出,为了保证估计的每个分量都是非负的,文献中有如下方法可用,即求如下的最小值问题的解,min{t=2n(yt2-a0-a1yt-12-a2yt-22-…-apyt-p2):a0>0,a10,…,ap0;(yt2-a0-a1yt-12-a2yt-22-…-apyt-p2)0}=t=2n(yt2-a0+-a1+yt-12-…-ap+yt-p2)2,(2.9)取估计*=(a0+,a1+,…,ap+)。这里叙述此方法的目的有三点可言,其一,这是最有效的保证估计的分量都是非负的;其二,有多种方法可获得ARCH模型系数的估计;其三,除了最小二乘估计,都不易计算,比如(2.9)式的求解问题,就是典型的优化求解问题,其计算的复杂性可想而知.2.2.2.极大似然估计对于序列y1,y2,…,yn,如果它们的联合分布的形式已知,其中只有有限个参数未知,那么,寻求合适的参数值,使得其分布在这些观测值y1,y2,…,yn处达到最大值,称其为极大概率估计方法.其合理性是不言而喻的.相对其它方法,可算是精细些.当然,其前提是联合分布的形式已知的.进而言之,如果已知联合分布密度函数时,使用上述的极大概率估计方法,应改为寻求合适的参数值,使得其分布密度函数在这些观测值y1,y2,…,yn处达到最大值,称其为极大概率密度估计方法.此情况有更广的应用背景,ARCH模型数估计就属于此情况.再进一步,如果已知联合分布密度函数呈现指数形式,改为寻求合适的参数值,使得其分布密度函数的对数函数(此函数被称为似然函数),在这些观测值y1,y2,…,yn处达到极大值,称其为极大似然估计方法.用极大化似然函数代替分布密度函数,只是讨论和应用时有方便之处,并无本质区别.极大化似然方法是统计学中熟知的,重要的方法.依上所述,使用极大似然估计方法,有两个关键步骤:一是,找出y1,y2,…,yn的联合分布密度函数,它仅依赖有限个未知参数,由此易得其似然函数;二是,寻找使似然函数达到极大值的参数,即参数的极大似然估计.一般说来,第一步仅是细心的推理,第二步是精心的计算,而且常常要使用近似的迭代算法.以下介绍ARCH模型参数的极大似然估计,就要对此两步作具体叙述.第一步:根据ARCH模型的假定,再使用条件概率密度的公式可得知,y1,y2,…,yn的联合分布密度函数f(y1,y2,…,yn)=f(Y),Y=(y1,y2,…,yn),有以下表达式f(Y)=f(y1,y2,…,yn)=f(yn|y1,y2,…,yn-1)f(y1,y2,…,yn-1)(依条件密度公式)=(2hn)-1/2exp{-yn2/2hn}f(y1,y2,…,yn-1)(依(1.5)式和tN(0,1),且n与{yn-1,yn-2,…}独立)=(0+1yn-12+…+pyn-p2)-1/2exp{-yn2/2(0+1yn-12+…+pyn-p2)f(y1,y2,…,yn-1)(依(1.6)式)=t=p+1n(0+1yt-12+…+pyt-p2)-1/2exp{-yt2/2(0+1yt-12+…+pyt-p2)}f(y1,y2,…,yp)(2)-(n-p)/2.(依反复递推)(2.10)记其对数函数为L()=logf(y1,y2,…,yn)=-(1/2)t=p+1n{log(0+1yt-12+…+pyt-p2)+yt2/(0+1yt-12+…+pyt-p2)}+logf(y1,y2,…,yp)+log(2)-(n-p)/2.(2.11)忽略上式中的常数项和常数因子-(1/2),再记l()=t=p+1nlt()-logf(y1,y2,…,yp),(2.12)其中lt()={log(0+1yt-12+…+pyt-p2)+yt2/(0+1yt-12+…+pyt-p2)}.显然,l()与L()只相差常数加项,所以,求解L()的最大值解,等价于求解I()的最小值解.以后我们总是考虑后者.在上述诸式中,f(y1,y2,…,yp)是y1,y2,…,yp的联合分布密度函数,为了使用极大似然估计方法,也应当将它表达成依赖于y1,y2,…,yp和0,1,…,p的明确形式.这不是一件容易的事情!仅以p=1为例,即可说明其难点所在.此时只须求出y0和y1所满足的共同的分布密度数,并使得它们满足关系式y1=(0+1y02)1/21,tN(0,1),且t与y0独立.此问题看似简单,但是很难解答.举此例的目的,还在于提请注意,当tN(0,1)时,由它驱动生成的平稳AR序列也是正态分布的,但是,由它驱动生成的平稳ARCH序列不是正态分布的.因为,在AR模型中,t以加项形式出现,在ARCH模型中,t以乘积因子形式出现(见(1.5)式).然而,在(2.10)式中,人们容易误以为f(y1,y2,…,yn)是正态分布的连乘积形式,所以它是多元正态密度.其实,因子f(y1,y2,…,yp)不是正态的密度函数(这一点容易被忽视),所以f(y1,y2,…,yn)也不是正态的.第二步:寻求极大似然估计,就是寻求使(2.11)式中L()取最大值的,等价于寻求使(2.12)式中l()取最小值的.很明显,此极大似然估计没有如同(2.8)式的显示表达式,于是只能寻找近似的数值解法.由于l()有很好的解性质,求(2.12)式中l()的最小值,可求l()/=0的解.即使如此,也还很难求解.进一步还要使用其它近似手段,即将未知项logf(y1,y2,…,yp)从(2.12)式中忽略掉,寻求t=p+1nlt()/=0(2.13)的解,其中lt()/=log(0+1yt-12+…+pyt-p2)/+{yt2/(0+1yt-12+…+pyt-p2)}/.(2.14)请注意,lt()/是向量,所以(2.13)式是(p+1)元代数方程组,利用(2.14)式很容易求得lt()/的表达式,而且lt()/有很简单的形式,但是,它是非线性的,所以(2.13)式是(p+1)元非线性代数方程组.求(2.13)式的数值解法,是计算数学中的简单问题,即可使用已有的软件,亦可自行编程计算,这里从略.2.3.GARCH模型的参数估计2.3.1.极大似然估计在这一小节,先介绍极大似然估计,因为这与前面联系紧密.如前所述,GARCH模型的参数估计,要比ARCH模型复杂.其复杂性表现在:GARCH模型的参数估计不仅没有显示的表达式,而且,其似然函数也没有显示的表达式,只有迭代 计算公式 六西格玛计算公式下载结构力学静力计算公式下载重复性计算公式下载六西格玛计算公式下载年假计算公式 .这一特点,对求解极大似然估计的算法,不带来实质困难,但是在叙述它时,会繁琐些.现在叙述GARCH模型似然函数.仿照(2.10)式可得f(Y)=f(y1,y2,…,yn)=f(yn|yn-1,..,yn-p;hn,…,hn-q+1)f(yn-1,..,yn-p;hn,…,hn-q+1)=(2hn)-1/2exp{-yn2/2hn}f(yn-1,..,yn-p;hn,…,hn-q+1)=t=1n(2ht)-1/2exp{-yt2/2ht}f(y0,..,y-p+1;h1,…,h-q+2).(2.15)仿照(2.12)式又有I()=t=p+1nlt()+logf(y0,..,y-p+1;h1,…,h-q+2),(2.16)其中lt()=loght+yt2/ht,=(0,1,…,p;1,…,q).再仿照(2.13)式和(2.14)式,在求解GARCH模型参数的极大似然估计时,近似为求解如下的方程组之解,即t=p+1nlt()/=0(2.17)的解,其中lt()/=loght/+(yt2/ht)/=ht-1(ht/)(1-yt2/ht).(2.18)在以上各式中,虽然都是明确的表达式,但是,(ht/)尚未被表达出来,实际上无法用显式表达它.幸运的是,它有递推关系式可利用.在设计求解方程(2.17)式时,有递推关系式也足够了.记zt=(1,yt-12,yt-22,…,yt-p2;ht-1,ht-2,…,ht-q),于是可得出(ht/)的递推关系式如下ht/=zt+k=1qk(ht-k/).(2.19)虽然有(2.19)式可用,但是,此迭代公式的初始值h1/,h0/,…,h-q+2/仍然未有明显表达式.在实际应用时,常用零值作为它们的近似值使用,于是可求得近似的极大似然估计值.当然,求解过程又常用迭代算法,这里从略.2.3.2.最小二乘估计对GARCH模型参数使用最小二乘估计方法,也同样遇到像极大似然估计类似的麻烦.在此,我们推荐使用平稳的ARMA模型参数的矩估计方法.细节可参看有关著作.尽管如此,当q值较大时,其算法也不比极大似然估计更方便.所以,最多使用的仍是极大似然估计方法.模型检验根据观测数据y1,y2,…,yn,判断所要拟合的模型是否适用,称为模型检验.在为数据y1,y2,…,yn建立模型时,一般都应当进行模型检验.对于GARCH模型也不例外.所谓模型检验,有在建立模型前进行的,有在之后进行的.对于GARCH模型来说,在为数据y1,y2,…,yn建立GARCH模型前,首先应当判断有没有必要.如前言所说到,平稳序列的条件方差S(yt,yt-1,…)可能是常数值,此时就不必建立GARCH模型.于是判断条件方差S(yt,yt-1,…)是否为常数,就应当在建模前完成.即使经判断后,条件方差不是常数,它也未必满足GARCH模型.然而目前GARCH模型是比较熟知的条件异方差模型,所以常用它来近似拟合观测数据.那么,在建模后还应当对所得到的模型进行检验,以判断其是否可接受.在建模前和后所进行的模型检验,其方法不一定相同.建模后使用的模型检验方法,还可作为确定GARCH模型阶数的辅助手段.以下分别介绍.3.1.条件异方差性检验在这一小节里,我们仍考虑{yt}为鞅差序列的情况,也就是(0.12)式中的(yt-1,yt-2,…)=0的情况,即yt=S(yt-1,yt-2,…)t,(3.1)其中{t}为标准化的鞅差序列,即E(t|yt-1,yt-2,…)=0,E(t2|yt-1,yt-2,…)=1.考查(3.1)式两边平方的模型yt2=S2(yt-1,yt-2,…)t2,当S(yt-1,yt-2,…)为常数时,不妨记为,即E{yt2|yt-1,yt-2,…}=S2(yt-1,yt-2,…)=2.(3.2)此时(3.1)式可写成yt2=2t2.(3.3)于是又有yt2-2=2(t2-1).(3.4)此时我们还发现E{(yt2-2)|yt-1,yt-2,…}=E{yt2|yt-1,yt-2,…}-2=2-2=0,这说明{yt2-2}也是鞅差序列.还容易看出,如果{yt}是任意一个鞅差序列,且Eyt2=2,{yt2-2}未必是鞅差序列.但是从上式不难看到,当且仅当(3.2)式成立时,{yt2-2}才是鞅差序列.此事实是进行条件异方差性检验的理论依据.以下介绍具体检验方法.计算n2=(1/n)t=1nyt2.计算n=(1/n)t=1n(yt-1,…,yt-m),2=(1/n)t=1n{(yt-1,…,yt-m)-n}2,其中是m元标准正态分布的密度函数,m<0,i0,i=1,2,…,p.其中{t}为i.i.d.的序列,tN(0,1),且t与{et-1,et-2,…}独立.在下文中,我们将叙述如何为联合模型(4.1)-(4.3)(4.4),以及联合模型(4.2)-(4.3)(4.4)建模.前者称为:线性回归---ARCH模型,后者称为:自回归---ARCH模型,并简记为AR---ARCH模型.为以上联合模型建模时,一种简便的方法是两步法.对联合模型(4.1)-(4.3)(4.4)式而言,在已获得观测数据(yt,xt1,xt2,…,xts),t=1,2,…,n,为了对(4.1)-(4.3)(4.4)式建模,分两步完成:第一步,按照线性回归分析方法,给出回归系数估计a0*,a1*,…,as*,(常用最小二乘估计)再給出拟合残差序列et*=yt-(a0*+a1*xt1+a2*xt2+…+as*xts),t=1,2,…,n.第二步,将序列e1*,e2*,…,en*看做观测e1,e2,…,en的近似值.按照前面为y1,y2,…,yn建立ARCH模型的方法,用于对e1,e2,…,en的分析,并获得模型(4.3)(4.4)式的估计.必要时,还可作模型检验.以上就是为(4.1)-(4.3)(4.4)式建模的两步法.当然,我们也可以使用熟知的极大似然估计方法为(4.1)-(4.3)(4.4)式建模.但是,其复杂程度和计算难度明显高于上述方法.只有当模型阶数不太高时才会使用.按照前面叙述的极大似然方法的要领,除了繁琐外并无实质性困难.但是,无法再利用前述的GARCH模型建模方法.对于联合模型(4.2)-(4.3)(4.4)式而言,在已获得观测数据y1,y2,…,yn,为(4.2)-(4.3)(4.4)式建模也分两步完成:第一步,按照线性自回归分析,给出自回归系数估计a0*,a1*,…,as*,(常用最小二乘估计)再給出拟合残差序列et*=yt-(a0*+a1*yt-1+…+as*yt-s),t=s+1,s+2,…,n.第二步,将序列es+1*,es+2*,…,en*看做观测es+1,es+1…,en的近似值.按照前面为y1,y2,…,yn建立ARCH模型的方法,用于对es+1,es+2,…,en的分析,并获得模型(4.3)(4.4)式的估计.需要时,还可作模型检验.以上就是为(4.2)-(4.3)(4.4)式建模的两步法.当然,对于(4.2)-(4.3)(4.4)式建模也可使用极大似然估计方法,其情况与(4.1)-(4.3)(4.4)式类似.最后,按照类似的推广方法,也可以讨论(自)回归与GARCH模型的联合模型.除了增加复杂性以外,并无本质困难,这里不再介绍了.4.2.在区间预报中的应用在自回归分析中,预报是重要的应用之一.除了数值预报外,有时还需要給出预报值的区间,称为区间预报.在经典自回归分析中,其预报值与预报步数和历史数据有关,但是,其预报区间只与预报步数有关,和历史数据无关.其根原在于,未考虑序列的条件异方差性.近年来,由于条件异方差模型的问世,比如GARCH模型,对区间预报提供了改进的方法.现在就介绍用GARCH模型给出更合理的区间预报方法.先回顾一下经典的区间预报方法.不妨只考查一步预报.以平稳AR(s)序列{yt}为例,它满足如下模型:yt=a1yt-1+a2yt-2+…+asyt-s+et,(4.5)其中{et}为i.i.d.序列,且et~N(0,2),et与{yt-1,yt-2,…}独立.其一步预报值为yn+1|n=a1yn+a2yn-1+…+asyn-s+1,(4.6)其预报误差及其方差分别为yn+1-yn+1|n=a1yn+a2yn-1+…+asyn-s+1+en+1-(a1yn+a2yn-1+…+asyn-s+1)=en+1,(4.
本文档为【GARCH模型与应用简介】,请使用软件OFFICE或WPS软件打开。作品中的文字与图均可以修改和编辑, 图片更改请在作品中右键图片并更换,文字修改请直接点击文字进行修改,也可以新增和删除文档中的内容。
该文档来自用户分享,如有侵权行为请发邮件ishare@vip.sina.com联系网站客服,我们会及时删除。
[版权声明] 本站所有资料为用户分享产生,若发现您的权利被侵害,请联系客服邮件isharekefu@iask.cn,我们尽快处理。
本作品所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用。
网站提供的党政主题相关内容(国旗、国徽、党徽..)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
下载需要: 免费 已有0 人下载
最新资料
资料动态
专题动态
个人认证用户
言言无悔一生
暂无简介~
格式:doc
大小:70KB
软件:Word
页数:49
分类:
上传时间:2022-01-21
浏览量:1