Lecture04一元时间序列分析方法

《Lecture04一元时间序列分析方法》由会员分享，可在线阅读，更多相关《Lecture04一元时间序列分析方法（46页珍藏版）》请在装配图网上搜索。

1、第四章第四章 一元时间序列分析方法一元时间序列分析方法l 学习目标：学习目标：了解平稳性和白噪声过程；熟悉随机序列模型；熟悉ARIMA过程；掌握时间序列的平稳性和单位根检验。第四章第四章 一元时间序列分析方法一元时间序列分析方法l第一节第一节 时间序列的相关概念时间序列的相关概念l第二节第二节 随机序列模型随机序列模型l第三节第三节 单整自回归移动平均模型单整自回归移动平均模型l第四节第四节 平稳性与单位根检验平稳性与单位根检验时间序列的相关概念l第一节 时间序列的相关概念一、平稳性l平稳性是时间序列分析的基础。判断一个序列平稳与否非常重要，因为一个序列是否平稳会对它的行为及其性质产生重要的影

2、响。在时间序列平稳性，一般包括下列两类平稳过程：l1、严格平稳过程（Strictly Stationary Process）如果对所有的t，任意正整数n和任意n个正整数()，（）的联合分布与()的联合分布是相同的,即：1,ntt1,nttyy1,ntmtmyy1111,tt nntmtn mnP ybybP ybyb时间序列的相关概念l2、弱平稳性过程（Weakly Stationary Process）如果一个时间序列 的均值，方差在时间过程上保持是常数，并且在任何两时期之间的协方差值仅依赖于该两时期间的距离或滞后，而不依赖于计算这个协方差的实际时间，则称时间序列 是弱平稳的。弱平稳的时间序

3、列有如下性质：可见，如果一个时间序列概率分布的所有阶矩都不随时间变化，那它就是严格平稳的；而如果仅仅是一阶矩和二阶矩（即均值和方差）不随时间变化，那它就是弱平稳的。ty ty2()()ttE yuyu()tE yuuE yytttt()()122112,tt时间序列的相关概念l二、自协方差（auto-covariance）决定 是如何与它自身的先前值相关的，对于一个平稳的时间序列，它只依赖于 与 之差。其中，被称为自协方差函数。另一种更为简洁的方法使用自相关系数来描述他们之间的关系。考虑弱平稳时间序列 ，当 与它的过去值 线性相关时，可以把相关系数的概念推广到自相关系数，与 的相关系数称为 的

4、间隔为 的自相关系数，通常记为 ，在弱平稳性的假定下它只是的函数，定义tyty1ty()(),0,1,2ttt st ssE yE yE yE ys tytyt lytyt lytylll(,)()()tt ltt lCov y yVar y Var y0(,)()tt lltCov y yVar y=时间序列的相关概念l三、白噪声过程如果时间序列 是一个有有限均值和有限方差的、独立同分布的随机变量序列，则称时间序列 为白噪声。特别的，若时间序列还服从均值为0，方差为 的正态分布，则这个序列称为高斯白噪声。它是其它各类型时间序列的重要组成部分，在金融市场效率理论中具有重要的意义。对于白噪声序列

5、，自相关系数为零。在实际应用中，如果所有样本的自相关函数接近为零，则认为这个序列为白噪声序列。若一个随机过程满足：则我们称之为白噪声过程（white noise process）。tyty2()tE y2()tVar y2tr0trt r若若随机序列模型l第二节 随机序列模型若对每一个固定的t，是一个随机变量，则 ，为随机时间序列。而揭示随机时间序列自身变化规律和相关关系的数学表达式就是时间序列分析模型。随机时间序列分析模型分为三类：自回归模型（auto-regressive model,AR）、移动平均模型（moving-average model,MA）和自回归移动平均模型（auto-re

