面板大数据地F检验,固定效应检验

1、面板数据模型(PANEL DATA)F检验，固定效应检验1面板数据定义。时间序列数据或截面数据都是一维数据。例如时间序列数据是变量按时间得到的数据；截面数据是变量在截面空间上的数据。面板数据(panel data )也称时间序列截面数据(timeseries and cross section data )或混合数据(pool data )。面板数据是同时在时间和截 面空间上取得的二维数据。面板数据示意图见图1。面板数据从横截面(cross section )上看，是由若干个体(en tity, un it, in dividual )在某一时刻构成的截面观测值，从纵剖 面(Iongitudi

2、nal section)上看是一个时间序列。面板数据用双下标变量表示。例如yi t , i = 1,2,N t = 1,2,TN表示面板数据中含有N个个体。T表示时间序列的最大长度。若固定 t不变，y，( i = 1,2,N)是横截面上的 N个随机变量；若固定i不变，y. t , ( t = 1,2,T)是纵剖面 上的一个时间序列(个体)。图1N=7 , T=50的面板数据示意图例如1990-2000年30个省份的农业总产值数据。 固定在某一年份上， 它是由30个农业总产 总值数字组成的截面数据；固定在某一省份上，它是由 11年农业总产值数据组成的 一个时 间序列。面板数据由 30个个体组成。

3、共有 330个观测值。对于面板数据 yt , i = 1,2,N; t = 1,2,T来说，如果从横截面上看，每个变量都有观测值，从纵剖面上看，每一期都有观测值，则称此面板数据为平衡面板数据(bala need pa nel data)。若在面板数据中丢失若干个观测值，则称此面板数据为非平衡 面板数据(unbalaneed panel data )。注意：EViwes 3.1、4.1、5.0既允许用平衡面板数据也允许用非平衡面板数据估计模型。例1(file:panel02 )： 1996-2002年中国东北、华北、华东 15个省级地区的居民家庭人均 消费(不变价格)和人均收入数据见表1和表2。

4、数据是7年的，每一年都有15个数据,共105组观测值。人均消费和收入两个面板数据都是平衡面板数据，各有15个个体。人均消费和收入的面板数据从纵剖面观察分别见图2和图3。从横截面观察分别见图4和图5。横截面数据散点图的表现与观测值顺序有关。图4和图5中人均消费和收入观测值顺序是按地区名的汉语拼音字母顺序排序的。表11999-2002年中国东北、华北、华东15个省级地区的居民家庭人均消费数据(不变价格)地区人均消费1996199719981999200020012002CP-AH（）3282.4663646.1503777.4103989.5814203.5554495.1744784.364CP

5、-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.96142

6、81.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-NMGQ2572.3422901.7223127.6333475.9423877.3454170.5964850.180CP-SD()3440.6843930.5744168.974454

7、6.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-200

8、3。表21999-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.

9、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.

10、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

11、.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:4pa nel02）图4 15个省级地区的人均消费散点图图5 15个省级地区的人均收入散点图（7个横截面叠加）（每条连线表示同一年度 15个地区的消费值）（每条连线表示同一年度 15个地区的收入值）用 CP表

12、示消费，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给出和1

13、996-2002年消费对收入散点图。从图中可以看出，无论是从收入还是从消费看的水平都低于市。2002年的收入与消费规模还不如市1996年的大。图9给出该15个省级地区1996和2002年的消费对收入散点图。可见6年之后15个地区的消费和收入都有了相应的提高。图8 和1996-2002年消费对收入时序图图91996和2002年15个地区的消费对收入散点图2 面板数据的估计。用面板数据建立的模型通常有3种。即混合估计模型、固定效应模型和随机效应模型。2.1混合估计模型。如果从时间上看，不同个体之间不存在显著性差异；从截面上看，不同截面之间也不存在显著性差异，那么就可以直接把面板数据混合在一起用普通

14、最小二乘法（OLS估计参数。如果从时间和截面看模型截距都不为零，且是一个相同的常数， 以二变量模型为例， 则建立如下模型，yit = a +b1 Xit +eit, i = 1,2, ,N; t = 1,2,T（1）a和b1不随i, t变化。称模型（1）为混合估计模型。以例1中15个地区1996和2002年数据建立关于消费的混合估计模型，得结果如下:图10EViwes估计方法：在打开工作文件窗口的基础上，点击主功能菜单中的Objects键，选NewObject功能，从而打开 NewObject （新对象）选择窗。在 Type of Object选择区选择 Pool（混合数据库），点击 OK键，

