时间序列分析实验指导

PAGE/NUMPAGES时间序列分析实验指导实验一EVIEWS的基本操作【实验目的】熟悉Eviews的操作：菜单方式，命令方式；练习并掌握与时间序列分析相关的函数操作。【实验内容】一、EViews软件的常用菜单方式和命令方式；二、各种常用差分函数 表 关于同志近三年现实表现材料材料类招标技术评分表图表与交易pdf视力表打印pdf用图表说话 pdf 达式；三、时间序列的自相关和偏自相关图与函数；【实验步骤】一、EViews软件的常用菜单方式和命令方式；㈠创建工作文件⒈菜单方式启动EViews软件之后，进入EViews主窗口在主菜单上依次点击File/New/Workfile，即选择新建对象的类型为工作文件，将弹出一个对话框，由用户选择数据的时间频率（frequency）、起始期和终止期。选择时间频率为Annual（年度），再分别点击起始期栏（Startdate）和终止期栏（Enddate），输入相应的日期，然后点击OK按钮，将在EViews软件的主显示窗口显示相应的工作文件窗口。工作文件窗口是EViews的子窗口，工作文件一开始其中就包含了两个对象，一个是系数向量C（保存估计系数用），另一个是残差序列RESID（实际值与拟合值之差）。⒉命令方式在EViews软件的命令窗口中直接键入CREATE命令，也可以建立工作文件。命令格式为：CREATE时间频率类型起始期终止期则菜单方式过程可写为：CREATEA19851998㈡输入Y、X的数据⒈DATA命令方式在EViews软件的命令窗口键入DATA命令，命令格式为：DATA<序列名1><序列名2>…<序列名n>本例中可在命令窗口键入如下命令：DATAYX⒉鼠标图形界面方式在EViews软件主窗口或工作文件窗口点击Objects/NewObject，对象类型选择Series，并给定序列名，一次只能创建一个新序列。再从工作文件目录中选取并双击所创建的新序列就可以展示该对象，选择Edit＋/－，进入编辑状态，输入数据。㈢生成log（Y）、log（X）、X^2、1/X、时间变量T等序列在命令窗口中依次键入以下命令即可：GENRLOGY=LOG(Y)GENRLOGX=LOG(X)GENRX1=X^2GENRX2=1/XGENRT=@TREND(84)㈣选择若干变量构成数组，在数组中增加变量。在工作文件窗口中单击所要选择的变量，按住Ctrl键不放，继续用鼠标选择要展示的变量，选择完以后，单击鼠标右键，在弹出的快捷菜单中点击Open/asGroup，则会弹出数组窗口，其中变量从左至右按在工作文件窗口中选择变量的顺序来排列。在数组窗口点击Edit＋/－，进入全屏幕编辑状态，选择一个空列，点击标题栏，在编辑窗口输入变量名，再点击屏幕任意位置，即可增加一个新变量增加变量后，即可输入数据。点击要删除的变量列的标题栏，在编辑窗口输入新变量名，再点击屏幕任意位置，弹出RENAME对话框，点击YES按钮即可。㈤在工作文件窗口中删除、更名变量。⒈在工作文件窗口中选取所要删除或更名的变量并单击鼠标右键，在弹出的快捷菜单中选择Delete（删除）或Rename（更名）即可⒉在工作文件窗口中选取所要删除或更名的变量，点击工作文件窗口菜单栏中的Objects/Deleteselected…（Renameselected…），即可删除（更名）变量⒊在工作文件窗口中选取所要删除的变量，点击工作文件窗口菜单栏中的Delete按钮即可删除变量。三、图形分析与描述统计分析㈠利用PLOT命令绘制趋势图在命令窗口中键入：PLOTY也可以利用PLOT命令将多个变量的变化趋势描绘在同一张图中，例如键入以下命令，可以观察变量Y、X的变化趋势PLOTYX㈡利用SCAT命令绘制X、Y的散点图在命令窗口中键入：SCATXY则可以初步观察变量之间的相关程度与相关类型二、各种常用差分函数表达式表1-1:1949年1月---1960年12月数据1949年1950年1951年1952年1953年1954年1955年1956年1957年1958年1959年1960年111211514517119620424228431534036041721181261501801961882332773013183423913132141178193236235267317356362406419412913516318123522726931334834839646151211251721832292342703183553634204726135149178218243264315374422435472535714817019923026430236441346549154862281481701992422722933474054675055596069136158184209237259312355404404463508101191331621912112292743063473594074611110411414617218020323727130531036239012118140166194201229278306306337405432（一）利用D(x)命令系列对时间序列进行差分(x为表1-1中的数据)。1、在命令窗口中键入：genrdx=D(x)则生成的新序列为序列x的一阶差分序列2、在命令窗口中键入：genrdxn=D(x,n)则生成的新序列为序列x的n阶差分。3、在命令窗口中键入：genrdxs=D(x,0,s)则生成的新序列为序列x的对周期长度为s一阶季节差分。4、在命令窗口中键入：genrdxsn=D(x,n,s)则生成的新序列为对周期长度为s的时间序列x取一阶季节差分后的序列再取n阶差分。5、在命令窗口中键入：genrdlx=Dlog(x)则生成的新序列为x取自然对数后，再取一阶差分。6、在命令窗口中键入：genrdlxsn=Dlog(x,n,s)则生成的新序列为周期长度为s的时间序列x先取自然对数，再取一阶季节差分，然后再对序列取n阶差分。在EVIEWS中操作的图形分别为：三、时间序列的自相关和偏自相关图与函数；（一）观察时间序列的自相关图。命令方式：（1）在命令行输入命令：Identx(x为序列名称)；（2）然后在出现的对话框中输入滞后时期数。（可取默认数）菜单方式：（1）双击序列图标。菜单操作方式：View—>Correlogram，在出现的对话框中输入滞后数。（可取默认数）（二）练习：观察一些文件中的序列自相关函数Autocorrelation，偏自相关函数Partialautocorrelation的特征练习1：操作文件：Stpoor~1.wf1(美国S&P500工业股票价格指数1980年1月~1996年2月)步骤：（1）打开该文件。（2）观察序列stpoorr的趋势图，自相关图（自相关函数，偏自相关函数）的特征。（3）对序列取一阶差分，生成新序列dsp：genrdsp=d(stpoor)，并观察其趋势图，自相关图（同上，下略）的特征。（4）对该序列的自然对数取一阶差分，生成新的序列dlnsp：genrdlnsp=dlog(stpoor)，并观察其趋势图，自相关图。练习2：操作文件：usagnp.wf1(美国1947年第一季度~1970年第四季度GNP数据)步骤：（1）打开该文件。（2）观察序列usagdp的趋势图的特征，自相关图的特征。（3）对该序列取一阶差分，生新的序列dgdp：Genrdgdp=d(usagdp)。观察其趋势图，自相关图。（4）对该序列的自然对数取一阶差分，生成新的序列dlngdp：Genrdlngdp=dlog(gdp)。观察其趋势图，自相关图。（5）对序列一阶季节差分，生成新序列dsgdp=d(usagdp,0,4)观察其趋势图，自相关图的特征。（6）对该序列的自然对数取一阶季节差分，生成新的序列：dslngdp=dlog(usagdp,0,4)，观察其趋势图、自相关图。实验二传统时间序列分析【实验目的】认识Eviews输出的时间序列自相关图的内容及含义：自相关函数、偏自相关函数、95%置信限、Q-statistic。学会通过自相关图的Q统计量判断序列是否为白噪声。通过观察序列的趋势图及自相关图判断序列是否为平稳序列。【实验内容】一、本次练习主要操作文件为ar1.wf1，ar2.wf1，ma1.wf1，ma2.wf1，arma11.wf1，arma21.wf1，各文件中包含的序列都是模拟生成的零均值平稳序列。二、总结各种过程自相关函数，偏自相关函数的特征。三、观察其他文件中的序列，看其是否平稳，若不平稳，试通过适当的差分变换、方差平稳化变换(取对数，平方根等)使其转化为平稳序列，然后观察序列的自相关函数，偏自相关函数的特征，并与自已总结的各种过程的特征对照。【实验步骤】练习1.操作文件：ar1.wf1说明：该文件中含有三个序列：at为模拟生成的正态白噪声序列；x、y均是模拟生成的ar(1)过程，其参数各不相同。文件中有两个模型：EQX、EQY分别是对x、y的估计结果。操作内容：（1）观察序列at的自相关图，看其是否为白噪声序列，为什么？（2）观察序列x的自相关图：样本自相关函数（SACF）呈指数衰减，样本偏自相关函数（SPACF）滞后一阶截尾。（3）观察序列y的自相关图：样本自相关函数呈正负交替的指数衰减，样本偏自相关函数滞后一阶截尾。（4）分别打开EQX、EQY，试写出对x、y的估计结果。练习2：操作文件：ar2.wf1说明：该文件中含有四个序列：at为模拟生成的白噪声序列；x，y，z均为模拟生成的AR(2)过程，且其参数各不相同。文件中有三个模型：分别是对x、y、z的估计结果。操作内容：（1）分别观察序列x，y，z的自相关图，看其样本自相关函数，样本偏自相关函数各有什么特征。（提示：其样本自相关函数分别呈混合指数衰减、正负交替的混合指数衰减、阻尼正弦波衰减；样本偏自相关函数均滞后二阶截尾）。（2）分别打开EQX、EQY、EQZ，写出对x、y、z的估计结果。练习3：操作方件：ma1.wf1说明：文件中的序列x、y分别为模拟生成的ma(1)过程，其参数各不相同。文件中的模型EQX、EQY为对x、y的估计结果。操作内容：（1）分别观察序列x，y的自相关图，看其样本自相关图，偏自相关图各有什么特征。（提示：其样本自相关函数均呈滞后一阶截尾，样本偏自相关函数分别呈指数衰减、正负交替的指数衰减）。（2）分别打开EQX、EQY、写出对x、y的估计结果。练习4：操作文件：ma2.wf2说明：文件中的序列分别为模拟生成的MA(2)过程，其参数各不相同。操作内容：（1）分别观察序列x，y的自相关图，看其样本自相关图，偏自相关图各有什么特征。（提示：各序列的样本自相关函数均滞后二阶截尾，样本偏自相关函数分别呈混合指数衰减、正负交替的混合指数衰减，阻尼正弦波衰减）。（2）分别打开EQX、EQY、写出对x、y的估计结果。练习5：操作文件：ARMA11.wf1说明：文件中的序列x，y，z分别为模拟生成的不同参数的ARMA(1,1)过程，EQX、EQY、EQZ分别为对各序列估计的结果。操作内容：（1）分别观察序列x，y的自相关图，看其样本自相关图，偏自相关图各有什么特征。（提示：各序列的自相关函数，偏自相关函数都呈指数衰减）。（2）写出各模型的估计结果。练习6：操作文件：ARMA21.wf1操作内容：（1）分别观察序列x，y的自相关图，看其样本自相关图，偏自相关图各有什么特征。（提示：各序列的自相关函数，偏自相关函数都呈指数衰减）。（2）写出各模型的估计结果。实验三AR模型的应用【实验目的】观察具有实际背景的经济数据，判断其是否平稳、是否含有季节性，均值是否为零。能运用合适的方法如差分、季节差分、取对数、平方根等，使序列变为平稳序列；平稳序列减去其均值，使其零均值化。【实验内容】一、判断序列的平稳性和可逆性，给出相应判断依据，并写出模型形式。二、找出自己感兴趣的数据，判断数据是否平稳，是否具有季节性，均值是否为零等。【实验步骤】练习一操作文件：ar1.wf1，ar2.wf1，ma1.wf1，ma2.wf1，arma11.wf1，arma21.wf1操作内容：一、（1）打开文件ar1.wf1，(2)依据EQX，写出关于序列x的模型形式：Xt=0.68Xt-1+at（3）写出用B算子表示的模型形式：(1－0.68B)Xt=at(4)判断模型是否平稳？说明原因。(5)写出该模型的传递形式。二、（1）打开文件ar2.wf1（2）依据EQX写出序列x的模型形式为：Xt=0.49Xt-1+0.25Xt-2+at(3)写出用B算子表示的形式：(4)判断模型是否平稳?说明原因。(5)试推导模型的传递形式。并写出其前5个格林函数。三、（1）打开文件ma1.wf1(2)依据EQX写出序列X的模型形式：Xt=at－0.82at-1(3)写出用B算子表示的形式：Xt=(1－0.82B)at(4)判断模型是否可逆？说明原因。（5）写出该模型的逆转形式。四、（1）打开文件arma1.wf1(2)依据EQX写出序列X的模型形式：Xt=0.92Xt-1+at－0.57at-1(3)写出用B算子表示的形式：(1－0.92B)Xt=(1－0.57B)at(4)判断模型是否平稳？是否平稳？说明原因。（5）试推该模型的传递函数形式。五、打开ma2.wf1，写出各序列模型形式及用B算子表示的形式，判断序列是否可逆，试推导其逆转形式。打开ARMA21.wf1，写出各序列模型形式及用B算子表示的形式，判断序列是否平稳，是否可逆，试推导其传递函数形式，逆转形式。练习二操作文件：zl1.wf1~zl20.wf1，gdp.wf1，gdpindex.wf1，stpoor.wf1，usagnp.wf1等。文件说明：（1）zl1wf1~zl20.wf1各文件是教材后附录III所列资料，各数据背景参见附录。（2）gdp.wf1为我国1978~2001各年GDP数据。Gdpindex.wf1为我国1953~2001各年GDP指数，即各年GDP发展速度数据。（3）stpoor.wf1，usagnp.wf1文件说明见第一次上机实习内容说明。判断是否平稳、是否具有季节性的方法：（1）通过序列的趋势图粗略的判断。（2）通过序列的自相关图判断。若序列自相关函数衰减缓慢，滞后较长时期仍不为零，则可初步断定序列非平稳。若序列的自相关函数周期性的显著不为零（如月度数据的滞后12期，24期，36期等自相关函数显著不为零；季度数据的滞后4，8，12，16各期自相关函数显著不为零）则可判断序列含有季节性。使序列平稳化的方法：（1）若数据方差非平稳，应先通过对数变换、平方根变换等方法，使序列方差平稳。（2）先通过差分消除序列的长期趋势（如果有的话）。（3）再通过季节差分消除序列的季节性（如果有的话）。差分函数的使用可见前两次上机实习内容。使平稳序列零均值化的方法：HYPERLINK"mailto:在Eviews中可通过函数@mean()求序列的均值"在Eviews中可通过函数@mean()求序列的均值。如要求平稳序列x的均值，并对序列x零均值化，则可用如下命令：ScalarHYPERLINK"mailto:m=@mean(x"m=@mean(x)Genry=x－m其中：Scalar命令在Eviews中表示生成标量数据（均值只是一个数，而不是序列）。Y为对x零均值化后的序列。当然，上述命令也可简化为：Genry=x－@mean(x)习题三：用自相关分析图识别1990年1月至1997年12月我国工业总产值的月度时间序列及其自然对数的平稳性，并说明理由。若不平稳试绘制自然对数序列的一阶逐期差分和一阶季节差分后的我国工业总产值序列的相关分析图。1990年1月至1997年12月我国工业总产值（单位：亿元）年月数据199011421.421367.431719.741759.651795.761848.171637.381670.991760.1101789.5111888.6121981.4199111757.821485.731893.941969.852033.76210371836.381914.792022.2102045.1112069.2122136199211984.221812.432274.742328.952373.162515.8722888232192441.1102502.6112608.8122823.8199312179.122408.732869.442916.753022.163274.572862.982864.292908102911.8113101.3123664.3199412903.322513.83340943499.553642.663871.47337383463.493663.74103753.38113973.17124469.02199512996.722740.333580.943746.353817.964046.673483.983510.693703.1103810.7114091124650.799199613476.622970.333942.644067.654746.89964417.29973806.883746.394011.1104129.6114372.899124991.5199713843.8423181.2634404.4944520.1854638.9964969.9374146.89984198.794536.839104718.91115034.939125545.74实验四ARMA模型的应用【实验目的】熟悉对零均值平稳序列建立ARMA模型的前三个阶段：模型识别、模型参数估计、诊断检验。根据时间序列自相关图对零均值平稳序列进行初步的模型识别。运用Eviews软件估计ARMA模型参数。对所建立的模型是否为适应性模型进行诊断检验。【实验内容】模型识别根据零均值平稳化后的序列的自相关函数和偏自相关函数表现出的特征，对序列进行初步的模型识别（注：这种方法并不总是有效）。模型参数估计Eviews建立ARMA模型的命令用到AR、MA、SAR，SMA等参数项。其中SAR、SMA两参数在建立季节性时间序列模型时要用到。例如：对一个零均值的平稳序列x建立ARMA(2,1)模型，（1）命令操作方式为：lsxar(1)ar(2)ma(1)（2）菜单操作方式：Quick---Estimateequation，输入：xar(1)ar(2)ma(1)，OK。以上述操作方式建模时，Eviews自动采用非线性最小二乘法估计模型参数。模型的诊断检验：判断模型是否为适应性模型判断模型是否为序列的适应性模型，主要根据模型残差是否为白噪声来判断，若残差是白噪声，则可认为此模型是序列的适应性模型，否则，不是。Eviews操作：在模型窗口，View----Residualtests----Correlogram—Qstatistics根据输出的残差的Q统计量判断残差是否为白噪声序列。模型中各项的取舍若建立的模型为适应性模型，还要看输出项中各变量是否显著(通过输出结果中的t统计量值及相应的P值)，对不显著的项，要剔除，然后重新建模。模型的选择（定阶）对于同一个序列来说，可能有多个适应性模型，要从这多个适应性模型中选择，通常根据多个模型输出项中的赤池信息准则(AIC，Akaikeinfocriterion)和施瓦茨准则(SBC，SchwartzBayescriterion)进行比较，一般认为这两个统计量值越小的模型越好。4．模型平稳性和可逆性的判断判断模型是适应性模型后，还应判断模型是否平稳和可逆，判断方法如下。模型输出结果最下方输出的两项，ARinvertedroot(如果有的话)和MAinvertedroot（如果有的话），其含义分别为：invertedARroot：为模型自回归AR部分所对应的差分方程的特征方程的特征根。若特征根的绝对值都小于1，则说明模型是平稳的；若其中有大于或等于1的，说明模型非平稳；若有等于1或很接近于1的，说明原序列为单位根过程，需要先对序列进行差分平稳化变换(有几个单位根，作几阶差分变换)，然后建模。invertedMAroot：为模型移动平均MA部分所对应的差分方程的特征方程的特征根。若特征根绝对值都小于1，则说明模型是可逆的；若有大于或等于1的，说明模型不可逆；若有等于1或很接近于1的，则很有可能在数据处理过程中，对原序列过度差分了，这时需要减少对序列差分的阶数，再重新建模。；【实验步骤】练习一操作文件：参见上机3练习一操作内容：打开一个文件，如arma2.wf1选取一个序列，如x，判断序列是否为平稳，均值是否零均值平稳序列（本例略）。观察该序列自相关图，根据自相关函数滞后二阶截尾，偏自相关函数表现为拖尾，初步判断模型阶数AR(2)。建模：lsxar(1)ar(2)诊断检验：模型是否为序列的适应性模型：检验模型残差是否为白噪声。模型中各项是否显著：用各变量的t检验值及相应的p值。模型选择：先记下拟合的ar(2)模型的AIC和SBC。再拟合其它模型如：ARMA(2,1)，记下输出的AIC和SBC。比较上述结果，看哪个更小。判断模型是否平稳：看invertedARroot是否小于1.练习二：操作文件：参见上机3练习二操作内容：打开一个文件选取一个序列，判断序列是否平稳，均值是否为零，若否，应先将序列转化为零均值平稳序列。转化方法见以前上机实习内容。（2）~（4）同上。实验五时间序列平稳性检验【实验目的】了解并掌握判断非平稳性的单位根检验和非参数检验方法；【实验内容】一、序列非平稳性的单位根检验；二、序列非平稳性的非参检验三、两种方法的比较【实验步骤】序列非平稳性的单位根检验操作文件：zl11.wf1（1）打开zl11.wf1（2）观察序列gjhy的折线图，发现图形有曲折上升的趋势，可以初步判断该序列存在季节性，同时不平稳。在进行单位根检验之前应先进行季节调整。进行季节调整后的图形如图9-2图9-1图9-2（3）进行单位根检验：操作：在VIEW窗口，选unitroottest菜单，在出现的对话框中，按照图示操作，选择ADF，增广的单位根检验。图9-3（4）判断。由于单位根检验是单侧检验，如果统计量的值在临界值右侧，则表示接受原假设，存在一个单位根，序列非平稳。反之序列平稳。由下图结果可以判断，该序列存在一个单位根。图9-4序列非平稳性的非参检验操作文件：zl11.wf1（1）打开zl11.wf1（2）观察序列gjhy的折线图，发现图形有曲折上升的趋势，可以初步判断该序列存在季节性，同时不平稳。在进行单位根检验之前应先进行季节调整。进行季节调整后的图形如图9-2。（3）进行单位根检验：操作：在VIEW窗口，选unitroottest菜单，在出现的对话框中，按照图示操作，选择PP，菲力普-配荣检验。图9-5（4）判断。结果如下：PP检验是非参数检验，判断依据也是根据统计量值和临界值做比较，根据比较结果，如果统计量的值在临界值右侧，则表示接受原假设，存在一个单位根，序列非平稳。反之序列平稳。从图9-6可以看出，采用非参数检验，得到的结果和ADF检验结果恰好相反，该结果是不存在单位根，那么数据的生成过程就可能是由时间趋势项产生的。图9-6两种方法结果的比较这个题目练习了单位根的参数和非参数检验两种方法，从这两种方法可以看出，采用不同的方法最后导致的结果可能是不相同的。但在多数情况下，两种检验方法的结果还是较为一致的。实验六模型预测【实验目的】：（1）进一步熟悉ARMA模型建模过程。（2）利用ARMA模型进行预测。预测说明：Eviews中有两种不同的预测处理方式：Dynamic（动态）和Static（静态）。熟悉对零均值平稳序列建立ARMA模型的前三个阶段：模型识别、模型参数估计、诊断检验。【实验内容】平稳时间序列模型预测非平稳时间序列模型的预测【实验步骤】平稳时间序列模型预测操作文件：zl14.wf1（1）打开zl14.wf1（2）对序列mlpm建立AR(2)模型操作命令：lsmlpmcar(1)ar(2)（3）进行追溯预测：操作：在Equation窗口，选Forecast菜单，在出现的对话框中，选static，将预测结果存入mlpmf1序列中，单击OK。观察输出结果mlpmf1。说明：static为一步超前预测。（4）进行向前多步预测。操作命令：expand1259smpl251259然后在Equation窗口，选Forecast菜单，在出现的对话框中，选Dynamic，并将预测结果保存在mlpmf2序列中，单击OK。观察输出结果mlpmf2。说明：Dynamic为动态预测。注：S.E用于存放预测的估计标准误差，便于计算置信区间。非平稳时间序列预测（操作文件：gdp.wf1）操作步骤：(1)打开gdp.wf1，(2)对序列dlog(gdp)建立ar(2)模型操作命令：lsdlog(gdp)ar(1)ar(2)(3)进行追溯预测：打开forecast对话框，选forecastofgdp，选static，预测结果保存在gdpf1中，单击OK。(4)进行向前多步预测操作命令：expand19782005smpl20022005打开forecast对话框，选forecastofgdp，选dynamic，预测结果保存在gdpf2中，单击OK。观察输出结果。友情提示： 方案 气瓶 现场处置方案 .pdf气瓶 现场处置方案 .doc见习基地管理方案.doc关于群访事件的化解方案建筑工地扬尘治理专项方案下载 范本是 经验 班主任工作经验交流宣传工作经验交流材料优秀班主任经验交流小学课改经验典型材料房地产总经理管理经验 性极强的领域，本范文无法思考和涵盖全面，供参考!最好找专业人士起草或审核后使用。