ARMA模型的eviews的建立-时间序列分析实验指导1_第1页
ARMA模型的eviews的建立-时间序列分析实验指导1_第2页
ARMA模型的eviews的建立-时间序列分析实验指导1_第3页
ARMA模型的eviews的建立-时间序列分析实验指导1_第4页
ARMA模型的eviews的建立-时间序列分析实验指导1_第5页
已阅读5页,还剩30页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

时间序列分析实验指导统计与应用数学学院前言随着计算机技术的飞跃发展以与应用软件的普与,对高等院校的实验教学提出了越来越高的要求。为实现教育思想与教学理念的不断更新,在教学中必须注重对大学生动手能力的培训和创新思维的培养,注重学生知识、能力、素质的综合协调发展。为此,我们组织统计与应用数学学院的部分教师编写了系列实验教学指导书。这套实验教学指导书具有以下特点:①理论与实践相结合,书中的大量经济案例紧密联系我国的经济发展实际,有利于提高学生分析问题解决问题的能力。②理论教学与应用软件相结合,我们根据不同的课程分别介绍了SPSS、SAS、MATLAB、EVIEWS等软件的使用方法,有利于提高学生建立数学模型并能正确求解的能力。这套实验教学指导书在编写的过程中始终得到安徽财经大学教务处、实验室管理处以与统计与应用数学学院的关心、帮助和大力支持,对此我们表示衷心的感谢!限于我们的水平,欢迎各方面对教材存在的错误和不当之处予以批评指正。统计与数学模型分析实验中心2007年2月目录TOC\o"1-1"\h\z\u实验一EVIEWS中时间序列相关函数操作 -1-实验二确定性时间序列建模方法 -9-实验三时间序列随机性和平稳性检验 -18-实验四时间序列季节性、可逆性检验 -21-实验五ARMA模型的建立、识别、检验 -27-实验六ARMA模型的诊断性检验 -30-实验七ARMA模型的预测 -31-实验八复习ARMA建模过程 -33-实验九时间序列非平稳性检验 -35-实验一EVIEWS中时间序列相关函数操作【实验目的】熟悉Eviews的操作:菜单方式,命令方式;练习并掌握与时间序列分析相关的函数操作。【实验内容】一、EViews软件的常用菜单方式和命令方式;二、各种常用差分函数表达式;三、时间序列的自相关和偏自相关图与函数;【实验步骤】一、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),观察其趋势图、自相关图。实验二确定性时间序列建模方法【实验目的】熟悉确定性时间序列模型的建模原理;掌握确定性时间序列建立模型的几种常用方法。【实验内容】一、多项式模型和加权最小二乘法的建立;二、单参数和双参数指数平滑法进行预测的操作练习;三、二次曲线和对数曲线趋势模型建立与预测;【实验步骤】