15、从而打开Pool （混合数据）窗口。在窗口中输入 15个地区 标识AH（）、BJ （）、ZJ （）。工具栏中点击 Sheet键，从而打开 Series List （列写 序列名）窗口，定义变量CP刑IP?，点击OKS, Pool （混合或合并数据库）窗口显示面板数据。在Pool窗口的工具栏中点击Estimate键，打开Pooled Estimation （混合估计）窗口如下图。图11在Dependent Variable （相依变量）选择窗填入CP?在Commoicoefficients（系数相同）选择窗填入ip? ； Cross section specific coefficients（截

16、面系数不同）选择窗保持空白；在Intercept（截距项）选择窗点击Common在 Weighting （权数）选择窗点击Noweighting。点击Pooled Estimation （混合估计）窗口中的 OK键。得输出结果如图10。相应表达式是=129.6313+0.7587 IPit（2.0）（79.7）R = 0.98, SSE = 4824588,t°.°5（ =1.9915个省级地区的人均支出平均占收入的76%如果从时间和截面上看模型截距都为零，就可以建立不含截距项的（a = 0 ）的混合估计模型。以二变量模型为例，建立混合估计模型如下，yit = b1 xit

17、 +eit,i = 1,2,N; t = 1,2,对于本例，因为上式中的截距项有显著性(t = 2.0 > t 0.05 (103) = 1.99 )，所以建立截距项为零的混合估计模型是不合适的。EViwes估计方法：在 Pooled Estimation(混合估计)对话框中Intercept(截距项)选择窗中选None,其余选项同上。2.2固定效应模型。在面板数据散点图中，如果对于不同的截面或不同的时间序列，模型的截距是不同的，则可以采用在模型中加虚拟变量的方法估计回归参数，称此种模型为固定效应模型(fixedeffects regressi on model )。固定效应模型分为 3

18、种类型，即个体固定效应模型(en tity fixed effects regressionmodel)、时刻固定效应模型(time fixed effects regression model)和时刻个体固定效 应模型(time and entity fixed effects regression model)。下面分另U介绍。(1) 个体固定效应模型。个体固定效应模型就是对于不同的个体有不同截距的模型。如果对于不同的时间序列(个体)截距是不同的，但是对于不同的横截面，模型的截距没有显著性变化，那么就应该建立个体固定效应模型，表示如下，yit = b1 xit +g1 W + g W +

19、+ gN WN +©t,t = 1,2,T(3)其中W=e*, i = 1,2,N t = 1,2,T,表示随机误差项。yj xj i = 1,2,N t =1,2,T分别表示被解释变量和解释变量。模型(3)或者表示为yit = gi + bi Xit + eit,y2t = g2 + bi X2t +e2 ti = 1 (对于第1个个体，或时间序列)，t = 1,2,i = 2 (对于第2个个体，或时间序列)，t = i, 2,yN t = gN +bi xn t +e n t, i = N (对于第 N个个体，或时间序列)，t = i, 2,T写成矩阵形式，yi = ( iXi)

20、 + ei = gi + xi b + eiyN =( ixn) + eN = gN + x n b +eN上式中yi,gi,e,x i都是N，i阶列向量。b为标量。当模型中含有k个解释变量时，b为k '1 阶列向量。进一步写成矩阵形式，=+ b +上式中的元素i, 0都是1阶列向量。面板数据模型用 OLS方法估计时应满足如下 5个假定条件：(1) E(eit|Xii, Xi2,，x iT, ai) = 0。以 Xii, Xi2,Xi, ai 为条件的 e 的期望等于零。(2) (Xii, Xi2,，XiT), ( yii, yi2,，y r), i = i, 2,N 分别来自于同一个

21、联合分 布总体，并相互独立。(3) (xit, eit)具有非零的有限值 4阶矩。(4) 解释变量之间不存在完全共线性。（5）Cov（et eis|Xit,Xis, a） = 0, t 1 s。在固定效应模型中随机误差项eit在时间上是非自相关的。其中Xit代表一个或多个解释变量。对模型（1）进行OLSf古计，全部参数估计量都是无偏的和一致的。模型的自由度是 NT - 1 - N当模型含有k个解释变量，且 N很大，相对较小时，因为模型中含有k+ N个被估参数，一般软件执行 OLS运算很困难。在计量经济学软件中是采用一种特殊处理方式进行OLS估计。估计原理是，先用每个变量减其组均值，把数据中心化

