高铁梅老师的EVIEWS教学第二十二章状态空间模型和卡尔曼滤波课件_第1页
高铁梅老师的EVIEWS教学第二十二章状态空间模型和卡尔曼滤波课件_第2页
高铁梅老师的EVIEWS教学第二十二章状态空间模型和卡尔曼滤波课件_第3页
高铁梅老师的EVIEWS教学第二十二章状态空间模型和卡尔曼滤波课件_第4页
高铁梅老师的EVIEWS教学第二十二章状态空间模型和卡尔曼滤波课件_第5页
已阅读5页,还剩50页未读 继续免费阅读

下载本文档

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

文档简介

第二十二章状态空间模型和卡尔曼滤波StateSpaceModelsandKalmanFilter

60年代初,由于工程控制领域的需要,产生了卡尔曼滤波(KalmanFiltering)。进入70年代初,人们明确提出了状态空间模型的标准形式,并开始将其应用到经济领域。80年代以后,状态空间模型已成为一种有力的建模工具。许多时间序列模型,包括典型的线性回归模型和ARIMA模型都能作为特例写成状态空间的形式,并估计参数值。在计量经济学文献中,状态空间模型被用来估计不可观测的时间变量:理性预期,测量误差,长期收入,不可观测因素(趋势和循环要素)。状态空间模型在经济计量学领域其他方面的大量应用请参见Hamilton(1994)和Harvey(1989)。在一般的统计模型中出现的变量都是可以观测到的,这些模型以反映过去经济变动的时间序列数据为基础,利用回归分析或时间序列分析等方法估计参数,进而预测未来的值。状态空间模型的特点是提出了“状态”这一概念。而实际上,无论是工程控制问题中出现的某些状态(如导弹轨迹的控制问题)还是经济系统所存在的某些状态都是一种不可观测的变量,正是这种观测不到的变量反映了系统所具有的真实状态,所以被称为状态向量。这种含有不可观测变量的模型被称为UC模型(UnobservableComponentModel),UC模型通过通常的回归方程式来估计是不可能的,必须利用状态空间模型来求解。状态空间模型建立了可观测变量和系统内部状态之间的关系,从而可以通过估计各种不同的状态向量达到分析和观测的目的。

EViews状态空间对象对单方程或多方程动态系统提供了一个直接的、易于使用的界面来建立、估计及分析方程结果。它提供了大量的建立、平滑、滤波及预测工具,帮助我们利用状态空间形式来分析动态系统。利用状态空间形式表示动态系统主要有两个优点:第一,状态空间模型将不可观测的变量(状态变量)并入可观测模型并与其一起得到估计结果;其次,状态空间模型是利用强有效的递归算法——卡尔曼滤波来估计的。卡尔曼滤波可以用来估计单变量和多变量的ARMA模型、MIMIC(多指标和多因果)模型、马尔可夫转换模型以及变参数模型。

§22.1

状态空间模型理论及方法

在本节中,我们仅就如何定义并预测一个线性状态空间模型做以简要的讨论,更为详细的内容可以查询Hamilton(1994)、Harvey(1993)。

一、模型表示维向量的动态线性状态空间表示可通过下面的方程组给出:(22.1)(22.2)式中,为维不可观测的状态向量,是服从于零均值正态分布的扰动向量。不可观测的状态向量假定服从于一阶向量自回归过程。我们将第一个方程称为“信号”或“量测”方程,第二个方程称为“状态”或“转移”方程。扰动向量的同一时刻的协方差矩阵为:

(22.3)

[例1]

一阶移动平均模型MA(1)

(22.4)通过定义状态向量可以写成状态空间形式

量测方程

(22.5)

状态方程

(22.6)这种形式的特点是不存在量测方程噪声。