6、gressive moving average model,ARMA）。对于任一个时间序列，怎样判断它是遵循纯AR过程（若是的话，阶数p取什么值），纯MA过程，（若是的话，阶数q取什么值）或是ARMA模型，此时p和q各取多少。我们将遵循以下四个步骤对这三个模型做一详细介绍：ty1y2yty随机序列模型l步骤一：识别。步骤一：识别。就是找出适当的p和q值。我们即将说明怎样借助相关图和偏相关图来解决此类问题。l步骤二：估计。步骤二：估计。一旦辨别适当的p和q值，下一步便是估计模型中所含自回归和移动平均项的参数。l步骤三：诊断。步骤三：诊断。选定模型并估计其参数之后，下一步就要看所选的模型对数据拟合

7、的是否够好。对所选模型的一个简单的检验，是看从该模型估计出来的残差是不是白噪声；如果是，就可接受这个具体的拟合；如果不是，我们必须重新在做。l步骤四：预测。步骤四：预测。ARMA建模方法之所以得以普及，理由之一是它在预测方面的成功。有许多事例用这个方法做出的预测比用传统的计量经济建模方法做出的预测更为可靠，特别是在短期预测方面。随机序列模型l一、自回归模型（AR）若一个时间序列可表示为 （4.12）其中，为白噪声，则称 为一阶自回归过程，或简称为 。自回归模型是时间序列 表示为它的先前值与一个误差项 的线性函数。在p阶自回归中，、，是自回归参数，它表明每改变一个单位时间值时，对 所产生的影响，

8、它是根据样本观测值来估计的参数。011tttyy t()0tE2var()ttr (1)ARtyt12pty随机序列模型l2、AR模型阶的识别在实际应用中，一个AR时间序列的p阶是未知的，必须根据实际情况来决定。这个问题叫做AR模型的阶的决定。一般可以通过两种方法：第一种方法是利用偏自相关函数（partial autocorrelation function,PACF），第二种方法是用某个信息准则函数。（1）偏自相关函数（PACF）偏自相关就是 和 之间的，除去居中的诸 （即 ）的影响后的相关。其相关程度可用偏自相关系数 度量。进行回归对一个 模型，间隔为的样本偏自相关系数不应为零，而对所有

9、，应接近零，我们利用这一性质来决定p阶。tyt lyy12,1,ttt lyyy,l l0111(1),1,2,ttltll ltltyyyytT()AR pjp,j j随机序列模型l（2）采用信息准则法判别模型阶数在实际应用中，很难利用自相关函数来确定模型的合理阶数。较为简便的方法是，所选定的阶数应使得信息准则的数值达到最小。对于信息准则，一般应用赤池（Akaike）准则信息准则（AIC）和许瓦兹（Schwarz）贝叶斯信息准则(SBIC)。随机序列模型l3 3、参数估计、参数估计对一个 模型，我们常用条件最小二乘法来估计其参数，条件最小二乘是从第 个观测值开始的。l4、模型验证、模型验证对

10、实际数据所时拟合的模型，要仔细地验证它的合理性。若模型是合理的，其残差序列应该是白噪声。残差的样本自相关函数和Ljung-Box统计量可用来检验 与一个白噪声的接近程度。对 模型，Ljung-Box统计量 渐进服从自由度为m-p的 分布。如果所拟合的模型经经验验证是不合理的，那么就需要对它进行修正。()AR p1p t()AR p()Q m2随机序列模型l5、预测、预测预测是时间序列分析的一个重要应用。向前一步预测向前两步预测向前多步预测随机序列模型l6、判定预测是否精确、判定预测是否精确l在实际中应用中，通常是对整个样本外的区间进行预测，然后将其与实际值比较，把他们之间的差异用某种方法加总。

11、对第i个观测值的预测误差定义为其实际值和预测值之间的差值，再求其平方或取其绝对值使各项为正后进行加总。随机序列模型l案例说明案例说明4-1上证指数收益率的上证指数收益率的AR建模建模l本案例数据来自高铁梅（2006）计量经济分析方法与建模，数据选取了上证收盘指数(1991年1月2003年3月)的月度时间序列S作为研究对象,用AR(1)模型描述其变化规律。l在此，对其做变化率 ，这样便得到了变化率序列。一般来讲，股价指数序列并不是一个平稳的序列，而通过变化后的变化率数据，是一个平稳序列，可以作为我们研究、建模的对象。l对上证收益率数据拟合。在此，记上证股价指数变化率序列为sr，建立如下模型：11

