《时间序列分析 (2)2幻灯片.ppt》由会员分享,可在线阅读,更多相关《时间序列分析 (2)2幻灯片.ppt(55页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、时间序列分析第1页,共55页,编辑于2022年,星期六第一节第一节 时间序列分析的一般问题时间序列分析的一般问题时间序列的含义时间序列是指被观察到的以时间为序排列的数据序列。时间序列可以以表格的形式或图形的形式表现。例:第2页,共55页,编辑于2022年,星期六上海180指数某时间段的变化第3页,共55页,编辑于2022年,星期六国际航运乘客资料(单位:千人)第4页,共55页,编辑于2022年,星期六19461970美国各季生产者耐用品支出(单位:十亿美元)第5页,共55页,编辑于2022年,星期六1952年1994年我国社会消费品零售总额(单位:亿元)第6页,共55页,编辑于2022年,星期
2、六时间序列分析的目的:(1)预测序列未来的发展方向。(2)分析序列的基本趋势、季节和随机项的构成。(3)分析待定的数据集合、模拟理论模型,尤其是建立数学模型,进而进行模型的结构分析和实证研究。第7页,共55页,编辑于2022年,星期六4.按序列的分布规律来分:高斯型时间序列和非高斯型时间序列。时间序列的主要分类:1.按所研究的对象的多少分:一元时间序列和多元时间序列。2.按时间的连续性分:离散时间序列和连续时间序列。3.按序列的统计特性分:平稳时间序列和非平稳时间序列。第8页,共55页,编辑于2022年,星期六第二节第二节 时间序列的建立时间序列的建立我们把获取时间序列以及对其进行检查、整理和
3、预处理等工作,称为时间序列的建立。第9页,共55页,编辑于2022年,星期六时间序列数据的采集相应于时间的连续性,系统在不同的时刻上的响应常常是时间t的连续函数。为了数字计算处理上的方便,往往只按照一定的时间间隔对所研究系统的响应进行记录和观察,我们称之为采样。相应地把记录和观察时间间隔称为采样间隔。通常采样采用等间隔采样。第10页,共55页,编辑于2022年,星期六离群点(Outlier)离群点(Outlier)是指一个时间序列中,远离序列一般水平的极端大值和极端小值。对时间序列离群点分析的方法,有时也被称作稳健估计(RobustEstimation),该方法最早由Box和Anderson于
4、1955年提出。第11页,共55页,编辑于2022年,星期六1.离群点(Outlier)产生的原因:(1)采样误差;(2)系统各种偶然非正常因素影响。2.离群点的数理描述:(1)它们是既定分布中的极端点(extemepoint),它们虽与数据主体来自同一分布,但本身应以极小的概率出现;(2)这种点与数据集的主体并非采自同一分布,而是在采集数据过程中受到其他分布的“污染”,致使现有数据集掺入不应有的“杂质”。第12页,共55页,编辑于2022年,星期六3.离群点(Outlier)的类型:(1)加性离群点(AdditiveOutlier),造成这种离群点的干扰,只影响该干扰发生的那一个时刻T上的序
5、列值,而不影响该时刻以后的序列值。(2)更新离群点(InnovationalOutlier),造成离群点的干扰不仅作用于XT,而且影响T时刻以后序列的所有观察值。第13页,共55页,编辑于2022年,星期六(3)水 平 移 位 离 群 点(Level ShiftOutlier),造成这种离群点的干扰是在某一时刻T,系统的结构发生了变化,并持续影响T时刻以后的所有行为,在数列上往往表现出T时刻前后的序列均值发生水平位移。(4)暂时变更离群点(TemporaryChangeOutlier),造成这种离群点的干扰是在T时刻干扰发生时具有一定初始效应,以后随时间根据衰减因子的大小呈指数衰减。(请大家观
6、看四种离群点的演示试验)第14页,共55页,编辑于2022年,星期六(1)直接进行剔除;(2)对数据模型进行修正处理分析。4.离群点(Outlier)的处理方法:第15页,共55页,编辑于2022年,星期六缺损值(Missingvalue)的补足依据系统运动轨迹或变化趋势,运用一定的方法对缺损值进行估计、推测,以补足缺损的数值。第16页,共55页,编辑于2022年,星期六第三节第三节 确定性时序分析方法概述确定性时序分析方法概述一个时间序列往往是以下几类变化形式的叠加或耦合:长期趋势变动(T)季节变动(S)循环变动(C)不规则变动(R)第17页,共55页,编辑于2022年,星期六常见的确定性时
7、间序列模型:()混合模型()加法模型()乘法模型第18页,共55页,编辑于2022年,星期六确定性时序分析方法:1.移动平均法 设观测序列为,正整数NT。一次移动平均值计算公式为:第19页,共55页,编辑于2022年,星期六二次移动平均值计算公式为:第20页,共55页,编辑于2022年,星期六移动平均模型的特点:(1)等加权;(2)与RandomWalk相比,对反转的反应滞后,平均滞后(N+1)/2期;(3)无现成统计理论用于预测区间的推导,要根据样本的数据计算标准差做经验估计。移动平均模型在数据处理中常用作预处理,用于消除周期波动(取N为周期长度)和减弱随机干扰的影响。第21页,共55页,编
8、辑于2022年,星期六2.指数平滑法,为加权系数,01,设观测序列p次指数平滑公式为:一次指数平滑公式为:二次指数平滑公式为:第22页,共55页,编辑于2022年,星期六一次指数平滑模型(SimpleExponentialSmoothing,SES)的特点:非等加权,距离越近权数越大;平滑指数为连续变量,可以通过最小均方误(meansquarederror)计算出最最佳的平滑系数;对反转平均滞后1/;预测区间比RandomWalk窄;ARIMA模型的特例,ARIMA(0,1,1),MA的系数为1-。第23页,共55页,编辑于2022年,星期六二次指数平滑模型(Linear(double)Exp
9、onentialSmoothing,LES)常用于含线性趋势数据;三次指数平滑模型(Quadratic(triple)SmoothingModle)常用于含曲线趋势的数据。第24页,共55页,编辑于2022年,星期六3.时间回归法(1)线性方程:(2)二次曲线:(3)指数曲线:(4)修正指数曲线:第25页,共55页,编辑于2022年,星期六(5)Gompertz(龚帕兹)曲线:(6)Logistic(逻辑斯谛)曲线:(7)振动曲线:f(t)为多项式(请同学们看一个实例)第26页,共55页,编辑于2022年,星期六4.季节周期预测法(1)乘法型季节模型(2)加法型季节模型第27页,共55页,编辑
10、于2022年,星期六第四节 随机时间序列分析的几个基本概念随机过程(StochasticProcess)随机过程是一簇随机变量Xt,tT,其中T表示时间t的变动范围,对每个固定的时刻t而言,Xt是一普通的随机变量,这些随机变量的全体就构成一个随机过程。当t=0,1,2,时,即时刻t只取整数时,随机过程Xt,tT可写成如下形式,Xt,t=0,1,2,。此类随机过程Xt是离散时间t的随机函数,称它为随机序列或时间序列。对于一个连续时间的随机过程的等间隔采样序列,即Xt,t=0,1,2,就是一个离散随机序列。第28页,共55页,编辑于2022年,星期六时间序列的概率分布和数值特征1.时间序列的概率分
11、布 一个时间序列的概率分布可以用它有限维分布簇来描述时间序列所有的一维分布是:,F-1(),F0(),F1(),所有二维分布是:Fij(,),i,j=0,1,2,(ij)一个时间序列的所有有限维分布簇的全体,称为该序列的有限维分布簇。第29页,共55页,编辑于2022年,星期六2.时间序列的均值函数一个时间序列的均值函数是指:其中:EXt表示在t固定时对随机变量Xt的求均值,它只一维分布簇中的分布函数Ft()有关。第30页,共55页,编辑于2022年,星期六3.时间序列的协方差函数与自相关函数协方差函数:其中Ft,s(X,Y)为(Xt,Xs)的二维联合分布。自相关函数:第31页,共55页,编辑
12、于2022年,星期六时间序列的自协方差函数有以下性质:对称性:非负定性:为对称非负定矩阵。自相关函数同样也具有上述性质且有(t,t)=1。对任意正整数m和任意m个整数k1,k2,km,方阵第32页,共55页,编辑于2022年,星期六平稳随机过程严平稳:如果对于时间t的任意n个值和任意实数,随机过程的n维分布满足关系式:则称为严平稳过程。第33页,共55页,编辑于2022年,星期六宽平稳:若随机过程的均值(一阶矩)和协方差存在,且满足则称为宽平稳随机过程。第34页,共55页,编辑于2022年,星期六严平稳与宽平稳的关系:1.严平稳不等于宽平稳;2.宽平稳不等于严平稳;3.对于严平稳序列,如果其二
13、阶距存在,其必为宽平稳,反之则一般不成立;4.对于高斯序列,严平稳与宽平稳是等价的。第35页,共55页,编辑于2022年,星期六平稳时间序列自协方差函数和自相关函数设平稳时间序列的均值为零,即。自协方差函数:自相关函数:第36页,共55页,编辑于2022年,星期六平稳序列的自协方差函数有以下性质:对称性:非负定性:为非负定矩阵。自相关函数同样也具有上述性质。m阶自协方差阵有界性:第37页,共55页,编辑于2022年,星期六平稳序列的样本统计量样本均值:时间序列无法获得多重实现,常用时间均值代替总体均值。上式的估计是无偏的。第38页,共55页,编辑于2022年,星期六样本自协方差函数第一式是有偏
14、估计,第二式是无偏估计,但有效性不如第一式。第39页,共55页,编辑于2022年,星期六几类特殊的随机过程(序列):1.纯随机过程:随机过程如果是由一个不相关的随机变量的序列构成的,则称其为纯随机过程。2.白噪声序列(Whitenoise):如果时间序列满足以下性质:(1)(2)式中,当ts时,简称白噪声。称此序列为白噪声序列,第40页,共55页,编辑于2022年,星期六3.独立同分布序列:如果时间序列中的随机变量Xt,t=0,1,2,为相互独立的随机变量,而且具有相同的分布,称这样的时间序列为独立同分布序列。4.二阶矩过程:若随机过程对每个 二阶矩过程。的均值和方差存在,则称之为5.正态过程
15、:若正态分布,则称的有限维分布都是为正态随机过程。第41页,共55页,编辑于2022年,星期六第二章第二章 平稳时间序列模型平稳时间序列模型(单变量)单变量)选择单变量时间序列的原因一是对相关序列间的可能关系还缺乏可靠的先验信息。单变量时间序列建模是一个有用的简化手段,并可以进行有效的短期预测;例如例如,金融市场里面的股票走势、利率变化等等。金融市场里面的股票走势、利率变化等等。随机时间序列分析模型,就是要通过序列过去的变化特随机时间序列分析模型,就是要通过序列过去的变化特征来预测未来的变化趋势征来预测未来的变化趋势。第42页,共55页,编辑于2022年,星期六二是如果有关经济结构的理论已相当
16、成熟,则这个结构的某种表现将是对该结构中的每个内生变量都给出一个与单变量模型方程相同的方程形式。例如,例如,对于如下最简单的宏观经济模型:这里,Ct、It、Yt分别表示消费、投资与国民收入。Ct与与Yt作作为为内内生生变变量量,它它们们的的运运动动是是由由作作为为外外生生变变量量的投资的投资It的运动及随机扰动项的运动及随机扰动项 t的变化决定的。的变化决定的。tttCYCaaaa+=-12110第43页,共55页,编辑于2022年,星期六上述模型可作变形如下:两个方程等式右边除去第一项外的剩余部分可看成一个综合性的随机扰动项,其特征依赖于投资项It的行为。如如果果It是是一一个个白白噪噪声声
17、,则消费序列Ct就成为一个1阶阶自自回回归归过过程程AR(1),而收入序列Yt就成为一个(1,1)阶阶的自回归移动平均过程的自回归移动平均过程ARMA(1,1)。ttttICCaaaaaaaa1111011211111-+-+-+-=-tttttIIYYaaaaaaaaa11121101121111111-+-+-+-=-第44页,共55页,编辑于2022年,星期六第一节自回归模型一阶自回归模型AR(1)其中Xt零均值平稳序列,t为随机扰动。AR(1)的特点:Xt对Xt-1有线性相关关系t为独立正态同分布序列第45页,共55页,编辑于2022年,星期六AR(1)与普通一元线性回归一元线性回归一
18、阶自回归两个变量,Y为随机变量,X为确定性变量;一个变量,为随机变量;为白噪声序列,还可假定 为正态分布。第46页,共55页,编辑于2022年,星期六随机游动模型模型的特性:系统具有极强的一期记忆性;在t-1时刻,系统的一步超前预测可表示成无限独立随机变量和的形式:第47页,共55页,编辑于2022年,星期六一般自回归模型AR(n)其中:为白噪声,第48页,共55页,编辑于2022年,星期六第二节第二节 移动平均模型移动平均模型一阶移动平均模型MA(1)其中:为白噪声。MA(1)模型的基本假设为:(1)Xt仅与其前一时刻进入系统的扰动t-1有一定的依存关系;(2)t为白噪声。第49页,共55页
19、,编辑于2022年,星期六一般移动平均模型MA(m)其中:(1)Xt仅与有关,而与(j=m+1,m+2,)无关;(2)t为白噪声。第50页,共55页,编辑于2022年,星期六AR模型和MA模型的比较:AR模型可以写成无穷阶MA模型的形式,体现出AR模型对随机干扰的长效记忆性。对于不同系数的AR(n)序列,受到随机干扰后变化方向是不同的,对于不同系数的MA(m)序列,随机冲击并不改变其变化方向,也就是说随机冲击对于MA模型来说并不是一个“严重问题”。第51页,共55页,编辑于2022年,星期六第三节第三节 自回归移动平均自回归移动平均(ARMA)模型模型ARMA(n,m)模型模型其假设条件与其假设条件与AR模型相同模型相同第52页,共55页,编辑于2022年,星期六补充练习学会EVIEWS软件的基本操作,结合本章所学的的模型,练习用软件模拟相关模型。第53页,共55页,编辑于2022年,星期六第54页,共55页,编辑于2022年,星期六时间序列分析主讲教师李志强山东工商学院统计学专业第55页,共55页,编辑于2022年,星期六
限制150内