对于任何特殊的统计模型,状态向量的定义是由结构确定的。它的元素一般包含具有实际解释意义的成分,例如趋势或季节要素。状态空间模型的目标是,所建立的状态向量包含了系统在时刻t的所有有关信息,同时又使用尽可能少的元素。所以如果状态空间模型的状态向量具有最小维数,则称为最小实现(MinimalRealization)。对一个好的状态空间模型,最小实现是一个基本准则。然而对于任一特殊问题的状态空间模型的表示形式却不是惟一的,这一点很容易验证。考虑通过定义一个任意的非奇异矩阵B,得到,为新的状态向量。用B矩阵左乘状态方程(22.2),得到

(22.7)式中,,。相应的量测方程是

(22.8)式中。

[例2]

对二阶自回归模型AR(2)

(22.9)考虑两个可能的状态空间形式()是

(22.10)

(22.11)换一种形式

(22.12)

[例3]

由于各种各样的外界冲击和政策变化等因素的影响,经济结构不断发生变化,用OLS等固定参数模型表现不出来这种经济结构的变化,因此,需要考虑采用变参数模型(Time-varyingParameterModel)。下面利用状态空间模型来构造变参数模型。

二、滤波(Filtering)考虑状态向量t在时刻s的条件分布,我们可以定义条件分布的均值和方差矩阵:

(22.16)

(22.17)

期望算子的下标表示条件分布期望作用的期间。令则可得到一个重要的条件分布,即可以得到状态向量t的向前一步均值和向前一步方差。当给定和时,t的条件分布的均值和估计误差的协方差矩阵由下式给定,即:(22.18)(22.19)上面的两个方程称为预测方程(PredictionEquations)。当得到新的观测值,就能够修正的估计,更新方程(UpdatingEquations)是:

(22.20)

(22.21)

(22.22)上述的(22.18)~(22.22)一起构成Kalman滤波的公式。t的条件分布的均值是t在最小均方误差意义下的最优估计量。给出一步向前状态条件均值,我们还可以得到yt的一步向前(线性)最小均方误差估计:

(22.23)一步向前预测误差可以通过下面的公式得到:

(22.24)预测误差的方差被定义为;

(22.25)卡尔曼滤波是在给出新的信息条件下更新状态向量的向前一步预测均值和方差的连续递归算法,具体的递归步骤如上所示。只要给出状态向量初值0=a0和协方差矩阵的初值P0,以及系统矩阵Zt

,Tt

,Ht

,Qt,ct

,dt的值和yt的观测值,就可以利用卡尔曼滤波计算出状态向量和最小均方误差矩阵的估计,同期的或“滤波”的状态向量和方差为{at,Pt},向前一步预测、预测误差、预测误差方差为。

(22.29)

(22.30)n步向前预测:

(22.31)相应的n步向前预测的最小均方误差矩阵为:

(22.32)对于n=1,2,…,同前面一样,在时刻t的可利用信息的基础上,仍然可以被解释为的最小均方误差估计,是误差协方差矩阵的最小均方估计值。n步向前预测可以通过少量地改变卡尔曼递归算法计算出来(Harvey

1989)。对于在s=t+n期间的预测,使用时刻t的信息,利用状态向量和状态协方差阵的预测值只需初始化时刻t+1的卡尔曼滤波,对每一期s预测样本进行重复操作,s=t

+1,…,t+n*。2、动态预测动态预测的概念和其他EViews估计对象的概念相似。在动态预测中,我们从时刻t的预测样本开始,对于预测间隔的每一期n=1,…,n*计算完整的一组n期向前预测。这样,如果我们希望从第t期开始动态地预测到第t+n*期,则对第t+1期我们计算一步向前预测,对第t+2期计算2步向前预测,持续下去,直到对第t+n*期计算n*步向前预测。这表明对n步向前预测,我们只需初始化时刻t+1的卡尔曼滤波,并且运用期内信息对期外进行滤波。然而对于动态预测,要求计算所有的预测值,因为从预测期开始信息系统没有被更新。

3、平滑预测相应地,我们可以计算平滑预测,就是使用预测样本的所有可利用信息,计算平滑值(例如)。这些向前预测值可以通过初始化预测期间的初始状态求得,使用所有相关的信号数据对整个预测期间进行卡尔曼平滑滤波。这个技术对于设置信号的信息被用来对全部预测样本插值时是很有用的。