12、100/(1,2,),ttttsrSSStT1tttsrcsru 1,2,tT案例说明4-1上证指数收益率的AR建模图图4-2：AR(1)回归结果回归结果案例说明4-1上证指数收益率的AR建模-4004080120160200919293949596979899000102SRSR F图图4-3：上证指数收益率序列及其拟合值：上证指数收益率序列及其拟合值l在图4-3中，实线是上证指数变化率序列，虚线是AR(1)模型的拟合值。从该图可以看出我国上证股价指数变化率序列在1991-1994年之间变化很大，而后逐渐趋于平稳。近年来，波动平缓，并且大多在3%下面波动。拟合曲线基本代表了这一时期的均值。随

13、机序列模型l案例说明4-2应用AR（1）进行预测下面，我们利用建立的AR(1)模型进行预测。我们选取2000年1月至2006年6月的我国广义货币供应量（M2）月度数据的时间序列，进行AR（1）建模并预测。1000001500002000002500003000003500004000002000200120022003200420052006M2FForecast:M2FActual:M2Forecast sample:2000M01 2006M07Adjusted sample:2000M02 2006M07Included observations:78Root Mean Squared

14、Error 3494.959Mean Absolute Error 2842.679Mean Abs.Percent Error 1.555453Theil Inequality Coefficient 0.008101 Bias Proportion 0.467290 Variance Proportion 0.048130 Covariance Proportion 0.484580图图4-5:利用利用AR模型进行预测模型进行预测随机序列模型l二、移动平均模型（二、移动平均模型（MA）若一个随机过程 可为下面形式：（4.40）则称方程式（4.40）表示的是q阶的移动平均过程(moving

15、average)，表示为 。在 模型中，为参数，为白噪声过程。最简单的移动平均过程是 ，可表达为：tx011tttqt qyc ()MA q()MA q12,q t(1)MA1(1)ttyL()MA q随机序列模型l1、MA模型阶的识别自相关函数是识别MA模型的阶的有用工具。一个时间序列 具有自相关函数 ，若 但对 有 ，则 服从一个 模型。l2、MA模型估计模型估计l估计MA模型通常用最大似然法。有两种方法求MA模型的似然函数。第一种是条件似然法，即假定初始的扰动（即 ，）都是0；这样 ，计算似然函数所需要的抖动可以递推得到。第二种方法是把初始抖动 ，当作模型的附加参数与其它参数一起估计出来

16、。tyl0llp0lty()MA qt0t 110yc2201 1,yc t0t 随机序列模型随机序列模型l3、MA模型预测模型预测l由于MA模型有有限记忆性，它的点预测很快可以打到序列的均值。设预测原点为 ,对MA（1）过程的向前1步预测，模型为l l取条件期望，我们有l向前1步预测误差的方差为 。h1011hhhyc 1101(1)(,)hhhhhrEyyyc 11(1)(1)hhhheyy2(1)hVar e随机序列模型l三、ARMA模型自回归模型和移动平均模型是时间序列中最基本的两种模型类别，将这两种基本的模型类别结合起来，就产生了自回归移动平均模型（ARMA）。若一个时间序列 可表示

17、为：(4.51)或者表达为：(4.52)则称时间序列模型为自回归移动平均模型，表示为 。在模型中，和 分别表示为滞后之后p和q阶的表达式，并称其为自回归算子和移动平均算子。tx1111ttptpttqt qyyy 11(1)(1)pqptqLLyLL()()ttL yL(,)ARMA p q()L()L随机序列模型l案例说明4-3应用Eviews建立ARMA模型的实例以中国联通（600050）为例使用的数据为联通公司股票的日股价序列，期限为2003-1-2日至2006年9月15日，共886个样本观测值。该模型涉及三个步骤：识别、估计和诊断性检验。首先，通过观察自相关系数，对数据结构加以识别。(

18、1)估计高达阶的自相关系数(2)采用信息准则法判别模型阶数随机序列模型l案例4-4在Eviews中运用ARMA模型进行预测 一旦选定了模型的阶数并利用一定的数据完成了模型估计后,就可以利用该模型对序列的未来值进行预测。估计出所需模型并在Eviews中打开输出结果窗口后，点击Forecast图标。Eviews使用两种预测方法：动态预测和静态预测。动态选项从预测第一项开始，计算多步预测值；而静态选项则计算一系列向前一步的预测值，此时每生成一个预测值，就将样本范围向前移动一个观测值，以便将真实值而非预测值作为滞后因变量。在此，我们应用了上证指数1990年12月19日至2006年8月31日时间区间，共