22、（en tity-demea ned），然后用变换的数据先估计个体固定效应模型的回归系数（不包括截距项），然后利用组均值等式计算截距项。这种方法计算起来速度快。具体分3步如下。（1）首先把变量中心化（en tity-demea ned )。仍以单解释变量模型（3）为例，则有=gi+ b1 + ,i=1,2,N(4)其中=,=,=,i=1,2,N=公式（1）、相减得，(屮-)=b( xit -)+ ( eit-)(5)令(屮-)=,(Xit -)=,-)=，上式写为=b1+用OLS法估计（1 ）、（ 6）式中的b1,结果是一样的，但是用（6）式估计，可以减少被估参数个数。(2 )用OLS法估计回

23、归参数(不包括截距项，即固定效应)。在k个解释变量条件下，把用向量形式 表示，则利用中心化数据,公式计算个体固定效应模型中回归参数估计量的方差协方差矩阵估计式如下，按OLS法估计()=(')-1其中=，是相对于的残差向量。(3)计算回归模型截距项，即固定效应参数gio(8)以例1 (file:panel02)为例得到的个体固定效应模型估计结果如下：注意：个体固定效应模型的EViwes输出结果中没有公共截距项。图12EViwes 估计方法：在 EViwes 的 Pooled Estimation对话框中 Intercepteffects 。其余选项同上。选项中选Fixed(1) 个体固定

24、效应模型的EViwes输出结果中没有公共截距项。(2) EViwes输出结果中没有给出描述个体效应的截距项相应的标准差和t值。不认为截距项是模型中的重要参数。(3) 当对个体固定效应模型选择加权估计时，输出结果将给出加权估计和非加权估计两种 统计量评价结果。(4) 输出结果的联立方程组形式可以通过点击View选Representations功能获得。(5) 点击View选 Wald Coefficient Tests功能可以对模型的斜率进行Wald检验。(6) 点击 View 选 Residuals/Table, Graphs, Covarianee Matrix, Correlation M

25、atrix功能可以分别得到按个体计算的残差序列表，残差序列图，残差序列的方差协方差矩阵，残差序列的相关系数矩阵。进一步点击Solve(7) 点击Procs选MakeModel功能，将会出现估计结果的联立方程形式,键，在随后出现的对话框中可以进行动态和静态预测。输出结果的方程形式是=+xit = 479.3 + 0.70 xit(55.0)=+ X2t = 1053.2 + 0.70 X2t(55.0)=+ Xi5t = 714.2 + 0.70 Xi5t(55.0)氏=0.99, SSE= 2270386,10.05 =1.98从结果看，、是消费函数截距(自发消费)最大的 3个地区。相对于混合

26、估计模型来说，是否有必要建立个体固定效应模型可以通过F检验来完成。原假设不同个体的模型截距项相同(建立混合估计模型)。备择假设H:不同个体的模型截距项不同（建立个体固定效应模型）。F统计量定义为：F= =（9）其中SSEE, SSEU分别表示约束模型（混合估计模型）和非约束模型（个体固定效应模型）的 残差平方和。非约束模型比约束模型多了N-1个被估参数。（混合估计模型给出公共截距项。）注意：当模型中含有 k个解释变量时，F统计量的分母自由度是 NF2k。用上例计算，已知 SSE= 4824588，SSEU = 2270386，F= = = = 7.15Fo.05（14, 89）= 1.81因为

27、F= 7.15> F 0.05（14, 89） = 1.81 ，所以，拒绝原假设。结论是应该建立个体固定效应模型。（2）时刻固定效应模型。时刻固定效应模型就是对于不同的截面（时刻点）有不同截距的模型。如果确知对于不同的截面，模型的截距显著不同，但是对于不同的时间序列（个体）截距是相同的，那么应该建 立时刻固定效应模型，表示如下，yit = bi Xit +ai + a2 D2 + + aT DT +eit, i = 1,2,N（10）其中D =et, i = 1,2,N; t = 1,2,，表示随机误差项。 y t , x it, i = 1,2,N; t=1,2,T分别表示被解释变量和

28、解释变量。模型（10）也可表示为yi i = a1 +b1 Xi 1 + e 1,Nyi 2 = ( ai + a2)+ b1 Xi 2 + ei2,t = 1,（对于第1个截面），i = 1, 2,t = 2 ,（对于第2个截面），i = 1,2,yiT = (ai + aT) + b XiT + eNt = T,（对于第T个截面），i = 1,2,N如果满足上述模型假定条件，对模型（2）进行OLS估计，全部参数估计量都具有无偏性和一致性。模型的自由度是NT - T-1。图13EViwes 估计方法：在 Pooled Estimation（混合估计）窗口中的 Dependent Variab

29、le （相依变量）选择窗填入CP?;在Commoncoefficients （系数相同）选择窗填入 IP?和虚拟变量 D1997, D1998, D1999, D2000, D2001, D2002;在 Cross section specific coefficients（截面系数不同）选择窗保持空白；在Intercept （截距项）选择窗点击Common在 Weighting（权数）选择窗点击No weighting。点击Pooled Estimation（混合估计）窗口中的OK键。以例1为例得到的时刻固定效应模型估计结果如下：=1996 + Xi 1 = 108.5057 + 0.778

30、9 Xi 1(1.5)(74.6)=1997 + Xi 2 = 108.5057 +28.1273 + 0.7789 Xi2(1.5)(0.4)(74.6)=2002 + Xi 7 = 108.5057 -199.8213 + 0.7789Xi 7(1.5)(0.4)(74.6)氏=0.9867,SSEE = 4028843, 仏。5 =1.98相对于混合估计模型来说，是否有必要建立时刻固定效应模型可以通过F检验来完成。H):对于不同横截面模型截距项相同(建立混合估计模型)。H:对于不同横截面模型的截距项不同(建立时刻固定效应模型)。F统计量定义为：F= =(11)其中SSEE，SSEU分别表

