eviews關(guān)于面板數(shù)據(jù)模型截距系數(shù)固定效應(yīng)還是隨機效應(yīng)的選取得檢驗方法及具體事例

1、面板數(shù)據(jù)模型1面板數(shù)據(jù)定義。時間序列數(shù)據(jù)或截面數(shù)據(jù)都是一維數(shù)據(jù)。例如時間序列數(shù)據(jù)是變量按時間得到的數(shù)據(jù)；截面數(shù)據(jù)是變量在截面空間上的數(shù)據(jù)。面板數(shù)據(jù)（panel data）也稱時間序列截面數(shù)據(jù)（time series and cross section data）或混合數(shù)據(jù)（pool data）。面板數(shù)據(jù)是同時在時間和截面空間上取得的二維數(shù)據(jù)。面板數(shù)據(jù)示意圖見圖1。面板數(shù)據(jù)從橫截面（cross section）上看，是由若干個體（entity, unit, individual）在某一時刻構(gòu)成的截面觀測值，從縱剖面（longitudinal section）上看是一個時間序列。面板數(shù)據(jù)用雙下標(biāo)變

2、量表示。例如yi t, i = 1, 2, , n; t = 1, 2, , tn表示面板數(shù)據(jù)中含有n個個體。t表示時間序列的最大長度。若固定t不變，yi ., ( i = 1, 2, , n)是橫截面上的n個隨機變量；若固定i不變，y. t, (t = 1, 2, , t)是縱剖面上的一個時間序列（個體）。圖1 n=7，t=50的面板數(shù)據(jù)示意圖例如1990-2000年30個省份的農(nóng)業(yè)總產(chǎn)值數(shù)據(jù)。固定在某一年份上，它是由30個農(nóng)業(yè)總產(chǎn)總值數(shù)字組成的截面數(shù)據(jù)；固定在某一省份上，它是由11年農(nóng)業(yè)總產(chǎn)值數(shù)據(jù)組成的一個時間序列。面板數(shù)據(jù)由30個個體組成。共有330個觀測值。對于面板數(shù)據(jù)yi t, i

3、= 1, 2, , n; t = 1, 2, , t來說，如果從橫截面上看，每個變量都有觀測值，從縱剖面上看，每一期都有觀測值，則稱此面板數(shù)據(jù)為平衡面板數(shù)據(jù)（balanced panel data）。若在面板數(shù)據(jù)中丟失若干個觀測值，則稱此面板數(shù)據(jù)為非平衡面板數(shù)據(jù)（unbalanced panel data）。注意：eviwes 3.1、4.1、5.0既允許用平衡面板數(shù)據(jù)也允許用非平衡面板數(shù)據(jù)估計模型。例1（file:panel02）：1996-2002年中國東北、華北、華東15個省級地區(qū)的居民家庭人均消費（不變價格）和人均收入數(shù)據(jù)見表1和表2。數(shù)據(jù)是7年的，每一年都有15個數(shù)據(jù)，共105組觀測

4、值。人均消費和收入兩個面板數(shù)據(jù)都是平衡面板數(shù)據(jù)，各有15個個體。人均消費和收入的面板數(shù)據(jù)從縱剖面觀察分別見圖2和圖3。從橫截面觀察分別見圖4和圖5。橫截面數(shù)據(jù)散點圖的表現(xiàn)與觀測值順序有關(guān)。圖4和圖5中人均消費和收入觀測值順序是按地區(qū)名的漢語拼音字母順序排序的。表1 1999-2002年中國東北、華北、華東15個省級地區(qū)的居民家庭人均消費數(shù)據(jù)（不變價格）地區(qū)人均消費1996199719981999200020012002cp-ah（安徽） 3282.466 3646.150 3777.410 3989.581 4203.555 4495.174 4784.364cp-bj（北京） 5133.97

