1、2002CP-AH(安徽)3282.4663646.1503777.4103989.5814203.5554495.1744784.364CP-BJ(北京)5133.9786203.0486807.4517453.7578206.2718654.43310473.12CP-FJ(福建)4011.7754853.4415197.0415314.5215522.7626094.3366665.005CP-HB(河北)3197.3393868.3193896.7784104.2814361.5554457.4635120.485CP-HLJ(黑龙江)2904.6873077.9893289.9903
2、596.8393890.5804159.0874493.535CP-JL(吉林)2833.3213286.4323477.5603736.4084077.9614281.5604998.874CP-JS(江苏)3712.2604457.7884918.9445076.9105317.8625488.8296091.331CP-JX(江西)2714.1243136.8733234.4653531.7753612.7223914.0804544.775CP-LN(辽宁)3237.2753608.0603918.1674046.5824360.4204654.4205402.063CP-NMG(内蒙
3、古)2572.3422901.7223127.6333475.9423877.3454170.5964850.180CP-SD(山东)3440.6843930.5744168.9744546.8785011.9765159.5385635.770CP-SH(上海)6193.3336634.1836866.4108125.8038651.8939336.10010411.94CP-SX(山西)2813.3363131.6293314.0973507.0083793.9084131.2734787.561CP-TJ(天津)4293.2205047.6725498.5035916.6136145.6
4、226904.3687220.843CP-ZJ(浙江)5342.2346002.0826236.6406600.7496950.7137968.3278792.210资料来源:中国统计年鉴1997-2003。表2 1999-2002年中国东北、华北、华东15个省级地区的居民家庭人均收入数据(不变价格)地区人均收入IP-AH(安徽)4106.2514540.2474770.4705178.5285256.7535640.5976093.333IP-BJ(北京)6569.9017419.9058273.4189127.9929999.70011229.6612692.38IP-FJ(福建)4884
5、.7316040.9446505.1456922.1097279.3938422.5739235.538IP-HB(河北)4148.2824790.9865167.3175468.9405678.1955955.0456747.152IP-HLJ(黑龙江)3518.4973918.3144251.4944747.0454997.8435382.8086143.565IP-JL(吉林)3549.9354041.0614240.5654571.4394878.2965271.9256291.618IP-JS(江苏)4744.5475668.8306054.1756624.3166793.43773
6、16.5678243.589IP-JX(江西)3487.2693991.4904209.3274787.6065088.3155533.6886329.311IP-LN(辽宁)3899.1944382.2504649.7894968.1645363.1535797.0106597.088IP-NMG(内蒙古)3189.4143774.8044383.7064780.0905063.2285502.8736038.922IP-SD(山东)4461.9345049.4075412.5555849.9096477.0166975.5217668.036IP-SH(上海)7489.4518209.03
7、78773.10010770.0911432.2012883.4613183.88IP-SX(山西)3431.5943869.9524156.9274360.0504546.7855401.8546335.732IP-TJ(天津)5474.9636409.6907146.2717734.9148173.1938852.4709375.060IP-ZJ(浙江)6446.5157158.2887860.3418530.3149187.28710485.6411822.00 图2 15个省级地区的人均消费序列(纵剖面) 图3 15个省级地区的人均收入序列(file:4panel02) 图4 15个省
8、级地区的人均消费散点图 图5 15个省级地区的人均收入散点图(7个横截面叠加)(每条连线表示同一年度15个地区的消费值) (每条连线表示同一年度15个地区的收入值)用CP表示消费,IP表示收入。AH, BJ, FJ, HB, HLJ, JL, JS, JX, LN, NMG, SD, SH, SX, TJ, ZJ分别表示安徽省、北京市、福建省、河北省、黑龙江省、吉林省、江苏省、江西省、辽宁省、内蒙古自治区、山东省、上海市、山西省、天津市、浙江省。15个地区7年人均消费对收入的面板数据散点图见图6和图7。图6中每一种符号代表一个省级地区的7个观测点组成的时间序列。相当于观察15个时间序列。图7中
9、每一种符号代表一个年度的截面散点图(共7个截面)。相当于观察7个截面散点图的叠加。图6 用15个时间序列表示的人均消费对收入的面板数据图7 用7个截面表示的人均消费对收入的面板数据(7个截面叠加) 为了观察得更清楚一些,图8给出北京和内蒙古1996-2002年消费对收入散点图。从图中可以看出,无论是从收入还是从消费看内蒙古的水平都低于北京市。内蒙古2002年的收入与消费规模还不如北京市1996年的大。图9给出该15个省级地区1996和2002年的消费对收入散点图。可见6年之后15个地区的消费和收入都有了相应的提高。图8 北京和内蒙古1996-2002年消费对收入时序图 图9 1996和2002
10、年15个地区的消费对收入散点图2面板数据的估计。用面板数据建立的模型通常有3种。即混合估计模型、固定效应模型和随机效应模型。2.1 混合估计模型。如果从时间上看,不同个体之间不存在显著性差异;从截面上看,不同截面之间也不存在显著性差异,那么就可以直接把面板数据混合在一起用普通最小二乘法(OLS)估计参数。如果从时间和截面看模型截距都不为零,且是一个相同的常数,以二变量模型为例,则建立如下模型, yit = a +b1 xit +eit, t = 1, 2, , T (1)a 和b1不随i,t变化。称模型(1)为混合估计模型。以例1中15个地区1996和2002年数据建立关于消费的混合估计模型,
11、得结果如下:图10 EViwes估计方法:在打开工作文件窗口的基础上,点击主功能菜单中的Objects键,选New Object功能,从而打开New Object(新对象)选择窗。在Type of Object选择区选择Pool(混合数据库),点击OK键,从而打开Pool(混合数据)窗口。在窗口中输入15个地区标识AH(安徽)、BJ(北京)、ZJ(浙江)。工具栏中点击Sheet键,从而打开Series List(列写序列名)窗口,定义变量CP?和IP?,点击OK键,Pool(混合或合并数据库)窗口显示面板数据。在Pool窗口的工具栏中点击Estimate键,打开Pooled Estimatio
12、n(混合估计)窗口如下图。图11在Dependent Variable(相依变量)选择窗填入CP?;在Common coefficients(系数相同)选择窗填入IP?Cross section specific coefficients(截面系数不同)选择窗保持空白;在Intercept(截距项)选择窗点击Common;在Weighting(权数)选择窗点击No weighting。点击Pooled Estimation(混合估计)窗口中的OK键。得输出结果如图10。相应表达式是= 129.6313 +0.7587 IPit (2.0) (79.7) R2 = 0.98, SSEr = 48
13、24588, t0.05 (103) = 1.9915个省级地区的人均支出平均占收入的76%。如果从时间和截面上看模型截距都为零,就可以建立不含截距项的(a = 0)的混合估计模型。以二变量模型为例,建立混合估计模型如下,yit = b1 xit +eit, (2)对于本例,因为上式中的截距项有显著性(t = 2.0 t0.05 (103) = 1.99),所以建立截距项为零的混合估计模型是不合适的。在Pooled Estimation(混合估计)对话框中Intercept(截距项)选择窗中选None,其余选项同上。2.2 固定效应模型。在面板数据散点图中,如果对于不同的截面或不同的时间序列,
14、模型的截距是不同的,则可以采用在模型中加虚拟变量的方法估计回归参数,称此种模型为固定效应模型(fixed effects regression model)。固定效应模型分为3种类型,即个体固定效应模型(entity fixed effects regression model)、时刻固定效应模型(time fixed effects regression model)和时刻个体固定效应模型(time and entity fixed effects regression model)。下面分别介绍。(1)个体固定效应模型。个体固定效应模型就是对于不同的个体有不同截距的模型。如果对于不同的时间
15、序列(个体)截距是不同的,但是对于不同的横截面,模型的截距没有显著性变化,那么就应该建立个体固定效应模型,表示如下, yit = b1 xit +g1 W1 + g2 W2 + +gN WN +eit, (3)其中Wi = eit, i = 1, 2, , N; t = 1, 2, , T,表示随机误差项。yit, xit, i = 1, 2, , N; t = 1, 2, , T分别表示被解释变量和解释变量。模型(3)或者表示为 y1t = g1 +b1 x1t +e1t, i = 1(对于第1个个体,或时间序列),t = 1, 2, , T y2t = g2 +b1 x2t +e2 t,
16、i = 2(对于第2个个体,或时间序列),t = 1, 2, , T yN t = gN +b1 xN t +e N t, i = N(对于第N个个体,或时间序列),t = 1, 2, , T写成矩阵形式,y1 = (1 x1) +e1 = g1 + x1 b +e1yN = (1 xN) +eN = gN + xN b +eN上式中yi,gi,ei,xi都是N1阶列向量。b为标量。当模型中含有k个解释变量时,b为k进一步写成矩阵形式,= + b + 上式中的元素1,0都是T 面板数据模型用OLS方法估计时应满足如下5个假定条件:(1)E(eit|xi1, xi2, , xiT, ai) =
17、0。以xi1, xi2, , xiT, ai为条件的eit的期望等于零。(2)(xi1, xi2, , xiT), ( yi1, yi2, , yiT), i = 1, 2, , N分别来自于同一个联合分布总体,并相互独立。(3)(xit, eit)具有非零的有限值4阶矩。(4)解释变量之间不存在完全共线性。(5)Cov(eit eis|xit,xis, ai) = 0, t s。在固定效应模型中随机误差项eit在时间上是非自相关的。其中xit代表一个或多个解释变量。对模型(1)进行OLS估计,全部参数估计量都是无偏的和一致的。模型的自由度是N T 1N。 当模型含有k个解释变量,且N很大,相
18、对较小时,因为模型中含有k + N个被估参数,一般软件执行OLS运算很困难。在计量经济学软件中是采用一种特殊处理方式进行OLS估计。估计原理是,先用每个变量减其组内均值,把数据中心化(entity-demeaned),然后用变换的数据先估计个体固定效应模型的回归系数(不包括截距项),然后利用组内均值等式计算截距项。这种方法计算起来速度快。具体分3步如下。(1)首先把变量中心化(entity-demeaned)。仍以单解释变量模型(3)为例,则有 = gi + b1 + , i = 1, 2, , N (4)其中 = , = , = , i = 1, 2, , N。公式(1)、(4)相减得,(y
19、it - ) = b1(xit - ) + (eit - )(5)令(yit - ) = ,(xit - ) = ,(eit - ) = ,上式写为= b1 + (6)用OLS法估计(1)、(6)式中的b1,结果是一样的,但是用(6)式估计,可以减少被估参数个数。(2)用OLS法估计回归参数(不包括截距项,即固定效应)。 在k个解释变量条件下,把 用向量形式 表示,则利用中心化数据,按OLS法估计公式计算个体固定效应模型中回归参数估计量的方差协方差矩阵估计式如下,( ) = ( )-1(7)其中 = , 是相对于 的残差向量。(3)计算回归模型截距项,即固定效应参数gi。= - (8)以例1(
20、file:panel02)为例得到的个体固定效应模型估计结果如下:个体固定效应模型的EViwes输出结果中没有公共截距项。图12在EViwes的Pooled Estimation对话框中Intercept选项中选Fixed effects。其余选项同上。(1)个体固定效应模型的EViwes输出结果中没有公共截距项。(2)EViwes输出结果中没有给出描述个体效应的截距项相应的标准差和t值。不认为截距项是模型中的重要参数。(3)当对个体固定效应模型选择加权估计时,输出结果将给出加权估计和非加权估计两种统计量评价结果。(4)输出结果的联立方程组形式可以通过点击View选Representation
21、s功能获得。(5)点击View选Wald Coefficient Tests功能可以对模型的斜率进行Wald检验。(6)点击View选Residuals/Table, Graphs, Covariance Matrix, Correlation Matrix功能可以分别得到按个体计算的残差序列表,残差序列图,残差序列的方差协方差矩阵,残差序列的相关系数矩阵。(7)点击Procs选Make Model功能,将会出现估计结果的联立方程形式,进一步点击Solve键,在随后出现的对话框中可以进行动态和静态预测。输出结果的方程形式是 = 安徽+ x1t = 479.3 + 0.70 x1t (55.0)
22、 = 北京+ x2t = 1053.2 + 0.70 x2t (55.0) = 浙江+ x15t = 714.2 + 0.70 x15t (55.0)R2 = 0.99, SSEr = 2270386, t0.05 (88) = 1.98从结果看,北京、上海、浙江是消费函数截距(自发消费)最大的3个地区。相对于混合估计模型来说,是否有必要建立个体固定效应模型可以通过F检验来完成。原假设H0:不同个体的模型截距项相同(建立混合估计模型)。备择假设H1:不同个体的模型截距项不同(建立个体固定效应模型)。F统计量定义为:F= = (9)其中SSEr,SSEu分别表示约束模型(混合估计模型)和非约束模
23、型(个体固定效应模型)的残差平方和。非约束模型比约束模型多了N-1个被估参数。(混合估计模型给出公共截距项。)当模型中含有k个解释变量时,F统计量的分母自由度是NT-N-k。 用上例计算,已知SSEr = 4824588,SSEu = 2270386,F= = = = 7.15F0.05(14, 89) = 1.81因为F= 7.15 F0.05(14, 89) = 1.81,所以,拒绝原假设。结论是应该建立个体固定效应模型。(2)时刻固定效应模型。时刻固定效应模型就是对于不同的截面(时刻点)有不同截距的模型。如果确知对于不同的截面,模型的截距显著不同,但是对于不同的时间序列(个体)截距是相同的,那么应该建立时刻固定效应模型,表示如下, yit = b1 xit +a1 + a2 D2 + +aT DT +eit, (10)Dt = yi t, xit, i = 1, 2, , N; t = 1,