31、示约束模型(混合估计模型的)和非约束模型(时刻固定效应模型的) 的残差平方和。非约束模型比约束模型多了T-1个被估参数。注意：当模型中含有 k个解释变量时，F统计量的分母自由度是NTT- k。用上例计算，已知 SSE= 4824588，SSE= 4028843，F= = = = 3.19Fo.O5(6, 87)= 2.2因为F= 3.19> F 0.05(14, 89) = 2.2 ，拒绝原假设，结论是应该建立时刻固定效应模型。(3) 时刻个体固定效应模型。时刻个体固定效应模型就是对于不同的截面（时刻点）、不同的时间序列（个体）都有不同 截距的模型。如果确知对于不同的截面、不同的时间序列

32、（个体）模型的截距都显著地不相同，那么应该建立时刻个体效应模型，表示如下，yit = bi Xit +ai+azD2 + +aT Dr+giWg2W/+gN WJ+er, i =1,2,，N, t = 1,2,T（12）其中虚拟变量D =（注意不是从1开始）W=（注意是从1开始）et, i = 1,2, ,N t = 1,2, ,T,表示随机误差项。yt , x it, （ i = 1,2,N;t = 1,2,T）分别表示被解释变量和解释变量。模型也可表示为yNi =a + gN +b1 xni + eNi,t = 1 , i = N （对于第1个截面、第N个个体）y12 =（ a1 + a2

33、）+ g1 +b X12 + e12,t = 2 , i = 1 （对于第2个截面、第1个个体）yn =a1 + g1 +b1 X11 +en,t = 1 , i=1 （对于第1个截面、第1个个体）y21 =a1 + g2 +b1 X21 +e21,t = 1 , i=2 （对于第1个截面、第2个个体）y22 = ( ai +a2)+ g2 +b X22 + e22,t = 2 , i = 2 （对于第2个截面、第2y订=(ai +aT)+ gi +b X12 +t = T, i = i （对于第T个截面、第i个个体）yN? = ( ai + a2)+ gN+bi xn2 +eN2,t = 2

34、 , i=N （对于第2个截面、第N个个体）个个体）y2T =( ai +aT)+ g2 +b x + s.t = T, i = 2 （对于第T个截面、第2个个体）yNT =(ai + aT) + gN +bi xnt + eNT,t = T, i = N （对于第T个截面、第 N个个体）如果满足上述模型假定条件，对模型（i2）进行OLS估计，全部参数估计量都是无偏的和一致的。模型的自由度是 NT- N-T。注意：当模型中含有k个解释变量时，F统计量的分母自由度是NT- N - T- k+i。以例i为例得到的截面、时刻固定效应模型估计结果如下：图i4EViwes 估计方法：在 Pooled E

