时间序列分析培训讲义dpvu.pptx

《时间序列分析培训讲义dpvu.pptx》由会员分享，可在线阅读，更多相关《时间序列分析培训讲义dpvu.pptx（180页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、时间序列分析西安交通大学经济与金融学院统计系赵春艳本课程内容体系：第一章：平稳时间序列分析导论第二章：平稳时间序列分析的基础知识第三章：平稳时间序列模型的建立第四章：协整理论导论第五章：单位根过程第六章：单位根过程的假设检验第七章：协整理论参考书目参考书目:1 1、陆懋祖，高等时间序列经济计量学，上海人民出版、陆懋祖，高等时间序列经济计量学，上海人民出版社，社，19991999年版；年版；2 2、王振龙主编，时间序列分析，中国统计出版社，、王振龙主编，时间序列分析，中国统计出版社，20002000；3 3、王耀东等编，经济时间序列分析，上海财经大学出、王耀东等编，经济时间序列分析，上海财经大学

2、出版社，版社，19961996；4 4、马薇，协整理论与应用，南开大学出版社，、马薇，协整理论与应用，南开大学出版社，20042004；5 5、王少平，宏观计量的若干前沿理论与应用，南开大、王少平，宏观计量的若干前沿理论与应用，南开大学出版社，学出版社，20032003。第一章 平稳时间序列分析导论一、时间序列1、含义：指被观察到的依时间为序排列的数据序列。2、特点：（1）现实的、真实的一组数据，而不是数理统计中做实验得到的。既然是真实的，它就是反映某一现象的统计指标，因而，时间序列背后是某一现象的变化规律。（2）动态数据。二、时间序列分析1、时间序列分析：是一种根据动态数据揭示系统动态结构和

3、规律的统计方法。其基本思基本思想想：根据系统的有限长度的运行记录（观察数据），建立能够比较精确地反映序列中所包含的动态依存关系的数学模型，并借以对系统的未来进行预报（王振龙）2、计量经济学中的建模方法和思想3、理论依据：尽管影响现象发展的因素无法探求，但其结果之间却存在着一定的联系，可以用相应的模型表示出来，尤其在随机性现象中。三、确定性时间序列分析与随机性时间序列分析时间序列依据其特征，有以下几种表现形式，并产生与之相适应的分析方法：（1）长期趋势变化 受某种基本因素的影响，数据依时间变化时表现为一种确定倾向，它按某种规则稳步地增长或下降。使用的分析方法有：移动平均法、指数平滑法、模型拟和法

4、等；（2）季节性周期变化 受季节更替等因素影响，序列依一固定周期规则性的变化，又称商业循环。采用的方法：季节指数；（3）循环变化 周期不固定的波动变化。(4)(4)随机性变化随机性变化由许多不确定因素引起的序列变化。它所使用的分析由许多不确定因素引起的序列变化。它所使用的分析方法就是我们要讲的时间序列分析。方法就是我们要讲的时间序列分析。确定性变化分析确定性变化分析 趋势变化分析趋势变化分析 周期变化分析周期变化分析 循环变化分析循环变化分析时间序列分析时间序列分析 随机性变化分析随机性变化分析 AR AR、MAMA、ARMAARMA模型模型四、发展历史1、时间序列分析奠基人：20世纪40年代

5、分别由Norbort Wiener 和Andrei Kolemogoner 独立给出的，他们对发展时间序列的参数模型拟和和推断过程作出了贡献，提供了与此相关的重要文献，促进了时间序列分析在工程领域的应用。2、时间序列分析在经济领域的应用 20世纪70年代，G.P.Box 和G.M.Jenkins发表专著时间序列分析：预测和控制，使时间序列分析的应用成为可能。3、现代时间序列分析的发展趋势（1）单位根检验（2）协整检验2003年度诺贝尔经济学奖的获得者是美国经济学家罗伯特.恩格尔和英国经济学家克莱夫.格兰杰。获奖原因：“今年的获得者发明了处理许多经济时间序列两个关键特性的统计方法：时间变化的变更