5、8 6203.048 6807.451 7453.757 8206.271 8654.433 10473.12cp-fj（福建） 4011.775 4853.441 5197.041 5314.521 5522.762 6094.336 6665.005cp-hb（河北） 3197.339 3868.319 3896.778 4104.281 4361.555 4457.463 5120.485cp-hlj（黑龍江） 2904.687 3077.989 3289.990 3596.839 3890.580 4159.087 4493.535cp-jl（吉林） 2833.321 3286.432

6、 3477.560 3736.408 4077.961 4281.560 4998.874cp-js（江蘇） 3712.260 4457.788 4918.944 5076.910 5317.862 5488.829 6091.331cp-jx（江西） 2714.124 3136.873 3234.465 3531.775 3612.722 3914.080 4544.775cp-ln（遼寧） 3237.275 3608.060 3918.167 4046.582 4360.420 4654.420 5402.063cp-nmg（內(nèi)蒙古） 2572.342 2901.722 3127.633

7、3475.942 3877.345 4170.596 4850.180cp-sd（山東） 3440.684 3930.574 4168.974 4546.878 5011.976 5159.538 5635.770cp-sh（上海） 6193.333 6634.183 6866.410 8125.803 8651.893 9336.100 10411.94cp-sx（山西） 2813.336 3131.629 3314.097 3507.008 3793.908 4131.273 4787.561cp-tj（天津） 4293.220 5047.672 5498.503 5916.613 614

8、5.622 6904.368 7220.843cp-zj（浙江） 5342.234 6002.082 6236.640 6600.749 6950.713 7968.327 8792.210資料來源：中國統(tǒng)計年鑒1997-2003。表2 1999-2002年中國東北、華北、華東15個省級地區(qū)的居民家庭人均收入數(shù)據(jù)（不變價格）地區(qū)人均收入1996199719981999200020012002ip-ah（安徽） 4106.251 4540.247 4770.470 5178.528 5256.753 5640.597 6093.333ip-bj（北京） 6569.901 7419.905 827

9、3.418 9127.992 9999.700 11229.66 12692.38ip-fj（福建） 4884.731 6040.944 6505.145 6922.109 7279.393 8422.573 9235.538ip-hb（河北） 4148.282 4790.986 5167.317 5468.940 5678.195 5955.045 6747.152ip-hlj（黑龍江） 3518.497 3918.314 4251.494 4747.045 4997.843 5382.808 6143.565ip-jl（吉林） 3549.935 4041.061 4240.565 4571

10、.439 4878.296 5271.925 6291.618ip-js（江蘇） 4744.547 5668.830 6054.175 6624.316 6793.437 7316.567 8243.589ip-jx（江西） 3487.269 3991.490 4209.327 4787.606 5088.315 5533.688 6329.311ip-ln（遼寧） 3899.194 4382.250 4649.789 4968.164 5363.153 5797.010 6597.088ip-nmg（內(nèi)蒙古） 3189.414 3774.804 4383.706 4780.090 5063.

11、228 5502.873 6038.922ip-sd（山東） 4461.934 5049.407 5412.555 5849.909 6477.016 6975.521 7668.036ip-sh（上海） 7489.451 8209.037 8773.100 10770.09 11432.20 12883.46 13183.88ip-sx（山西） 3431.594 3869.952 4156.927 4360.050 4546.785 5401.854 6335.732ip-tj（天津） 5474.963 6409.690 7146.271 7734.914 8173.193 8852.470

12、 9375.060ip-zj（浙江） 6446.515 7158.288 7860.341 8530.314 9187.287 10485.64 11822.00資料來源：中國統(tǒng)計年鑒1997-2003。 圖2 15個省級地區(qū)的人均消費序列（縱剖面） 圖3 15個省級地區(qū)的人均收入序列（file:4panel02） 圖4 15個省級地區(qū)的人均消費散點圖 圖5 15個省級地區(qū)的人均收入散點圖（7個橫截面疊加）(每條連線表示同一年度15個地區(qū)的消費值) (每條連線表示同一年度15個地區(qū)的收入值)用cp表示消費，ip表示收入。ah, bj, fj, hb, hlj, jl, js, jx, ln,

