时间序列初探—平稳性分析及R实现

1、1基本概念时间序列的平稳性假定某个时间序列是由某一随机过程 <stochastic process ）生 成的，即假定时间序列Xt<t=1,2,）的每一个数值都是从一个概率分布中随机得到，如果满足下列条件：b5E2RGbCAP1 ）均值E（Xt>= 是与时间t无关的常数；2 ）方差Var（Xt>= 2是与时间t无关的常数；3 ）协方差Cov（Xt,Xt+k>= k是只与时期间隔k有关，与时 间t无关的常数；则称该随机时间序列是平稳的 stationary，而该随机过程是一 平后急随机过程 <stationary stochastic process ） 。

2、plEanqFDPw 时间序列的非平稳性平稳时间序列的均值为常数，自协方差函数与起点无关，而非平 稳时间序列则不满足这两条要求。常见的非平稳类型有趋势和突变 DXDiTa9E3d趋势趋势是指变量随时间持续长期的运动，时间序列变量围绕其趋势波动。可以用线性趋势、二次趋势、季节性均值趋势和余弦趋势来估计一般的非常数均值趋势模型的参数。RTCrpUDGiT突变突变来自总体回归系数在某一特定日期上的离散变化或来自系数 在长时期内的渐变。平稳性判断图示判断?给出一个随机时间序列，首先可通过该序列的时间路径图来粗 略地判断它是否是平稳的。? 一个平稳的时间序列在图形上往往表现出一种围绕其均值不断 波动的过

3、程；?而非平稳序列则往往表现出在不同的时间段具有不同的均值<如持续上升或持续下降）。函数1:时间序列及趋势绘制 参数1:时间序列 功能：绘制时间序列绘制时间序列的趋势函数返回值：无 单位根检验单位根检验<unit root test ）是针对宏观经济数据序列、货币金融数据序列中是否具有某种统计特性而提出的一种平稳性检验的特殊方法，单位根检验的方法有很多种，包括 DF检验、ADF检验、PP检验、NP检验等。5PCzVD7HxA单位根检验时间序列的单位根研究是时间序列分析的一个热点问题。时间序列特性的时变行为实际上反映了时间序列的非平稳性质。对非平稳时间序列的处理方法一般是将其转变为平

4、稳序列，这样就可以应用有关平稳时间序列的方法来进行相应得研究。对时间序列单位根的检验就是对时间序列平稳性的检验，非平稳时间序列如果存在单位根，则一般可以通过差分的方法来消除单位根，得到平稳序列。对于存在单位根的时间序列，一般都显示出明显的记忆性和波动的持续性，因此单位根检验是有关协整关系存在性检验和序列波动持续性讨论的基础。jLBHrnAILg自相关函数（ACF冽断平稳时间序列的自相关函数（ACF懑么是截尾的，要么是拖尾的。因此我们可以根据这个特性来判断时间序列是否为平稳序列。若时间序列具有上升或下降的趋势，那么对于所有短时滞来说，自相 关系 数大 且为正 ，而 且随 着时滞 k 的增 加而缓

5、慢地下 降。xHAQX74J0X若序列无趋势，但是具有季节那么于按月采集的数据，时滞12, 24, 36的自相关系数达到最大（如果数据是按季度采集，则最大自相关系数出现在 4， 8， 12，>，并且随着时滞的增加变得较小。LDAYtRyKfE平稳时间序列模型自回归AR模型由于经济系统惯性的作用，经济时间序列往往存在着前后依存关系。最简单的一种前后依存关系就是变量当前的取值主要与其前一时期的取值状况有关。用数学模型来描述这种关系就是如下的一阶自回归模型：Zzz6ZB2LtkXt=（|）Xt-1+ et<2.1.1）常记作AR（1>。其中Xt为零均值 <即已中心化处理）平稳

6、序列，0为Xt对Xt 1的依赖程度，£ t为随机扰动项序列 <外部冲击）。 dvzfvkwMI1若0 =1则Xt包含了一个随机性趋势，是非平稳的。若 0的绝对值 <1 则其是平稳的。如果 Xt 与过去时期直到Xt-p 的取值相关，则需要使用包含Xt 1 ,Xt-p在内的p阶自回归模型来加以刻画。P阶自回归模型的一般形式为：rqyn14ZNXIXt= 0 1 Xt 1+0 2 Xt 2+ d p Xt -p+£ t <2.1.2）为了简便运算和行文方便，我们引入滞后算子来简记模型。设 B 为 滞 后 算 子 ， 即 BXt=Xt-1, 则 B（Bk-1Xt&