我们解释一下上面描述的预测方法。对于传统的n步向前预测和动态预测,在预测窗口的开始处,通常利用一步向前预测的状态向量和方差初始化状态向量。对于平滑预测,一般使用状态向量和方差的相应的平滑值进行预测的初始化。对某些情况,预测滤波和平滑可以选择初始值的不同设置。EViews的预测程序提供了可控的初始设置(在“状态空间”过程中叙述)。然而,如果选择了不同的设置,根据可利用信息,预测的解释将发生变化。

五、估计未知参数在应用卡尔曼滤波时,我们必须先把系统矩阵中的未知元素用它们的估计值代替。在和服从正态分布条件下,样本的对数似然函数为:

(22.33)其中,。

EViews利用数值微分和标准迭代技术求解带有未知参数的似然函数。§22.2在EViews中定义状态空间模型EViews可以处理大量的单方程和多方程状态空间模型,提供了指定系统方程、协方差矩阵和初始条件控制的详细方法。在定义和估计一个状态空间模型时,第一步是创建一个状态空间对象。从主菜单中选择Objects/NewObject/Sspace,或在命令窗口键入命令sspace。EViews将创建一个状态空间对象,并打开一个空的状态空间说明窗口。有两种方法定义一个状态空间模型,最简单的方法就是利用EViews中的“自动指定”功能引导状态空间模型的标准形式。这种方式只需在状态空间过程Procs中选择DefineStateSpace功能

,就可以弹出定义对话框,指导创建一个状态空间的过程。这一方式的详细介绍见“自动定义”一节。描述状态空间模型的更一般方法是使用关键字和文本来描述量测方程、状态方程、误差结构、初始条件和待估参数的初值。下面来介绍描述状态空间对象的一般语法。

一、模型指定的语句

1、量测方程作为缺省,如果一个方程通过关键字“@SIGNAL”来明确定义,或没有用关键字,EViews将把其作为量测方程处理。要注意以下几点:(1)量测方程的因变量可以包含表达式。(2)量测方程中不能包含量测变量的当期和未来值,在量测方程中任何滞后量测变量都被看作多步向前预测的预测值看待。(3)量测方程必须是同期状态向量的线性方程。状态向量的非线性或存在超前或滞后状态变量将导致错误的信息。(4)量测方程中可以包含外生变量和未知参数,也可以是这些元素的非线性形式。量测方程可以包含误差或误差方差指定的选项,如果方程中不包含误差或误差方差,方程是确定性的。状态空间模型中误差指定的详细内容参看后面的“误差和方差”。

例子:下面是有效的量测方程的定义(注:下面量测方程中的sv1,sv2,sv3,sv4是状态向量)@signaly=sv1+sv2*x1+sv3*x2+sv4*y(-1)+[var=exp(c(1))]

log(p)=sv1+c(1)+c(3)*x+sv2*yz=c(1)+sv1+sv2*x1+sv3*x2+[var=exp(c(2))]下面是不正确的方程的指定:

@signaly=sv1*sv2*x1+[var=exp(c(1))]

log(p)=c(1)+c(3)*x+sv1(-1)z=sv1+sv2*x1+c(3)*z(1)+c(1)+[var=exp(c(2))]因为它们至少违背了上面描述条件中的一个条件(其顺序是:状态向量的非线性、状态向量的滞后、量测向量的超前)。

例子:下面两个状态方程定义了一个服从AR(2)过程的不可观测误差:@statesv1=c(2)*sv1(-1)+c(3)*sv2(-1)+[var=exp(c(5))]@statesv2=sv1(-1)第一个关于sv1的方程,根据AR(1)的系数c(2),和AR(2)的系数c(3),确定AR(2)模型的参数。误差方差的指定在方框中给出。sv2的状态方程定义为变量sv1的一步滞后,所以sv2(-1)表示sv1的两步滞后。下面是不正确的状态方程:@stateexp(sv1)=sv1(-1)+[var=exp(c(3))]@statesv2=log(sv2(-1))+[var=exp(c(3))]@statesv3=c(1)+c(2)*sv3(-2)+[var=exp(c(3))]因为它们至少违背了上面描述条件中的一个条件(其次序是:状态方程因变量是表示式,状态变量是非线性的,出现状态变量的多期滞后)。