13、nmg, sd, sh, sx, tj, zj分別表示安徽省、北京市、福建省、河北省、黑龍江省、吉林省、江蘇省、江西省、遼寧省、內(nèi)蒙古自治區(qū)、山東省、上海市、山西省、天津市、浙江省。15個地區(qū)7年人均消費對收入的面板數(shù)據(jù)散點圖見圖6和圖7。圖6中每一種符號代表一個省級地區(qū)的7個觀測點組成的時間序列。相當(dāng)于觀察15個時間序列。圖7中每一種符號代表一個年度的截面散點圖（共7個截面）。相當(dāng)于觀察7個截面散點圖的疊加。圖6 用15個時間序列表示的人均消費對收入的面板數(shù)據(jù)圖7 用7個截面表示的人均消費對收入的面板數(shù)據(jù)（7個截面疊加） 為了觀察得更清楚一些，圖8給出北京和內(nèi)蒙古1996-2002年消費對收

14、入散點圖。從圖中可以看出，無論是從收入還是從消費看內(nèi)蒙古的水平都低于北京市。內(nèi)蒙古2002年的收入與消費規(guī)模還不如北京市1996年的大。圖9給出該15個省級地區(qū)1996和2002年的消費對收入散點圖。可見6年之后15個地區(qū)的消費和收入都有了相應(yīng)的提高。 圖8 北京和內(nèi)蒙古1996-2002年消費對收入時序圖 圖9 1996和2002年15個地區(qū)的消費對收入散點圖2面板數(shù)據(jù)的估計。用面板數(shù)據(jù)建立的模型通常有3種。即混合估計模型、固定效應(yīng)模型和隨機效應(yīng)模型。2.1 混合估計模型。如果從時間上看，不同個體之間不存在顯著性差異；從截面上看，不同截面之間也不存在顯著性差異，那么就可以直接把面板數(shù)據(jù)混合在

15、一起用普通最小二乘法（ols）估計參數(shù)。如果從時間和截面看模型截距都不為零，且是一個相同的常數(shù)，以二變量模型為例，則建立如下模型， yit = a +b1 xit +eit, i = 1, 2, , n; t = 1, 2, , t (1)a 和b1不隨i，t變化。稱模型(1)為混合估計模型。以例1中15個地區(qū)1996和2002年數(shù)據(jù)建立關(guān)于消費的混合估計模型，得結(jié)果如下：圖10 eviwes估計方法：在打開工作文件窗口的基礎(chǔ)上，點擊主功能菜單中的objects鍵，選new object功能，從而打開new object（新對象）選擇窗。在type of object選擇區(qū)選擇pool（混合數(shù)

16、據(jù)庫），點擊ok鍵，從而打開pool（混合數(shù)據(jù)）窗口。在窗口中輸入15個地區(qū)標(biāo)識ah（安徽）、bj（北京）、zj（浙江）。工具欄中點擊sheet鍵，從而打開series list（列寫序列名）窗口，定義變量cp?和ip?，點擊ok鍵，pool（混合或合并數(shù)據(jù)庫）窗口顯示面板數(shù)據(jù)。在pool窗口的工具欄中點擊estimate鍵，打開pooled estimation（混合估計）窗口如下圖。圖11在dependent variable（相依變量）選擇窗填入cp?；在common coefficients（系數(shù)相同）選擇窗填入ip?；cross section specific coefficien

17、ts（截面系數(shù)不同）選擇窗保持空白；在intercept（截距項）選擇窗點擊common；在weighting（權(quán)數(shù)）選擇窗點擊no weighting。點擊pooled estimation（混合估計）窗口中的ok鍵。得輸出結(jié)果如圖10。相應(yīng)表達(dá)式是= 129.6313 +0.7587 ipit (2.0) (79.7) r2 = 0.98, sser = 4824588, t0.05 (103) = 1.9915個省級地區(qū)的人均支出平均占收入的76%。如果從時間和截面上看模型截距都為零，就可以建立不含截距項的（a = 0）的混合估計模型。以二變量模型為例，建立混合估計模型如下， yit =