35、stimation （混合估计）窗口中的 Dependent Variable （相 依变量）选择窗填入 CP?;在Commoncoefficients （系数相同）选择窗填入 IP?和虚拟变量 D1997, D1998, D1999, D2000, D2001, D2002;在 Cross section specific coefficients（截面系数不同）选择窗保持空白；在In tercept（截距项）选择窗中选Fixed effects ;在 Weighting （权数）选择窗点击 No weighting。点击Pooled Estimation（混合估计）窗口中的OK键。（1）

36、对于第1个截面（t=1） EViwes输出结果中把（a1 + g） （ i = 1,2,N估计在一 起。（2） 对于第2,T个截面（t=1） EViwes输出结果中分别把（a + 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年市）

37、2002 + +X15,7 = （183.3882 +870.4197） + 0.6712X15,1,（ 2002 年省）氏=0.9932, SSE = 2045670, toa =1.98F检验来完成。相对于混合估计模型来说，是否有必要建立时刻个体固定效应模型可以通过H）:对于不同横截面，不同序列，模型截距项都相同（建立混合估计模型）。H:不同横截面，不同序列，模型截距项各不相同（建立时刻个体固定效应模型）。F统计量定义为：F=（13）其中SSEE, SSEU分别表示约束模型（混合估计模型的）和非约束模型（时刻个体固定效应模型的）的残差平方和。非约束模型比约束模型多了N+T个被估参数。注意：

38、当模型中含有 k个解释变量时，F统计量的分母自由度是 NT2T- k+1。用上例计算，已知 SSEE= 4824588，SSEU= 2045670,F= = = = 5.6Fo.O5（2O, 81）= 1.64因为F= 5.6> Fo.O5（14, 89） = 1.64 ，拒绝原假设，结论是应该建立时刻个体固定效应模型。（4）随机效应模型在固定效应模型中采用虚拟变量的原因是解释被解释变量的信息不够完整。也可以通过对误差项的分解来描述这种信息的缺失。yit = a + b1 Xit +et其中误差项在时间上和截面上都是相关的，用3个分量表示如下。et = u i + v t +w(15)其

39、中Ui N(0, si)表示截面随机误差分量；Vt N(0, Sv)表示时间随机误差分量；Wt N(0,Sw2)表示混和随机误差分量。同时还假定Ui , Vt，w之间互不相关，各自分别不存在截面自相关、时间自相关和混和自相关。上述模型称为随机效应模型。随机效应模型和固定效应模型比较，相当于把固定效应模型中的截距项看成两个随机变量。一个是截面随机误差项(Ui), 个是时间随机误差项(vt)。如果这两个随机误差项都服从正态分布，对模型估计时就能够节省自由度，因为此条件下只需要估计两个随机误差项的均值和方差。假定固定效应模型中的截距项包括了截面随机误差项和时间随机误差项的平均效应，而且对均值的离差分

40、别是 Ui和Vt，固定效应模型就变成了随机效应模型。为了容易理解，先假定模型中只存在截面随机误差项Ui，不存在时间随机误差分量(vt)，yit = a + bi Xit + ( wt + u i) = a + bi Xit+eit(16)截面随机误差项Ui是属于第个个体的随机波动分量，并在整个时间围(t = 1,2,T)保持不变。随机误差项 Ui, wt应满足如下条件：E(Ui) =0,E(wt ) = 02 2E(w/t ) = s w ,2 2E(Ui )= Su ,E(wt)=0,包括所有的i , t, j。E(wt Wjs) =0, i 1 j , t 1 sE(ui u j) =0,

