09联立方程.docx
- 文档编号:8195377
- 上传时间:2023-01-29
- 格式:DOCX
- 页数:21
- 大小:407KB
09联立方程.docx
《09联立方程.docx》由会员分享,可在线阅读,更多相关《09联立方程.docx(21页珍藏版)》请在冰豆网上搜索。
09联立方程
1.5.8联立方程模型(simultaneous-equationsmodel)
1.联立方程模型的概念
有时由于两个变量之间存在双向因果关系,用单一方程模型就不能完整的描述这两个变量之间的关系。
有时为全面描述一项经济活动只用单一方程模型是不够的。
这时应该用多个方程的组合来描述整个经济活动。
从而引出联立方程模型的概念。
联立方程模型:
对于实际经济问题,描述变量间联立依存性的方程体系。
联立方程模型的最大问题是E(X'u)0,当用OLS法估计模型中的方程参数时会产生联立方程偏倚,即所得参数的OLS估计量
是有偏的、不一致的。
给出三个定义:
内生变量(endogenousvariable):
由模型内变量所决定的变量。
外生变量(exogenousvariable):
由模型外变量所决定的变量。
前定变量(predeterminedvariable):
包括外生变量、外生滞后变量、内生滞后变量。
例如:
yt=0+1yt-1+0xt+1xt-1+ut
yt为内生变量;xt为外生变量;yt-1,xt,xt-1为前定变量。
联立方程模型必须是完整的。
所谓完整即“方程个数内生变量个数”。
否则联立方程模型是无法估计的。
2.联立方程模型的分类(结构模型,简化型模型,递归模型)
⑴结构模型(structuralmodel):
把内生变量表述为其他内生变量、前定变量与随机误差项的方程体系。
例:
如下凯恩斯模型(为简化问题,对数据进行中心化处理,从而不出现截距项)
ct=1yt+ut1消费函数,行为方程(behaviorequation)
It=1yt+2yt-1+ut2投资函数,行为方程
yt=ct+It+Gt国民收入等式,定义方程(definitionalequation)
(1)
其中,ct消费;yt国民收入;It投资;Gt政府支出。
1,1,2称为结构参数。
模型中内生变量有三个ct,yt,It。
外生变量有一个Gt。
内生滞后变量有一个yt-1。
Gt,yt-1又称为前定变量。
因模型中包括三个内生变量,含有三个方程,所以是一个完整的联立模型。
内生变量与外生变量的划分不是绝对的,随着新的行为方程的加入,外生变量可以转化为内生变量;随着行为方程的减少,内生变量也可以转化为外生变量。
⑵简化型模型(reduced-formequations):
把内生变量只表示为前定变量与随机误差项函数的联立模型。
仍以凯恩斯模型为例其简化型模型为,
ct=11yt-1+12Gt+vt1
It=21yt-1+22Gt+vt2
yt=31yt-1+32Gt+vt3
(2)
或
=
+
,
其中ct,yt,It为内生变量,yt-1,Gt为前定变量,ij,(i=1,2,3,j=1,2),为简化型参数。
用如下矩阵符号表示上式
Y=X+v(3)
显然结构模型参数与简化型模型参数之间存在函数关系。
把结构模型
(1)中的内生变量全部移到方程等式的左边得
ct-1yt=ut1
It-1yt=2yt-1+ut2
-ct-It+yt=Gt(4)
用矩阵形式表达
=
+
用如下矩阵符号表示上式
Y=X+u(5)
则
Y=-1X+-1u(6)
比较联立方程模型(3)和(6),结构参数和简化型参数有如下关系存在,
=-1
=
=
其中,A-1=
。
A=
=
。
adj(A)=
=
。
的伴随矩阵是的代数余子式组成的矩阵的转置。
v=-1u
=
⑶递归模型(recursivesystem):
在结构方程体系中每个内生变量只是前定变量和比其序号低的内生变量的函数。
y1=11x1+…+1kxk+u1
y2=21x1+…+2kxk+21y1+u2
y3=31x1+…+3kxk+31y1+32y2+u3
…..
ym=m1x1+…+mkxk+m1y1+m2y1+…+mm-1ym-1+um(7)
其中yi和xj分别表示内生变量和外生变量。
其随机误差项应满足
E(u1u2)=E(u1u3)=…=E(u2u3)=…=E(um-1um)=0
3.联立方程模型的识别(identification)
例:
关于粮食的需求供给模型如下,
Dt=0+1Pt+u1(需求函数)
St=0+1Pt+u2(供给函数)
St=Dt(平衡条件)(8)
其中Dt需求量,St供给量,Pt价格,ui,(i=1,2)随机项。
当供给与需求在市场上达到平衡时,Dt=St=Qt(产量),当用收集到的Qt,Pt样本值,而无其他信息估计回归参数时,则无法区别估计值是对0,1的估计还是对0,1的估计。
从而引出联立方程模型的识别问题。
也许有人认为若样本显示的是负斜率,则为需求函数;若是正斜率,则为供给函数。
其实样本点所代表的只是不同需求与供给曲线的交点而已。
显然为区别需求与供给曲线应进一步获得其他信息。
例如收入和偏好的变化会影响需求曲线随时间变化产生位移,而对供给曲线不会产生什么影响。
所以这些观测点就会描绘出供给曲线的位置。
也就是说供给曲线是可识别的。
同理耕种面积、气候条件等因素只会影响供给曲线,不会对需求曲线产生影响。
需求曲线就是可识别的。
可见一个方程的可识别性取决于它是否排除了联立模型中其他方程所包含的一个或几个变量。
QtQt
需求曲线
供给曲线供给曲线,耕地面积不同
需求曲线,收入水平不同
PtPt
在模型(8)的需求函数和供给函数中分别加入收入变量It和滞后价格变量Pt-1,
Dt=0+1Pt+2It+u1(需求函数)
St=0+1Pt+2Pt-1+u2(供给函数)
St=Dt(平衡条件)
于是行为方程成为可识别方程。
从代数意义上讲,当与上述结构方程参数相与对应的简化型方程参数有一一对应关系时,结构模型是恰好识别的。
举例说明。
上模型写为,
Qt=0+1Pt+2It+u1
Qt=0+1Pt+2Pt-1+u2
有6个结构参数。
相应简化型模型为
Qt=10+11It+12Pt-1+vt1
Pt=20+21It+22Pt-1+vt2
当简化型参数多于结构参数时,结构模型是过度识别的。
当简化型参数少于结构参数时,结构模型是不可识别的。
识别问题是完整的联立方程模型所特有的问题。
只有行为方程才存在识别问题,对于定义方程或恒等式不存在识别问题。
识别问题不是参数估计问题,但是估计的前提。
不可识别的模型则不可估计。
识别依赖于对联立方程模型中每个方程的识别。
若有一个方程是不可识别的,则整个联立方程模型是不可识别的。
可识别性分为恰好识别和过度识别。
不可识别
模型的识别恰好识别
可识别
过度识别
识别方法:
1阶条件(ordercondition)
不包含在待识别方程中的变量(被斥变量)个数(联立方程模型中的方程个数–1)
阶条件是必要条件但不充分,即不满足阶条件是不可识别的,但满足了阶条件也不一定是可识别的。
2秩条件(rankcondition)
待识别方程的被斥变量系数矩阵的秩=(联立方程模型中方程个数–1)
秩条件是充分必要条件。
满足秩条件能保证联立方程模型内每个方程都有别于其他方程。
识别的一般过程是
(1)先考查阶条件,因为阶条件比秩条件判别起来简单。
若不满足阶条件,识别到此为止。
说明待识别方程不可识别。
若满足阶条件,则进一步检查秩条件。
(2)若不满足秩条件,说明待识别方程不可识别。
若满足秩条件,说明待识别方程可识别,但不能判别可识别方程是属于恰好识别还是过度识别。
对此还要返回来利用阶条件作判断。
(3)若阶条件中的等式(被斥变量个数=方程个数–1)成立,则方程为恰好识别;若阶条件中的不等式(被斥变量个数>方程个数–1)成立,则方程为过度识别。
例:
某结构模型为,
y1=12y2+11x1+12x2+u1(恰好识别)
y2=23y3+23x3+u2(过度识别)
y3=31y1+32y2+33x3+u3(不可识别)(9)
试考查第二个方程的可识性。
由于结构模型有3个方程,3个内生变量,所以是完整的联立方程模型。
对于第2个方程,被斥变量有3个y1,x1,x2,(方程个数–1)=2。
所以满足阶条件。
结构模型的系数矩阵是,
(10)
从系数阵中划掉第2个方程的变量y2,y3,x3的系数所在的相应行和列,得第2个方程被斥变量的系数阵如下,
(11)
因为
0,
0,(12)
被斥变量系数阵的秩=2,已知(方程个数)-1=2,所以第2个方程是可识别的。
下面用阶条件判断第2个方程的恰好识别性或过度识别性。
因为被斥变量个数是3>2,所以第2个方程是过度识别的。
现考查第3个方程的可识性。
对于第3个方程,被斥变量有2个x1,x2,(方程个数–1)=2。
所以满足阶条件。
从系数阵中划掉第3个方程的变量y1,y2,y3,x3的系数所在的相应行和列,得第3个方程的被斥变量系数阵如下
因为
=0
被斥变量系数阵的秩=1,已知(方程个数)-1=2,所以第3个方程是不可识别的。
4.联立方程模型的估计方法
y1=11x1+…+1kxk+u1
y2=21x1+…+2kxk+21y1+u2
y3=31x1+…+3kxk+31y1+32y2+u3
…..
递归模型的估计方法是OLS法。
解释如下。
首先看第一个方程。
由于等号右边只含有外生变量和随机项,外生变量和随机项不相关,符合假定条件,所以可用OLS法估计参数。
对于第二个方程,由于等号右边只含有一个内生变量y1,以及外生变量和随机项。
根据假定u1和u2不相关,所以y1和u2不相关。
对于y2来说,y1是一个前定变量。
因此可以用OLS法估计第2个方程。
以此类推可以用OLS法估计递归模型中的每一个方程。
参数估计量具有无偏性和一致性。
简化型模型可用OLS法估计参数。
由于简化型模型一般是由结构模型对应而来,每个方程只含有一个内生变量且为被解释变量。
它是前定变量和随机项的唯一函数。
方程中解释变量都是前定变量,自然与随机项无关。
所以用OLS法得到的参数估计量为一致估计量。
对于结构模型有两种估计方法。
一种为单一方程估计法,即有限信息估计法;另一种为方程组估计法,系统估计法,即完全信息估计法。
前者只考虑被估计方程的参数约束问题,而不过多地考虑方程组中其他方程所施加的参数约束,因此称为有限信息估计方法。
后者在估计模型中的所有方程的同时,要考虑由于略去或缺少某些变量而对每个方程所施加的参数约束。
因此称为完全信息估计法。
显然对于联立方程模型,理想的估计方法应当是完全信息估计法,例如完全信息极大似然法(FIML)。
然而这种方法并不常用。
因为①这种方法计算工作量太大,②将导致在高度非线性的情况下确定问题的解,这常常是很困难的,③若模型中某个方程存在设定误差,这种误差将传播到其他方程中去。
所以对于联立方程模型常用的估计方法是单一方程估计法。
常用的单一方程估计法有①间接最小二乘法(ILS),②工具变量法(IV),③两段最小二乘法(2SLS),④有限信息极大似然法(LIML)。
工具变量法与2SLS法一起介绍。
有限信息极大似然法放在第2章介绍。
ILS法只适用于恰好识别模型。
具体估计步骤是先写出与结构模型相对应的简化型模型,然后利用OLS法估计简化型模型参数。
因为简化型模型参数与结构模型参数存在一一对应关系,利用=-1可得到结构参数的唯一估计值。
ILS估计量是有偏的,但具有一致性和渐近有效性。
当结构方程为过度识别时,其相应简化型方程参数的OLS估计量是有偏的,不一致的。
采用ILS法时,简化型模型的随机项必须满足OLS法的假定条件。
viN(0,2),cov(vi,vj)=0,cov(xi,vj)=0。
当不满足上述条件时,简化型参数的估计误差就会传播到结构参数中去。
2SLS法。
对于恰好识别和过度识别的结构模型可采用2SLS法估计参数。
2SLS法即连续两次使用OLS法。
使用2SLS法的前提是结构模型中的随机项和简化型模型中的随机项必须满足通常的假定条件,前定变量之间不存在多重共线性。
以如下模型为例作具体说明。
y1=1y2+1x1+u1(13)
y2=2y1+2x2+u2(14)
其中uiN(0,i2),i=1,2;plimT-1(xiuj)=0,(i,j=1,2);E(u1u2)=0。
第一步,作如下回归,
y2=
x1+
x2+
(15)
因为
=
x1+
x2是x1和x2的线性组合,而x1,x2与u1,u2无关,所以
也与u1,u2无关。
是y2的OLS估计量,自然与y2高度相关。
所以可用
作为y2的工具变量。
第二步,用
代替方程(13)中的y2,得
y1=1
+1x1+u1
用OLS法估计上式。
定义W=(
x1),则
=(W'W)-1(W'y1)
为2SLS估计量。
是有偏的、无效的、一致估计量。
可以证明当结构模型为恰好识别时,2SLS估计值与ILS估计值相同。
例1:
河南省国民收入计量模型(1952-1982年数据,递归模型,OLS法估计参数)
⑴Y1=-21.0982+0.0486X1+0.033X4+20.5486D1(农业生产函数)
(7.63)(9.99)(9.04)R2=0.9845,F=572.9,DW=2.20
⑵LnY2=0.0876+0.2184LnX2+0.6545LnX5+0.3503D2(重工业生产函数)
(1.54)(5.19)(2.45)R2=0.8165,F=38.54,DW=1.27
⑶LnY3=0.5946+0.3728LnX3+0.7798LnX6(轻工业生产函数)
(5.10)(6.86)R2=0.7939,F=51.98,DW=2.12
⑷Y4=Y2+Y3(定义方程)
⑸Y5=2.1586+0.4271Y1+0.5854Y4+16.8646D3(国民收入函数)
(4.34)(10.37)(5.26)R2=0.9874,F=709.1,DW=1.34
变量定义:
Y1,农业总产值(亿元)
X1,农业劳动力人数(万人)
Y2,重工业总产值(亿元)
X2,重工业劳动力人数(万人)
Y3,轻工业总产值(亿元)
X3,轻工业劳动力人数(万人)
Y4,工业总产值(亿元)
X4,农机总动力(万马力)
Y5,国民收入(亿元)
X5,重工业固定资产原值(亿元)
X6,轻工业固定资产原值(亿元)
D1,D2,D3,虚拟变量(区别经济困难时期)
(1)在河南省国民收入计量模型中若删去1号方程,则Y1变为外生变量。
(2)若在模型中加入方程
X4=f(可灌溉亩数,农机台数,副业产值),
则X4由外生变量转化为内生变量。
(3)若在5号方程中加入交通运输业变量Y6,则Y6为外生变量。
若加入方程
Y6=f(货运量,铁路运营公里数,公路运营公里数),
则Y6由外生变量转化为内生变量。
例2:
美国电力需求模型
(摘自ReviewofEconometricsandStatisticsVol.57,p12-18,1975)
电销量,电边际价格,人均年收入,天然气价格,取暖天数,7月平均气温,农村人口比率,家庭人口
LnQ=-0.21-1.15LnP+0.51LnY+0.04LnG-0.02LnD+0.54LnJ+0.21LnR-0.24LnH
(-38.3)(8.5)(4.0)(1.0)(4.5)(10.5)(2.0)R2=0.91
电边际价格,电销量,劳动力成本,上市发电比率,电成本,农村人口比率,工民电销比,时间
LnP=-0.57-0.60LnQ+0.24LnL-0.02LnK+0.01LnF+0.03LnR-0.12LnI+0.004LnT
(-20.0)(6.0)(2.0)(3.3)(3.0)(12.0)(1.3)R2=0.97
其中,
Q:
民用电年平均销售量。
P:
民用电边际价格。
Y:
人均年收入。
G:
民用天然气价格。
D:
取暖天数。
J:
:
7月份平均气温。
R:
农村人口比率。
H:
平均家庭人口。
L:
劳动力成本。
K:
上市电力企业发电比重。
F:
每度电平均成本。
I:
工业用电与民用电销量比。
T:
时间。
上模型中内生变量是Q和P。
并互做解释变量。
因为每个方程中各有5个区别于另外方程的外生变量,所以上模型为过度识别模型。
2SLS估计的步骤是
(1)用模型中每个内生变量对模型中全部外生变量进行最小二乘回归,
(2)用得到的Q和P的估计值替代结构方程右侧的相应内生变量,并进行最小二乘估计,从而得到上述结果。
用的是1961-1969年美国48个州的时序与截面混合数据。
实际分析:
从第一个方程看,与电销售量对其他变量的弹性系数值相比,只有电销量的价格弹性系数值(绝对值)最大。
这说明近年来,居民用电量的增长主要是因为电价下降的结果。
例3:
中国宏观经济的联立方程模型(file:
simu4)
消费方程:
Ct=0+1Yt+2Ct-1+u1t
投资方程:
It=0+1Yt-1+u2t
收入方程;Yt=Ct+It+Gt
其中:
Ct消费;Yt国民收入;It投资;Gt政府支出
用中国1978-2000数据得估计结果如下,
估计结果表达式是,
消费方程:
Ct=362.0544+0.3618Yt+0.2467Ct-1+
(3.5)(17.0)(4.9)R2=0.9995
投资方程:
It=625.9373+0.4095Yt-1+
(1.0)(26.0)R2=0.9713
收入方程;Yt=Ct+It+Gt
联立方程模型的两段最小二乘估计(Eviews)
在打开工作文件窗口的基础上,点击主功能菜单上的Objects键,选NewObject功能,
从而打开NewObject(新对象)选择窗。
选择System,并在NameofObject处为联立方程模型起名(图中显示为Untitled)。
然后点击OK键。
从而打开System(系统)窗口。
在System(系统)窗口中键入联立方程模型。
消费方程:
Ct=0+1Yt+2Ct-1+u1t
投资方程:
It=0+1Yt-1+u2t
收入方程;Yt=Ct+It+Gt
在Eviews命令中若用Cons表示Ct,用gdp表示Yt,用Inv表示It,用Gov表示Gt,选Ct-1,Yt-1,Gt为工具变量,则模型的Eviews表达是
把如上的方程式键入System(系统)窗口如下图。
点击系统(System)窗口上的估计(estimation)键,立刻弹出系统估计方法窗口(见下图)。
共有9种估计方法可供选择。
他们是OLS,WLS,SUR(SeeminglyUnrelatedRegression),2SLS,WTSLS,3SLS,FIML,GMM(White协方差矩阵,用于截面数据),GMM(HAC协方差矩阵,用于时间序列数据)。
选择2SLS估计,点击”OK”键,得估计结果如前。
附数据如下:
obs
GDP
CONS
INV
GOV
1978
3605.6
1759.1
1377.9
480
1979
4074
2005.4
1474.2
614
1980
4551.3
2317.1
1590
659
1981
4901.4
2604.1
1581
705
1982
5489.2
2867.9
1760.2
770
1983
6076.3
3182.5
2005
838
1984
7164.4
3674.5
2468.6
1020
1985
8792.1
4589
3386
1184
1986
10132.8
5175
3846
1367
1987
11784.7
5961.2
4322
1490
1988
14704
7633.1
5495
1727
1989
16466
8523.5
6095
2033
1990
18319.5
9113.2
5444
2252
1991
21280.4
10315.9
7617
2830
1992
25863.7
12459.8
9636
3492.3
1993
34500.7
15682.4
14998
4499.7
1994
46690.7
20809.8
19260.6
5986.2
1995
58510.5
26944.5
23877
6690.5
1996
68330.4
32152.3
26867.2
7851.6
1997
74894.2
34854.6
28457.6
8724.8
1998
79003.3
36921.1
29545.9
9484.8
1999
82673.1
39334.4
30701.6
10388.3
2000
89112.5
42911.9
32255
11705.3
例4:
1999年度中国宏观经济计量模型框图(原书1~56页)
原始资料来源:
《中国社会科学院数量经济与技术经济研究所经济模型集》,汪同三、沈利生
主编,社会科学文献出版社,2001,第4页。
本人有修改。
1999年度中国宏观经济计量模型分为8个模块(蓝色区域),共174个方程。
含174个内生变量,37个外生变量。
其中
1.
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 09 联立方程