18、 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 固定效應(yīng)模型。在面板數(shù)據(jù)散點圖中，如果對于不同的截面或不同的時間序列，模型的截距是不同的，則可以采用在模型中加虛擬變量的方法估計回歸參數(shù)，稱此種模型為固定效應(yīng)模型（fixed effects regression mode

19、l）。固定效應(yīng)模型分為3種類型，即個體固定效應(yīng)模型（entity fixed effects regression model）、時刻固定效應(yīng)模型（time fixed effects regression model）和時刻個體固定效應(yīng)模型（time and entity fixed effects regression model）。下面分別介紹。（1）個體固定效應(yīng)模型。個體固定效應(yīng)模型就是對于不同的個體有不同截距的模型。如果對于不同的時間序列（個體）截距是不同的，但是對于不同的橫截面，模型的截距沒有顯著性變化，那么就應(yīng)該建立個體固定效應(yīng)模型，表示如下， yit = b1 xit +g1

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

21、 = 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為標(biāo)量。當(dāng)模型中含有k個解釋變量時，b為k1階列向量。進(jìn)一步寫成矩陣形式，= +b + 上式中的元素1，0都是t1階列向量。 面板數(shù)據(jù)模型用ols方法估計時應(yīng)滿足如下5個假定條件：（1）e(eit|xi1, xi2, , xit, ai) = 0。以xi1, xi2, , xit, ai為條件的eit的期望等于零。（2）(xi1, xi2, , xi

22、t), ( yi1, yi2, , yit), i = 1, 2, , n分別來自于同一個聯(lián)合分布總體，并相互獨立。（3）(xit, eit)具有非零的有限值4階矩。（4）解釋變量之間不存在完全共線性。（5）cov(eit eis|xit,xis, ai) = 0, t s。在固定效應(yīng)模型中隨機誤差項eit在時間上是非自相關(guān)的。其中xit代表一個或多個解釋變量。對模型（1）進(jìn)行ols估計，全部參數(shù)估計量都是無偏的和一致的。模型的自由度是n t 1n。 當(dāng)模型含有k個解釋變量，且n很大，相對較小時，因為模型中含有k + n個被估參數(shù)，一般軟件執(zhí)行ols運算很困難。在計量經(jīng)濟(jì)學(xué)軟件中是采用一種特殊

23、處理方式進(jìn)行ols估計。估計原理是，先用每個變量減其組內(nèi)均值，把數(shù)據(jù)中心化（entity-demeaned），然后用變換的數(shù)據(jù)先估計個體固定效應(yīng)模型的回歸系數(shù)（不包括截距項），然后利用組內(nèi)均值等式計算截距項。這種方法計算起來速度快。具體分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 -) =

24、，上式寫為 = b1+ (6)用ols法估計（1）、（6）式中的b1，結(jié)果是一樣的，但是用（6）式估計，可以減少被估參數(shù)個數(shù)。（2）用ols法估計回歸參數(shù)（不包括截距項，即固定效應(yīng)）。 在k個解釋變量條件下，把用向量形式表示，則利用中心化數(shù)據(jù)，按ols法估計公式計算個體固定效應(yīng)模型中回歸參數(shù)估計量的方差協(xié)方差矩陣估計式如下，() = ()-1 (7)其中=，是相對于的殘差向量。（3）計算回歸模型截距項，即固定效應(yīng)參數(shù)gi。=- (8)以例1（file:panel02）為例得到的個體固定效應(yīng)模型估計結(jié)果如下：注意：個體固定效應(yīng)模型的eviwes輸出結(jié)果中沒有公共截距項。圖12eviwes估計方法