3、误差与方差在误差项的处理中,状态空间对象方程的指定在某种程度上是唯一的。EViews总是把一个隐含的误差项加到一个方程或系统对象的各个方程中去。但如不特殊指定,状态空间量测或状态方程中不能包含误差项。误差项必须被加到(在方括号中)指定方程的后面。把一个误差项加到状态空间方程中去的最简单的方法是指定误差项的方差。即加一个误差表达式到已存在的方程中去。误差表达式由关键字“var”和一个赋值语句组成(用方括号括起)。@signaly=c(1)+sv1+sv2+[var=1]@statesv1=sv1(-1)+[var=exp(c(2))]@statesv2=c(3)+c(4)*sv2(-1)+[var=exp(c(2)*x)]指定的方差可以是已知常数值,也可以是包含待估计未知参数的表达式。还可以在方差中使用序列表达式建立时变参数模型。这种方差的直接指定方法不允许不同方程的误差之间存在相关关系。作为默认,EViews假定误差项之间的协方差为零。如果指定误差项间存在相关关系,需要使用“命名误差”方法指定它们间的关系。“命名误差”方法包括两部分:

(1)首先,必须通过加一个由关键字“ename”后接等号和变量名的误差表达式为方程中的残差序列命名。y=c(1)+sv1*x1+[ename=e1]@statesv1=sv1(-1)+[ename=e2]

(2)其次,需要键入由关键字“@evar”后接一个误差的方差或两个误差之间的协方差的赋值语句。@evarcov(e1,e2)=c(2)@evarvar(e1)=exp(c(3))@evarvar(e2)=exp(c(4))*x

二、模型指定的例子下面的例子描述一个解释变量带有随机系数sr1和递归系数sr2的模型,sr1和因变量cp的误差项间存在相关关系(22-1/cp_y):

例1量测方程:状态方程:~在方程中,cp(实际消费),y(实际收入),r(利率)。

cp=c(1)+sr1*y+sr2*r(-2)+[ename=e1,var=exp(c(2))]@statesr1=sr1(-1)+[ename=e2,var=exp(c(3))]@statesr2=sr2(-1)@evarcov(e1,e2)=c(4)模型的方差、和协方差g由参数c(2)、c(3)和c(4)确定,cp和sr1的方差被限制为参数的非负函数。

1、指定未知参数(超参数)的初始值上例中c(1),c(2),c(3),c(4)的初值分别为3.88,0.165,0.001,0.28。除非另外指定,否则EViews将用相应系数向量的当前值初始化所有参数。可以通过在状态空间形式指定中使用PARAM或@PARAM语句来明确指定合适的参数值。例:paramc(1)3.88c(2)0.165c(3)0.001c(4)0.28

2、指定状态向量和方差矩阵的初始条件缺省时,EViews将自动处置初始条件。对一些平稳模型,其稳定状态条件使我们能够解出0和P0的值。当不能解出初始条件时,EViews将把初始值处理为扩散的,设置0=0,给P0一个任意大的正数乘单位矩阵代表其值的不确定性。可能有0和P0的先验信息,这样,可以使用关键字@mprior或@vprior创建一个包含适当值的向量或矩阵。设置初始状态值时,键入“@mprior”,后接向量对象名。向量对象的长度必须与状态空间的维数相匹配。其元素的顺序要与指定窗口中状态向量的顺序相一致。@mpriorv1@vpriorm1要设置状态向量方差矩阵P的初始值,键入“@vprior”后接一个标志对象名(表示这是一个标志对象,而不是一个普通的矩阵对象)。标志的维数必须与状态空间的维数相匹配,其顺序要与指定的状态向量的顺序相一致。如果你想设置一个元素为扩散的,只需把该元素设置为“NA”缺省值。EViews将重新设置所有相应的方差和协方差是扩散的。相比较,会发现使用命令创建和初始化向量和矩阵更容易一些。可以在命令窗口中键入下面的命令:vector(2)a0a0.fill1,0matrix(2,2)p0p0.fill(b=c)1,0.5,0.5,2然后在编辑状态空间模型指定窗口,可以把下面两行加到状态空间对象中去:@mpriora0@vpriorp0同样,也可以在命令窗口键入下面的命令:ss1.append@mpriora0ss1.append@vpriorp0关于矩阵对象与填加“fill”或追加“append”命令的详细说明,参见命令与程序指南“CommandandProgrammingReference”。