19、3856交易日的样本数。采用Eviews软件的Froecast预测模型，动态预测和静态预测的结果分别如下：随机序列模型-.06-.04-.02.00.02.04.06.08.10500100015002000250030003500Forecast:SZ_RETURNFActual:SZ_RETURNForecast sample:1 3856Adjusted sample:2 3856Included observations:3854Root Mean Squared Error 0.026716Mean Absolute Error 0.014488Mean Abs.Percent Er

20、ror 118.9587Theil Inequality Coefficient 0.962871 Bias Proportion 0.000001 Variance Proportion 0.949712 Covariance Proportion 0.050287(a)动态预测随机序列模型-.08-.04.00.04.08.12500100015002000250030003500Forecast:SZ_RETURNFActual:SZ_RETURNForecast sample:1 3856Adjusted sample:2 3856Included observations:3854R

21、oot Mean Squared Error 0.026621Mean Absolute Error 0.014431Mean Abs.Percent Error 133.6966Theil Inequality Coefficient 0.909112 Bias Proportion 0.000000 Variance Proportion 0.834198 Covariance Proportion 0.165801（b）静态预测单整自回归移动平均模型单整自回归移动平均模型l第三节 单整自回归移动平均模型单整自回归移动平均模型（autoregressive integrate moving

22、 average models，ARIMA）最先由博克斯（Boy）和詹金斯（Jenkins）在1976年提出的。该模型是指将非平稳时间序列转化为平稳时间序列，然后将因变量仅对它的滞后值及其随机误差项的现值和滞后值进行回归所建立的模型。目前，该模型已经在众多领域和研究中得到应用，并证明了其较强的解释力和适应性。单整自回归移动平均模型l一、ARIMA模型介绍假定存在一个随机过程含有 个单位根，则经 次差分后就变成一个平稳过程，这样的性质称为齐次非平稳性。即若 是平稳时间序列，则称 是d阶齐次非平稳序列，这里 表示d阶差分。考虑如下形式的模型：（4.56）其中，是平稳的自回归算子，为可逆的移动平均算

23、子。而 是对序列 进行d阶差分之后的序列，并且得到的该序列具有平稳性特征。若用 替代 ，则（4.56）式就可以表示为：（4.57）则该表达式与前面所属的ARMA模型的表达式相同。而方程（4.57）则表示的是一个ARIMA模型。dddttyx txd()()dttLxL()L()Ldtxtxtydtx()()ttL yL 单整自回归移动平均模型l二、ARIMA模型的确定以上证指数为例为说明 模型参数的相关参数的确定，在此我们选取上证指数为例进行解释。在此，选取上证指数2000年2月1日至2006年6月30日为观测区间。首先，确定ARIMA(p,d,q)模型中的d值。其次，对ARIMA(p,1,q

24、)模型中的p和q数值进行确定。从收益率的自相关系数和偏自相关系数图中我们可以看到,它们都是拖尾的，因此可设定为 过程。收益率的自相关函数第1阶是显著的，从第2阶开始下降很大，数值也不太显著，因此先设定q值为1。收益率的偏相关函数也是第1阶很显著，从第2阶开始下降很大，因此设定p值为1，于是初步建立 模型。(,)ARIMA p d qARMA(1,1,1)ARIMA单整自回归移动平均模型l三、ARIMA过程应用和结果解释l为分析ARIMA过程，在此，我们选取上证指数1990年12月19日至2006年8月31日为观测区间，共3856个观测样本。应用的软件为SAS8.0。单整自回归移动平均模型l四、