25、：在eviwes的pooled estimation對話框中intercept選項中選fixed effects。其余選項同上。注意：（1）個體固定效應(yīng)模型的eviwes輸出結(jié)果中沒有公共截距項。（2）eviwes輸出結(jié)果中沒有給出描述個體效應(yīng)的截距項相應(yīng)的標(biāo)準(zhǔn)差和t值。不認(rèn)為截距項是模型中的重要參數(shù)。（3）當(dāng)對個體固定效應(yīng)模型選擇加權(quán)估計時，輸出結(jié)果將給出加權(quán)估計和非加權(quán)估計兩種統(tǒng)計量評價結(jié)果。（4）輸出結(jié)果的聯(lián)立方程組形式可以通過點擊view選representations功能獲得。（5）點擊view選wald coefficient tests功能可以對模型的斜率進(jìn)行wald檢驗。（6）

26、點擊view選residuals/table, graphs, covariance matrix, correlation matrix功能可以分別得到按個體計算的殘差序列表，殘差序列圖，殘差序列的方差協(xié)方差矩陣，殘差序列的相關(guān)系數(shù)矩陣。（7）點擊procs選make model功能，將會出現(xiàn)估計結(jié)果的聯(lián)立方程形式，進(jìn)一步點擊solve鍵，在隨后出現(xiàn)的對話框中可以進(jìn)行動態(tài)和靜態(tài)預(yù)測。輸出結(jié)果的方程形式是 = 安徽+ x1t = 479.3 + 0.70 x1t (55.0) = 北京+x2t = 1053.2 + 0.70 x2t (55.0) = 浙江+x15t = 714.2 + 0.7

27、0 x15t (55.0) r2 = 0.99, sser = 2270386, t0.05 (88) = 1.98從結(jié)果看，北京、上海、浙江是消費函數(shù)截距（自發(fā)消費）最大的3個地區(qū)。相對于混合估計模型來說，是否有必要建立個體固定效應(yīng)模型可以通過f檢驗來完成。原假設(shè)h0：不同個體的模型截距項相同（建立混合估計模型）。備擇假設(shè)h1：不同個體的模型截距項不同（建立個體固定效應(yīng)模型）。f統(tǒng)計量定義為：f= (9)其中sser，sseu分別表示約束模型（混合估計模型）和非約束模型（個體固定效應(yīng)模型）的殘差平方和。非約束模型比約束模型多了n-1個被估參數(shù)。（混合估計模型給出公共截距項。）注意：當(dāng)模型中含

28、有k個解釋變量時，f統(tǒng)計量的分母自由度是nt-n-k。 用上例計算，已知sser = 4824588，sseu = 2270386，f= 7.15f0.05(14, 89) = 1.81因為f= 7.15 f0.05(14, 89) = 1.81，所以，拒絕原假設(shè)。結(jié)論是應(yīng)該建立個體固定效應(yīng)模型。（2）時刻固定效應(yīng)模型。時刻固定效應(yīng)模型就是對于不同的截面（時刻點）有不同截距的模型。如果確知對于不同的截面，模型的截距顯著不同，但是對于不同的時間序列（個體）截距是相同的，那么應(yīng)該建立時刻固定效應(yīng)模型，表示如下， yit = b1 xit +a1 + a2 d2 + +at dt +eit, i =

29、 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, ,

30、 n如果滿足上述模型假定條件，對模型（2）進(jìn)行ols估計，全部參數(shù)估計量都具有無偏性和一致性。模型的自由度是n t t-1。圖13eviwes估計方法：在pooled estimation（混合估計）窗口中的dependent variable（相依變量）選擇窗填入cp?；在common coefficients（系數(shù)相同）選擇窗填入ip? 和虛擬變量d1997, d1998, d1999, d2000, d2001, d2002；在cross section specific coefficients（截面系數(shù)不同）選擇窗保持空白；在intercept（截距項）選擇窗點擊common；在we