6、率和非平稳性。”两人是时间序列经济学的奠基人。时间变化的变更率指方差随时间变化而变化的频率，这主要是指恩格尔在1982年发表的条件异方差模型（ARCH），最初主要用于研究英国的通货膨胀问题，后来广泛用作金融分析的高级工具；传统的计量经济学研究中，通常假定经济数据和产生这些数据的随机过程是平稳的。格兰杰的贡献主要是在非平稳过程假定下所进行的严格计量模型的建立。（协整检验）第二章第二章 平稳时间序列分析的基础知识平稳时间序列分析的基础知识 第一节第一节 随机序列随机序列一、随机过程一、随机过程 1、定义：在数学上，随机过程被定义为一组随机变量，即，其中，T表示时间t的变动范围，对每个固定的时刻 t

7、而言，Zt是一随机变量，这些随机变量的全体就构成一个随机过程。2、特征（1）随机过程是随机变量的集合（2）构成随机过程的随机变量是随时间产生的，在任意时刻，总有随机变量与之相对应。二、随机序列（时间序列）二、随机序列（时间序列）1、当 时，即时刻t只取整数时，随机过程 可写成此类随机过程 称为随机序列，也成时间序列。可见（1）随机序列是随机过程的一种，是将连续时间的随机过程等间隔采样后得到的序列；（2）随机序列也是随机变量的集合，只是与这些随机变量联系的时间不是连续的、而是离散的。三、时间序列的分布、均值、协方差三、时间序列的分布、均值、协方差 函数函数1、分布函数(1)一维分布函数:随机序列

8、中每个随机变量的分布函数.F1(z),F2(z),Ft-1(z),Ft(z)(2)二维分布函数:随机序列中任意两个随机变量的联合分布函数Fi,j(zi,zj).i,j=,-2,-1,0,1,2,(3)柯尔莫哥洛夫定理与有限维概率分布柯尔莫哥洛夫定理表明,一个随机序列的特征,可以用它的有限维分布表示出来。2 2、均值函数、均值函数对随机序列中的任一随机变量取期望。对随机序列中的任一随机变量取期望。当当t t取遍所有可能整数时，就形成了离散时间的函数取遍所有可能整数时，就形成了离散时间的函数u ut t称称u ut t 为时间序列的均值函数。为时间序列的均值函数。3 3、自协方差函数和自相关函数、

9、自协方差函数和自相关函数 自相关函数：当t,s取遍所有可能的整数时，就形成了时间序列的自相关函数，它描述了序列的自相关结构。它的本质等同于相关系数。第二节 平稳时间序列一、平稳时间序列一、平稳时间序列1、定义：时间序列zt是平稳的。如果zt有有穷的二阶中心矩，而且满足：（1）ut=Ezt=c;（2）r(t,s)=E(zt-c)(zs-c)=r(t-s,0)则称zt是平稳的。含义：a a有穷二阶矩意味着期望和自协方差存在；有穷二阶矩意味着期望和自协方差存在；b b平稳时间序列任意时刻所对应的随机变量的均值相平稳时间序列任意时刻所对应的随机变量的均值相等；等；c c自协方差函数只与时间间隔有关，而

10、与时间起点无自协方差函数只与时间间隔有关，而与时间起点无关。关。二、平稳时间序列的均值、自协方差和自相关函数1、均值函数：平稳时间序列均值为常数，为分析方便，假定E zt=0，当均值不为零时，给每个值减去均值后再求均值，即等于0。2 2、自协方差函数：平稳时间序列的自协方差仅与时间间隔有关，而与具体时刻无关，所以，自协方差函数仅表明时间间隔即可。3、自相关函数k 平稳时间序列自协方差仅与时间隔有关，当间隔为 零时，自协方差应相等：4、自协方差与自相关函数的性质 (1)rk=r-k k=-k k、k仅是时间先后顺序上的差异，它们代表的间隔是相同的。(2)三、偏自相关函数（PACF)1、偏自相关函

11、数用来考察扣除zt 和zt+k之间zt+1，zt+2，zt+k-1影响之后的zt 和zt+k之间的相关性。2 2、偏自相关函数的定义、偏自相关函数的定义设设zzt t 为零均值平稳序列，为零均值平稳序列，z zt+1 t+1，z zt+2t+2，z zt+k-1t+k-1对对z zt t 和和z zt+k t+k 的线性估计为：的线性估计为：kk表示偏自相关函数，则：3、PACF的涵义设有zt+1，zt+2，zt+34、pacf的推导四、随机序列的特征描述（1）样本均值（2）样本自协方差函数（3）样本自相关函数（4）样本偏自相关函数）样本偏自相关函数例1、设动态数据16，12，15，10，9，