7、gt;=BkXt=Xt-k B（C>=C（C 为 常 数 >。 利 用 这 些 记 号 ， <2.1.2 ） 式 可 化 为 ： EmxvxOtOcoXt= d 1BXt+ d 2B2Xt+ d 3B3Xt+ © pBpXt+ £ t从而有：<1- j 1B-2B2-j pBp) Xt= & t记算子多项式0<B) = <1-0 1B-0 2B2-0pBP)，则模型可以表示成0<B) Xt= £ t(2.1.3>例如，二阶自回归模型Xt=0.7Xt-1+0.3Xt-2+0.3Xt-3+£ t可写成&

8、lt;1-0.7B-0.3B2 ) Xt= & tSixE2yXPq5若AR(P消一个等于1的根，则称序列有一个单位自回归根或称为单位根，从而也说明它包含了随机性趋势，是非平稳的。当且仅当AR 特 征方 程的每一个根 绝对 值大于 1 ， 时间 序列是 平稳 的。6ewMyirQFL滑动平均模型<MA)有时，序列Xt 的记忆是关于过去外部冲击值的记忆，在这种情况下 ， Xt 可 以表示成过去 冲击 值和现 在冲 击值 的线性 组合 ，即kavU42VRUsXt= £ t- 0 1et-1- 0 2et-2- -0 q £ t-q (2.1.4>此模型常称

9、为序列Xt 的滑动平均模型，记为MA(q>，其中 q 为滑动平均的阶数，0 1,0 2- 0 q为参滑动平均的权数。相应的序列Xt 称为滑动平均序列。y6v3ALoS89使用滞后算子记号，<2.1.4 )可写成Xt=<1- 0 1B- 0 2B2- -0 qBq) qt= 0 (B> £ t (2.1.5>自回归滑动平均模型<ARM) A如果序列 Xt的当前值不仅与自身的过去值有关，而且还与其以前进入系统的外部冲击存在一定依存关系，则在用模型刻画这种动态特征时，模型中既包括自身的滞后项，也包括过去的外部冲击，这种模型叫做自回归滑动平均模型，其一般结

10、构为：M2ub6vSTnPXt= 0 1Xt-1+ 0 2Xt- 2+ + 0 pXt-p+ £ t- 0 1 £ t-1- 0 2 £ t-2- 0 q e t-q (2.1.6>0YujCfmUCw简记为ARMA(p, q>。利用滞后算子，此模型可写为0 <B) Xt= 0 (B> e t<2.1.7)R中实现判断时间序列的平稳性例一> x=rnorm(500># 生成 500 个服从正太分布的数> y=cumsum(x> #累加x 的数对应得到y绘制时序图> plot.ts(x>Time>

11、; plot.ts(y>Time从两个图的不同可以看出 x时间序列趋势不随时间的变化而变化, 其随机性比较强。而y序列则有明显的时间趋势。eUts8ZQVRd ADF.test 检验install.packages("tseries"># 安装时间序列包library("tseries",lib.loc="e:/ProgramFiles/R/R-2.15.2/library"># 载入时间序列包sQsAEJkW5T> adf.test(x>Augmented Dickey-Fuller Testdata:

12、 xDickey-Fuller = -8.0878, Lag order = 7, p-value0.01GMsIasNXkAalternative hypothesis: stationary结论： p-value = 0.01 拒绝原假设<原假设认为时间序列是非平稳的)，即可认为x 是平稳的。> adf.test(y>Augmented Dickey-Fuller Testdata: yDickey-Fuller = -1.1291, Lag order = 7, p-value0.9179TIrRGchYzgalternative hypothesis: station

13、ary结论： p-value = 0.9179 不能拒绝原假设，所以认为y 是非平稳的。函数2: ADF检验时间序列的平稳性：ADFTEST参数 1：时间序列P 临界值，默认值为0.05返回结果：用框架来组织返回结果结论<1：平稳，0：不平稳)adf.test 函数的返回值PP检验> pp.test(x>Phillips-Perron Unit Root Testdata: xDickey-Fuller Z(alpha> = -510.4566, Truncation lagparameter = 5, p-value = 0.017EqZcWLZNXalternati

14、ve hypothesis: stationary警告信息：In pp.test(x> : p-value smaller than printed p-valuelzq7IGf02E结论：p-value = 0.01拒绝非平稳性假设，即认为 x是平稳的。> pp.test(y>Phillips-Perron Unit Root Testdata: yDickey-Fuller Z(alpha> = -3.9888, Truncation lag parameter= 5, p-value = 0.8872zvpgeqJ1hkalternative hypothesis

15、: stationary结论：p-value = 0.8872 不能拒绝原假设y是非平稳的，所以认为 y 是非平稳的。函数3: PP检验时间序列的平稳性：PPTEST参数1：时间序列P 临界值，默认值为0.05返回结果：用框架来组织返回结果结论<1：平稳，0：不平稳)R语言pp检验函数的返回值ACF自相关函数判断> modelx=lm(xtime(x>>> summary(modelx>Call:lm(formula = x time(x>>Residuals:Min 1Q Median 3Q Max-2.87920 -0.75003 0.011

16、03 0.70595 3.15625Coefficients:Estimate Std. Error t value Pr(>|t|>(Intercept> 0.1524849 0.0915359 1.6660.0964 .NrpoJac3v1time(x> -0.0005077 0.0003166 -1.603 0.1095-Signif. codes: 0 * 0.001 * 0.01 * 0.05 . 0.1 1 1nowfTG4KIResidual standard error: 1.022 on 498 degrees offreedomfjnFLDa5ZoM

17、ultiple R-squared: 0.005136, Adjusted R-squared: 0.003138 tfnNhnE6e5F-statistic: 2.571 on 1 and 498 DF, p-value: 0.1095HbmVN777sL > acf(rstudent(modelx>,main=' 关于 x 的 acf 自相关系数'>V7l4jRB8Hs从图中可以看出其 K阶滞后自相关系数都非常小呈截现象，因此判断时间系列为平稳性是合理的。函数4: ACF检验函数参数：时间序列检验p值，默认为0.05图形保存路径，默认为空返回值：以框架形式

18、线性回归函数各个系数的检验p 值ACF 函数的返回值> modely=lm(ytime(y>>> summary(modely>Call:lm(formula = y time(y>>Residuals:Min 1Q Median 3Q Max-13.2206 -5.6292 -0.6742 6.3185 13.6971Coefficients:Estimate Std. Error t value Pr(>|t|>(Intercept> 8.659376 0.620538 13.96 <2e-16*83lcPA59W9time

19、(y> 0.032966 0.002146 15.36 <2e-16*mZkklkzaaP-Signif. codes: 0 * 0.001 * 0.01 * 0.05 . 0.1 1 AVktR43bpwResidual standard error: 6.927 on 498 degrees of freedomORjBnOwcEdMultiple R-squared: 0.3214,Adjusted R-squared: 0.32012MiJTy0dTTF-statistic: 235.9 on 1 and 498 DF, p-value: < 2.2e-16gIiSp

20、iue7A> acf(rstudent(modely>,main=' 关于 y 的 acf 自相关系数'>从图中可以看出ACF随着k的增大而缓慢下降，自相关系数大且为正因此判断y序列为非平稳时间序列是合理的。uEh0U1Yfmh例二 以TSA自带的数据tempdub为例验证数据的平稳性检验绘图>library("TSA",lib.loc="e:/ProgramFiles/R/R- 2.15.2/library”>IAg9qLsgBX从图中可以看出此时间序列具有非常明显的周期性趋势WwghWvVhPEAdf检验>

21、adf.test(tempdub>Augmented Dickey-Fuller Testdata: tempdubDickey-Fuller = -11.0773, Lag order = 5, p-value =0.01asfpsfpi4kalternative hypothesis: stationary结论：p-value = 0.01 所以tempdub时间序列是平稳的。Pp检验> pp.test(tempdub>Phillips-Perron Unit Root Testdata: tempdubDickey-Fuller Z(alpha> = -51.0795, Truncation lag parameter= 4, p-value = 0.01ooeyYZTjj1 alternative hypothesis: stationary警告信息：In pp.test(tempdub> : p-value smaller than printed p- valueBkeGuInkxI结论： p-v