时间序列实验习题汇总

时间序列实验习题汇总实验二时间序列纯随机性检验和平稳性检验3、Eviews操作方法打开要检验的序列，单击View\Correlogram二、序列平稳性检验主要有时序图检验，自相关图检验，单位根检验。（一）时序图检验根据平稳时间序列均值、方差为常数的性质，平稳序列的时序图应该显示出该序列始终在一个常数值附近随机波动，而且波动的范围有界的特点，如果该序列的时序图有明显的趋势性和周期性，那么该序列通常不是平稳序列。如下图所示可知，x，y两个序列均为非平稳序列。01002003004005006001965197019751980198519901995X5006007008009001,0001962196419661968197019721974Y（二）自相关图检验平稳序列的自相关函数很快衰减为0；非平稳序列的自相关函数衰减为0的速度通常比较慢，这是我们通过利用自相关图进行平稳性判断的标准。下面左图我们发现序列的自相关函数衰减到0的速度相当缓慢，在很长的一段时间里，自相关系数为正，然后一直为负，在自相关图上显示出明显的三角对称性，这是具有单调趋势的非平稳序列的一种典型的自相关图的形式。这与上面x的时序图显示的显著的单调递增性是一致的。下面右图中自相关图呈现出明显的正弦波动规律，这是具有周期变化规律的非平稳序列的典型特征。这与上面y的时序图显示的显著的单调递增性是一致的。（三）单位根检验1、ADF检验法（AugmentedDickey—FullerTest）（1）ADF检验法是由迪基（Dickey）和福勒（Fuller）在1979年提出的，是DF方法的推广。DF假定{εt}是独立同分布序列，ADF假定随机扰动项{μt}是稳定过程。（2）原假设0H：有一个单位根。2、迪基---福勒(DF)检验法（1）DF检验法是由Dickey、Fuller在20世纪70、80年代的一系列文章中建立起来的。（2）原假设0H：有一个单位根3、Eviews操作方法打开要检验的序列，单击View\unitroottest【实验目的】认识Eviews输出的时间序列自相关图的内容及含义：自相关函数、偏自相关函数、95%置信限、Q-statistic。学会通过自相关图的Q统计量判断序列是否为纯随机序列（白噪声序列）。通过观察序列的趋势图及单位根检验结果判断序列是否为平稳序列。【实验内容】序列纯随机性检验和平稳性检验【实验步骤】案例1：对1949-1998年北京市最高气温做纯随机序列检验打开序列x，单击view，出现如下的对话框，单击correlogram，出现如下对话框，要对原序列水平做纯随机序列检验，所以在correlogramof复选框下面选择level，而不是1stdifference（原序列一阶差分），2nddifferenc（原序列二阶差分），滞后阶数为24（一般是默认值）。点击ok。出现下面结果。从自相关图、偏自相关图（均在2倍标准差以内）可以看出，能以95%的水平保证自相关函数、偏自相关函数为0，所以，该序列为纯随机序列；另外，对于显著性水平0.05，Q-Stat的P值均大于0.05，所以，该序列为纯随机序列。北京市最高气温的变动属于纯随机波动，我们很难根据历史信息预测未来年份的最高气温。给序列的分析到此结束。实验三ARMA模型的识别、建立、检验【实验目的】熟悉对零均值平稳序列建立ARMA模型的前三个阶段：模型识别、模型参数估计、诊断检验。（1）根据时间序列自相关图对零均值平稳序列进行初步的模型识别。（2）运用Eviews软件估计ARMA模型参数。（3）对所建立的模型是否为适应性模型进行诊断检验。【实验内容】一、模型识别根据零均值平稳化后的序列的自相关函数和偏自相关函数表现出的特征（即拖尾或者截尾），对序列进行初步的模型识别（注：这种方法并不总是有效）。二、模型参数估计（最小二乘估计）Eviews建立ARMA模型的命令用到AR、MA等参数项。例如：对一个零均值的平稳序列x建立ARMA(2,1)模型，方法一：命令操作方式，即在命令编辑窗口输入：lsxar(1)ar(2)ma(1)；方法二：菜单操作方式：Quick---Estimateequation，出现下面的对话框，输入：xar(1)ar(2)ma(1)，点击确定。以上述操作方式建模时，Eviews自动采用非线性最小二乘法估计模型参数。三、模型的诊断检验：（一）模型的适应性检验检验残差序列t是否是一白噪声序列，关键是t的独立性检验。主要的检验方法为德宾—沃森统计量（Durbin-Watsonstat）检验、相关图和Q统计量检验（已介绍）、Breush-GodfreyLM（LagrangMultiplierTest，拉格朗日乘子检验）检验。下面逐一做介绍。（1）德宾—沃森统计量（Durbin-Watsonstat）检验检验规则：在D—W小于等于2时，存在正的一阶序列自相关，特别是接近于0时，存在严重的序列自相关;在D—W大于2时,存在负的一阶序列自相关。由上面输出结果可以显示，残差序列存在一阶序列正相关。对模型进行修正之后的输出结果为我们可以认为残差序列不存在一阶序列相关性。但是并不能下结论：残差序列不存在序列相关性。因为德宾—沃森统计量（D—W）检验仅仅能检验一阶序列相关，而“相关图和Q统计量检验”与“Breush-GodfreyLM（LagrangMultiplierTest）检验”则可以检验一般的序列相关形式，即可以检验残差序列是否存在高阶序列相关性。（2）相关图和Q统计量检验主要根据模型残差是否为白噪声来判断，若残差是白噪声，则可认为此模型是序列的适应性模型，否则，不是。Eviews操作：在模型窗口，View----Residualtests----Correlogram—Qstatistics根据输出的残差的Q统计量判断残差是否为白噪声序列。(ap，不拒绝原假设，该系列为白噪声序列。)（3）Breush-GodfreyLM（拉格朗日乘子）检验Breush-GodfreyLM（LagrangMultiplierTest，拉格朗日乘子检验）检验是由Breush-Godfrey提出的，主要用来检验回归方程中残差序列是否存在高阶序列相关性。该检验的原假设H0：残差序列不存在序列相关性。Eviews操作：在模型窗口，View----Residualtests----SerialCorrelationLMTest然后，出现下面的关于滞后阶数的对话框，默认阶数为2，我们可以改为3，4，或者更高阶。这里我们选择检验的阶数为4阶。点击ok，出现下面的结果检验统计量为Obs*R-squred，其值为3.402357，而p值为0.4929。如果给定的显著性水平为0.05，（ap）则接受原假设，即认为残差序列不存在序列相关性，残差检验通过。（二）模型中各项的取舍若建立的模型为适应性模型，还要看输出项中各变量是否显著(通过输出结果中的t统计量值及相应的P值)，对不显著的项，要剔除，然后重新建模。（三）模型的选择（定阶）对于同一个序列来说，可能有多个适应性模型，要从这多个适应性模型中选择，通常根据多个模型输出项中的赤池信息准则(AIC，Akaikeinfocriterion)和施瓦茨准则(SBC，SchwartzBayescriterion)（在Eviews输出项中为SC）进行比较，一般认为这两个统计量值越小的模型越好。（四）模型平稳性和可逆性的判断判断模型是适应性模型后，还应判断模型是否平稳和可逆，判断方法如下。模型输出结果最下方输出的两项，ARinvertedroot(如果有的话)和MAinvertedroot（如果有的话），其含义分别为：invertedARroot：为模型自回归AR部分所对应的差分方程的特征方程的特征根。若特征根的绝对值都小于1，则说明模型是平稳的；若其中有大于或等于1的，说明模型非平稳；若有等于1或很接近于1的，说明原序列为单位根过程，需要先对序列进行差分平稳化变换(有几个单位根，作几阶差分变换)，然后建模。invertedMAroot：为模型移动平均MA部分所对应的差分方程的特征方程的特征根。若特征根绝对值都小于1，则说明模型是可逆的；若有大于或等于1的，说明模型不可逆；若有等于1或很接近于1的，则很有可能在数据处理过程中，对原序列过度差分了，这时需要减少对序列差分的阶数，再重新建模。实验四ARMA模型的预测【实验目的】：（1）进一步熟悉ARMA模型建模过程。（2）利用ARMA模型进行预测。【实验内容】平稳时间序列模型预测、非平稳时间序列模型的预测【实验步骤】平稳时间序列模型预测操作文件：ar（2）.wf1（1）打开ar（2）.wf1（2）对序列x建立AR(2)模型操作命令：lsxcar(1)ar(2)（3）进行向前多步预测。Equation窗口点击estimate，样本范围改为：sampl1160然后在Equation窗口，选Forecast菜单，在出现的对话框中，选Dynamic，预测范围为forecastsample：160162，并将预测结果保存在xf序列中，单击OK。观察输出结果xf。说明：Dynamic为动态预测。注：S.E用于存放预测的估计标准误差，便于计算置信区间。非平稳时间序列预测（操作文件：北京市民用车拥有量.wf1）操作步骤：(1)打开北京市民用车拥有量.wf1，(2)对序列dlog(x)建立ar（1）模型操作命令：lsdlog(x)car(1)(3)进行向前多步预测Equation窗口点击estimate，样本范围改为：sampl19501997然后在Equation窗口，选Forecast菜单，在出现的对话框中，选Dynamic，预测范围为forecastsample：19971999，并将预测结果保存在xf序列中，单击OK。观察输出结果xf。实验五复习ARMA建模过程【实验目的】复习利用Eviews对时间序列建立ARMA模型的过程【实验内容】ARMA模型建模前的准备：判断序列是否平稳a.通过序列自相关图、趋势图等进行判断b.若序列不平稳：均值非平稳序列通过差分变换转换为平稳方差非平稳序列通过对数变换等转化为平稳序列c.模型平稳化以后，将序列零均值化（1）模型识别主要通过序列的自相关函数、偏自相关函数表现的特征，进行初步的模型识别（2）模型参数估计a.在Eviews中估计ARMA模型的方法b.估计模型以后要能写出模型的形式(差分方程形式和用B算子表示的形式)（3）模型的诊断检验a.根据模型残差是不是白噪声来判断模型是否为适应性模型b.能根据输出结果判断模型是否平稳，是否可逆c.若有多个序列是模型的适应性模型，会用合适的方法从这些模型中进行选择，如比较模型的残差方差，AIC，SC等。（4）模型应用a.掌握预测的操作方法