时间序列的平稳性及其检验.pptx
《时间序列的平稳性及其检验.pptx》由会员分享,可在线阅读,更多相关《时间序列的平稳性及其检验.pptx(75页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、一、问题的引出:非平稳变量与经一、问题的引出:非平稳变量与经典回归模型典回归模型第1页/共75页常见的数据类型常见的数据类型到目前为止,经典计量经济模型常用到的数据有:时间序列数据(time-series data);截面数据(cross-sectional data)平行/面板数据(panel data/time-series cross-section data)时间序列数据是最常见,也是最常用到的数据。第2页/共75页经典回归模型与数据的平稳性经典回归模型与数据的平稳性经典回归分析暗含着一个重要假设:数据是平稳的。数据非平稳,大样本下的统计推断基础“一致性”要求被破怀。经典回归分析的假设
2、之一:解释变量X是非随机变量放宽该假设:X是随机变量,则需进一步要求:(1)X与随机扰动项 不相关Cov(X,)=0依概率收敛:(2)第3页/共75页 第(2)条是为了满足统计推断中大样本下的“一致性”特性:第(1)条是OLS估计的需要如果X是非平稳数据(如表现出向上的趋势),则(2)不成立,回归估计量不满足“一致性”,基于大样本的统计推断也就遇到麻烦。因此:注意:在双变量模型中:第4页/共75页 表现在:两个本来没有任何因果关系的变量,却有很高的相关性(有较高的R2):例如:如果有两列时间序列数据表现出一致的变化趋势(非平稳的),即使它们没有任何有意义的关系,但进行回归也可表现出较高的可决系
3、数。在现实经济生活中:情况往往是实际的时间序列数据是非平稳的,而且主要的经济变量如消费、收入、价格往往表现为一致的上升或下降。这样,仍然通过经典的因果关系模型进行分析,一般不会得到有意义的结果。数据非平稳,往往导致出现数据非平稳,往往导致出现“虚假回归虚假回归”问题问题第5页/共75页 时间序列分析模型方法就是在这样的情况下,以通过揭示时间序列自身的变化规律为主线而发展起来的全新的计量经济学方法论。时间序列分析已组成现代计量经济学的重要内容,并广泛应用于经济分析与预测当中。第6页/共75页二、时间序列数据的平稳性二、时间序列数据的平稳性第7页/共75页 时间序列分析中首先遇到的问题是关于时间序
4、列数据的平稳性问题。假定某个时间序列是由某一随机过程(stochastic process)生成的,即假定时间序列Xt(t=1,2,)的每一个数值都是从一个概率分布中随机得到,如果满足下列条件:1)均值E(XE(Xt t)=)=是与时间t 无关的常数;2)方差Var(XVar(Xt t)=)=2 2是与时间t 无关的常数;3)协方差Cov(XCov(Xt t,X,Xt+kt+k)=)=k k 是只与时期间隔k有关,与时间t 无关的常数;则称该随机时间序列是平稳的(stationary),而该随机过程是一平稳随机过程(stationary stochastic process)。第8页/共75页
5、 例4.1一个最简单的随机时间序列是一具有零均值同方差的独立分布序列:Xt=t ,tN(0,2)例4.2另一个简单的随机时间列序被称为随机游走(random walk),该序列由如下随机过程生成:Xt=Xt-1+t这里,t是一个白噪声。该序列常被称为是一个白噪声(white noise)。由于Xt具有相同的均值与方差,且协方差为零,由定义,一个白噪声序列是平稳的。第9页/共75页 为了检验该序列是否具有相同的方差,可假设Xt的初值为X0,则易知 X1=X0+1 X2=X1+2=X0+1+2 X Xt t=X=X0 0+1+2+t 由于X0为常数,t是一个白噪声,因此Var(Xt)=t2 即Xt
6、的方差与时间t t有关而非常数,它是一非平稳序列。容易知道该序列有相同的均值:E(Xt)=E(Xt-1)第10页/共75页然而,对X取一阶差分(first difference):Xt=Xt-Xt-1=t由于t是一个白噪声,则序列Xt是平稳的。后面将会看到:如果一个时间序列是非平稳的,它常常可通过取差分的方法而形成平稳序列。事实上,随机游走过程是下面我们称之为1 1阶自回归AR(1)AR(1)过程的特例 X Xt t=X Xt-1t-1+t 不难验证:1)|1|1时,该随机过程生成的时间序列是发散的,表现为持续上升(1)1)或持续下降(-1)1),因此是非平稳的;第11页/共75页 第二节中将
7、证明:只有当-1-1 10,样本自相关系数近似地服从以0为均值,1/n 为方差的正态分布,其中n为样本数。也可检验对所有k0k0,自相关系数都为0 0的联合假设,这可通过如下Q QLBLB统计量进行:第18页/共75页 该统计量近似地服从自由度为m的2分布(m为滞后长度)。因此:如果计算的Q Q值大于显著性水平为 的临界值,则有1-1-的把握拒绝所有 k k(k0)(k0)同时为0 0的假设。例例4.3:4.3:4.3:4.3:表4.14.1序列Random1Random1是通过一随机过程(随机函数)生成的样本容量为1919的随机时间序列。第19页/共75页表 4.1 4.1 一个纯随机序列与
8、随机游走序列的检验 序号 Random1 自相关系数 kr(k=0,1,17)LBQ Random2 自相关系数 kr(k=0,1,17)LBQ 1-0.031 K=0,1.000 -0.031 1.000 2 0.188 K=1,-0.051 0.059 0.157 0.4800.480 5.116 3 0.108 K=2,-0.393 3.679 0.264 0.018 5.123 4-0.455 K=3,-0.147 4.216-0.191-0.069 5.241 5-0.426 K=4,0.280 6.300-0.616 0.028 5.261 6 0.387 K=5,0.187 7.
9、297-0.229-0.016 5.269 7-0.156 K=6,-0.363 11.332-0.385-0.219 6.745 8 0.204 K=7,-0.148 12.058-0.181-0.063 6.876 9-0.340 K=8,0.315 15.646-0.521 0.126 7.454 10 0.157 K=9,0.194 17.153-0.364 0.024 7.477 11 0.228 K=10,-0.139 18.010-0.136-0.249 10.229 12-0.315 K=11,-0.297 22.414-0.451-0.404 18.389 13-0.377
10、K=12,0.034 22.481-0.828-0.284 22.994 14-0.056 K=13,0.165 24.288-0.884-0.088 23.514 15 0.478 K=14,-0.105 25.162-0.406-0.066 23.866 16 0.244 K=15,-0.094 26.036-0.162 0.037 24.004 17-0.215 K=16,0.039 26.240-0.377 0.105 25.483 18 0.141 K=17,0.027 26.381-0.236 0.093 27.198 19 0.236 0.000 第20页/共75页容易验证:该样
11、本序列的均值为0 0,方差为0.07890.0789。从图形看:它在其样本均值0 0(a)附近上下波动,且样本自相关系数(b)迅速下降到0 0,随后在0 0附近波动且逐渐收敛于0 0。第21页/共75页 由于该序列由一随机过程生成,可以认为不存在序列相关性,因此该序列为一白噪声。根据Bartlett的理论:kN(0,1/19)因此任一rk(k0)的95%的置信区间都将是 可以看出:k0k0时,r rk k的值确实落在了该区间内,因此可以接受 k k(k0)k0)为0 0的假设。同样地,从Q QLBLB统计量的计算值看,滞后1717期的计算值为26.3826.38,未超过5%5%显著性水平的临界
12、值27.5827.58,因此,可以接受所有的自相关系数 k k(k0)k0)都为0 0的假设。因此,该随机过程是一个平稳过程。第22页/共75页 序列Random2是由一随机游走过程 Xt=Xt-1+t 生成的一随机游走时间序列样本。其中,第0项取值为0,t是由Random1表示的白噪声。第23页/共75页 样本自相关系数显示:r1=0.48,落在了区间-0.4497,0.4497之外,因此在5%的显著性水平上拒绝1的真值为0的假设。该随机游走序列是非平稳的。图形表示出:该序列具有相同的均值,但从样本自相关图看,虽然自相关系数迅速下降到0,但随着时间的推移,则在0附近波动且呈发散趋势。第24页
13、/共75页例4.44.4 检验中国支出法GDP时间序列的平稳性。表4.2 4.2 1978200019782000年中国支出法GDPGDP(单位:亿元)年份GDP年份GDP年份GDP19783605.6198610132.8199446690.719794073.9198711784199558510.519804551.3198814704199668330.419814901.4198916466199774894.219825489.2199018319.5199879003.319836076.3199121280.4199982673.119847164.4199225863.620
14、0089112.519858792.1199334500.6第25页/共75页 图形:表现出了一个持续上升的过程,可初步判断是非平稳的。样本自相关系数:缓慢下降,再次表明它的非平稳性。图 4.54.5 1978197820002000 年中国 GDPGDP 时间序列及其样本自相关图 -0.4-0.20.00.20.40.60.81.01.2246810121416182022GDPACF020000400006000080000100000788082848688909294969800GDP第26页/共75页 拒绝:该时间序列的自相关系数在滞后1期之后的值全部为0的假设。结论:1978200
15、0年间中国GDP时间序列是非平稳序列。从滞后18期的QLB统计量看:QLB(18)=57.1828.86=20.05第27页/共75页例4.54.5 检验3.10中关于人均居民消费与人均国内生产总值这两时间序列的平稳性。原图 样本自相关图 第28页/共75页从图形上看:人均居民消费(CPC)与人均国内生产总值(GDPPC)是非平稳的。从滞后1414期的QLB统计量看:CPC与GDPPC序列的统计量计算值均为57.18,超过了显著性水平为5%时的临界值23.68。再次表明它们的非平稳性。就此来说,运用传统的回归方法建立它们的回归方程是无实际意义的。不过,第三节中将看到,如果两个非平稳时间序列是协
16、整的,则传统的回归结果却是有意义的,而这两时间序列恰是协整的。第29页/共75页四、平稳性的单位根检验四、平稳性的单位根检验第30页/共75页 对时间序列的平稳性除了通过自相关系数及其图形直观判断外,运用统计量进行统计检验则是更为准确与重要的。单位根检验(unit root test)是统计检验中普遍应用的一种检验方法。1 1、DFDF检验检验已知道,随机游走序列 Xt=Xt-1+t是非平稳的,其中t是白噪声。而该序列可看成是随机模型 Xt=Xt-1+t中参数=1时的情形。第31页/共75页也就是说,我们对式 Xt=Xt-1+t (*)做回归,如果确实发现=1,就说随机变量XtXt有一个单位根
17、。(*)式可变形式成差分形式:Xt=(1-)Xt-1+t =Xt-1+t (*)检验(*)式是否存在单位根=1,也可通过(*)式判断是否有=0。第32页/共75页 一般地:检验一个时间序列XtXt的平稳性,可通过检验带有截距项的一阶自回归模型 X Xt t=+X Xt-1t-1+t t (*)中的参数 是否小于1 1。或者:检验其等价变形式 X Xt t=+X Xt-1t-1+t t (*)中的参数 是否小于0 0。在第二节中将证明,(*)式中的参数 11或=1=1时,时间序列是非平稳的;对应于(*)式,则是 00或 =0。第33页/共75页因此,针对式 X Xt t=+X Xt-1t-1+t
18、 t 我们关心的检验为:零假设 H0:=0。备择假设 H1:0 上述检验可通过OLS法下的t检验完成。然而,在原假设(序列非平稳)下,即使在大样本下t统计量也是有偏误的(向下偏倚),通常的t 检验无法使用。Dicky和Fuller于1976年提出了这一情形下t统计量服从的分布(这时的t统计量称为 统计量),即DF分布(见表3.3)。由于t统计量的向下偏倚性,它呈现围绕小于零值的偏态分布。第34页/共75页 因此,可通过OLS法估计 X Xt t=+X Xt-1t-1+t t 并计算t统计量的值,与DF分布表中给定显著性水平下的临界值比较:如果:t临界值,则拒绝零假设H0:=0,认为时间序列不存
19、在单位根,是平稳的。表 4.3 DF 分布临界值表 样 本 容 量 显著性水平 25 50 100 500 t分布临界值(n=)0.01-3.75-3.58-3.51-3.44-3.43-2.33 0.05-3.00-2.93-2.89-2.87-2.86-1.65 0.10-2.63-2.60-2.58-2.57-2.57-1.28 第35页/共75页注意:在不同的教科书上有不同的描述,但是结果是相同的。例如:“如果计算得到的t统计量的绝对值大于临界值的绝对值,则拒绝=0”的假设,原序列不存在单位根,为平稳序列。第36页/共75页 进一步的问题:在上述使用 X Xt t=+X Xt-1t-1
20、+t t对时间序列进行平稳性检验中,实际上假定了时间序列是由具有白噪声随机误差项的一阶自回归过程AR(1)生成的。但在实际检验中,时间序列可能由更高阶的自回归过程生成的,或者随机误差项并非是白噪声,这样用OLS法 进 行 估 计 均 会 表 现 出 随 机 误 差 项 出 现 自 相 关(autocorrelation),导致DF检验无效。另外,如果时间序列包含有明显的随时间变化的某种趋势(如上升或下降),则也容易导致上述检验中的自相关随机误差项问题。为了保证DF检验中随机误差项的白噪声特性,Dicky和 Fuller对 DF检 验 进 行 了 扩 充,形 成 了ADF(Augment Dic
21、key-Fuller)检验。2 2、ADFADF检验检验第37页/共75页ADF检验是通过下面三个模型完成的:模型3 中的t是时间变量,代表了时间序列随时间变化的某种趋势(如果有的话)。检验的假设都是:针对H1:500-2.58-2.23-1.95-1.6125-3.75-3.33-3.00-2.6250-3.58-3.22-2.93-2.60100-3.51-3.17-2.89-2.58250-3.46-3.14-2.88-2.57500-3.44-3.13-2.87-2.57dt500-3.43-3.12-2.86-2.57253.412.972.612.20503.282.892.562
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 时间 序列 平稳 及其 检验
限制150内