VAR模型与向量VECM模型.docx
- 文档编号:10887931
- 上传时间:2023-05-28
- 格式:DOCX
- 页数:47
- 大小:55.84KB
VAR模型与向量VECM模型.docx
《VAR模型与向量VECM模型.docx》由会员分享,可在线阅读,更多相关《VAR模型与向量VECM模型.docx(47页珍藏版)》请在冰点文库上搜索。
VAR模型与向量VECM模型
第7章向量自回归模型(VAR与向量误差修正模型(VEC
§向量自回归模型(VAR(p))
传统的经济计量学联立方程模型建摸方法,是以经济理论为基础来描述经济变量之间的结构关系,采
用的是结构方法来建立模型,所建立的就是联立方程结构式模型。
这种模型其优点是具有明显的经济理论含义。
但是,从计量经济学建摸理论而言,也存在许多弊端而受到质疑。
一是在模型建立之处,首先需要明确哪些是内生变量,哪些是外生变量,尽管可以根据研究问题和目
的来确定,但有时也并不容易;
二是所设定的模型,每一结构方程都含有内生多个内生变量,当将某一内生变量作为被解释变量出现
在方程左边时,右边将会含有多个其余内生变量,由于它们与扰动项相关,从而使模型参数估计变得十
分复杂,在未估计前,就需要讨论识别性;
三是结构式模型不能很好地反映出变量间的动态联系。
为了解决这一问题,经过一些现代计量经济学家门的研究,就给出了一种非结构性建立经济变量之间
关系模型的方法,这就是所谓向量自回归模型(VectorAutoregressionModel)。
VAR莫型最早是1980年,
由引入到计量经济学中,它实质上是多元AR莫型在经济计量学中的应用,
VAR莫型不是以经济理论为基础描述经济变量之间的结构关系来建立模型的,它是以数据统计性质为
基础,把某一经济系统中的每一变量作为所有变量的滞后变量的函数来构造模型的。
它是一种处理具有相
关关系的多变量的分析和预测、随机扰动对系统的动态冲击的最方便的方法。
而且在一定条件下,多元MA
模型、ARM模型,也可化为VAR莫型来处理,这为研究具有相关关系的多变量的分析和预测带来很大方便。
7.1.1VAR模型的一般形式
1、非限制性VAR模型(高斯VAR模型),或简化式非限制性VARI型
(y1ty2t..jkt)为
k维随机时间序列,p为滞后阶数,
Ut(u1t山…也)为一k维随机扰动的
且有结构关系
(1)
(1)
(1)
⑵
(2)
⑵
y
a11y1t1
a12y2t1
■■-a1kykt1
a11y1t2
a12y2t2■
a1kykt2■■■
」p)
(p)
a(p)
a11y1tpa
112y2tp■‘
…a1kyktp
U1t
a
(1)
a
(1)
a
(1)
a
(2)
a
(2)
a
(2)
y2t
a21y1t1
a22y2t1
…a2kykt1
a21y1t2
a22y2t2
…a2kykt2■■■
_(p)
a(p)
」p)
a21y2t
pa12y2t
p■■■a2k
YktpU2t
(1)
(1)
(1)
⑵
(2)
(2)
ykt
ak1y1t1
ak2y2t1
■■■akkYkt1
ak1y1t2
a12y2t2
■■■a1kykt2■■■
„(p)
(p)
」p)
ak1y1t
pak2y2t
a
pkk
yktpUkt
时间序列,
设y
t12…,T
(7.1.1)
若引入矩阵符号,记
可写成
A
a⑴a⑴
2122
a⑴
2k-
i
1,2,...,p
a⑴a⑴
k1k2
a⑴
kk
ytAyt1
A2yt2...
ApytpUt,t
1,2,...,T
(7.1.2)
(i)(i)(i)
a11a12...a你
进一步,若引入滞后算子L,则又可表示成
A(L)%ut,t1,2,...,T
其中:
A(L)IkALAL2...ApLP,为滞后算子多项式.
如果模型满足的条件:
①参数阵Ap0,p0;
(7.1.3)
②特征方程det[A(L)]lkALA.L2...ApLp
0的根全在单位园外;
③Ut~iidN(0,),t1,2,...,T,即Ut相互独立,同服从以E(ut)0为期望向量、
Cov(ut)E(qut)为方差协方差阵的k维正态分布。
这时,Ut是k维白噪声向量序列,由于Ut没有
结构性经济含义,也被称为冲击向量;
Cov(5Xtj)E(UtXtj)0,j1,2,...,即ut与xt及各滞后期不相
关。
则称上述模型为非限制性VARI型(高斯VARI型),或简化式非限制性VAR莫型。
2、受限制性VAR莫型,或简化式受限制性VAR莫型
如果将y(y1ty2t...ykt)做为一k维内生的随机时间序列,
受d维外生的时间序列xt(x1tx2t..xdt)
影响(限制),贝UVAR莫型为
ytAyt1A?
yt2
ApytpDxt
Ut,t1,2,...,T
(7•1.4)
或利用滞后算子表示成
A(L)yt
DxtUt,t1,2,...,T
(7.1.5)
其中:
dt1dt2...dtd
d21d22...d2d
dk1dk2...dkd
此时称该模型为受限制性VAR莫型,
简化式受限制性VAR莫型。
对于受限制性VAR莫型,可通过
yt(y1ty2t...ykt)对x佃作ols^归,得到残差估计
%yt?
,从而将%变换成
(15.1.2)或()形式的非限制性VARI型,即
...Ap%pUt,t1,2,...,T
(7.1.6)
A(L)%Ut,t1,2,...,T
(7.1.7)
这说明受限制性VARI型可化为非限制性VAR莫型。
简化式非限制、受限制VAR莫型,皆简记为VAR(p)。
3、结构式非限制性VAR模型
如果y(yity2t...ykt)中的每一分量受其它分量当期影响,无d维外生的时间序列为(刈x2t..Xdt)
影响(限制),则模型化为
A)yt
Ayt1A稣2...ApytpUt,t
1,2,...,T
(7.1.8)
或利用滞后算子表示成
A(L)ytUt,
t1,2,...
T
(7.1.9)
1a(o)12
a(°)
…a1k
(o)
(o)
其中:
A21
a
2k,这时的A(L)
AoAL
A.L2...ApLp
a化2
此时称该模型为结构式非限制性VAR莫型。
如果A可逆,既逆阵A10存在,则结构式非限制性VAR莫型可化为简化式非限制性VAR模型
1111ytAoAyt1AoAyt2...A°ApytpAoUt,t1,2,...,T(7.1.10)
或利用滞后算子表示成
1
A(L)ytAoUt,t1,2,...,T(7.1.11)
这时,其中的A(L)IA1oA1LA1oA2L2...A;ApLP
4、结构式受限制性VAR模型
如果将yt(y1ty2t...ykt)做为一k维内生的随机时间序列,其中每一分量受其它分量当期影响,且还
受d维外生的时间序列xt(x1tx2t..Xdt)影响(限制),则VAR模型为
(7.1.12)
A0%A』t1A?
%2…Ap%pDxtUt,t1,2,...,T
或利用滞后算子表示成
A(L)ytDxtUt,t1,2,...,T(7.1.13)
此时称该模型为结构式受限制性VAR莫型。
如果A可逆,既逆阵A10存在,则结构式受限制性VAR莫型可化为简化式受限制性VAR模型
11111ytAoAyt1aoAyt2...AoApytpA°DXtAoUt,t1,2,...,T(7.1.14)
或利用滞后算子表示成
11
A(L)ytAoDxtA。
山,t1,2,...,T(7.1.15)
这时,其中的A(L)IA1oA1LA1oA2L2...A1oApLp
结构式非限制、受限制VAR莫型,皆简记为SVAR(p)。
7.1.2简化式VAR莫型的参数估计
VAR模型参数估计,简化式VARI型比较简单可采用Yule-Walker估计、OLS古计、极大似然估计法等进行估计,且可获得具有良好统计性质的估计量。
结构式VAR莫型参数估计比较复杂,可有两种途径:
一种是
化成简化式,直接估计简化式模型参数,然后再通过简化式模型参数与结构式模型参数的关系,求得结构式模型参数估计,但这存在一个问题是否可行,什么情况下可行,这与结构式模型的识别性有关。
另一种途径是直接对结构式模型参数进行估计,但这也存在一个问题,上述方法不可应用,原因是每一方程含有众多内生的与扰动项相关变量,那么,如何估计?
这也与结构式模型的识别性有关。
对于简化式VARI型(15.1.1)—(),在冲击向量满足假设ut~iidN(0,),t1,2,…,T,即ut相互独立,同服从以E(ut)0为期望向量、Cov(ut)E(utut)为方差协方差阵的k维正态分布。
这时,
ut是k维白噪声向量序列的条件下,模型参数阵A,Ae,...,Ap及也可采用Yule-Walker估计、OLS古计、
极大似然估计。
设yt(y1ty2t...ykt),t1,2,…,t为长度为t的样本向量
1、Yule-Walker估计
在T充分大时,
首先估计自协方差阵
T
?
・h
ytyt
h/T
(7.1.16)
th
1
?
?
?
?
A
0
"1■■■
・p1
・1
1
令?
?
?
・0・・・
?
p2
?
?
2
A
&
・・・
M
M
?
p1
?
p2
...?
)
?
・p
AP
则可得模型参数阵的Yule-Walker估计(矩估计)为
?
■o
?
K?
1
?
?
■1
M
..
?
■p1
2、OLS估计
模型参数阵A,A?
.••,Ap的OLS古计,
即求使
?
1■■■
?
1
■p1
?
?
?
?
■0…
■p2
■2
(7.1.17)
M,
?
■p2
?
0
?
■p
pp
Ajy「)(ytjtj)
j1j1
q(A,A2,…,Ap)7T(yt
Tjp1
min
下的兄局,…,Ap作为Ai,A2,.・・,Ap估计。
T
记
?
ytyt
tp1
h/T
(7.1.18)
由此可推得
A
?
■o
?
"1■■■
?
■p1
1?
■1
A
?
?
?
?
AA2
?
1?
1
■0…
■p2
■2
(7.1.19)
M
...
M,
ap
?
■p1
?
■p2
?
0
?
■p
由此可见,模型参数阵AA,...,Ap的OLS古计(7.1.15)
与Yule-Walker估计形式相同
但式中的彳的计算不同.但是,当T充分大时,(7.1.16)与相差很小,这时与相差也很小,这时二者的估计及估计量的性质等价。
因此,在T充分大时,可直接采用Yule-Walker估计比较简单方便。
1T
而的估计为?
?
)A?
A—utut(7.1.20)
Tt1
其中:
?
yt%t1定yt2...Apytp
3、极大似然估计
可证明,模型参数阵A,A2,...,Ap的极大似然估计与OLS古计完全等价。
除此之外,还有递推估计法(参见:
马树才,《经济时序分析》,辽宁大学出版社,)这里不在赘述。
7.1.3简化式VAR莫型的预测
在已知yt1,yt2,...时,对yt的一步线性预测
(7.1.21)
其一步预测误差为
%yt
?
心)
e
一步预测误差的方差阵为
E%%
Eetet
S的估计为
S?
(1
kP)1
T
p
(?
0&?
)
i1
?
t1
(1)A1Yt1A2yt2...Apytp
在已知yt1,yt2,...时,如果利用模型参数的估计量
(7.1.22)
R,A2,...,Ap,对yt进行一步线性预测,则
yt的实际一步线性预测为
其一步预测误差为
1
(1)
AYt1
A2yt2
Apytp
(7.1.23)
%yt?
t1
(1)
(AA1)Y1(A2
心%2...(Ap/?
p)ytpet
步预测误差的方差阵为E%%Eq©D的估计为
D?
(1kP)(1kP)1(?
oAj?
)(7.1.24)
TTii
7.1.4VAR模型阶数p的确定
VAR模型的定阶是一个矛盾过程,阶数p的确定,既不能太大,又不能太小,必须兼顾。
因为,一方
面,希望滞后阶数p要大一些,以便使模型能更好地反映出动态特征,但另一方面,又不希望太大,否则,阶数p太大,会造成需要估计的模型参数过多,而使模型自由度减少。
因此,在定阶时需要综合考虑,以既要有足够大的滞后项,又能有足够大的自由度为原则确定阶数。
VAR莫型的定阶方法有多种:
1、FPE准则(最小最终预测误差准则)
FPE准则(最小最终预测误差准则),即利用一步预测误差方差进行定阶。
因为,如果模型阶数合适,则模型对实际数据拟合优度必然会高,其一步预测误差方差也必然会小;反之,则相反。
设给定时间序列向量长度为T的样本向量为y(y1ty2t-..ykt),t1.2.....T,则其一步预测误差方差阵的估计量为(7.1.24)式,它是一个kk阶阵,因此可定义其最终预测误差为
FPEk(p)detD?
(1®)k(1kp)kdet(?
0卩A?
?
)(7.1.25)
TTi1
显然,FPEk(p)是p的函数。
所谓最小最终预测误差准则,就是分别取p=1,2,…,M,来计算FPEk(p),使FPEk(p)min值
所对应的p,为模型合适阶数。
相应的模型参数估计A1,A,,...,&p为最佳模型参数估计。
其中,M为预先选定的阶数上界,一般取MT/10k~T/5k之间。
在实际计算过程中,可如下判断:
1如果FPEk(p)的值,随着p从1开始逐渐增大就一直上升,则可判定p=1;
2如果FPEk(p)的值,随着p从1开始逐渐增大就一直下降,则可判定该随机时间序列不能用AR(p)
模型来描述;
3如果FPEk(p)的值,在某一p值下降很快,而后又缓慢下降,则可判定该p值为所确定的阶数;
4如果FPEk(p)的值,随着p从1开始逐渐增大而上下剧烈跳动,难以找到最小值,这可能由于样本数据长度T太小造成的,应增大样本长度,重新进行定阶、估计模型参数,建立模型。
利用FPE言息准则还可以用来检验模型的建立是否可由部分分量,比如前r(rk)个分量y1ty2t...yrt,
t1,2,...,T来进行,方法如下:
pp
记(7.1.21)式中的kk阶矩阵(?
0A?
)的左上角r阶子方阵为代A?
)rr,则前r个分
i1i1
量y1ty2t...yrt,t1,2,...,T的最终预测误差为
kp
FPEr(p)detD?
r(1kP)r(1kP)rdet(
P
?
A?
)
■0i1rr
i1
(7.1.26)
当rk时,(7.1.26)为式。
如果,minFPE「(p)minFPEk(p),则可认为仅用前r个分量y1ty2t...yrt,t1,2,...,T建立模型
即可,没有必要采用k维随机时间序列yt(yity2f-ykt)建立模型,因为从最小最终预测误差准则角度,
1,2,...,T建立模型,
用k维随机时间序列yt(yity2t...ykt)建立模型比仅采前r个分量y2t-..yrt,t
带来拟合优度的显著改善;反之,则相反。
2、AIC(AkaikeInformationCriterion)
与SC(BayesInformationCriterion
)信息准则
AIC、SC言息准则,也称最小信息准则,定义
AIC
2l/T2n/T,SC21/TnlnT/T
(7.1.27)
其中:
ITk(1
2
ln2)
T
ln
2
n为模型需要估计参数个数,对(7.1.1
2
),npk;对于
nk(dpk);对于(
n(P
22
1)k;对于(),nk(dpk)k。
所谓最小信息准则,就是分别取p=1,2,…,来计算aic或者sc,使aic或scmin值所对应的p,为
模型合适阶数。
相应的模型参数估计A局,...,AP为最佳模型参数估计。
3、似然比检验法(
LikelihoodRatio,LR
检验):
由于ut~iidN(0,
)t12T
/?
kI,A-,...,I?
即Ut相互独立,同服从以E(ut)0为期望向量、
Cov(ut)
E(UtUt)
为方差协方差阵的
k维正态分布。
因此,
条件分布为
yt1
yM,A
ytp
AA2LAp,则在给yt1,yt2,...,yp1的条件下,乂(%t丫玄…丫减)的
ytyt1,Yt2,...,yP1~N(AY,)
是,在给yt1,yt2,...,y
p1的条件下,y1,y2,...,y的联合分布密度,
即似然函数为
L(A,)(2
Tk/2
T/2
exp{(
T
[(ytAYt)
1(yt
AYt)]}
对数似然函数为lnL(A,
Tkln(2
2
)—ln
2
T
[(yt
AYt)
1
(ytAYt)]
将参数估计代入,则有
lnL(A,)
Tin
2
(Ut?
1ut
1T
又?
_?
?
Tt1
因此,有
Tkt
lnL(A,)ln
(2)—ln
22
Tk
2
(7.1.28)
现在,欲检验假设H°:
样本数据是由滞后阶数为
p的VAR莫型生成;H1:
样本数据是由滞后阶数为
p1的VAR莫型生成
取似然比统计量为
LR2[lnL(A,?
pjlnL(A,?
p)]T(ln
在给定的显著性水平下,当LR
22
(k),则拒绝
In
22
):
(k)分布
(7.1.29)
H°,表明增加滞后阶数,可显著增大似然函数值;
否则,则相反。
LR检验在小样本下,可取似然比统计量为
LR(Tm)(ln
In
22
):
(k)分布
(7.1.30)
其中,mdkp.
7.1.5VAR模型的Granger因果关系检验
Granger因果关系,这也是建
VAR莫型的另一重要应用是可用来检验一个变量与另一变量间是否存在立VAR莫型所需要的。
1、Granger因果关系的涵义
设yt(y1ty2t)为一2维随机时间序列,如果在给定、血的滞后值下的条件分布与仅在给定的y1t的滞后值下y1t的条件分布相同,即
f(y1ty1t1,y1t2,...皿pM心2,.
..,y2tp)f(y1ty1t1,%2,...,%p)
则称y2t对y1t存在Granger非因果性关系,否则,y2t对y1t存在Granger因果性关系。
y2t的滞后值,并不对只由yn的滞
Granger因果性关系涵义的另一表述:
在其条件不变下,如果加上
后值下对y1t进行预测有显著改善,则称y2t对y1t存在Granger非因果性关系,否则,y2t对y1t存在Granger
因果性关系。
2、Granger因果关系检验
设yt(y1ty2t)为一2维随机时间序列,p为滞后阶数,ut
(U1tU2t)
为一2维随机扰动的时间序
列,则有2元VAR莫型为
(1)
(1)⑵
y1ta11y1t1a12y2t1a11y1t2
(1)
(1)⑵
y2ta21y1t1a22y2t1a21yit2
(2)
a12y2t2
(2)
a22y2t2
a(p)nY1t
a(p)21y2tp
(p)
pa12y2tpU1t
a(p)ya12y2tp
U2t
(7•
1.31)
显然,欲检验y2t对y1t是否存在Granger非因果性关系,等价地,
检验假设H°:
a⑴12a⑵12...a(p)120;H1:
a⑴伐,a⑵仁,...a(p)12中至少有一个不为0。
其用于检验的统计量为
F(P,T
2p1)
(SS0SSR^/p
SSFU/(T2p1)
其中,SSFyi,y2为模型(7.1.31)中第1方程残差平方和,SS&为模型()中第1方程去掉y各期滞后项后拟合残差平方和。
在给定的显著性水平下,当FF(p,T2p1)时,拒绝H。
。
如果模型(7.1.31)满足ut~iidN(0,),t1,2,…,T,即ut相互独立,同服从以E(uJ0为
期望向量、Cov(ut)E(utut)
为方差协方差阵的k维正态分布条件,则
也可采用如下统计量进行检验
2Z(S^~2(p)(…)
在给定的显著性水平下,当22(p)时,拒绝H0,
上述Granger因果性关系检验,可推广到对任意k维VAR莫型以及SVAf模型中的某一或某几个随机时
间序列(包括内生、外生变量)是否对另一时间序列具有Granger因果性的检验上去。
§VAR(p)模型的脉冲响应函数与方差分解
在实际应用中,由于通常所设定的VAR莫型都是非经济理论性的简化式模型,出它无需对变量作任何
先验性约束,因此,在分析应用中,往往并不利用VAf模型去分析某一变量的变化对另一变量的影响如何,
而是分析当某一扰动项发生变化,或者说模型受到某种冲击时,对系统的动态影响,这钟分析方法称为脉
冲响应函数方法(Imp
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- VAR 模型 向量 VECM