-平稳时间序列预测法.ppt

7 7 平稳时间序列预测法平稳时间序列预测法7.1 概述概述7.2 时间序列的自相关分析时间序列的自相关分析7.3 单位根检验和协整检验单位根检验和协整检验7.4 ARMA模型的建模模型的建模.7.1 概概 述述时间序列取自某一个随机过程，则称：一、平稳时间序列一、平稳时间序列过程是平稳的过程是平稳的随机过程的随机特征不随时间变化而变化过程是非平稳的过程是非平稳的随机过程的随机特征随时间变化而变化.宽平稳时间序列的定义：宽平稳时间序列的定义：设时间序列，对于任意的t,k和m，满足：则称 宽平稳。严平稳时间序列的定义：严平稳时间序列的定义：所有的统计特性不随时间的平移而变化.qBox-Jenkins基本思想：用数学模型描述时间序列自身的相关性，并假定这种自相关性一直延续，用该模型预测未来的值。qARMA模型是描述平稳随机序列的最常用的一种模型。qBox-Jenkins方法提供了对时间序列进行分析、预测，以及对ARMA模型识别、估计和诊断的系统方法。.ARMA模型的三种基本形式：q 自回归模型（AR：Auto-regressive）；q 移动平均模型（MA：Moving-Average）；q 混合模型（ARMA：Auto-regressiveMoving-Average）。.如果时间序列满足其中是独立同分布的随机变量序列,且满足：则称时间序列 服从p阶自回归模型。二、自回归模型二、自回归模型.滞后算子多项式的根均在单位圆外，即的根大于1。自回归模型的平稳条件：自回归模型的平稳条件：.例例1 AR(1)模型的平稳性条件。模型的平稳性条件。对1阶自回归模型AR(1)方程两边平方再求数学期望，得到Xt的方差由于Xt仅与t相关，因此，E(Xt-1t)=0。如果该模型稳定，则有E(Xt2)=E(Xt-12)，从而上式可变换为：在稳定条件下，该方差是一非负的常数，从而有|1。.而AR(1)的特征方程的根为 z=1/AR(1)稳定，即|1，意味着特征根大于1。对高阶自回模型对高阶自回模型AR(p)来说来说，多数情况下没有必要直接计算其特征方程的特征根，但有一些有一些有用的规则可用来检验高阶自回归模型的稳定性用的规则可用来检验高阶自回归模型的稳定性：(1)AR(p)模型稳定的必要条件是模型稳定的必要条件是：1+2+p1 (2)(2)由于i(i=1,2,p)可正可负，AR(p)模型模型稳定的充分条件是：稳定的充分条件是：|1|+|2|+|p|p，t与t-k间的偏自相关系数偏自相关系数为零。.样本的偏自相关函数的计算样本的偏自相关函数的计算.其中：.1、时间序列的随随机机性性，是指时间序列各项之间没有相关关系的特征。使用自相关分析图判断时间序列的随机性，一般给出如下准则：q若时间序列的自相关函数基本上都落入置信区间，则该时间序列具有随机性；q若较多自相关函数落在置信区间之外，则认为该时间序列不具有随机性。时间序列特性分析时间序列特性分析注：注：在B-J方法中，测定时间序列的随机性，多用于模型残差，以评价模型优劣。.2、判断时间序列是否平稳是否平稳，是一项很重要的工作。运用自相关分析图判定时间序列平稳性的准则是：q若时间序列的自相关函数在k3时都落入置 信区间 ，且逐渐趋于零，则该时间序列具有平稳性；q若时间序列的自相关函数更多地落在置信区间外面，则该时间序列不具有平稳性。.注：注：在B-J方法中，只有平稳的时间序列平稳的时间序列才能建立ARMA模型，否则必须经过适当的处理使序列满足平稳性要求。例对某种趋势的时间序列进行差分处理。但很多序列不能通过差分达到平稳，而且差分虽然消除了序列的趋势易于建模，但也消除了序列的长期特征，实际的经济序列差分一般不超过两次。.3、时间序列的季节性季节性判定准则：q月度数据，考察k=12,24,36,时的自相关系数是否与0有显著差异；q季度数据，考察k=4,8,12,时的自相关系数是否与0有显著差异。注注1 1：实际问题中常遇到季节性和趋势性同时存在的情况，应先剔除序列趋势性，在识别季节性。注注2 2：包含季节性的时间序列也不能直接建模，应先进行季节差分消除，季节差分一般不超过一阶。.三、三、ARMA模型的自相关分析模型的自相关分析 qAR(p)模型的偏自相关函数是以p步截尾的，自相关函数拖尾；qMA(q)模型的自相关函数具有q步截尾性，偏自相关函数拖尾；（可用以上两个性质来识别（可用以上两个性质来识别AR和和MA模型的阶数）模型的阶数）qARMA(p,q)模型的自相关函数和偏相关函数都 是拖尾的。.图图 ARMA(p,q)模型的模型的 ACF 与与 PACF理论模式理论模式 ACF PACF 模型模型1：tttXXe+=-17.00.00.20.40.60.812345678ACF10.00.20.40.60.812345678PACF1.7.4 ARMA模型的建模模型的建模 一、一、模型阶数的确定模型阶数的确定（1）基于自相关函数和偏相关函数的定阶方法基于自相关函数和偏相关函数的定阶方法 对于ARMA(p,q)模型，可以利用其样本的自相关函数和样本偏自相关函数的截尾性判定模型的阶数。q如果样本的偏自相关函数是以p步截尾的，模型为AR(p)；q如果样本的自相关函数具有q步截尾性，模型为MA(q)；q如果样本的自相关函数和偏相关函数都是拖尾的，模型为ARMA(p,q)。.（1）自相关函数的截尾性统计检验：）自相关函数的截尾性统计检验：q对于每一个q,计算.（M 取为或者 ），对于MA(q)模型，当kq时，q q考察其中满足 的个数是否占M个的68.3%或者95.5%以上。.（2）偏自相关函数的截尾性统计检验：）偏自相关函数的截尾性统计检验：q对于每一个p,计算（M 取为或者 ），对于AR(p)模型，当kp时，q q考察其中满足 的个数是否占M个的68.3%或者95.5%以上。.q 如果对于序列和截尾，即不存在上述的来说，均不和判定平稳时间序列，则可以为ARMA模型。.一般地，对ARMA模型 它们均值为0，可递推得到残量估计现作假设检验：是来自白噪声的样本 令（3）残差项的白噪声检验：（）残差项的白噪声检验：（Q统计量检验）统计量检验）.其中取 左右。则当成立时，服从的分布。对给定的显著性水平，若，则拒绝，即模型与原随机序列之间拟合得不好，则认为模型与原随机序列之间拟合需重新考虑得较好，模型检验被通过。建模；若自由度为注：上机操作时，一般看Q统计检验的相伴概率.(1)用AR（1）拟合时间序列，考察其残差样本的自相关函数 是否q1 1步截尾，则模型为ARMA(1,q1 1),否则；(2)用AR（2）拟合时间序列，考察其残差样本的自相关函数 是否q2 2步截尾，则模型为ARMA(1,q2 2)，否则；(3)继续增大p，重复上述做法，直至残差序列的样本自相关 函数截尾为止（4）Tasy和TiaoARMA模型定阶法.1950年-1998年北京城乡居民定期储蓄比例选 择 合 适的ARMA模型拟合可 以 考 虑拟 合 模 型为AR(1),ARMA(1,3).连续读取70个化学反应数据可以尝试使用AR(1)，MA(1)和ARMA(1,1)模型拟合该序列.（2 2）基于）基于F F 检验确定阶数检验确定阶数（3 3）利用信息准则法定阶（）利用信息准则法定阶（AICAIC准则和准则和BICBIC准则）准则）此外，常用的方法还有：.1967年，瑞典控制论专家K.J.Astrm教授将F检验准则用于对时间序列模型的定阶。原理(模型阶数简约原则 parsimony principle)：设yt(1tn)是零均值平稳序列，用模型AR模型拟合检验统计量：结论若FF，则拒绝原假设，认为AR(p)合适；若FF，则拒绝原假设，模型阶数仍有上升的可能；若F1时，ARMA(1,1)预测值也是由如下差分方程决定的。(3)向前L步预测公式(L2).三、三、预测误差预测误差.由于预测只能建立在到t时刻为止的可用信息的基础上，因此，根据最小均方误预测的第二个准则，以及平稳可逆序列可以表示成传递函数形式的论断，可以将预测值 表示成能够估计的项t，t-1，的加权和的形式：.由上得以t为原点，向前L步的预测误差为：由于t是白噪声，故有：.误差方差为：误差方差为：.注：预测误差的估算是注：预测误差的估算是1，p算和算和1，q估计都为估计都为正确的假设，实际中参数通过估计得到的，且估计量是随机正确的假设，实际中参数通过估计得到的，且估计量是随机变量，有均值和方差，因而实际误差大于理论估计误差。变量，有均值和方差，因而实际误差大于理论估计误差。.五、预测误差的置信区间五、预测误差的置信区间对于正态过程，预测误差的分布为所以：对所以：对yt+l预测的预测的95%的置信区间为：的置信区间为：因此：.设为一AR(2)序列，其中求的自协方差函数。例1解答：解答：Yule-WalkerYule-Walker方程为：方程为：且：联合上面三个方程，解出：.例例 2考虑如下AR(2)序列：若已知观测值（1）试预报（2）给出（1）预报的置信度为95%的预报区间。,.解答：（1）（2）预报的置信度为95%的预报区间分别为：.7.3 单位根检验和协整检验单位根检验和协整检验 一、平稳性的检验一、平稳性的检验引言：引言：前面我们讨论的是平稳时间序列的建模和预测方法，即所讨论的时间序列都是宽平稳的。一个宽平稳的时间序列的均值和方差均值和方差都是常数，并且它的协方差有时间上协方差有时间上的不变性。但是许多经济领域产生的时间序列都是非平稳的。呈现出明显得趋势性和周期性，序列不平稳，导致预测无效，产生谬误回归谬误回归等问题。.1、通过时间序列的趋势图来判断这种方法通过观察时间序列的趋势图来判断时间序列是否存在趋势性或周期性。优点：简便、直观。对于那些明显为非平稳的时间序列，可以采用这种方法。缺点：对于一般的时间序列是否平稳，不易用这种方法判断出来。.2、通过自相关函数(ACF)判断平稳时间序列的自相关函数(ACF)要么是截尾的，要么是拖尾的。因此我们可以根据这个特性来判断时间序列是否为平稳序列。若时间序列具有上升或下降的趋势若时间序列具有上升或下降的趋势，那么对于所有短时滞来说，自相关系数大且为正，而且随着时滞k的增加而缓慢地下降。.若序列无趋势若序列无趋势，但是具有季节性但是具有季节性，那末对于按月采集的数据，时滞12，24，36的自相关系数达到最大(如果数据是按季度采集，则最大自相关系数出现在4，8，12，)，并且随着时滞的增加变得较小。若序列是有趋势的，且具有季节性若序列是有趋势的，且具有季节性，其自相关函数特性类似于有趋势序列，但它们是摆动的，对于按月数据，在时滞12，24，36，等处具有峰态；如果时间序列数据是按季节的，则峰出现在时滞4，8，12，等处。.3、随机游走的单位根检验(Unit root test)随机游走是一种非平稳过程，其实随机游走一种特殊的齐次非平稳过程。检验序列是否为随机游走，通常利用David Dickey和Wayne Fuller的单位根检验。单位根的含义和检验原理如下：.（1）单位根的含义单位根的含义.（2）单位根的单位根的检验检验用Eviews进行单位根检验时给出了上述选项。如果DF检验统计量比给定显著水平临界值大，不能拒绝原假设，认为序列存在单位根，是非平稳的。.ADF检验检验在在DF检验中，常常因为序列存在高阶滞后相关，使检验中，常常因为序列存在高阶滞后相关，使得随机扰动不符合白噪声假设，得随机扰动不符合白噪声假设，ADF检验修正了检验修正了DF检检验中的自相关问题。验中的自相关问题。此外还有：此外还有：PP检验：检验具有一般形式的单位根过程检验：检验具有一般形式的单位根过程DFGLS检验：检验：DF及及ADF检验对含有时间趋势的退势平检验对含有时间趋势的退势平稳时间序列的检验失效稳时间序列的检验失效古典的回归方法：只能对平稳的时间序列进行回归分古典的回归方法：只能对平稳的时间序列进行回归分析，或将非平稳的序列转化为平稳序列再做回归析，或将非平稳的序列转化为平稳序列再做回归.非平稳时间序列分析逐期差分后平稳，建立求和自回归移动平均模型,记为ARIMA(p,d,q).非平稳时间序列分析非平稳时间序列分析逐期差分+季节差分后平稳，建立乘积季节模型，即随机季节模型与ARIMA模型结合，记为ARIMA(p,d,q)(P,D,Q).二、协整检验二、协整检验q如果两个或多个非平稳的时间序列，其某个线性组合后的序列呈平稳性，这样的时间序列间就被称为有协整关系存在。.协整检验一般用于误差修正模型，用于对回归模型协整检验一般用于误差修正模型，用于对回归模型进行修正，针对一元和多元回归，常用的协整检验进行修正，针对一元和多元回归，常用的协整检验两种两种qEngle-Granger两步协整检验法.q误差修正模型误差修正模型.多元回归，常用的协整检验为多元回归，常用的协整检验为Johansen协整检验。协整检验。.