31、ighting（權(quán)數(shù)）選擇窗點擊no weighting。點擊pooled estimation（混合估計）窗口中的ok鍵。以例1為例得到的時刻固定效應(yīng)模型估計結(jié)果如下： = 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 (9

32、7) = 1.98相對于混合估計模型來說，是否有必要建立時刻固定效應(yīng)模型可以通過f檢驗來完成。h0：對于不同橫截面模型截距項相同（建立混合估計模型）。h1：對于不同橫截面模型的截距項不同（建立時刻固定效應(yīng)模型）。f統(tǒng)計量定義為：f= (11)其中sser，sseu分別表示約束模型（混合估計模型的）和非約束模型（時刻固定效應(yīng)模型的）的殘差平方和。非約束模型比約束模型多了t-1個被估參數(shù)。注意：當(dāng)模型中含有k個解釋變量時，f統(tǒng)計量的分母自由度是nt-t- k。 用上例計算，已知sser= 4824588，sseu= 4028843，f= 3.19f0.05(6, 87) = 2.2因為f= 3.1

33、9 f0.05(14, 89) = 2.2，拒絕原假設(shè)，結(jié)論是應(yīng)該建立時刻固定效應(yīng)模型。（3）時刻個體固定效應(yīng)模型。時刻個體固定效應(yīng)模型就是對于不同的截面（時刻點）、不同的時間序列（個體）都有不同截距的模型。如果確知對于不同的截面、不同的時間序列（個體）模型的截距都顯著地不相同，那么應(yīng)該建立時刻個體效應(yīng)模型，表示如下，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,

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

35、個個體） 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個

36、截面、第n個個體）如果滿足上述模型假定條件，對模型（12）進(jìn)行ols估計，全部參數(shù)估計量都是無偏的和一致的。模型的自由度是n t nt。注意：當(dāng)模型中含有k個解釋變量時，f統(tǒng)計量的分母自由度是nt n -t- k+1。以例1為例得到的截面、時刻固定效應(yīng)模型估計結(jié)果如下：圖14eviwes估計方法：在pooled estimation（混合估計）窗口中的dependent variable（相依變量）選擇窗填入cp?；在common coefficients（系數(shù)相同）選擇窗填入ip? 和虛擬變量d1997, d1998, d1999, d2000, d2001, d2002；在cross se

37、ction specific coefficients（截面系數(shù)不同）選擇窗保持空白；在intercept（截距項）選擇窗中選fixed effects；在weighting（權(quán)數(shù)）選擇窗點擊no weighting。點擊pooled estimation（混合估計）窗口中的ok鍵。注意：（1）對于第1個截面（t=1）eviwes輸出結(jié)果中把(a1 +gi), (i = 1, 2, , n)估計在一起。（2）對于第2, , t個截面（t=1）eviwes輸出結(jié)果中分別把(a1 +at), (t = 2, , t)估計在一起。輸出結(jié)果如下： = 1996 + x11 = 537.9627 + 0

38、.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年浙江?。﹔2 = 0.9932, sser = 2045670, t0.05 (83) = 1.98相對于混合估計模型來說，是否有

39、必要建立時刻個體固定效應(yīng)模型可以通過f檢驗來完成。h0：對于不同橫截面，不同序列，模型截距項都相同（建立混合估計模型）。h1：不同橫截面，不同序列，模型截距項各不相同（建立時刻個體固定效應(yīng)模型）。f統(tǒng)計量定義為：f= (13)其中sser，sseu分別表示約束模型（混合估計模型的）和非約束模型（時刻個體固定效應(yīng)模型的）的殘差平方和。非約束模型比約束模型多了n+t個被估參數(shù)。注意：當(dāng)模型中含有k個解釋變量時，f統(tǒng)計量的分母自由度是nt-n-t- k+1。 用上例計算，已知sser= 4824588，sseu= 2045670，f= 5.6f0.05(20, 81) = 1.64因為f= 5.6