一、多项式模型和加权最小二乘法的建立;1、我国1974—1994年的发电量资料列于表中,已知1995年的发电量为10077.26亿千瓦小时,试以表1.1中的资料为样本:据拟合优度和外推检验的结果建立最合适的多项式模型。采用加权最小二乘法估计我国工业发电量的线性趋势,并与普通最小二乘法估计的线性模型进行比较,列出OLS方法预测值和W=0.6,W=0.7时1992到1995年预测值以与相对误差。74-7879-8384-8889-9394-9516682820377058489281195830064107621210077.26203130934495677522343277497375392566351454528395操作过程:建立WORKFILE:CREATEA19741995生成新序列Y:datay生成新的时间趋势序列t:genrt=@trend(1973)建立系列方程:smpl19741994lsyctlsyctt^2lsyctt^2t^3通过拟合优度和外推检验的结果发现一元三次多项式模型效果最好。首先生成权数序列:genrm=sqr(0.6^(21-t))加权最小二乘法的命令方式:ls(w=m)yct普通最小二乘法命令方式:lsyct进行预测:打开对应的方程窗口,点forecast按纽,将出现对话框,修改对话框samplerangeforforecast中的时间期限的截止日期为预测期.相对误差的计算公式为:(实际值-预测值)/实际值二、单参数和双参数指数平滑法进行预测的操作练习2、某地区1996~2003年的人口数据如表1.2,运用二次指数平滑法预测该镇2004年底的人口数(单位:人)。19961997199819992000200120022003114333115823117171118517119850121121122389123626建立WORKFILE:createU19962004建立新序列Y和T:datay然后输入数值。genrt=@trend(1995)打开y序列,点击exponentialsmoothing按纽,出现如图所示对话框按照图示选项点击确定即可。3、某地区1996—2003年农村用电量数据见表1.3,试利用Holt双参数指数平滑法预测该地区2004年该地区农村用电量(单位:千瓦时)。19961997199819992000200120022003844.5963.21106.91244.81473.91655.71812.71980.1建立WORKFILE:createU19962004建立新序列Y和T:datay然后输入数值。genrt=@trend(1995)打开y序列,点击exponentialsmoothing按纽,出现如图所示对话框按照图示选项点击确定即可。三、二次曲线和对数曲线趋势模型建立与预测;4、我国民航客运量数据的季节调整。有关数据如表1.4,对序列进行季节调整。(1指1993年10月,54指1998年3月)并对调整后序列建立二次曲线和对数曲线趋势模型,得到两个方程的民航客运量趋势估计值,并进行季节调整,求出两个趋势方程建立的季节模型预测值。(选做)12345673282632512412493163441112131415161738436840136333636633121222324252627397.31463509474508458.944123132333435363744748343951455048953441424344454647416451486.2507458.9949356251525354398442404.55428实验三时间序列随机性和平稳性检验【实验目的】认识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)写出各模型的估计结果。实验四时间序列季节性、可逆性检验【实验目的】观察具有实际背景的经济数据,判断其是否平稳、是否含有季节性,均值是否为零。能运用合适的方法如差分、季节差分、取对数、平方根等,使序列变为平稳序列;平稳序列减去其均值,使其零均值化。【实验内容】一、判断序列的平稳性和可逆性,给出相应判断依据,并写出模型形式。二、找出自己感兴趣的数据,判断数据是否平稳,是否具有季节性,均值是否为零等。【实验步骤】练习一操作文件: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)再通过季节差分消除序列的季节性(如果有的话)。差分函数的使用可见前两次上机实习内容。使平稳序列零均值化的方法:在Eviews中可通过函数@mean()求序列的均值。如要求平稳序列x的均值,并对序列x零均值化,则可用如下命令:Scalarm=@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)菜单操作方式:QuickEstimateequation,输入:xar(1)ar(2)ma(1),OK。以上述操作方式建模时,Eviews自动采用非线性最小二乘法估计模型参数。模型的诊断检验:判断模型是否为适应性模型判断模型是否为序列的适应性模型,主要根据模型残差是否为白噪声来判断,若残差是白噪声,则可认为此模型是序列的适应性模型,否则,不是。Eviews操作:在模型窗口,ViewResidualtestsCorrelogram—Qstatistics根据输出的残差的Q统计量判断残差是否为白噪声序列。怎么根据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)同上。实验六ARMA模型的诊断性检验【实验目的】通过练习,进一步熟悉建模步骤:模型识别,参数估计,诊断检验(适应性检验、模型定阶等)。【实验内容】(1)三个模型是否都为适应性模型?(2)哪个模型更佳?(3)三个模型中均包含了常数项,其与序列均值有何关系?(4)各个模型的估计中,实际用到的观察值的个数分别为多少?【实验步骤】操作文件:zl1.wf1~zl20.wf1与其它具有实际背景的数据。练习一zl14.wf1磨轮剖面数据,见附录。操作步骤:1、判断序列mlpm是否平稳,均值是否为零。2、根据自相关图,进行模型识别。3、建立模型:lsmlpmcar(1)ar(2)4、模型诊断检验:看此模型是否合适。5.再分别建立两个模型,且重复上述步骤。lsmlpmcma(1)ma(2)lsmlpmcar(1)ma(1)练习二、依据上述思路,请分别对zl17.wf1,zl18.wf1,gdpindex.wf1等文件进行操作。(该组文件中的数据本身为平稳数据)练习三、操作文件:zl1.wf1,zl3.wf1,zl4.wf1,zl5.wf1,zl19.wf1,gdp.wf1,stpoor.wf1。(该组文件中的数据均非平稳,建模前需先作适当变换)练习四、操作文件:zl2.wf1,zl11.wf1,zl20.wf1,usagnp.wf1等。(该组文件中的数据均含有季节性,建模前需作适当变换)实验七ARMA模型的预测【实验目的】:(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。观察输出结果。实验八复习ARMA建模过程【实验目的】复习利用Eviews对时间序列建立ARMA模型的过程【实验内容】ARMA模型建模前的准备:判断序列是否平稳a.通过序列自相关图、趋势图等进行判断b.若序列不平稳:均值非平稳序列通过差分变换转换为平稳方差非平稳序列通过对数变换等转化为平稳序列c.模型平稳化以后,将序列零均值化模型识别主要通过序列的自相关函数、偏自相关函数表现的特征,进行初步的模型识别模型参数估计在Eviews中估计ARMA模型的方法估计模型以后要能写出模型的形式(差分方程形式和用B算子表示的形式)模型的诊断检验根据模型残差是不是白噪声来判断模型是否为适应性模型能根据输出结果判断模型是否平稳,是否可逆若有多个序列是模型的适应性模型,会用合

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论