25、ARIMA过程的SAS程序模拟除了利用SAS程序进行ARIMA过程分析之外，我们还可以通过这一程序进行ARIMA过程模拟。在此，我们应用了朱世武（2004）所著的基于SAS系统的金融计算中的一个例子进行说明。对ARIMA(0,1,1)进行SAS模拟实现。假定初始数值 =0.01，产生1000个来自 ，的随机时间序列。1r110.8ttttrr1(0,1)tiidN单整自回归移动平均模型图图4-17：ARIMA的的SAS模拟模拟 平稳性与单位根检验l第四节 平稳性与单位根检验l一、非平稳性检验的必要性从前几节知，当时间序列含有单位根时，它就是一个非平稳时间序列。而非平稳时间序列恰好具有这种齐次非

26、平稳特征，即通过足够次数的差分就可以转换为一个平稳的时间序列。l1.单整性的定义若一个非平稳时间序列 必须经过d次差分后才能变换成一个平稳的、可逆的 时间序列，则称 具有d阶单整性，用 表示。tyARMAty()tyI d平稳性与单位根检验l2.伪回归问题如果对非平稳性数据进行回归，则在回归结果中，我们可能会发现R2很高，t值也极高，这似乎表示变量之间存在着很好的拟合关系。但是，同时会发现杜宾-沃森d值偏低。这时，则可能存在伪回归（spurious regressions）现象发生。即回归结果是不正确的。Granger和 Newbold曾经提出一个良好的经验规则：当 时，所估计的回归就有谬误之

27、嫌。有时候时间序列的高度相关仅仅是因为两者同时随时间有或上或下变动的趋势，并没有真正的联系。这种情况就称为伪回归。2Rd平稳性与单位根检验l二、两种类型的平稳性通常，有两种类型被用来描述非平稳性，它们是带漂移的随机游走模型和趋势平稳过程。其中，带漂移的随机游走模型表达为：趋势平稳过程是因其在线性趋势附近而得名，此过程表达为：在上述情况下，是白噪声扰动项。1tttyyttytt平稳性与单位根检验图图4-18 随机游走与带漂移的随机游走时间序列图随机游走与带漂移的随机游走时间序列图平稳性与单位根检验l 三、单位根检验 l1.ADF检验l检验经济时间序列是否平稳，需要先检验单位根的存在。常用测验单位

28、根的方法是由Dickey 和Fuller(Fuller,1976;Dickey and Fuller,1979)提出的Dickey-Fuller(DF)检验，即单位根检验。l开始模型：l其中，是随机误差项。1tttyytu平稳性与单位根检验l2.ADF检验模型的确定lADF检验模型是一般形式，然而是否应该包含常数项 ，是否包含时间趋势项 ，以及如何确定最优滞后阶数p，这是一个需要解决的现实问题。l首先，我们来看如何判断检验模型是否应该包含常数项与时间趋势项。l其次，我们来看如何确定检验模型的最优滞后阶数。l3.菲利普斯-配荣（Phillips-Perron,PP）检验lPP检验针对的是回归模型

29、的干扰项 存在异方差或序列相关的现象。回归模型的三种形式及检验规则与DF检验相同。但PP检验下，这两个统计量的计算相对复杂，是在对应DF统计量的形式上加以修正。但PP检验比照的临界值分布表和DF检验下的三种回归形式下的临界值分布表相同。12tt平稳性与单位根检验l案例说明4-5运用EVIEWS进行上证指数单位根检验l对上证指数1990.12.19-2006-8-31的上证指数日收盘序列 进行单位根检验。l根据模型所示，计算得到的ADF统计量0.022985，大于临界值，故不能拒绝被检验的指数序列非平稳的零假设。对其一阶差分序列进行ADF检验，此时的统计量为-59.79900，小于相应临界值，故

30、拒绝指数差分序列非平稳的假设。综合上两个结果，得出指数序列为一阶单整的结论。,1,2,3856ty t 平稳性与单位根检验l案例说明4-6股市与经济增长的关联性检验l李志刚（2005）选取了1992年第1季度至2005年第2季度的时间序列，共53组观察值。本文首先由各季度的累计GDP求得每季度的GDP增加值，并采用Y-11方法进行季节调整（记为QGDP，单位：亿元人民币），而后取自然对数，即economy=lnQGDP。上证指数的季度数据，以每季度最后一个交易日的收盘价作为该季度的指数值（记为QSHY），而后取自然对数，即 stock=lnQSHY。平稳性与单位根检验5678910111992

