时间序列分析——最经典的

【时间简“识”】 说明：本文摘自于经管之家(原人大经济论坛) 作者：胖胖小龟宝。原版请到经管之家(原人大经济论坛) 查看。 1.带你看看时间序列的简史 现在前面的话——

时间序列作为一门统计学，经济学相结合的学科，在我们论坛，特别是五区计量经济学中是热门讨论话题。本月楼主推出新的系

列专题——时间简“识”，旨在对时间序列方面进行知识扫盲（扫盲，仅仅扫盲而已……），同时也想借此吸引一些专业人士能够协

助讨论和帮助大家解疑答惑。 在统计学的必修课里，时间序列估计是遭吐槽的重点科目了，其理论性强，虽然应用领域十分广泛，但往往在实际操作中会遇到很多“令人发指”的问题。所以本帖就从基础开始，为大家絮叨絮叨那些关于“时间”的故事！ Long long ago，有多long？估计大概7000年前吧，古埃及人把尼罗河涨落的情况逐天记录下来，这一记录也就被我们称作所谓的时间序列。记录这个河流涨落有什么意义？当时的人们并不是随手一记，而是对这个时间序列进行了长期的观察。结果，他们发现尼罗河的涨落非常有规律。掌握了尼罗河泛滥的规律，这帮助了古埃及对农耕和居所有了规划，使农业迅速发展，从而创建了埃及灿烂的史前文明。

好~~从上面那个故事我们看到了 1、时间序列的定义——按照时间的顺序把随机事件变化发展的过程记录下来就构成了一个时间序列。 2、时间序列分析的定义——对时间序列进行观察、研究，找寻它变化发展的规律，预测它将来的走势就是时间序列分析。

既然有了序列，那怎么拿来分析呢？ 时间序列分析方法分为描述性时序分析和统计时序分析。 1、描述性时序分析——通过直观的数据比较或绘图观测，寻找 序列中蕴含的发展规律，这种分析方法 就称为描述性时序分析

• 描述性时序分析方法具有操作简单、直观有效的特点，它通常是人们进行统计 时序分析的第一步。 2、统计时序分析

（1）频域分析方法 • 原理：假设任何一种无趋势的时间序列都可以分解成若干不同频率的周期波动 • 发展过程： 1）早期的频域分析方法借助富里埃分析从频率的角度揭示时间 序列的规律 2）后来借助了傅里叶变换，用正弦、余弦项之和来逼近某个函数 3）20世纪60年代，引入最大熵谱估计理论，进入现代谱分析阶段 • 特点：非常有用的动态数据分析方法，但是由于分析方法复杂，结 果抽象，有一定的使用局限性 （2）时域分析方法

• 原理：事件的发展通常都具有一定的惯性，这种惯性用统 计的语言来描述就是序列值之间存在着一定的相关 关系，这种相关关系通常具有某种统计规律。

• 目的：寻找出序列值之间相关关系的统计规律，并拟合出 适当的数学模型来描述这种规律，进而利用这个拟 合模型预测序列未来的走势

• 特点：理论基础扎实，操作步骤规范，分析结果易于解 释，是时间序列分析的主流方法 楼主，说了半天，你终于到正题了，时域分析才是我们经常接触的，你赶紧说说怎么做吧？ ★时域分析方法的分析步骤： • 考察观察值序列的特征 • 根据序列的特征选择适当的拟合模型 • 根据序列的观察数据确定模型的口径 • 检验模型，优化模型 • 利用拟合好的模型来推断序列其它的统 计性质或预测序列将来的发展

时域分析方法的发展过程 • 基础阶段—— G.U.Yule：1927年，AR模型 G.T.Walker：1931年，MA模型，ARMA模型 • 核心阶段——G.E.P.Box和G.M.Jenkins 1970年，出版《Time Series Analysis Forecasting and Control》 提出ARIMA模型（Box—Jenkins 模型） Box—Jenkins模型实际上是主要运用于单变量、同方差场合的线性模型 • 完善阶段—— 异方差场合：Robert F.Engle，1982年，ARCH模型 Bollerslov，1985年GARCH模型 多变量场合：C.Granger ，1987年，提出了协整（co- integration）理论 非线性场合：汤家豪等，1980年，门限自回归模型

用哪些软件可以做时间序列分析呢？ S-plus，Matlab，Gauss，TSP，Eviews 和SAS 上述软件楼主觉得Eviews是基础版，Gauss是小众版，Matlab&S-pluss是正常小青年~~SAS，万能的软件BOSS啊~~~ 下一辑——时间序列的预处理！敬请关注！ 【时间简“识”】2.那些必不可少的预处理 - 计量经济学与统计软件 - 经管之家(原人大经济论坛) bbs.pinggu.org 2012-7-27 本帖最后由

经管之家(原人大经济论坛)