40、f0.05(14, 89) = 1.64，拒絕原假設(shè)，結(jié)論是應(yīng)該建立時刻個體固定效應(yīng)模型。（4）隨機效應(yīng)模型在固定效應(yīng)模型中采用虛擬變量的原因是解釋被解釋變量的信息不夠完整。也可以通過對誤差項的分解來描述這種信息的缺失。 yit = a + b1 xit + eit (14)其中誤差項在時間上和截面上都是相關(guān)的，用3個分量表示如下。eit = ui + vt + wit (15)其中ui n(0, su2)表示截面隨機誤差分量；vt n(0, sv2)表示時間隨機誤差分量；wit n(0, sw2)表示混和隨機誤差分量。同時還假定ui，vt，wit之間互不相關(guān)，各自分別不存在截面自相關(guān)、時間自

41、相關(guān)和混和自相關(guān)。上述模型稱為隨機效應(yīng)模型。隨機效應(yīng)模型和固定效應(yīng)模型比較，相當(dāng)于把固定效應(yīng)模型中的截距項看成兩個隨機變量。一個是截面隨機誤差項（ui），一個是時間隨機誤差項（vt）。如果這兩個隨機誤差項都服從正態(tài)分布，對模型估計時就能夠節(jié)省自由度，因為此條件下只需要估計兩個隨機誤差項的均值和方差。假定固定效應(yīng)模型中的截距項包括了截面隨機誤差項和時間隨機誤差項的平均效應(yīng)，而且對均值的離差分別是ui和vt，固定效應(yīng)模型就變成了隨機效應(yīng)模型。為了容易理解，先假定模型中只存在截面隨機誤差項ui，不存在時間隨機誤差分量（vt）， yit = a + b1 xit + (wit+ ui) = a + b

42、1 xit +eit (16)截面隨機誤差項ui是屬于第個個體的隨機波動分量，并在整個時間范圍(t = 1,2, , t)保持不變。隨機誤差項ui, wit應(yīng)滿足如下條件：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因為根據(jù)上式有eit = wit+ ui所以這種隨機效應(yīng)模型又稱為誤差分量模型（error component model）。有結(jié)論，e(eit ) = e(wit +uj) = 0,(16)式，y

43、it = a + b1 xit + (wit+ ui)，也可以寫成yit = (a + ui) + b1 xit + wit。服從正態(tài)分布的截距項的均值效應(yīng)au被包含在回歸函數(shù)的常數(shù)項中。e(eit 2) = e(wit +uj)2 = sw2 +su2,e(eit eis) = e(wit+ ui)(wis+ ui) = e(wit wis + ui wis + wit ui + ui2) =su2, t s 令ei = (ei1, ei2, eit)則w = e(eiei) = =sw2 i(tt) +su2 1(t1) 1(t1) 其中i(tt)是(tt)階單位陣，1(t1)是(t1)階

44、列向量。因為第i期與j期觀測值是相互獨立的，所以nt個觀測值所對應(yīng)的隨機誤差項的方差與協(xié)方差矩陣v是v = = w = inn w其中inn表示由(t1)階列向量為元素構(gòu)成的單位陣，其中每一個元素1或0都是(t1)階列向量。表示科羅內(nèi)克積（kronecker product）。其運算規(guī)則是 ankb =檢驗個體隨機效應(yīng)的原假設(shè)與檢驗統(tǒng)計量是h0：su2 = 0。（混合估計模型）h1：su2 0。（個體隨機效應(yīng)模型）lm= =其中表示由個體隨機效應(yīng)模型計算的殘差平方和。表示由混合估計模型計算的殘差平方和。統(tǒng)計量lm服從1個自由度的c2分布?？梢詫﹄S機效應(yīng)模型進(jìn)行廣義最小二乘估計。以觀測值方差的倒