31、199419961998200020022004economystock图图 4-22 股指和经济增长运动的共同漂移股指和经济增长运动的共同漂移平稳性与单位根检验l如图所示，economy与stock这两个序列在一定程度上是一起漂移的，说明二者可能存在着协整关系。本文采用ADF（Augmented Dickey-Fuller）检验来考察各变量是否平稳。如表2所示，无论是否有时间趋势，各变量的原序列均无法一致拒绝存在单位根的零假设，而一阶差分序列则均拒绝该零假设。因此，economy、stock均为1 阶单整的时间序列。9、静夜四无邻，荒居旧业贫。22.7.2922.7.29Friday,Jul

32、y 29,202210、雨中黄叶树，灯下白头人。17:42:3617:42:3617:427/29/2022 5:42:36 PM11、以我独沈久，愧君相见频。22.7.2917:42:3617:42Jul-2229-Jul-2212、故人江海别，几度隔山川。17:42:3617:42:3617:42Friday,July 29,202213、乍见翻疑梦，相悲各问年。22.7.2922.7.2917:42:3617:42:36July 29,202214、他乡生白发，旧国见青山。2022年7月29日星期五下午5时42分36秒17:42:3622.7.2915、比不了得就不比，得不到的就不要。2

33、022年7月下午5时42分22.7.2917:42July 29,202216、行动出成果，工作出财富。2022年7月29日星期五17时42分36秒17:42:3629 July 202217、做前，能够环视四周；做时，你只能或者最好沿着以脚为起点的射线向前。下午5时42分36秒下午5时42分17:42:3622.7.299、没有失败，只有暂时停止成功！。22.7.2922.7.29Friday,July 29,202210、很多事情努力了未必有结果，但是不努力却什么改变也没有。17:42:3617:42:3617:427/29/2022 5:42:36 PM11、成功就是日复一日那一点点小小

34、努力的积累。22.7.2917:42:3617:42Jul-2229-Jul-2212、世间成事，不求其绝对圆满，留一份不足，可得无限完美。17:42:3617:42:3617:42Friday,July 29,202213、不知香积寺，数里入云峰。22.7.2922.7.2917:42:3617:42:36July 29,202214、意志坚强的人能把世界放在手中像泥块一样任意揉捏。2022年7月29日星期五下午5时42分36秒17:42:3622.7.2915、楚塞三湘接，荆门九派通。2022年7月下午5时42分22.7.2917:42July 29,202216、少年十五二十时，步行夺得

35、胡马骑。2022年7月29日星期五17时42分36秒17:42:3629 July 202217、空山新雨后，天气晚来秋。下午5时42分36秒下午5时42分17:42:3622.7.299、杨柳散和风，青山澹吾虑。22.7.2922.7.29Friday,July 29,202210、阅读一切好书如同和过去最杰出的人谈话。17:42:3617:42:3617:427/29/2022 5:42:36 PM11、越是没有本领的就越加自命不凡。22.7.2917:42:3617:42Jul-2229-Jul-2212、越是无能的人，越喜欢挑剔别人的错儿。17:42:3617:42:3617:42Fr

36、iday,July 29,202213、知人者智，自知者明。胜人者有力，自胜者强。22.7.2922.7.2917:42:3617:42:36July 29,202214、意志坚强的人能把世界放在手中像泥块一样任意揉捏。2022年7月29日星期五下午5时42分36秒17:42:3622.7.2915、最具挑战性的挑战莫过于提升自我。2022年7月下午5时42分22.7.2917:42July 29,202216、业余生活要有意义，不要越轨。2022年7月29日星期五17时42分36秒17:42:3629 July 202217、一个人即使已登上顶峰，也仍要自强不息。下午5时42分36秒下午5时42分17:42:3622.7.29MOMODA POWERPOINTLorem ipsum dolor sit,eleifend nulla ac,fringilla purus.Nulla iaculis tempor felis amet,consectetur adipiscing elit.Fusce id urna blanditut cursus.感 谢 您 的 下 载 观 看感 谢 您 的 下 载 观 看专家告诉