计量经济学实验八--李子奈

1实验八联立方程计量经济学模型一实验目的：掌握联立方程模型，理解掌握恰好识别方程与可识别方程的估计方法，熟悉Eviews的基本操作。二实验要求：应用教材P187习题8案例估计模型，采用工具变量法、间接最小二乘法与二阶段最小二乘法进行估计。三实验原理：工具变量法、间接最小二乘法、二阶段最小二乘法。四预备知识：可识别的概念、工具变量法、间接最小二乘法、二阶段最小二乘法。五实验内容：下列为一个完备的联立方程计量经济学模型：01121tttttYMCI0132ttttMYP其中，M为货币供给量，Y为国内生产总值，P为价格总指数。C，I分别为居民消费与投资。以如下中国的实际数据为资料，估计上述联立模型。要求恰好识别的方程按工具变量法与二阶段最小二乘法估计。年份货币与准货币2M(亿元)国内生产总值GDP(亿元)居民消费价格指数P(1978=100)居民消费CONS(亿元)固定资产投资I(亿元)199015293.419347.8165.29450.94517199119349.922577.4170.810730.65594.5199225402.227565.2181.713000.18080.1199334879.836938.1208.516412.113072.3199446923.550217.4258.721844.217042.1199560750.563216.9302.928369.720019.3199676094.974163.6328.133955.922913.5199790995.381658.5337.336921.524941.11998104498.586531.6334.639229.328406.21999119897.991125329.941920.429854.72000134610.498749331.245854.632917.72001158301.9108972.4333.549213.237213.52002185007120350.3330.952571.343499.92003221222.8136398.8334.856834.455566.62004254107160280.4347.963833.570477.42005298755.7188692.1354.271217.588773.62006345603.6221651.3359.580476.9109998.22007403442.2263242.5376.793317.2137323.9六实验步骤：26.1分析联立方程模型。题设模型为：0112101321,2,...tttttttttYMCItnMYP其结构参数矩阵为：1012113110100tttttYMCIPB易知，该结构式模型中内生变量个数为2g，先决变量个数为4k。对于第1个方程，有003()B，00()11RBg，且1111kkg，所以第1个结构方程为恰好识别的结构方程。对于第2个方程，有0012()B，有00()11RBg，且有22211kkg，所以第2个结果方程为过度识别方程。模型的简化式模型为：101112131202122232,1,2,...,ttttttttttYCIPtnMCIP6.2建立工作文件并录入数据，如图1所示。图16.3估计国内生产总值方程301121tttttYMCI6.3.1使用狭义的工具变量法估计国内生产总值方程选取国内生产总值方程中未包含的先决变量*0X作为内生解析变量的工具变量，得到结果参数的工具变量法估计量，利用公式进行估计：10**00000010''BXXYXXXY其中，1tYY，0tYM，*0tXP，0[1]ttXCI。(注意，这里估计的0B，0的含义已不同于上述结构式识别条件中的0B，0。)利用Matlab进行矩阵的计算，其部分过程及结果如下图2所示：图2根据Matlab计算出来的结果得到：1(1)Br，0(2)Br，1(3)Br，2(4)Br，其中(),1,2,3,4Brii为Matlab计算中Br矩阵中的第i个元素。于是得到参数的估计为：0173.585710.04939811.66929720.9407076.3.2使用间接最小二乘法估计国内生产总值方程有6.1的分析有国内生产总值方程中包含的内生变量的简化式方程为：101112131202122232,,ttttttttttYCIPMCIP其参数关系体系为：101200111211121222131230使用普通最小二乘法估计简化式方程，在Eviews中点击主界面菜单Quick\EstimateEquation，在弹出的对话框中输入YCTIP，点击确定，即可4得到第一个简化方程回归结果，如图3所示；同样的，点击主界面菜单Quick\EstimateEquation，在弹出的对话框中输入MCTIP，点击确定，可得到第二个简化方程回归结果，如图4所示。图3图4根据图3中的数据，可以得到：102152.238111.397255120.9267051323.40697根据图4中的数据，可以得到：2040055.10215.507125220.28345723473.8420于是，由参数关系体系计算得到结构参数间接最小二乘估计值为：13231/0.049398102001173.58571121111.6692971222210.9407076.3.3使用二阶段最小二乘法估计国内生产总值方程Step1：用普通最小二乘法估计内生变量的简化式方程，如图4所示，由图中的数据得到：40055.105.5071250.283457473.8420ttttMCIPStep2：据此方程计算tM，替换结果方程中的tM，再用普通最小二乘法估计变换了的结构式方程。点击主菜单中Object\GenerateSeries...，在弹出的对话框中输入：40055.105.507125*0.283457*473.8420MMtCTIP，产生序列MMt。点击主界面菜单Quick\EstimateEquation，在弹出的对话框中输入YCMMtCTI，点击确定即可得到回归结果，如图5所示。5图5由图5中数据，得到国内生产总值方程的二阶段最小二乘估计量为：0173.585710.04939811.66929720.940707比较上述国内生产总值方程的3种估计结果，说明这3种方法对于恰好识别的结构方程是等价的。6.4估计货币供给量方程0132ttttMYP由6.1的分析知，货币供给量方程为过度识别的结构方程，最能用二阶段最小二乘法进行参数的估计。同样得，仿照6.3.3的步骤有：Step1：用普通最小二乘法估计内生变量的简化式方程，如图3所示，由图中的数据得到：2152.2381.3972550.92670523.40697ttttYCIPStep2：据此方程计算tY，替换结果方程中的tY，再用普通最小二乘法估计变换了的结构式方程。点击主菜单中Object\GenerateSeries...，在弹出的对话框中输入：2152.2381.397255*0.926705*23.40697*YYtCTIP，产生序列YYt。点击主界面菜单Quick\EstimateEquation，在弹出的对话框中输入MCYYtP，点击确定即可得到回归结果，如图6所示。6图6由图6中的数据得到货币供给量方程的二阶段最小二乘估计量为：01986.22311.8096113146.9317至此，我们已经完成了该模型系统的估计，并完成了题目的要求。6.5模型的直接计算机估计EViews软件拥有强大的功能，在Eviews软件中有直接应用的二阶段最小二乘估计的函数，即我们可以不用自己一步步地进行二阶段最小二乘估计而直接输入结构方程要估计的参数跟结构式模型中的可决变量，Eviews软件就可以帮我们实现二阶段最小二乘法的估计。在Eviews主菜单中点击Quick\EstimateEquation，会打开方程定义对话框，在对话框项的下拉列表框中选择选项TSLS，屏幕会出现如图7所示的二阶段最小二乘估计对话框。图7图87TSLS估计对话框的设置。首先，在对话框上方EquationSpecification项编辑框中设定待估计方程的形式，对于第1个结构方程(国内生成总值方程)使用列表法输入：YCMCTI。然后再对对话框Instrumentlist项编辑框中输入结构式模型中的可决变量名，本例中输入：CCTIP。设定对话框后，点击确定，屏幕将输出如图8所示的方程2SLS的估计结果。根据图8中的数据，得到国内生成总值方程的2SLS估计为：173.58570.0493981.6692970.940707(0.190)(2.051)(24.373)(21.526)ttttYMCI20.999830R20.999794R..1.554243DW图9图10同样的，对于货币供给量方程有同样的操作，在TSLS估计对话框上方的EquationSpecification项编辑框中输入：MCYP，在Instrumentlist项编辑框中输入：CCTIP，如图9所示。点击确定即可得到如图10所示的方程2SLS的估计结果。根据图10中数据，得到货币供给量方程的2SLS估计为：1986.2221.809611146.9316(0.137)(28.668)(2.338)tttMYP20.992274R20.991244R..0.389408DW可以看到模型的直接计算机估计，跟自己分两步做的估计结果是一致的，估计量的差别只是很小的计算误差。