41、 i 1 j因为根据上式有©t = wt + u ierror comp onent model)。有结论,所以这种随机效应模型又称为误差分量模型（E(eit ) = E( wt +u) = 0,(16)式，yit= a +b1xit+ (wt+ui)，也可以写成yit= ( a + uj + bxit+wt。服从正态分布的截距项的均值效应au被包含在回归函数的常数项中。2 2 2 2E(eit ) = E( w +Uj) = s w +Su,2 2E(eit eis) = E(wt + 山)(wt+ u J = E(wt wis + u i wis + wit u i + u i

42、) = Su, t 1 s令ei = ( ei 1, ei2,er)'则W= E( ei e')=2 . 2 =Sw I (TZT) + Su1(T 1) 1(T 1)其中I（TT ）是（T'T）阶单位阵，1（1）是（T' 1）阶列向量。因为第i期与j期观测值是相互独立的, 所以NT个观测值所对应的随机误差项的方差与协方差矩阵V = = ? W = I N'N ?W其中Inn表示由（T' 1）阶列向量为元素构成的单位阵，其中每一个元素1或0都是（1）阶列向量。？表示科罗克积（Kronecker product ）。其运算规则是Ank?B =检验个

43、体随机效应的原假设与检验统计量是H）: Su2 = 0。（混合估计模型）H： Su2 1 0。（个体随机效应模型）LM=其中 表示由个体随机效应模型计算的残差平方和。表示由混合估计模型计算的残差平方和。统计量LM服从1个自由度的c2分布。可以对随机效应模型进行广义最小二乘估计。以观测值方差的倒数为权。为了求权数，必须采用两阶段最小二乘法估计。因为各随机误差分量的方差一般是未知的，第一阶段用普通最小二乘估计法对混合数据进行估计（采用固定效应模型）。用估计的残差计算随机误差分量的方差。第二步用这些估计的方差计算参数的广义最小二乘估计值。如果随机误差分量服从的是正态分布，模型的参数还可以用极大似然法

44、估计。仍以例1为例给出随机效应模型估计结果如下：图15注意：随机效应模型 EViwes输出结果中含有公共截距项。图16以例1为例，用个体随机效应模型和混合模型计算的统计量的值是2LM= = =8.75 ' (24.4)= 5209Fo.o5 =3.84因为F= 5209 > Fo.05（1）= 3.84，所以拒绝原假设，结论是应该建立个体随机效应模型。假定截面截距和时间截距都是随机的。分别服从均值为au和av,方差为su2和sv2的正态分布。随机误差项将由3部分组成，并有方差。2 2 2Var（ et） = Var（ Ui） + Var（ vt） + Var（ wt） = Su

45、+ Sv +Sw当su2和sv 2都等于零，随机效应模型退化为固定效应模型。随机效应模型和固定效应模型哪一个更好些？实际是各有优缺点。随机效应模型的好处是节省自由度。对于从时间序列和截面两方面上看都存在较大变化的数据，随机效应模型能明确地描述出误差来源的特征。固定效应模型的好处是很容易分析任意截面数据所对应的因变量 与全部截面数据对应的因变量均值的差异程度。此外，固定效应模型不要求误差项中的个体效应分量与模型中的解释变量不相关。当然，这一假定不成立时， 可能会引起模型参数估计的不一致性。（5 ）回归系数不同的面板数据模型当认为对于不同个体，解释变量的回归系数存在显著性差异时，还可以建立回归系数

46、不同的面板数据模型。EViwes 估计方法：在 Pooled Estimation（混合估计）窗口中的Dependent Variable（相依变量）选择窗填入 CP?;在Common coefficients（系数相同）选择窗保持空白（如果需要估计时刻固定效应也可输入虚拟变量D1997, D1998, D1999, D2000, D2001, D2002）;在 Cross section specific coefficients（截面系数不同）选择窗填入IP?;在 Intercept（截距项）选择窗中选Fixed effects（也可以做其他选择）；在 Weight ing （权数）选择窗

47、点击No weighting（也可以做其他选择）。点击 Pooled Estimation（混合估计）窗口中的OK键。图17=+ xit = 161.62 + 0.76xit（9.1）=+ X2t = 36.22 + 0.81 X2t（31.0）=+ X15t = 1328.26 + 0.63X15t（21.1）氏=0.995, SSE = 1409247用EViwes建立面板数据估计模型步骤。利用19962002年15个省级地区城镇居民家庭年人均消费性支出和年人均收入数据（不变价格数据）介绍面板数据模型估计步骤。（1）建立混合数据库（Pool）对象。首先建立工作文件。在打开工作文件窗口的基础上，点击EViwes主功能菜单上的 Objects键，选New Object功能（如图18），从而打开 New Object （新对象）选择窗。在 Type ofObj