自回归移动平均模型解析.doc
《自回归移动平均模型解析.doc》由会员分享,可在线阅读,更多相关《自回归移动平均模型解析.doc(29页珍藏版)》请在咨信网上搜索。
1、(完整word)自回归移动平均模型解析第二章 自回归移动平均模型一些金融时间序列的变动往往呈现出一定的平稳特征,由Box和Jenkins创立的ARMA模型就是借助时间序列的随机性来描述平稳序列的相关性信息,并由此对时间序列的变化进行建模和预测。第一节 ARMA模型的基本原理ARMA模型由三种基本的模型构成:自回归模型(AR,Auto-regressive Model),移动平均模型(MA,Moving Average Model)以及自回归移动平均模型(ARMA,Autoregressive Moving Average Model).2.1。1 自回归模型的基本原理1AR模型的基本形式AR模
2、型的一般形式如下:其中,c为常数项, 模型的系数,为白噪声序列。我们称上述方程为阶自回归模型,记为AR()。2AR模型的平稳性此处的平稳性是指宽平稳,即时间序列的均值,方差和自协方差均与时刻无关。即若时间序列是平稳的,即,。为了描述的方便,对式(2.1)的滞后项引入滞后算子。若,定义算子“L”,使得,L称为滞后算子。由此可知,。对于式子(2。1),可利用滞后算子改写为:移项整理,可得:AR()的平稳性条件为方程的解均位于单位圆外。3AR模型的统计性质(1)AR模型的均值。假设AR()模型是平稳的,对AR()模型两边取期望可得:根据平稳序列的定义知,,由于随即干扰项为白噪声序列,所以,因此上式可
3、化简为:所以,(2)AR模型的方差。直接计算AR()模型的方差较困难,这里引入Green函数。AR()模型可以改写成如下形式:设为平稳AR()模型的反特征根,则。进一步,其中,为常数,,称为Green函数,因为均在单位圆内,所以Green函数是呈负指数下降的.对上式两边取方差,可得:由于随机干扰项为白噪声序列,所以。因为Green函数是呈负指数下降,所以,这说明平稳时间序列方差有界,且等于常数.(3)自协方差函数。假设将原序列已经中心化,则,则对AR()模型等号两边同时乘以,两边取期望得: 因为当期的随机干扰项与过去的时间序列值无关,所以:。因此,上式可以化为:其中,表示阶自协方差。2.1。2
4、 移动平均模型的基本原理1MA模型的基本形式MA模型的一般形式如下:其中,为常数项,为模型的系数,为白噪声序列。我们称上述方程为阶移动平均模型,记为MA().2、MA模型的可逆性对于一个MA()模型:将其写成滞后算子的形式: 若方程的根全部落在单位圆外,则称MA模型是可逆的。可逆性可以保证MA模型可以改写成:即MA模型可以转化为AR模型,同时可以保证参数估计的唯一性。3、MA模型的数字特征(1)均值当时,对于一般的MA(q)模型:两边取期望,可得:即一般的MA(q)模型的期望值即为模型中的常数项.(2)方差对MA(q)模型,两边取方差:(3)协方差函数 化简可得:2.1。3 自回归移动平均模型
5、的基本原理1、ARMA模型的基本形式ARMA模型的一般形式如下:显然ARMA(p,q)模型可看成是AR(p)模型和MA(q)模型相结合的混合形式。2、ARMA模型的平稳性和可逆性对于一个ARMA(p,q)模型,将其写为滞后算子的形式:两边同时除以其中:由此可以看出,ARMA模型的平稳性完全取决于AR(p)模型的参数,与MA(q)模型的参数无关。类似地,ARMA模型的可逆性完全取决于MA(q)模型的参数,与AR(p)模型的参数无关.3、ARMA模型的数字特征(1)期望对于一个一般的ARMA(p,q)模型两边同时取期望,化简得:(2)自协方差函数 第二节 时间序列的相关性分析与平稳性2。2.1 时
6、间序列的自相关系数2.2。1.1 自相关函数(ACF)1、AR(p)的自相关函数在上一节中已经介绍了AR(p)模型的协方差函数满足下式:由于自相关系数,因此:该式表示自相关系数满足p阶差分方程。根据差分方程解的性质,上差分方程的通解可以写为:其中,为任意不全为0的常数,是滞后多项式的反特征根.根据平稳性的性质,.从自相关系数的一般形式可看出,始终不为0,但是随着滞后阶数的增加,自相关系数慢慢逼近0,在图形上表现出一定的拖尾性。2、MA模型的自相关函数根据上一节推导的MA模型的自协方差函数的表达式,MA模型的自相关函数表示为:因此,当kq时,自相关函数为0,也就是说MA(q)模型的自相关函数在q
7、步以后是截尾的。3、ARMA模型的自相关函数根据ARMA模型的自协方差函数,不难得到ARMA模型的自相关函数:由此可以看出,ARMA模型的自相关函数不具有截尾性.事实上,ARMA模型若满足可逆性,其形式相当于一个无穷阶的AR模型,因此自相关函数与AR模型一样具有拖尾性。2。2。1。2 偏自相关函数(PACF)1、偏自相关函数的定义自相关函数不能纯粹地表示与之间的相关性,两者的相关性还会受到、的间接影响,为了单纯地表示与之间的相关性,这里引入偏自相关函数。偏自相关函数表示在固定、的情况下与之间的相关性.下面介绍偏自相关函数的计算方法。设序列yt可由下回归方程估计:根据回归方程的性质,式中估计系数
8、即为偏自相关函数。为了估计回归系数,采用OLS方法,即达到最小。对L关于各回归系数求偏导,可得到以下方程组:该方程组称为Yule- Wolker方程。根据自相关系数,求解Y-W方程即可得到偏自相关系数.2、AR(p)的偏自相关函数对于AR(p)模型,时,由于与序列的滞后项无关,因此,且当由此,AR(p)模型的偏自相关函数在kp后等于0,即AR(p)模型的偏自相关函数具有截尾性。事实上,AR模型偏相关函数的截尾性也可直接从该模型的表达式看出。AR(p)模型实质上假设序列至多只与滞后p阶的值相关,因此偏自相关函数至多在p阶处非0。3、MA(q)和ARMA(p,q)的偏自相关函数由于MA(q)和AR
9、MA(p,q)相当于无穷阶的AR模型,因此这两个模型的偏自相关函数均不具有截尾性,而是拖尾性.2。2.1.3 ARMA模型自相关系数与偏自相关系数的估计与检验根据以上分析,不同ARMA模型自相关系数与偏自相关系数的表现存在明显的差异。表2。1给出了三类模型ACF与PACF的特征。 表2.1 ARMA类模型ACF与PACF的特征模型自相关系数偏自相关函数AR(p)拖尾p阶截尾MA(q)q阶截尾拖尾ARMA(p,q)拖尾拖尾因此,我们可以通过观察偏自相关函数来识别并确定AR模型的滞后阶数,通过自相关函数来识别并确定MA模型的滞后阶数q。那么对于给定的样本数据,如何估计样本ACF与PACF,并从统计
10、角度检验两者是否为0呢?下面分别介绍ACF与PACF的估计与检验。1、样本ACF与PACF的估计与实现对于给定样本,只需估计样本的自协方差与方差,将两者相除即可得到样本ACF.具体而言,样本自协方差表示为:其中表示样本均值。那么SACF=.对于PACF主要是利用Yule- Wolker方程求解.当滞后阶数较大时,YW方程直接计算较难,目前多采用递推算法来求解。2、样本ACF与PACF的显著性检验常用的检验方法主要包括两类:正态检验法和Portmanteau检验法。若序列满足独立性,则由统计渐进分布的有关定理可知,当样本个数充分大时,ACF和PACF均满足均值为0,方差1/T的正态分布,即,。因
11、此若,则可认为样本数据是独立的,即自相关系数和偏自相关系数均不显著异于0.该检验法即为正态检验法。Portmanteau检验法是联合检验法,即检验直到k阶的自相关系数是否同时为0.该检验法使用Q统计量进行检验。Q统计量具体形式为:其中T为样本容量,k为设定的滞后阶数。Q统计量服从分布.当Q统计量超过设定的临界值时,就拒绝原假设,即序列至少存在k阶以内的自相关性。2。2。2 时间序列平稳性检验建立ARMA的前提是序列是平稳的。检验平稳性常用的方法主要有三种:经验法、自/偏自相关系数法、单位根检验法。1、经验法经验法是通过观察图形的方式来初步判断时间序列是否平稳的。首先画出时间序列的图形,如果该图
12、形围绕某一直线上下以较小的幅度波动,则该序列一般是平稳的,否则是不平稳的.2、自/偏自相关系数法由于ARMA模型的自/偏自相关系数要么是截尾的,要么是拖尾的,因此可以观察时间序列的自/偏自相关图,如果时间序列的自/偏自相关系数从某个滞后期开始均与0无差异,可以认为该时间序列是平稳的;若自/偏自相关系数衰减很慢,且与0存在明显的差异,则时间序列是非平稳的。3、单位检验法常用的单位根检验法主要包括DF检验法和ADF检验法。(1)DF检验法DF检验包括三种形式:其中,c为常数项,表示线性趋势,随机干扰项独立同分布,且服从。根据平稳性的概念,若序列yt是不平稳的,则回归系数。一般较易识别。因此判断序列
13、yt是否平稳,主要是判断是否为1。如果,则说明序列存在单位根,是不平稳的,否则是平稳的。进一步,上述三个方程两边同时减去,得:其中,因此可以将DF检验的原假设和备择假设分别为:相应的统计量为:DF的形式与t统计量相似,但是该统计量并不服从t分布,Dickey和Fuller(1979)给出了利用蒙特卡罗模拟方法模拟的临界值,因此该检验称为DF检验。DF检验是左侧检验,且不同形式的方程临界值是不同的。注意DF检验只有当时间序列为AR(1)过程时才有效。如果存在高阶滞后相关,那么将违背随机干扰项独立同分布的假设。因此,DickeyFuller提出来ADF检验来弥补DF检验的不足。(2)ADF检验假设
14、时间序列存在p阶自相关,那么用p阶自回归方程来判断单位根,形式为:上式两边同时减去,通过整理可得:其中 ,,上述检验形式是在DF检验方程中加入了yt的高阶滞后项,因此可以看成是DF检验的增广形式,简称ADF检验。与DF检验类似,ADF检验也存在三种形式,不难看出,当p=1时,ADF检验就是DF检验,因此DF检验是ADF检验的特例。检验的原假设和备择假设为:即原假设为序列至少存在一个单位根,备择假设为序列不存在单位根。使用ADF检验时,应该注意如下几个问题:首先,要确定合理的滞后阶数。其次,因为检验统计量的临界值依赖于方程的形式,因此选择检验的方程形式很重要。再者,如果检验的结果是拒绝原假设,那
15、么原序列就不存在单位根,即原序列是平稳的;如果接受原假设,则序列是不平稳的,需要进行若干次差分,直到拒绝原假设,从而确定序列单整的阶数。例 2.1 以2003年2月到2010年4月上证国债(交易代码000012)月末收盘指数为原始数据进行分析。首先对指数序列取对数,然后对这个收益率序列进行单位根检验,结果如下: 表2.2 上证综指收益率单位根检验表Null Hypothesis: LSP has a unit rootExogenous: Constant, Linear TrendLag Length: 0 (Automatic based on SIC, MAXLAG=11)tStatis
16、ticProb。*Augmented Dickey-Fuller test statistic2。0279890。5777Test critical values:1 level4。0682905% level-3.46291210% level-3.157836从检验结果来看,应该接受原假设,即认为该序列存在单位根.第三节 季节性ARMA模型2。3。1 时间序列的季节性在现实中,存在一些经济时间序列具有明显的季节性,比如GDP、消费支出等.它们具有共同的特点:有规律,每年重复出现,表现为逐年同时期有相同的变化方向和大致相同的变动幅度。一般我们认为季节性是一种对数据的干扰,在数据分析中通常将其
17、过滤.那么该如何识别时间序列的季节性。常用的分析方法是自相关图分析法。观察时间序列的自相关图,如果是月度数据,可以观察滞后期为12、24、36的自相关系数;如果是季度数据,可以观察滞后期为4、8、12的自相关系数。如果自相关系数与0无差异,那么同月或者同季度之间不相关,即不存在季节性,否则该时间序列具有季节性.值得注意的是,若时间序列的趋势性较强时,这种趋势性会掩盖季节性,因此分析季节性时,应该先消除时间序列的趋势性。由于季节性是对数据的一种干扰,因此在数据分析中常常消除季节性。常用的处理季节性的方法主要包括三种:(1)移动平均法;(2)用X11方法去掉时间序列的季节性;(3)采用逐期差分和季
18、节差分,建立SARMA模型。下面先简单介绍一下移动平均法和X-11方法,再重点介绍SARMA模型。1、移动平均法移动平均法是将原时间序列的两个或多个时期的数据进行平均,用平均值代替原序列值,以此来平滑周期性。假设周期长度为s,则在时点t上的移动平均值可以表示为:2、X11方法X-11方法最初是由美国人口普查局在1965年研究开发出来的方法。它的基本原理是滑动平均法法.这种方法能适应各经济指标的性质,并根据季节调整的目的进行相应的调整及进行算法的选择。X11方法包括两种模型:加法模型和乘法模型。加法模型的基本形式为:乘法模型的基本形式为:其中Tt为趋势循环要素,St为季节要素,It为不规则变动要
19、素,Dt为周工作日变动要素。常用的是乘法模型,但注意乘法模型仅适用于时间序列数据为正的情况。 例 2。2 以江苏省2001年1月到2009年12月的每月的社会消费品零售总额为原始数据,进行季节性分析。数据来源于中经网统计数据库。数据的时间序列图如下。从图中可以看出,该序列存在一定的周期性季节性变动。采用X-11方法进行季节性调整,结果如下:从图中可以看出,经过调整,季节性因素与不规则因素得到了一定程度的消除.2.3。2 季节性ARMA模型SARMA模型是在回归项中加入季节自回归项(SAR)和季节移动平均项(SMA)来考虑季节性的.假设季节长度为S,则SMA(q)(m)S的一般形式为:其中q为一
20、般移动平均项的滞后阶数,m为季节性滞后阶数.SAR(p)(k)S的一般形式为:其中p为一个自回归项的滞后阶数,k为季节性滞后阶数.SARMA的一般形式只需将SAR和SMA合并即可,简记为:SARMA(p,q)(k,m)S。为了进一步理解季节模型的特点,以滞后两阶的AR模型为例,无季节性的AR(2)模型如下:用滞后算子表示:对于带有季节性的季度数据,若k=1,则模型就会变成下式:它等价于:因此该季节性AR模型实质上具有稀疏系数的更高阶AR模型。SARMA模型的自相关系数与偏自相关系数与普通的ARMA模型类似,只是ACF和PACF的计算需要计算滞后期分别为S、2S.。mS处的大小,性质与ARMA模
21、型相同,即季节自回归模型的季节偏相关函数按季节周期的增加而截尾;季节移动平均模型的季节自相关函数按季节周期增加而截尾;混合的季节自回归移动平均模型的季节自相关函数和偏相关函数均按指数衰减。例 2。3 引用例2.2的数据进行季节性分析。为了便于分析,首先将数据取自然对数,然后进行差分去除趋势性。下图是新序列的自相关图,结果如下:从图中可以看出,该序列在滞后6期、12期处的自相关函数和偏自相关函数显著异于0,因此可认为该序列存在一个周期为6的季节性变动。建立ARMA模型时,可在模型中加入周期为6的sar或sma项。第四节 ARMA模型的构建与Eviews实现2。4。1 ARMA模型的具体构建步骤A
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 回归 移动 平均 模型 解析
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【天****】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【天****】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。