三、模型定义视窗状态空间模型是比较复杂的。为了帮助检验模型定义,EViews提供了视窗功能,允许在交互方式下查看文本指定,在当前参数值下检验系统估计矩阵估计值。点击View菜单选择Specification…,不管状态空间模型是否被估计,下面的指定窗口都可以被使用。

(1)文本窗口

这是一个常见的模型指定的文本视窗。当创建或编辑状态空间模型指定时,可以使用这个窗口。文本窗口也可以通过点击状态空间工具栏的Spec按钮进入。

(2)系数描述

状态空间模型指定结构的文本描述。左边的变量yt和被表示为状态向量和残差项的线性函数。矩阵的元素是相应的系数。例如,例1模型的系数描述视图如下:

(3)协方差描述

状态空间模型协方差矩阵的文本描述。例如,例1模型有下面的协方差描述视图。

(4)系数值

用当前参数估计的量测方程和状态方程结构的数字描述。如果系统系数矩阵是时变的,EViews将提示对矩阵估计选择一个日期/观测值。

(5)协方差值

用当前参数估计的状态空间模型指定结构的数值描述。如果系统协方差矩阵是时变参数的,EViews将提示对矩阵估计选择日期/观测值。

四、自动指定为了帮助创建一个状态空间模型,EViews提供了一个“自动指定”工具栏,可以在对话框中为模型创建一个文本表示。如果模型是具有固定参数、递归参数、及不同的随机系数,或者误差项有一般ARMA结构的标准回归模型,这个工具是非常有用的。在状态空间过程procs中,选择Procs/DefineStateSpace…。EViews将打开一个三标签的对话框。第一个标签对话框BasicRegression被用来描述模型的基本回归部分。键入因变量和带有固定或递归系数的回归变量。在建立指定时EViews使用COEF对象代表未知参数。在底部,可以指定误差项一个ARMA结构。在这里,我们为上面的例子指定一个说明。第二个标签对话框StochasticRegressors被用来加带有随机系数的回归变量。在四个编辑区域中键入合适的回归变量。EViews允许定义具有如下四项组合的回归变量:固定均值系数、AR(1)系数、随机游动系数、带有漂移的随机游动系数。

最后,EViews允许选择状态空间模型的基本方差结构。点击第三个标签对话框VarianceSpecification,

为量测方程或状态方程选择方差矩阵类型:单位矩阵(Identity)、共同对角矩阵(CommonDiagonal,对角元素是共同的方差)、一般对角矩阵(Diagonal)、无限制矩阵(Unrestricted)。对话框还允许为量测方程和状态方程选择非零的误差协方差阵。

需要强调指出的是,状态空间模型可以不必被对话框提供的选择限制。如果发现自动指定对话框的限制了模型指定,可以简单地使用它建立一个基本的指定,然后利用更一般的文本工具描述模型。

五、估计状态空间模型一旦已经指定了一个状态空间模型,并且验证模型定义是正确的,打开估计对话框估计模型,点击工具菜单的Esimate按钮或者选择Procs/Estimate…。和其他估计对象一样,EViews允许选择估计样本区间,循环的最大次数,收敛值,估计算法,导数计算设置和是否显示初始值。对大部分问题,缺省设置提供一个好的初始设置。在进行模型估计时要注意下面两点:

(1)尽管EViews中卡尔曼滤波程序可以自动处理样本中的缺省值,但EViews要求估计样本必须是连续的,连续的观测值之间不能有缺口。

