应用时间序列分析实验手册2目录目录.........................................................................................................2第二章时间序列的预处理......................................................................3一、平稳性检验..................................................................................3二、纯随机性检验..............................................................................9第三章平稳时间序列建模实验教程....................................................10一、模型识别....................................................................................10二、模型参数估计(如何判断拟合的模型以及结果写法)........14三、模型的显著性检验....................................................................17四、模型优化....................................................................................18第四章非平稳时间序列的确定性分析................................................19一、趋势分析....................................................................................19二、季节效应分析............................................................................34三、综合分析....................................................................................38第五章非平稳序列的随机分析............................................................44一、差分法提取确定性信息............................................................44二、ARIMA模型.............................................................................58三、季节模型....................................................................................623第二章时间序列的预处理一、平稳性检验时序图检验和自相关图检验(一)时序图检验根据平稳时间序列均值、方差为常数的性质,平稳序列的时序图应该显示出该序列始终在一个常数值附近随机波动,而且波动的范围有界、无明显趋势及周期特征例2.1检验1964年——1999年中国纱年产量序列的平稳性1.在Eviews软件中打开案例数据图1:打开外来数据图2:打开数据文件夹中案例数据文件夹中数据4文件中序列的名称可以在打开的时候输入,或者在打开的数据中输入图3:打开过程中给序列命名图4:打开数据52.绘制时序图可以如下图所示选择序列然后点Quick选择Scatter或者XYline;绘制好后可以双击图片对其进行修饰,如颜色、线条、点等图1:绘制散点图图2:年份和产出的散点图6010020030040050060019601970198019902000YEAROUTPUT图3:年份和产出的散点图(二)自相关图检验例2.3导入数据,方式同上;在Quick菜单下选择自相关图,对Qiwen原列进行分析;可以看出自相关系数始终在零周围波动,判定该序列为平稳时间序列。图1:序列的相关分析7图2:输入序列名称图2:选择相关分析的对象图3:序列的相关分析结果:1.可以看出自相关系数始终在零周围波动,判定该序列为平稳时间序列2.看Q统计量的P值:该统计量的原假设为X的1期,2期……k期的自相关系数均等于0,备择假设为自相关系数中至少有一个不等于0,因此如图知,该P值都5%的显著性水平,所以接受原假设,即序列是纯随机序列,即白噪声序列(因为序列值之间彼此之间没有任何关联,所以说过去的行为对将来的发展没有丝毫影响,因此为纯随机序列,即白噪声序列.)有的题目平稳性描述可以模仿书本33页最后一段.(三)平稳性检验还可以用:8单位根检验:ADF,PP检验等;非参数检验:游程检验图1:序列的单位根检验图2:单位根检验的方法选择表示不包含截距项9图3:ADF检验的结果:如图,单位根统计量ADF=-0.016384都大于EVIEWS给出的显著性水平1%-10%的ADF临界值,所以接受原假设,该序列是非平稳的。二、纯随机性检验计算Q统计量,根据其取值判定是否为纯随机序列。例2.3的自相关图中有Q统计量,其P值在K=6、12的时候均比较大,不能拒绝原假设,认为该序列是白噪声序列。另外,小样本情况下,LB统计量检验纯随机性更准确。10第三章平稳时间序列建模实验教程一、模型识别1.打开数据图1:打开数据2.绘制趋势图并大致判断序列的特征图2:绘制序列散点图11图3:输入散点图的两个变量图4:序列的散点图123.绘制自相关和偏自相关图图1:在数据窗口下选择相关分析图2:选择变量图3:选择对象13图4:序列相关图4.根据自相关图和偏自相关图的性质确定模型类型和阶数如果样本(偏)自相关系数在最初的d阶明显大于两倍标准差范围,而后几乎95%的自相关系数都落在2倍标准差的范围以内,而且通常由非零自相关系数衰减为小值波动的过程非常突然。这时,通常视为(偏)自相关系数截尾。截尾阶数为d。本例:自相关图显示延迟3阶之后,自相关系数全部衰减到2倍标准差范围内波动,这表明序列明显地短期相关。但序列由显著非零的相关系数衰减为小值波动的过程相当连续,相当缓慢,该自相关系数可视为不截尾偏自相关图显示除了延迟1阶的偏自相关系数显著大于2倍标准差之外,其它的偏自相关系数都在2倍标准差范围内作小值随机波动,而且由非零相关系数衰减为小值波动的过程非常突然,所以该偏自相关系数可视为一阶截尾所以可以考虑拟合模型为AR(1)自相关系数偏相关系数模型定阶拖尾P阶截尾AR(p)模型Q阶截尾拖尾MA(q)模型拖尾拖尾ARMA(P,Q)模型具体判别什么模型看书58到62的图例。:就是常数项)。表示的是求出来的系数(其中模型中的模型:)(模型:模型:)1(MA)1(arB*)P(ARB*)2(ARB*)1(AR1B*)q(MAB*)2(MAB*)1(MA1ARMAB*)q(MAB*)2(MAB*)1(MA1MAB*)P(ARB*)2(ARB*)1(AR11ARtP2q2tXtq2tXtP2tX14二、模型参数估计根据相关图模型确定为AR(1),建立模型估计参数在ESTIMATE中按顺序输入变量cxccx(-1)或者cxcar(1)选择LS参数估计方法,查看输出结果,看参数显著性,该例中两个参数都显著。细心的同学可能发现两个模型的C取值不同,这是因为前一个模型的C为截距项;后者的C则为序列期望值,两个常数的含义不同。图1:建立模型15图2:输入模型中变量,选择参数估计方法图3:参数估计结果图4:建立模型16图5:输入模型中变量,选择参数估计方法图6:参数估计结果tB703332.01132034.81txAR模型:17三、模型的显著性检验检验内容:整个模型对信息的提取是否充分;参数的显著性检验,模型结构是否最简。图1:模型残差18图2:残差的平稳性和纯随机性检验对残差序列进行白噪声检验,可以看出ACF和PACF都没有显著异于零,Q统计量的P值都远远大于0.05,因此可以认为残差序列为白噪声序列,模型信息提取比较充分。常数和滞后一阶参数的P值都很小,参数显著;因此整个模型比较精简,模型较优。四、模型优化当一个拟合模型通过了检验,说明在一定的置信水平下,该模型能有效地拟合观察值序列的波动,但这种有效模型并不是唯一的。当几个模型都是模型有效参数显著的,此时需要选择一个更好的模型,即进行优化。优化的目的,选择相对最优模型。优化准则:最小信息量准则(AnInformationCriterion)指导思想似然函数值越大越好未知参数的个数越少越好AIC准则的缺陷在样本容量趋于无穷大时,由AIC准则选择的模型不收敛于真实模型,它通常比真实模型所含的未知参数个数要多但是本例中滞后二阶的参数不显著,不符合精简原则,不必进行深入判断。)(2)ˆln(2未知参数个数nAIC))(ln()ˆln(2未知参数nnSBC19第四章非平稳时间序列的确定性分析第三章介绍了平稳时间序列的分析方法,但是自然界中绝大多数序列都是非平稳的,因而对非平稳时间序列的分析跟普遍跟重要,人们创造的分析方法也更多。这些方法分为确定性时序分析和随机时序分析两大类,本章主要介绍确定性时序分析方法。一个序列在任意时刻的值能够被精确确定(或被预测),则该序列为确定性序列,如正弦序列、周期脉冲序列等。而某序列在某时刻的取值是随机的,不能给以精确预测,只知道取某一数值的概率,如白噪声序列等。Cramer分解定理说明每个序列都可以分成一个确定序列加一个随机序列,平稳序列的两个构成序列均平稳,非平稳时间序列则至少有一部分不平稳。本章先分析确定性序列不平稳的非平稳时间时间序列的分析方法。确定性序列不平稳通常显示出非常明显的规律性,如显著趋势或者固定变化周期,这种规律性信息比较容易提取,因而传统时间序列分析的重点在确定性信息的提取上。常用的确定性分析方法为因素分解。分析目的为:①克服其他因素的影响,单纯测度某一个确定性因素的影响;②推断出各种因素彼此之间作用关系及它们对序列的综合影响。一、趋势分析绘制序列的线图,观测序列的特征,如果有明显的长期趋势,我们就要测度其长期趋势,测度方法有:趋势拟合法、平滑法。(一)趋势拟合法1.线性趋势拟合例1:以澳大利亚政府1981-1990年每季度消费支出数据为例进行分析。图1:导入数据20图2:绘制线图,序列有明显的上升趋势长期趋势具备线性上升的趋势,所以进行序列对时间的线性回归分析。图3:序列支出(zc)对时间(t)进行线性回归分析21图4:回归参数估计和回归效果评价可以看出回归参数显著,模型显著,回归效果良好,序列具有明显线性趋势。图5:运用模型进行预测22图6:预测效果(偏差率、方差率等)图7:绘制原序列和预测序列的线图23图8:原序列和预测序列的线图图9:残差序列的曲线图可以看出残差序列具有平稳时间序列的特征,我们可以进一步检验剔除了长期趋势后的残差序列的平稳性,第三章知识这里不在叙述。242.曲线趋势拟合例2:对上海证券交易所1991.1-2001.10每月月末上正指数序列进行拟合。图1:导入数据图2:绘制曲线图可以看出序列不是线性上升,而是曲线上升,尝试用二次模型拟合序列的发展。25图3:模型参数估计和回归效果评价因为该模型中T的系数不显著,我们去掉该项再进行回归分析。图4:新模型参数估计和回归效果评价26图5:新模型的预测效果分析图6:原序列和预测序列值27图