非平稳时间序列模型

引言：前面我们讨论的是平稳时间序列的建模和预测方法，即所讨论的时间序列都是宽平稳的。一个宽平稳的时间序列的均值和方差都是常数，并且它的协方差有时间上的不变性。但是许多经济领域产生的时间序列都是非平稳的。对协方差过程，非平稳时间序列会出现各种情形，如它们具有非常数的均值μt，或非常数的二阶矩，如非常方差σt2，或同时具有这两种情形的非平稳序列。第七章非平稳时间序列模型第七章非平稳时间序列模型第一节非平稳时间序列模型的种类第二节非平稳性的检验第三节求和自回归滑动平均模型(ARIMA)第一节非平稳时间序列模型的种类一、均值非平稳过程二、方差和自协方差非平稳过程返回本节首页下一页上一页一、均值非平稳过程均值非平稳过程指随机过程的均值随均值函数的变化而变化。我们可以引进两种非常有用的均值非平稳过程：确定趋势模型和随机趋势模型。返回本节首页下一页上一页（一）确定趋势模型当非平稳过程均值函数可由一个特定的时间趋势表示时，一个标准的回归模型曲线可用来描述这种现象。.,::,,1010模型来描述前面介绍的可以用程是一个零均值的平稳过其中趋势模型表示如下则原序列可用确定的有服从线性趋势若均值例如ARMAyytxtttttttttyttxtt22102210:,原序列可用下式表示对二次均值函数此外，均值函数还可能是指数函数、正弦—余弦波函数等，这些模型都可以通过标准的回归分析处理。处理方法是先拟合出μt的具体形式，然后对残差序列yt={xt－μt}按平稳过程进行分析和建模。（二）随机趋势模型随机趋势模型又称齐次非平稳ARMA模型。为理解齐次非平稳ARMA模型，可先对ARMA模型的性质作一回顾。.1)(1)(:)()(:),(221221为白噪声序列其中模型如下假设有一个tqqpttaBBBBBBBBaBxBqpARMA.,,0)(.0)(:,就是非平稳的么那的根不都在单位圆外如果根都在单位圆外的则必须有为满足平稳性txBBttddaBxBBBBBdB)()1)((:)1)(()(:,,0)(于是原模型可写为则可令而其它根都在单位圆外个根落在单位圆上恰有现假设.)()()(:,)1(.,运算后可变为平稳序列差分次程经过若干次可见一个齐次非平稳过则令称为齐次性的阶为齐次非平稳过程这时我们就称daBwBxBwdxtttdtt可见我们所能分析处理的仅是一些特殊的非平稳序列，即齐次非平稳序列。由于齐次非平稳序列模型恰有d个特征根在单位圆上，即有d个单位根，因此齐次非平稳序列又称单位根过程。二、方差和自协方差非平稳过程一个均值平稳过程不一定是方差和自协方差平稳过程，同时一个均值非平稳过程也可能是方差和自协方差非平稳过程。不是所有的非平稳问题都可以用差分方法解决，还有期望平稳和方差非平稳序列，为了克服这个问题，我们需要适当进行方差平稳化变换。返回本节首页下一页上一页010ln:,)(tttxxx表示如下稳一般用幂变换使方差平这个变换最早由BOX和COX于1964年提出，因此称作BOX—COX变换。其中λ为变换参数。第二节非平稳性的检验一、通过时间序列的趋势图来判断二、通过自相关函数(ACF)判断三、特征根检验法四、用非参数检验方法判断序列的平稳性五、随机游走的单位根检验返回本节首页下一页上一页一、通过时间序列的趋势图来判断这种方法通过观察时间序列的趋势图来判断时间序列是否存在趋势性或周期性。优点：简便、直观。对于那些明显为非平稳的时间序列，可以采用这种方法。缺点：对于一般的时间序列是否平稳，不易用这种方法判断出来。返回本节首页下一页上一页二、通过自相关函数(ACF)判断平稳时间序列的自相关函数(ACF)要么是截尾的，要么是拖尾的。因此我们可以根据这个特性来判断时间序列是否为平稳序列。若时间序列具有上升或下降的趋势，那么对于所有短时滞来说，自相关系数大且为正，而且随着时滞k的增加而缓慢地下降。返回本节首页下一页上一页若序列无趋势，但是具有季节性，那末对于按月采集的数据，时滞12，24，36……的自相关系数达到最大(如果数据是按季度采集，则最大自相关系数出现在4，8，12，……)，并且随着时滞的增加变得较小。若序列是有趋势的，且具有季节性，其自相关函数特性类似于有趋势序列，但它们是摆动的，对于按月数据，在时滞12，24，36，……等处具有峰态；如果时间序列数据是按季节的，则峰出现在时滞4，8，12，……等处。三、特征根检验法（P146）.1;,1:,,,:则该序列就是非平稳的如果的则可以认为序列是平稳即件都满足平稳性条若所有的特征根特征根组成的特征方程的求由该适应模型的参数然后先拟合序列的适应模型基本思想返回本节首页下一页上一页根据拟合出的时序模型参数检验(P146)基本思想：时间序列模型的平稳性条件不仅可以用特征根来表示，也可以用模型的自回归参数表示，因此要检验一个序列是否平稳，可以先拟合适应的模型，然后再根据求出的自回归参数来检验序列是否平稳。检验方法：参见课本146.;,.1:),(111:,),2(1,),1(,.21212211作进一步的判定如果满足则需要列原序列肯定为非平稳序那么即如果上述条件不满足种便捷的途径是否平稳提供了一这为我们判断时间序列都有如下必要条件模型要满足平稳性一般的平稳性条件为模型平稳性条件为模型如例如时间序列模型较为方便归阶数较低的这种检验方法对于自回pqpARMAqARMAqARMA四、用非参数检验方法判断序列的平稳性（一）什么是参数检验和非参数检验？参数检验：参数检验是这样一种检验，它的模型对抽出研究样本的总体的分布作了限制性假定。。如果对总体的分布不知道或了解很少，则参数检验方法就不可靠，甚至会发生较大偏差。返回本节首页下一页上一页非参数检验：非参数检验是一种不依赖于总体分布知识的检验方法。由于非参数检验不对总体分布加以限制性假定，所以它也称为自由分布检验。非参数检验与参数检验相比有如下优点：a.检验条件比较宽松，适应性强。b.参数检验对样本容量的要求极低。c.检验方法灵活，用途更广泛。非参数检验主要用顺序统计量进行检验，因此它既可检验定距数据和定比数据，又可以检验定类数据和定序数据；而参数检验只能处理定距数据和定比数据。因为这些优点，非参数检验比参数检验应用更广泛。d.非参数检验计算相对简单，易于理解。非参数检验的缺点：如果参数统计模型的所有假设在数据中事实上都能满足，而且测量达到了所要求的水平(定距数据或定比数据)，那么用非参数检验就浪费了数据中的信息。也就是说此时非参数检验的功效不如参数检验高。(二)非参数检验方法在检验序列平稳性中的应用1.游程检验方法(1)什么是游程一个游程定义为一个具有相同符号的连续串，在它前后相接的是与其不同的符号或完全无符号。例如，观察的结果用加、减标志表示，得到一组这样的记录顺序：++---+----++-+这个样本的观察结果共有7个游程。（2）用游程检验方法检验时间序列平稳性的基本思想..,,,},{程数并可求出这个序列的游序列这样就形成了一个符号号的观察值记为大比号小的观察值记为对序列中比设其样本均值为对于一个时间序列xxxxt如果符号序列是随机的，那么“+”和“-”将随机出现，因此它的游程数既不会太多，又不会太少；反过来说如果符号序列的游程总数太少或太多，我们就可以认为时间序列存在某种趋势性或周期性。.)1,0()()(:)15()1()12(2)(12)(::,,212212121212121服布渐近服从有大于或在大样本情况下的期望和方差分别如下数游程总明，对于随机序列可以证总数为出现的次数，游程与为记号序列中分别和设序列长度为NrDrErZNNNNNNNNrDNNNNrErrNNNNNN(3)检验方法a.小样本情况零假设：H0：加号和减号以随机的方式出现检验方法：取显著性水平α(一般取0.05),查单样本游程检验表，得出抽样分布的临界值rL、rU判定：若rLrrU则不能拒绝零假设，即不能拒绝序列是平稳的；若rrL或rrU则拒绝零假设，序列是非平稳的。b.大样本情况零假设：H0：加号和减号以随机的方式出现检验方法：给定显著性水平α(一般取0.05)查标准正态分布表，得出抽样分布的临界值-zα,+zα。并计算统计量:)()(rDrErZ判定：若-zαz+zα,则不能拒绝零假设，即不能拒绝序列是平稳的;否则拒绝零假设，序列是非平稳的。（例见P151例6.5）非参数检验可以很方便的通过SPSS软件进行，游程检验可见操作。实例：用游程检验S&T数据的平稳性；步骤如下：1.打开SPSS输入数据2.依次单击Analyze—NonparmetricTests—Runs;打开Runs对话框。3.在源变量对话框中选择“stpoor”进入“TestVariablelist”栏内4.选中“cutpoint”栏中“mean”选项5.单击“OK”按纽，开始进行统计分析。RunsTest288.9746102921944-13.532.000TestValueaCasesTestValueCases=TestValueTotalCasesNumberofRunsZAsymp.Sig.(2-tailed)STPOORMeana.输出结果分析：因为P值（sig.）极小，所以拒绝零假设，故原序列是非平稳的。也可以通过其它的非参数检验方法来判断序列是否平稳，如Spearman等级相关系数，Kendallτ相关系数等。五、随机游走的单位根检验(Unitroottest)在第三章我们已经讲过，随机游走是一种非平稳过程，其实随机游走一种特殊的齐次非平稳过程。检验序列是否为随机游走，通常利用DavidDickey和WayneFuller的单位根检验。单位根的含义和检验原理如下：返回本节首页下一页上一页.,,1:..),,,(:11间序列序列就叫做随机游走时一个有位单位根的时间根有一个单位那么就说随机变量并确实发现如果我们作如下回归非平稳的情况它是一种的单位根问题这时我们就遇到了所谓无自相关恒定方差点零均值为白噪声其中可由下式描述假设序列tttttttttyuyyuuyyy)(.)(.,0ˆ.;,,0ˆ,ˆ)1(:11检验由此得名统计量进行检验算的计和这时需要根据分布不服从值并通常计算的下然而要在原假设序列是平稳的否则它是非平稳的是随机游走过程则如果检验出对上述模型作回归下等价形式随机游走模型可写成如DFtauFullerDickeyttyuyuyytttttt检验时：如果所计算的τ统计量的绝对值(即|τ|)大于DF分布表中临界值的绝对值，则拒绝δ=0的假设，原序列是平稳的；否则，如果它小于临界值，则时间序列是非平稳的。注：Eviews输出结果中直接计算出了τ统计量及其临界值。（所列出的是麦金农MacKinnon对DF分布表扩充后的临界值）.,0,).(),(int:)3()2()1(:,1121111存在一个单位根即原假设都是在上面每一种有形中为趋势项为常数项其中检验以下形式的回归作人们常用原因由于理论上和实践上的trendterceptuytyuyyuyyFullerDickeyttttttttt用Eviews进行单位根检验时给出了上述选项。.)().(,0.,,,:)3(,2111121检验该检验被称为是非平稳的即个单位根有一即时该检验的原假设仍是独立的要做到使上式中少个滞后的差分项模型中要包含多等如分项上式中包含有滞后的差式修改如下就把上述如果误差项是自相关的FullerDickeyAugmentedADFyyyyyyytytttttttmiititt1.单位根检验例可见操作。2.补充利用Ev