(2)如果模型定义中有未知系数,为用卡尔曼滤波估计状态空间模型,需要指定初值。

六、解释估计结果在选择方差选项并点击OK以后,EViews在状态空间窗口显示估计结果。例2建立变参数的IS模型:其中GDP/P是实际收入,r是利率,则自发支出0用状态变量sv1表示,1用状态变量sv2表示,1/1(即sv2的倒数)是IS曲线的斜率,随着时间的推移,可以观察IS曲线的变化,是越来越陡峭还是越来越平坦,从而看出货币政策的有效性。(22_1/is):

y=sv1+sv2*r(-3)+[var=abs(c(4))]@statesv1=sv1(-1)+[var=abs(c(2))]@statesv2=c(1)*sv2(-1)+[var=abs(c(3))]paramc(1)0.94c(2)3.5c(3)0.00001c(4)81

估计这个模型,EViews将打开估计输出视窗:

输出视窗与其他EViews估计对象相似。上面的信息描述估计的基本信息:状态空间对象名、估计方法、估计的时间和日期、样本区间、包含的样本数、收敛信息和系数估计。还显示最终的状态向量的一步向前预测值和误差协方差矩阵对角线元素的平方根。视窗的底部描述了最大对数似然估计值,估计参数的数目,以及相关的信息准则。

对于状态向量和协方差矩阵整个路径的设置,EViews提供了一系列视窗和过程来检验状态结果。§22.3

状态空间模型的视窗和过程EViews提供了一系列专门的工具用来指定和检验状态空间模型。与其他的估计对象相比较,状态空间对象提供了附加的视窗和过程来检验估计结果,处理推断和指定检验,并且提取结果到其他EViews对象中去。

一、状态空间视窗

1、一般视窗许多状态空间视窗与前面讨论的相似:

(1)模型定义视窗(Specification)查看模型文本定义,系数和协方差定义。仅看文本窗口还可以点击系统工具条的“Spec”。

(2)估计结果(EstimationOutput)还可以点击系统工具条的“Stats”,显示估计结果。

(3)梯度视窗(GradientsandDerivatives)和其他的估计对象视窗相似,如果状态空间包含待估参数,该视窗提供了被估计参数(已估计)的对数似然估计的梯度的简要可视信息或当期参数值。

(4)实际值、拟合值和残差(Actual,Predicted,ResidualGraph)用图表的方式显示量测方程因变量实际值和一步向前拟合值,和一步向前标准残差。

(5)估计系数协方差矩阵(CoefficientCovarianceMatrix)

(6)Wald检验(WaldCoefficientTests)

允许做估计系数的假设检验。

(7)Label视窗允许为状态空间对象做注释。注意,除了Label和模型定义(Specification)视窗之外,其余的视窗只有在状态空间模型被正确估计的情况下才可以使用。

2、量测视窗当点击View/SignalViews,EViews显示一个包含视窗选择的次级菜单。·ActualSignalTable和ActualSignalGraph显示量测方程因变量的表和图的形式。如果有多个量测方程,Eviews将按其顺序显示多个序列。在状态空间模型没有被估计的条件下,这两个选项也是可以利用的。

·GraphSignalSeries…,可以打开一个对话框,选择显示结果。对话框允许在下列选项中做出选择:量测变量一步向前预测,相应的一步预测残差,标准化的一步残差,平滑的量测变量,平滑的量测方程扰动项,或标准平滑的量测方程扰动项。加上2倍的标准误差的点线图。·STd.ResidualCorrelationMatrix和Std.ResidualCovarianceMatrix显示量测方程一步向前预测标准差的相关阵和协方差阵。

3、状态视窗为了检验不可观测的状态变量,点击View/StateViews显示状态方程子菜单。EViews允许检验状态变量的初值和终值,或者画状态向量的各种平滑和滤波序列图。在估计前后,视窗中有两个选项是可利用的:·InitialStateVector和InitialStateCovarianceMatrix显示状态向量的初始值,和协方差阵。如果未知参数已被估计,EViews将使用估计值计算初始条件。如果状态空间模型没有被估计,使用当期系数值来估计初始条件。

