面板数据的F检验-固定效应检验.doc

. 面板数据模型(PANEL DATA)F检验,固定效应检验 1面板数据定义。时间序列数据或截面数据都是一维数据。例如时间序列数据是变量按时间得到的数据；截面数据是变量在截面空间上的数据。面板数据（panel data）也称时间序列截面数据（time series and cross section data）或混合数据（pool data）。面板数据是同时在时间和截面空间上取得的二维数据。面板数据示意图见图1。面板数据从横截面（cross section）上看，是由若干个体（entity, unit, individual）在某一时刻构成的截面观测值，从纵剖面（longitudinal section）上看是一个时间序列。面板数据用双下标变量表示。例如yi t, i = 1, 2, , N; t = 1, 2, , TN表示面板数据中含有N个个体。T表示时间序列的最大长度。若固定t不变，yi ., ( i = 1, 2, , N)是横截面上的N个随机变量；若固定i不变，y. t, (t = 1, 2, , T)是纵剖面上的一个时间序列（个体）。图1 N=7，T=50的面板数据示意图例如1990-2000年30个省份的农业总产值数据。固定在某一年份上，它是由30个农业总产总值数字组成的截面数据；固定在某一省份上，它是由11年农业总产值数据组成的一个时间序列。面板数据由30个个体组成。共有330个观测值。对于面板数据yi t, i = 1, 2, , N; t = 1, 2, , T来说，如果从横截面上看，每个变量都有观测值，从纵剖面上看，每一期都有观测值，则称此面板数据为平衡面板数据（balanced panel data）。若在面板数据中丢失若干个观测值，则称此面板数据为非平衡面板数据（unbalanced panel data）。注意：EViwes 3.1、4.1、5.0既允许用平衡面板数据也允许用非平衡面板数据估计模型。例1（file:panel02）：1996-2002年中国东北、华北、华东15个省级地区的居民家庭人均消费（不变价格）和人均收入数据见表1和表2。数据是7年的，每一年都有15个数据，共105组观测值。人均消费和收入两个面板数据都是平衡面板数据，各有15个个体。人均消费和收入的面板数据从纵剖面观察分别见图2和图3。从横截面观察分别见图4和图5。横截面数据散点图的表现与观测值顺序有关。图4和图5中人均消费和收入观测值顺序是按地区名的汉语拼音字母顺序排序的。表1 1999-2002年中国东北、华北、华东15个省级地区的居民家庭人均消费数据（不变价格）地区人均消费1996199719981999200020012002CP-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.9903596.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（内蒙古）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.6226904.3687220.843CP-ZJ（浙江）5342.2346002.0826236.6406600.7496950.7137968.3278792.210资料来源：中国统计年鉴1997-2003。表2 1999-2002年中国东北、华北、华东15个省级地区的居民家庭人均收入数据（不变价格）地区人均收入1996199719981999200020012002IP-AH（安徽）4106.2514540.2474770.4705178.5285256.7535640.5976093.333IP-BJ（北京）6569.9017419.9058273.4189127.9929999.70011229.6612692.38IP-FJ（福建）4884.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.4377316.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.0378773.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资料来源：中国统计年鉴1997-2003。 图2 15个省级地区的人均消费序列（纵剖面） 图3 15个省级地区的人均收入序列（file:4panel02） 图4 15个省级地区的人均消费散点图 图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中每一种符号代表一个年度的截面散点图（共7个截面）。相当于观察7个截面散点图的叠加。图6 用15个时间序列表示的人均消费对收入的面板数据图7 用7个截面表示的人均消费对收入的面板数据（7个截面叠加） 为了观察得更清楚一些，图8给出北京和内蒙古1996-2002年消费对收入散点图。从图中可以看出，无论是从收入还是从消费看内蒙古的水平都低于北京市。内蒙古2002年的收入与消费规模还不如北京市1996年的大。图9给出该15个省级地区1996和2002年的消费对收入散点图。可见6年之后15个地区的消费和收入都有了相应的提高。 图8 北京和内蒙古1996-2002年消费对收入时序图 图9 1996和2002年15个地区的消费对收入散点图2面板数据的估计。用面板数据建立的模型通常有3种。即混合估计模型、固定效应模型和随机效应模型。2.1 混合估计模型。如果从时间上看，不同个体之间不存在显著性差异；从截面上看，不同截面之间也不存在显著性差异，那么就可以直接把面板数据混合在一起用普通最小二乘法（OLS）估计参数。如果从时间和截面看模型截距都不为零，且是一个相同的常数，以二变量模型为例，则建立如下模型， yit = a +b1 xit +eit, i = 1, 2, , N; t = 1, 2, , T (1)a 和b1不随i，t变化。称模型(1)为混合估计模型。以例1中15个地区1996和2002年数据建立关于消费的混合估计模型，得结果如下：图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 Estimation（混合估计）窗口如下图。图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 = 4824588, t0.05 (103) = 1.9915个省级地区的人均支出平均占收入的76%。如果从时间和截面上看模型截距都为零，就可以建立不含截距项的（a = 0）的混合估计模型。以二变量模型为例，建立混合估计模型如下， yit = b1 xit +eit, i = 1, 2, , N; t = 1, 2, , T (2)对于本例，因为上式中的截距项有显著性（t = 2.0 t0.05 (103) = 1.99），所以建立截距项为零的混合估计模型是不合适的。EViwes估计方法：在Pooled Estimation（混合估计）对话框中Intercept（截距项）选择窗中选None，其余选项同上。2.2 固定效应模型。在面板数据散点图中，如果对于不同的截面或不同的时间序列，模型的截距是不同的，则可以采用在模型中加虚拟变量的方法估计回归参数，称此种模型为固定效应模型（fixed effects regression model）。固定效应模型分为3种类型，即个体固定效应模型（entity fixed effects regression model）、时刻固定效应模型（time fixed effects regression model）和时刻个体固定效应模型（time and entity fixed effects regression model）。下面分别介绍。（1）个体固定效应模型。个体固定效应模型就是对于不同的个体有不同截距的模型。如果对于不同的时间序列（个体）截距是不同的，但是对于不同的横截面，模型的截距没有显著性变化，那么就应该建立个体固定效应模型，表示如下， yit = b1 xit +g1 W1 + g2 W2 + +gN WN +eit, t = 1, 2, , T (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, 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为k1阶列向量。进一步写成矩阵形式，= + b + 上式中的元素1，0都是T1阶列向量。 面板数据模型用OLS方法估计时应满足如下5个假定条件：（1）E(eit|xi1, xi2, , xiT, ai) = 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很大，相对较小时，因为模型中含有k + N个被估参数，一般软件执行OLS运算很困难。在计量经济学软件中是采用一种特殊处理方式进行OLS估计。估计原理是，先用每个变量减其组内均值，把数据中心化（entity-demeaned），然后用变换的数据先估计个体固定效应模型的回归系数（不包括截距项），然后利用组内均值等式计算截距项。这种方法计算起来速度快。具体分3步如下。（1）首先把变量中心化（entity-demeaned）。仍以单解释变量模型（3）为例，则有 = gi + b1 + , i = 1, 2, , N (4)其中 = ， = ， = , i = 1, 2, , N。公式(1)、(4)相减得， (yit - ) = b1(xit - ) + (eit - ) (5)令(yit - ) = ，(xit - ) = ，(eit - ) = ，上式写为 = b1 + (6)用OLS法估计（1）、（6）式中的b1，结果是一样的，但是用（6）式估计，可以减少被估参数个数。（2）用OLS法估计回归参数（不包括截距项，即固定效应）。 在k个解释变量条件下，把 用向量形式 表示，则利用中心化数据，按OLS法估计公式计算个体固定效应模型中回归参数估计量的方差协方差矩阵估计式如下，( ) = ( )-1 (7)其中 = ， 是相对于 的残差向量。（3）计算回归模型截距项，即固定效应参数gi。= - (8)以例1（file:panel02）为例得到的个体固定效应模型估计结果如下：注意：个体固定效应模型的EViwes输出结果中没有公共截距项。图12EViwes估计方法：在EViwes的Pooled Estimation对话框中Intercept选项中选Fixed effects。其余选项同上。注意：（1）个体固定效应模型的EViwes输出结果中没有公共截距项。（2）EViwes输出结果中没有给出描述个体效应的截距项相应的标准差和t值。不认为截距项是模型中的重要参数。（3）当对个体固定效应模型选择加权估计时，输出结果将给出加权估计和非加权估计两种统计量评价结果。（4）输出结果的联立方程组形式可以通过点击View选Representations功能获得。（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) = 北京+ 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分别表示约束模型（混合估计模型）和非约束模型（个体固定效应模型）的残差平方和。非约束模型比约束模型多了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, i = 1, 2, , N (10)其中Dt = eit, i = 1, 2, , N; t = 1, 2, , T，表示随机误差项。yi t, xit, i = 1, 2, , N; t = 1, 2, , T分别表示被解释变量和解释变量。模型（10）也可表示为 yi1 = a1 +b1 xi1 + ei1, t = 1，（对于第1个截面），i = 1, 2, , N yi2 = (a1 +a2) +b1 xi2 + ei2, t = 2，（对于第2个截面），i = 1, 2, , N yiT = (a1 +aT) +b1 xiT + eiT, t = T，（对于第T个截面），i = 1, 2, , N如果满足上述模型假定条件，对模型（2）进行OLS估计，全部参数估计量都具有无偏性和一致性。模型的自由度是N T T-1。图13EViwes估计方法：在Pooled Estimation（混合估计）窗口中的Dependent Variable（相依变量）选择窗填入CP?；在Common coefficients（系数相同）选择窗填入IP? 和虚拟变量D1997, D1998, D1999, D2000, D2001, D2002；在Cross section specific coefficients（截面系数不同）选择窗保持空白；在Intercept（截距项）选择窗点击Common；在Weighting（权数）选择窗点击No weighting。点击Pooled Estimation（混合估计）窗口中的OK键。以例1为例得到的时刻固定效应模型估计结果如下： = 1996 + xi1 = 108.5057 + 0.7789 xi1 (1.5) (74.6)= 1997 + xi2 = 108.5057 +28.1273 + 0.7789 xi2 (1.5) (0.4) (74.6) = 2002 + xi7 = 108.5057 -199.8213 + 0.7789 xi7 (1.5) (0.4) (74.6) R2 = 0.9867, SSEr = 4028843, t0.05 (97) = 1.98相对于混合估计模型来说，是否有必要建立时刻固定效应模型可以通过F检验来完成。H0：对于不同横截面模型截距项相同（建立混合估计模型）。H1：对于不同横截面模型的截距项不同（建立时刻固定效应模型）。F统计量定义为：F= = (11)其中SSEr，SSEu分别表示约束模型（混合估计模型的）和非约束模型（时刻固定效应模型的）的残差平方和。非约束模型比约束模型多了T-1个被估参数。注意：当模型中含有k个解释变量时，F统计量的分母自由度是NT-T- k。 用上例计算，已知SSEr= 4824588，SSEu= 4028843，F= = = = 3.19F0.05(6, 87) = 2.2因为F= 3.19 F0.05(14, 89) = 2.2，拒绝原假设，结论是应该建立时刻固定效应模型。（3）时刻个体固定效应模型。时刻个体固定效应模型就是对于不同的截面（时刻点）、不同的时间序列（个体）都有不同截距的模型。如果确知对于不同的截面、不同的时间序列（个体）模型的截距都显著地不相同，那么应该建立时刻个体效应模型，表示如下，yit = b1 xit +a1+a2D2 +aT DT +g1W1+g2W2 +gN WN+eit, i=1,2,N，t = 1, 2, , T (12)其中虚拟变量Dt = （注意不是从1开始）Wi = （注意是从1开始）eit, i = 1, 2, , N; t = 1, 2, , T，表示随机误差项。yi t, xit, (i = 1, 2, , N; t = 1, 2, , T)分别表示被解释变量和解释变量。模型也可表示为 y11 = a1 +g1 +b1 x11 + e11, t = 1，i = 1（对于第1个截面、第1个个体） y21 = a1 +g2 +b1 x21 + e21, t = 1，i = 2（对于第1个截面、第2个个体） yN1 = a1 +gN +b1 xN1 + eN1, t = 1，i = N（对于第1个截面、第N个个体） y12 = (a1 +a2) +g1 +b1 x12 + e12, t = 2，i = 1（对于第2个截面、第1个个体） y22 = (a1 +a2) +g2 +b1 x22 + e22, t = 2，i = 2（对于第2个截面、第2个个体） yN2 = (a1 +a2) +gN +b1 xN2 + eN2, t = 2，i = N（对于第2个截面、第N个个体） y1T = (a1 +aT) +g1 +b1 x12 + e1T, t = T，i = 1（对于第T个截面、第1个个体） y2T = (a1 +aT) +g2 +b1 x22 + e2T, t = T，i = 2（对于第T个截面、第2个个体） yNT = (a1 +aT) +gN +b1 xNT + eNT, t = T，i = N（对于第T个截面、第N个个体）如果满足上述模型假定条件，对模型（12）进行OLS估计，全部参数估计量都是无偏的和一致的。模型的自由度是N T NT。注意：当模型中含有k个解释变量时，F统计量的分母自由度是NT N -T- k+1。以例1为例得到的截面、时刻固定效应模型估计结果如下：图14EViwes估计方法：在Pooled Estimation（混合估计）窗口中的Dependent Variable（相依变量）选择窗填入CP?；在Common coefficients（系数相同）选择窗填入IP? 和虚拟变量D1997, D1998, D1999, D2000, D2001, D2002；在Cross section specific coefficients（截面系数不同）选择窗保持空白；在Intercept（截距项）选择窗中选Fixed effects；在Weighting（权数）选择窗点击No weighting。点击Pooled Estimation（混合估计）窗口中的OK键。注意：（1）对于第1个截面（t=1）EViwes输出结果中把(a1 +gi), (i = 1, 2, , N)估计在一起。（2）对于第2, , T个截面（t=1）EViwes输出结果中分别把(a1 +at), (t = 2, , T)估计在一起。输出结果如下： = 1996 + x11 = 537.9627 + 0.6712 x11， （1996年安徽省） = 1996 + x21 = 1223.758 + 0.6712x21， （1996年北京市） = 1997 + x11 = 98.91126 + 0.6712 x11， （1997年安徽省） = 1997 + x21 = 98.91126 +1223.758 + 0.6712x21， （1997年北京市） = 2002 + + x15,7 = (183.3882 +870.4197) + 0.6712 x15,1，（2002年浙江省）R2 = 0.9932, SSEr = 2045670, t0.05 (83) = 1.98相对于混合估计模型来说，是否有必要建立时刻个体固定效应模型可以通过F检验来完成。H0：对于不同横截面，不同序列，模型截距项都相同（建立混合估计模型）。H1：不同横截面，不同序列，模型截距项各不相同（建立时刻个体固定效应模型）。F统计量定义为：F= = (13)其中SSEr，SSEu分别表示约束模型（混合估计模型的）和非约束模型（时刻个体固定效应模型的）的残差平方和。非约束模型比约束模型多了N+T个被估参数。注意：当模型中含有k个解释变量时，F统计量的分母自由度是NT-N-T- k+1。 用上例计算，已知SSEr= 4824588，SSEu= 2045670，F= = = = 5.6F0.05(20, 81) = 1.64因为F= 5.6 F0.05(14, 89) = 1.64，拒绝原假设，结论是应该建立时刻个体固定效应模型。（4）随机效应模型在固定效应模型中采用虚拟变量的原因是解释被解释变量的信息不够完整。也可以通过对误差项的分解来描述这种信息的缺失。 yit = a + b1 xit + eit (14)其中误差项在时间上和截面上都是相关的，用3个分量表示如下。eit = ui + vt + wit (15)其中ui N(0, su2)表示截面随机误差分量；vt N(0, sv2)表示时间随机误差分量；wit N(0, sw2)表示混和随机误差分量。同时还假定ui，vt，wit之间互不相关，各自分别不存在截面自相关、时间自相关和混和自相关。上述模型称为随机效应模型。随机效应模型和固定效应模型比较，相当于把固定效应模型中的截距项看成两个随机变量。一个是截面随机误差项（ui），一个是时间随机误差项（vt）。如果这两个随机误差项都服从正态分布，对模型估计时就能够节省自由度，因为此条件下只需要估计两个随机误差项的均值和方差。假定固定效应模型中的截距项包括了截面随机误差项和时间随机误差项的平均效应，而且对均值的离差分别是ui和vt，固定效应模型就变成了随机效应模型。为了容易理解，先假定模型中只存在截面随机误差项ui，不存在时间随机误差分量（vt）， yit = a + b1 xit + (wit+ ui) = a + b1 xit +eit (16)截面随机误差项ui是属于第个个体的随机波动分量，并在整个时间范围(t = 1,2, , T)保持不变。随机误差项ui, wit应满足如下条件：E(ui) =0, E(wit) = 0E(wit 2) = sw2, E(ui 2)= su2,E(wit uj) =0, 包括所有的i, t, j。E(wit wjs) =0, i j, t sE(ui uj) =0, i j因为根据上式有eit = wit+ ui所以这种随机效应模型又称为误差分量模型（error component model）。有结论，E(eit ) = E(wit +uj) = 0,(16)式，yit = a + b1 xit + (wit+ ui