12、17，11，16，10，14，求样本均值、样本自相关函数（SACF）和偏自相关函数（SPACF）（各求前三项）第三节第三节 线性平稳时间序列模型线性平稳时间序列模型一、自回归过程一、自回归过程(A R(p)）1、2、AR(P)模型的模型的ACF、PACF特征特征以AR(1)为例例：例：k k1 12 23 34 45 56 67 78 89 91010 k k0.880.880.760.760.670.670.570.570.480.480.40.40.340.340.280.280.210.210.170.17 kkkk0.880.880.010.01-0.010.010.110.110.0

13、20.02-0.010.010.010.01-0.020.02-0.060.060.050.05计算结果表明，ACF逐渐衰减，但不等于零；PACF在k=1后，与零接近，是截尾的。结论：ACF呈指数衰减，是拖尾的；PACF在一步后为零，是截尾的。二、滑动平均模型（二、滑动平均模型（MA(q)）1、形如zt=at-1at-1-2at-2-qat-q模型为滑动平均模型，其中，简化形式zt=(B)at(B)=1-1B-2B2-qBq,满足(B)=0的根在单位圆外，即B，1此时该过程是可逆的。2、MA模型的模型的ACF及及PACF(3)PACF例：用zt=(1-0.5B)at模拟产生250个观察值，at

14、为白噪声序列，得到序列自相关和偏自相关函数如下：可见，ACF在一步后截尾，PACF是拖尾的。结论：MA(q)的ACF是截尾的，PACF是拖尾的。k k1 12 23 34 45 56 67 78 89 91010ACFACF-0.44-0.440 00.020.02-0.03-0.03-0.01-0.01-0.050.050.040.04-0.030.03-0.030.030.020.02PACFPACF-0.44-0.44-0.240.24-0.11-0.11-0.08-0.08-0.07-0.07-0.120.12-0.060.06-0.070.07-0.1-0.1-0.080.08三、自

15、回归滑动平均模型（三、自回归滑动平均模型（AR M A(p,q)）1、2、ARMA(p,q)的的ACF和和PACF(2)ACF、PACF均是拖尾的例：(1-0.9B)zt=(1-0.5B)at模拟产生250个观察值，ACF、PACF如下表所示：k k1 12 23 34 45 56 67 78 89 91010acfacf0.50.57 70.50.50.470.470.350.350.30.31 10.250.250.20.21 10.10.18 80.10.10.10.12 2pacfpacf0.50.57 70.20.26 60.180.18-0.03-0.030.00.01 1-0.0

16、10.010.00.01 10.00.01 1-0.080.080.00.05 5本节介绍了三类模型的形式、特性及自相关和偏自相关函数的特征，现绘表如下：AR(p)AR(p)MA(q)MA(q)ARMA(p,q)ARMA(p,q)模型方程模型方程(B)=a(B)=at tz zt t=(B)a(B)at t(B)z(B)zt t=(B)a(B)at t平稳性条件平稳性条件(B)=0(B)=0的根在的根在单位圆外单位圆外无无(B)=0(B)=0的根在单位圆外的根在单位圆外可逆性条件可逆性条件无无(B)=0(B)=0的根在的根在单位圆外单位圆外(B)=0(B)=0的根在单位圆外的根在单位圆外自相关

17、函数自相关函数拖尾拖尾Q Q步截尾步截尾拖尾拖尾偏自相关函数偏自相关函数P P步截尾步截尾拖尾拖尾拖尾拖尾第三章第三章 平稳时间序列模型的建立平稳时间序列模型的建立第一节第一节 模型识别与定阶模型识别与定阶一、模型识别一、模型识别1、含义：对一个观察序列，选择一个与其实际过程相吻合的模型结构。2、方法：利用序列的acf、pacf识别。判断截尾、拖尾的主观性较大，只是初步识别。二、模型定阶二、模型定阶二、模型定阶二、模型定阶（一）（一）（一）（一）a c fa c f、p a c fp a c f方法方法方法方法（1 1）M A(q):M A(q):BartlettBartlett公式：当公式：

18、当kqkq时，时，N N充分大，充分大，(2)AR（P）：（二）残差方差图：（二）残差方差图：（1）残差：在多元回归y=a1x1+a2x2+.+an x n+at,存在自变量x的选择问题。如果x选择不够，模型拟合不足，表现为y与 差异较大；若x选择多，则过度拟合，y与差异减小速度很慢。将(y-)称为残差，多元回归就是利用此确定模型的自变量，即新增或减少变量是否会显著影响残差。（2）将该思想应用到时间序列模型定阶上。(3)利用a2的变化规律，确定模型阶数。随着模型阶数的增大，分母减小；分子在不足拟合时，一直减小，速度较快；过拟合时，分子虽减小，但速度很慢，几乎不变。a2取决于分子、分母减小的速度

19、。在不足拟合时，a2一直减小；过拟合时，a2却增大。选择a2的最低点为模型的最优阶数。（三）（三）（三）（三）F F 检验定阶法：检验定阶法：检验定阶法：检验定阶法：（1 1）F F分布：分布：（2）用F分布检验两个回归模型是否有显著差异。(3)对于ARMA(p,q)模型定阶例如：在ARMA(p,q)和ARMA(p-1,q-1)选择。例：每隔20分钟进行一次观察的造纸过程入口开关调节器的观察值（第241页，18）1、series Mean S.D Max Min z 32.02 0.74 34 30.7令z1=z 32.022 2、1 2 3 4 5 6 7 8 9 10 1 2 3 4 5

20、6 7 8 9 10 acf 0.868 0.782 0.708 0.663 0.627 0.617 0.594 0.559 0.5 0.48 acf 0.868 0.782 0.708 0.663 0.627 0.617 0.594 0.559 0.5 0.48 pacf 0.868 0.115 0.028 0.099 0.055 0.122 0.01 0.04 -0.099 0.1 pacf 0.868 0.115 0.028 0.099 0.055 0.122 0.01 0.04 -0.099 0.1 3、定阶（1）acf、pacf：从 acf、pacf可知，acf拖尾，pacf截尾，初

21、步识别为AR模型。具体阶数：(2)残差方差：(3)F检验：（四）最佳准则函数定阶法1、基本思想：确定一个函数，该函数既要考虑用某一模型拟合原始数据的接近程度，同时又考虑模型中所含参数的个数。当该函数取最小值时，就是最合适的阶数。衡量模型拟合数据的接近程度的指标是残差方差。残差方差=2、最佳准则函数包括FPE、AIC、BIC准则。3、AIC准则(1)该准则既适合于AR，也适合于ARMA模型。关于ARMA模型的定阶1、ACF、PACF都呈现一定的拖尾性，试拟合ARMA模型。Pandit-Wu于1977年提出了不同于Box-Jenkins的系统建模方法。该方法认为，任一平稳序列总可以用一个ARMA(

22、n,n-1)表示，AR(n)、MA(m)、ARMA(n,m)都是ARMA(n,n-1)的特例。2、建模思想：逐渐增加模型阶数，直到剩余平方和不再减小为止。3、如何在不同模型之间取舍第四章 协整理论绪论一、协整理论产生的背景1、20世纪70年代以前的建模技术以时间序列平稳为前提设计的。2、理论假定与现实的矛盾。3、协整理论的产生-计量经济学方法研究的新阶段-Granger首先提出了伪回归问题（1974）；-1978年，EngleGranger发表论文“协整与误差修正”，正式提出“协整”（cointegration）概念二、与协整检验有关的两个问题：单位根和误差修正模型1、单位根：协整检验处理的是

23、非平稳时间序列，单位根检验就是要说明一个时间序列的平稳性。包括DF和ADF检验2、误差修正模型（Error Correction Model,ECM）：ECM由、Hendry、Srba于1978年提出的。三、本部分的体系单位根检验-协整检验-误差修正模型第五章 单位根过程第一节 单位根过程的定义一、随机游动过程的定义1、随机过程y t,t=1,2,若y t=yt-1+t，其中t为独立同分布序列，E（t）=0，D（t）=E（t 2）=2则称y t为随机游动过程。2、随机游动过程是一非平稳过程(1)y t=yt-1+t =yt-2+t-1+t =yt-3+t-2+t-1+t =.=y0+1+2+t

24、E(y t)=y0(2)D(yt)=E(yt-y0)2=E(1+2+t)2=t2二、单位根过程的定义1、随机过程y t,t=1,2,若y t=yt-1+t，其中=1，t 为稳定过程，E（t）=0，Cov（t，t-s）=s1时，1时，就是平稳过程。4、单位根过程与稳定过程的本质区别第二节 与单位根过程形式接近的几种模型一、带常数项的随机游动过程1、2、深圳股票综合指数 二、长期趋势1、形如 称为确定趋势模型。2、前两类模型的图形接近。3、判别单位根的必要性。yt=0.1 t+ut 生成的序列 图三、含随机趋势和确定性趋势的混合随机过程1、yt=0.1+0.1t+yt-1+ut生成的序列 图四、近

25、单位根过程1、第六章 单位根过程的假设检验第一节 迪基-福勒(DF)检验法一、DF检验法产生的背景1、DF检验法是由Dickey、Fuller在20世纪70、80年代的一系列文章中建立起来的。2、3、这种方法不能用来检验H0:=1，当零假设成立时，t T不再服从t分布，因而无法得到临界值。此时，只能用模拟方法得到临界值。DF检验中用到两个统计量：T（T-1）和t T，它们不存在小样本分布，只有当样本容量T足够大时，它们的极限分布才有实际的应用价值。二、情况一的DF检验1、假设数据由 产生，并在其中检验H0:=1;H1:12、适用于数据是非平稳且没有趋势的情况。3、例：利用1947年第二季度到1

26、989年第一季度的数据对美国财政部债券利息率作不带常数的一阶自回归如下：三、情况二的DF检验1、假设数据由 产生，在一般先检验=1，若接受H0，再检验=0。若=0，则为 ，若 0，则为2、情况二适用的数据图形是有趋势，但不稳定的情况。这时，就在随机性非平稳及有漂移趋势的非平稳之间选择。3、例：仍利用美国财政部债券利率数据，估计带常数项的一阶自回归模型：四、情况三的DF检验1、情况三的DF检验（1）假设数据是由带常数项的单位根过程（2）缺陷五、情况四的DF检验1、（2）适用于序列有趋势的情况3、例：美国1947年一季度至1989年第二季度GNP的实际值，对图中数据进行模型拟合。解：（1）图中数据

27、有明显的长期趋势；（2）这类图形可能适合的模型有：（3）六、DF检验小结第二节 增广的迪基-福勒(ADF)检验法一、ADF检验法（Augmented DickeyFuller Test）1、ADF检验法是由迪基（Dickey）和福勒（Fuller）在1979年提出的，是DF方法的推广。DF假定t是独立同分布序列，ADF假定随机扰动项t是稳定过程。2、原理：ADF假设数据服从有单位根的P阶自回归过程，即证明：二、情况二的ADF检验1、2、例：利用ADF检验法对美国财政部债券利率进行单位根检验。解：H0:=1;H1:0）查表6,=0.05，情况三，临界值为29.509,38.8529.509,拒绝

28、H0。(2)4、协整向量最大特征值 1=0.1105对应的特征向量就是协整向量，有第四节 误差修正模型(ECM)一、协整系统的表述2、假定yt的元素之间存在k个独立的协整关系，且协整向量为二、Granger表述定理设 yt是n维I(1)随机过程，若yt中有k个协整关系，即存在n*k阶矩阵A，r(A)=k,使得三、误差修正模型（ECM）1、ECM的主要形式是由Davidson、Hendry、Srba和Yeo于1978年提出。它将变量之间的短期与长期联系有机地结合在一起。2、ECM的形式（1）对于（1，1）阶自回归分布滞后模型式（2）为误差修正模型，为误差修正项。（2）对模型的理解ECM的被解释变量是yt，因此，实际上是一个短期模型，反映了yt的短期波动yt是如何被决定的。若y与z之间存在长期均衡关系，即y=az，在式（1）中，若则，不考虑常数项误差修正模型的含义：y t受另一变量的z t的影响；受（yt-1-）的影响。均值对一序列而言是恒定的，实际序列值与均值离差对y t的的影响，表示一种恒定力量对y t的作用，称为长期均衡影响。表明 是由 决定的，也就是说，y与z之间有长期均衡关系，它们的均值之间才会存在稳定关系。