在EViews正在利用系统矩阵的当前值求解初始条件时,这个信息是特别有意义的。在开始估计有困难的情况下,可以从任意初始参数值出发来估计初始条件。窗口中的其他选项,只对已成功估计的模型有效:·FinalStateVetor和FinalStateCovarianceMatrix,显示状态向量终值,和协方差矩阵终值,在对参数估计后估计得到。·GraphStateSeries…菜单,显示包含状态向量信息选项的对话框。可以画下列变量的线性图:状态向量的一步向前预测,经过滤波得到的同期状态向量,平滑的状态向量估计,平滑的状态扰动项的估计,标准的平滑状态扰动项。在每一线性图中,显示数据被包在其倍的标准差带中。

二、状态空间过程可以使用Eviews过程创建、估计、预测状态空间模型和从指定的状态空间模型生成数据。

1、自动指定状态空间模型(DefineStateSpace…)激活自动指定对话框。可以在交互的方式下指定状态空间模型。

2、估计(Estimate)估计指定模型的参数(22-1/IS)。上面两项功能在模型估计前后都可以使用。自动指定工具将代替存在的状态空间指定和清除任何结果。估计将替代已存在的结果。

3、预测如果已经对状态空间模型进行了估计,EViews提供生成数据的其他工具:预测允许利用选择的预测方法和初始化方法,来产生状态变量、量测变量和联合标准误差的预测。

(1)

选择预测方法

可以在动态预测,平滑预测和n期向前预测三种方法中选择其一。注意,任何在量测方程右边的延迟内生变量都被做为外生变量的看待。

(2)预测结果输出

EViews允许在工作文件中以序列形式存储预测输出结果。只须点击输出框,在相应的编辑区域指定序列名。可以指定一列变量名或一个通配符表达式。如果选择列变量名,变量名的数目必须与指定的量测变量的数目相匹配。如果输出序列的名字已经在工作文件中存在,EViews将全部覆盖序列的内容。如果使用一个通配符表达式,EViews将利用通配符表达式在适当的位置替代每一个量测变量的名字。例如,如果有一个具有量测变量y1、y2的模型,选择通配符“*F”存储一步预测的结果,EViews将使用序列名y1F和y2F存贮输出结果。对该功能有两点限制:一是不能使用通配符表达式“*”存储量测变量的结果,因为这将导致对原始量测数据的覆盖。二是当量测方程因变量通过表达式指定,或量测变量出现在多个方程中,不能使用通配符。对于这两点,EViews将不能产生新的序列,只产生错误信息。需要注意的是,如果量测方程的因变量是一个表达式,EViews将只提供表达式的预测,因此,如果量测变量是log(y)

,EViews将只预测y的对数。我们分别利用动态预测,平滑预测和1(n=1)期向前预测三种方法在期间1993:1至2000:4进行预测,分别得到预测结果YF、YF1、YF2。下面是3个预测值和实际值Y画在一起的图形,可以看出动态预测YFD和平滑预测YFS的结果重合了。1期向前预测YF1与实际值接近。

(3)设置初始条件

有几个选项是可以选用的。可以使用状态变量和状态协方差的One-stepahead预测或Smoothed预测估计值做为预测的初始值。这两种初始化方法在使用估计样本的信息数量上是不同的。一步向前预测使用直到预测期开始的信息,平滑预测使用全部估计期的信息。相应地,还可以使用EViews计算初始条件。在预测开始时,EViews将求解Riccati代数方程,获得状态向量和状态向量协方差的初始值。如果不能得到初始状态的解,EViews将使用扩散的初始值。最后,可以选择提供一个向量和包含预测初始值的特征对象。选择User并在编辑区域输入有效的EViews对象名。当计算n步向前预测时,EViews将调整预测期的起始时期,以便使用指定的方法可以得到每一期的初始条件。对于一步向前预测和平滑预测方法,意味着在预测期间最初n-1个观测值被做为估计样本,更

温馨提示

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

评论

0/150

提交评论