45、數(shù)為權(quán)。為了求權(quán)數(shù)，必須采用兩階段最小二乘法估計。因為各隨機誤差分量的方差一般是未知的，第一階段用普通最小二乘估計法對混合數(shù)據(jù)進(jìn)行估計（采用固定效應(yīng)模型）。用估計的殘差計算隨機誤差分量的方差。第二步用這些估計的方差計算參數(shù)的廣義最小二乘估計值。如果隨機誤差分量服從的是正態(tài)分布，模型的參數(shù)還可以用極大似然法估計。仍以例1為例給出隨機效應(yīng)模型估計結(jié)果如下：圖15注意：隨機效應(yīng)模型eviwes輸出結(jié)果中含有公共截距項。圖16 以例1為例，用個體隨機效應(yīng)模型和混合模型計算的統(tǒng)計量的值是lm =8.75(24.4)2 = 5209f0.05 (1) = 3.84因為f= 5209 f0.05 (1) =

46、 3.84，所以拒絕原假設(shè)，結(jié)論是應(yīng)該建立個體隨機效應(yīng)模型。假定截面截距和時間截距都是隨機的。分別服從均值為au和av，方差為su2和sv2的正態(tài)分布。隨機誤差項將由3部分組成，并有方差。var(eit) = var(ui) + var(vt) + var(wit) =su2 +sv2+sw2當(dāng)su2和sv 2都等于零，隨機效應(yīng)模型退化為固定效應(yīng)模型。隨機效應(yīng)模型和固定效應(yīng)模型哪一個更好些？實際是各有優(yōu)缺點。隨機效應(yīng)模型的好處是節(jié)省自由度。對于從時間序列和截面兩方面上看都存在較大變化的數(shù)據(jù)，隨機效應(yīng)模型能明確地描述出誤差來源的特征。固定效應(yīng)模型的好處是很容易分析任意截面數(shù)據(jù)所對應(yīng)的因變量與全部

47、截面數(shù)據(jù)對應(yīng)的因變量均值的差異程度。此外，固定效應(yīng)模型不要求誤差項中的個體效應(yīng)分量與模型中的解釋變量不相關(guān)。當(dāng)然，這一假定不成立時，可能會引起模型參數(shù)估計的不一致性。（5）回歸系數(shù)不同的面板數(shù)據(jù)模型當(dāng)認(rèn)為對于不同個體，解釋變量的回歸系數(shù)存在顯著性差異時，還可以建立回歸系數(shù)不同的面板數(shù)據(jù)模型。eviwes估計方法：在pooled estimation（混合估計）窗口中的dependent variable（相依變量）選擇窗填入cp?；在common coefficients（系數(shù)相同）選擇窗保持空白（如果需要估計時刻固定效應(yīng)也可輸入虛擬變量d1997, d1998, d1999, d2000,

48、d2001, d2002）；在cross section specific coefficients（截面系數(shù)不同）選擇窗填入ip?；在intercept（截距項）選擇窗中選fixed effects（也可以做其他選擇）；在weighting（權(quán)數(shù)）選擇窗點擊no weighting（也可以做其他選擇）。點擊pooled estimation（混合估計）窗口中的ok鍵。圖17 = 安徽+ x1t = 161.62 + 0.76 x1t (9.1) = 北京+x2t = 36.22 + 0.81 x2t (31.0) = 浙江+x15t = 1328.26 + 0.63 x15t (21.1)

49、r2 = 0.995, sser = 1409247用eviwes建立面板數(shù)據(jù)估計模型步驟。利用19962002年15個省級地區(qū)城鎮(zhèn)居民家庭年人均消費性支出和年人均收入數(shù)據(jù)（不變價格數(shù)據(jù)）介紹面板數(shù)據(jù)模型估計步驟。（1）建立混合數(shù)據(jù)庫（pool）對象。首先建立工作文件。在打開工作文件窗口的基礎(chǔ)上，點擊eviwes主功能菜單上的objects鍵，選new object功能（如圖18），從而打開new object（新對象）選擇窗。在type of object選擇區(qū)選擇pool（合并數(shù)據(jù)庫），并在name of object選擇區(qū)為混合數(shù)據(jù)庫起名pool01（初始顯示為untitled）。如圖19，點擊ok鍵，從而打開混合數(shù)據(jù)庫（