胖胖小龟宝 于 2014-12-12 09:12 编辑 上一辑预告说啦~~本期的主题是时间序列的预处理~~序列在建模前到底要做哪些预处理呢？首先，大伙都知道的平稳性检验是必须的！ 说到平稳，其实有两种平稳—— 宽平稳、严平稳 严平稳相较于宽平稳来说，条件更多更严格，而我们时常运用的时间序列，大多宽平稳就够了~~ 什么是严平稳： 是在固定时间和位置的概率分布与所有时间和位置的概率分布相同的随机过程。这样，数学期望和方差这些参数也不随时间和位置变化。（比如白噪声） 什么是宽平稳： 宽平稳是使用序列的特征统计量来定义的一种平稳性。它认为序列的统计性质主要由它的低阶矩决定，所以只要保证序列低阶矩平稳（二阶），就能保证序列的主要性质近似稳定。

两者关系： 一般关系： 严平稳条件比宽平稳条件苛刻，通常情况下，严平稳（低阶矩存在）能推出宽平稳成立，而宽平稳序列不能反推严平稳成立。 特例： 不存在低阶矩的严平稳序列不满足宽平稳条件，例如服从柯西分布的严平稳序列就不是宽平稳序列。当序列服从多元正态分布时，宽平稳可以推出严平稳。

如何判断序列是平稳的？ 咱们这次先从图形法上看（通常越是简单的方法，往往越能看到问题，图形给出的第一感觉也许就是真相哦~~~~） 时序图，例如(eviews画滴)： 分析：什么样的图不平稳，先说下什么是平稳，平稳就是围绕着一个常数上下波动。 看看上面这个图，很明显的增长趋势，不平稳。 我们还可以根据自相关和偏相关系数来查看： 还以上面的序列为例：用eviews得到自相关和偏相关图，Q统计量和伴随概率。 分析：平稳的序列的自相关图和偏相关图不是拖尾就是截尾。截尾就

是在某阶之后，系数都为 0 ，怎么理解呢，看上面偏相关的图，当阶数为 1 的时候，系数值还是很大， 0.914. 二阶长的时候突然就变成了 0.050. 后面的值都很小，认为是趋于 0 ，这种状况就是截尾。再就是拖尾，拖尾就是有一个衰减的趋势，但是不都为 0 。 自相关图既不是拖尾也不是截尾。以上的图的自相关是一个三角对称的形式，这种趋势是单调趋势的典型图形。

下面是通过自相关的其他功能 如果自相关是拖尾，偏相关截尾，则用 AR 算法 如果自相关截尾，偏相关拖尾，则用 MA 算法 如果自相关和偏相关都是拖尾，则用 ARMA 算法， ARIMA 是 ARMA 算法的扩展版，用法类似 。

一定有同学要问了：楼主检测出来不是平稳的怎么办啊？（楼主：当然要把它整平稳啦~~） 如果遇到数据检测出来不平稳，可以考虑使用差分这个最常用的办法（当然，还有好多种其他方法处理） 还是上面那个序列，两种方法都证明他是不靠谱的，不平稳的。确定不平稳后，依次进行1阶、2阶、3阶...差分，直到平稳位置。先来个一阶差分： 从图上看，一阶差分的效果不错，看着是平稳的。

在图形检验法中，我们能够较为直观的看到数据的一个大致变动趋势，如果它有周期或者上升等趋势，一般就不太平稳，需要做些处理，但图形始终是个主管判断为主的方法，这次，就来说说平稳检验的另一个方法：单位根检验（ADF检验）。 ADF检验简介： 检查序列平稳性的标准方法是单位根检验。有6种单位根检验方法：ADF检验、DFGLS检验、PP检验、KPSS检验、ERS检验和NP检验，本节将介绍DF检验、ADF检验。ADF检验和PP检验方法出现的比较早，在实际应用中较为常见，但是，由于这2种方法均需要对被检验序列作可能包含常数项和趋势变量项的假设，因此，应用起来带有一定的不便；其它几种方法克服了前2种方法带来的不便，在剔除原序列趋势的基础上，构造统计量检验序列是否存在单位根，应用起来较为方便。ADF检验是在Dickey-Fuller检验(DF检验)基础上发展而来的。因为DF检验只有当序列为AR(1)时才有效。如果序列存在高阶滞后相关，这就违背了扰动项是独立同分布的假设。在这种情况下，可以使用增广的DF检验方法（augmented Dickey-Fuller test ）来检验含有高阶序列相关的序列的单位根。

检验步骤（一般进行ADF检验要分3步）： 1 对原始时间序列进行检验，此时第二项选level，第三项选None.如果没通过检验，说明原始时间序列不平稳； 2 对原始时间序列进行一阶差分后再检验，即第二项选1st difference，第三项选intercept,若仍然未通过检验，则需要进行二次差分变换； 3 二次差分序列的检验，即第二项选择2nd difference ，第四项选择Trend and intercept.一般到此时间序列就平稳了！ tips: 在进行ADF检验时，必须注意以下两个实际问题： （1）必须为回归定义合理的滞后阶数，通常采用AIC准则来确定给定时间序列模型的滞后阶数。在实际应用中，还需要兼顾其他的因素，如系统的稳定性、模型的拟合优度等。 （2）可以选择常数和线性时间趋势，选择哪种形式很重要，因为检验显著性水平的 t 统计量在原假设下的渐近分布依赖于关于这些项的定义。 ① 若原序列中不存在单位根，则检验回归形式选择含有常数，意味着所检验的序列的均值不为0；若原序