(轉(zhuǎn)載)面板數(shù)據(jù)分析簡要步驟與注意事項(面板單位根—面板協(xié)整—回歸分析)

步驟一:分析數(shù)據(jù)的平穩(wěn)性(單位根檢驗)

按照正規(guī)程序,面板數(shù)據(jù)模型在回歸前需檢驗數(shù)據(jù)的平穩(wěn)性。李子奈曾指出,一些非平穩(wěn)的經(jīng)濟時間序列往往表現(xiàn)出共同的變化趨勢,而這些序列間本身不一定有直接的關(guān)聯(lián),此時,對這些數(shù)據(jù)進行回歸,盡管有較高的R平方,但其結(jié)果是沒有任何實際意義的。這種情況稱為稱為虛假回歸或偽回歸(spurious regression)。他認為平穩(wěn)的真正含義是:一個時間序列剔除了不變的均值(可視為截距)和時間趨勢以后,剩余的序列為零均值,同方差,即白噪聲。因此單位根檢驗時有三種檢驗?zāi)J剑杭扔汹厔萦钟薪鼐唷⒅挥薪鼐唷⒁陨隙紵o。

因此為了避免偽回歸,確保估計結(jié)果的有效性,我們必須對各面板序列的平穩(wěn)性進行檢驗。而檢驗數(shù)據(jù)平穩(wěn)性最常用的辦法就是單位根檢驗。首先,我們可以先對面板序列繪制時序圖,以粗略觀測時序圖中由各個觀測值描出代表變量的折線是否含有趨勢項和(或)截距項,從而為進一步的單位根檢驗的檢驗?zāi)J阶鰷蕚洹?/p>

單位根檢驗方法的文獻綜述:在非平穩(wěn)的面板數(shù)據(jù)漸進過程中,Levin andLin(1993) 很早就發(fā)現(xiàn)這些估計量的極限分布是高斯分布,這些結(jié)果也被應(yīng)用在有異方差的面板數(shù)據(jù)中,并建立了對面板單位根進行檢驗的早期版本。后來經(jīng)過Levin et al. (2002)的改進,提出了檢驗面板單位根的LLC 法。Levin et al. (2002) 指出,該方法允許不同截距和時間趨勢,異方差和高階序列相關(guān),適合于中等維度(時間序列介于25~250 之間,截面數(shù)介于10~250 之間) 的面板單位根檢驗。Im et al. (1997) 還提出了檢驗面板單位根的IPS 法,但Breitung(2000) 發(fā)現(xiàn)IPS 法對限定性趨勢的設(shè)定極為敏感,并提出了面板單位根檢驗的Breitung 法。Maddala and Wu(1999)又提出了ADF-Fisher和PP-Fisher面板單位根檢驗方法。

由上述綜述可知,可以使用LLC、IPS、Breintung、ADF-Fisher 和PP-Fisher5種方法進行面板單位根檢驗。

其中LLC-T 、BR-T、IPS-W 、ADF-FCS、PP-FCS 、H-Z 分別指Levin, Lin & Chu t* 統(tǒng)計量、Breitung t 統(tǒng)計量、lm Pesaran & Shin W 統(tǒng)計量、ADF- Fisher Chi-square統(tǒng)計量、PP-Fisher Chi-square統(tǒng)計量、Hadri Z統(tǒng)計量,并且Levin, Lin & Chu t* 統(tǒng)計量、Breitung t統(tǒng)計量的原假設(shè)為存在普通的單位根過程,lm Pesaran & Shin W 統(tǒng)計量、ADF- Fisher Chi-square統(tǒng)計量、PP-Fisher Chi-square統(tǒng)計量的原假設(shè)為存在有效的單位根過程, Hadri Z統(tǒng)計量的檢驗原假設(shè)為不存在普通的單位根過程。

有時,為了方便,只采用兩種面板數(shù)據(jù)單位根檢驗方法,即相同根單位根檢驗LLC(Levin-Lin-Chu)檢驗和不同根單位根檢驗Fisher-ADF檢驗(注:對普通序列(非面板序列)的單位根檢驗方法則常用ADF檢驗),如果在兩種檢驗中均拒絕存在單位根的原假設(shè)則我們說此序列是平穩(wěn)的,反之則不平穩(wěn)

如果我們以T(trend)代表序列含趨勢項,以I(intercept)代表序列含截距項,T&I代表兩項都含,N(none)代表兩項都不含,那么我們可以基于前面時序圖得出的結(jié)論,在單位根檢驗中選擇相應(yīng)檢驗?zāi)J健?/p>

但基于時序圖得出的結(jié)論畢竟是粗略的,嚴格來說,那些檢驗結(jié)構(gòu)均需一一檢驗。具體操作可以參照李子奈的說法:ADF檢驗是通過三個模型來完成,首先從含有截距和趨勢項的模型開始,再檢驗只含截距項的模型,最后檢驗二者都不含的模型。并且認為,只有三個模型的檢驗結(jié)果都不能拒絕原假設(shè)時,我們才認為時間序列是非平穩(wěn)的,而只要其中有一個模型的檢驗結(jié)果拒絕了零假設(shè),就可認為時間序列是平穩(wěn)的。

此外,單位根檢驗一般是先從水平(level)序列開始檢驗起,如果存在單位根,則對該序列進行一階差分后繼續(xù)檢驗,若仍存在單位根,則進行二階甚至高階差分后檢驗,直至序列平穩(wěn)為止。我們記I(0)為零階單整,I(1)為一階單整,依次類推,I(N)為N階單整。

步驟二:協(xié)整檢驗或模型修正

情況一:如果基于單位根檢驗的結(jié)果發(fā)現(xiàn)變量之間是同階單整的,那么我們可以進行協(xié)整檢驗。協(xié)整檢驗是考察變量間長期均衡關(guān)系的方法。所謂的協(xié)整是指若兩個或多個非平穩(wěn)的變量序列,其某個線性組合后的序列呈平穩(wěn)性。此時我們稱這些變量序列間有協(xié)整關(guān)系存在。因此協(xié)整的要求或前提是同階單整。

但也有如下的寬限說法:如果變量個數(shù)多于兩個,即解釋變量個數(shù)多于一個,被解釋變量的單整階數(shù)不能高于任何一個解釋變量的單整階數(shù)。另當解釋變量的單整階數(shù)高于被解釋變量的單整階數(shù)時,則必須至少有兩個解釋變量的單整階數(shù)高于被解釋變量的單整階數(shù)。如果只含有兩個解釋變量,則兩個變量的單整階數(shù)應(yīng)該相同。

也就是說,單整階數(shù)不同的兩個或以上的非平穩(wěn)序列如果一起進行協(xié)整檢驗,必然有某些低階單整的,即波動相對高階序列的波動甚微弱(有可能波動幅度也不同)的序列,對協(xié)整結(jié)果的影響不大,因此包不包含的重要性不大。而相對處于最高階序列,由于其波動較大,對回歸殘差的平穩(wěn)性帶來極大的影響,所以如果協(xié)整是包含有某些高階單整序列的話(但如果所有變量都是階數(shù)相同的高階,此時也被稱作同階單整,這樣的話另當別論),一定不能將其納入?yún)f(xié)整檢驗。

協(xié)整檢驗方法的文獻綜述:(1)Kao(1999)、Kao and Chiang(2000)利用推廣的DF和ADF檢驗提出了檢驗面板協(xié)整的方法,這種方法零假設(shè)是沒有協(xié)整關(guān)系,并且利用靜態(tài)面板回歸的殘差來構(gòu)建統(tǒng)計量。(2)Pedron(1999)在零假設(shè)是在動態(tài)多元面板回歸中沒有協(xié)整關(guān)系的條件下給出了七種基于殘差的面板協(xié)整檢驗方法。和Kao的方法不同的是,Pedroni的檢驗方法允許異質(zhì)面板的存在。(3)Larsson et al(2001)發(fā)展了基于Johansen(1995)向量自回歸的似然檢驗的面板協(xié)整檢驗方法,這種檢驗的方法是檢驗變量存在共同的協(xié)整的秩。

主要采用的是Pedroni、Kao、Johansen的方法。

通過了協(xié)整檢驗,說明變量之間存在著長期穩(wěn)定的均衡關(guān)系,其方程回歸殘差是平穩(wěn)的。因此可以在此基礎(chǔ)上直接對原方程進行回歸,此時的回歸結(jié)果是較精確的。

這時,我們或許還想進一步對面板數(shù)據(jù)做格蘭杰因果檢驗(因果檢驗的前提是變量協(xié)整)。但如果變量之間不是協(xié)整(即非同階單整)的話,是不能進行格蘭杰因果檢驗的,不過此時可以先對數(shù)據(jù)進行處理。引用張曉峒的原話,“如果y和x不同階,不能做格蘭杰因果檢驗,但可通過差分序列或其他處理得到同階單整序列,并且要看它們此時有無經(jīng)濟意義。”

下面簡要介紹一下因果檢驗的含義:這里的因果關(guān)系是從統(tǒng)計角度而言的,即是通過概率或者分布函數(shù)的角度體現(xiàn)出來的:在所有其它事件的發(fā)生情況固定不變的條件下,如果一個事件X的發(fā)生與不發(fā)生對于另一個事件Y的發(fā)生的概率(如果通過事件定義了隨機變量那么也可以說分布函數(shù))有影響,并且這兩個事件在時間上又有先后順序(A前B后),那么我們便可以說X是Y的原因。考慮最簡單的形式,Granger檢驗是運用F-統(tǒng)計量來檢驗X的滯后值是否顯著影響Y(在統(tǒng)計的意義下,且已經(jīng)綜合考慮了Y的滯后值;如果影響不顯著,那么稱X不是Y的“Granger原因”(Granger cause);如果影響顯著,那么稱X是Y的“Granger原因”。同樣,這也可以用于檢驗Y是X的“原因”,檢驗Y的滯后值是否影響X(已經(jīng)考慮了X的滯后對X自身的影響)。

Eviews好像沒有在POOL窗口中提供Granger causality test,而只有unit root test和cointegration test。說明Eviews是無法對面板數(shù)據(jù)序列做格蘭杰檢驗的,格蘭杰檢驗只能針對序列組做。也就是說格蘭杰因果檢驗在Eviews中是針對普通的序列對(pairwise)而言的。你如果想對面板數(shù)據(jù)中的某些合成序列做因果檢驗的話,不妨先導出相關(guān)序列到一個組中(POOL窗口中的Proc/Make Group),再來試試。

情況二:如果如果基于單位根檢驗的結(jié)果發(fā)現(xiàn)變量之間是非同階單整的,即面板數(shù)據(jù)中有些序列平穩(wěn)而有些序列不平穩(wěn),此時不能進行協(xié)整檢驗與直接對原序列進行回歸。但此時也不要著急,我們可以在保持變量經(jīng)濟意義的前提下,對我們前面提出的模型進行修正,以消除數(shù)據(jù)不平穩(wěn)對回歸造成的不利影響。如差分某些序列,將基于時間頻度的絕對數(shù)據(jù)變成時間頻度下的變動數(shù)據(jù)或增長率數(shù)據(jù)。此時的研究轉(zhuǎn)向新的模型,但要保證模型具有經(jīng)濟意義。因此一般不要對原序列進行二階差分,因為對變動數(shù)據(jù)或增長率數(shù)據(jù)再進行差分,我們不好對其冠以經(jīng)濟解釋。難道你稱其為變動率的變動率?

步驟三:面板模型的選擇與回歸

面板數(shù)據(jù)模型的選擇通常有三種形式:

一種是混合估計模型(Pooled Regression Model)。如果從時間上看,不同個體之間不存在顯著性差異;從截面上看,不同截面之間也不存在顯著性差異,那么就可以直接把面板數(shù)據(jù)混合在一起用普通最小二乘法(OLS)估計參數(shù)。一種是固定效應(yīng)模型(Fixed Effects Regression Model)。如果對于不同的截面或不同的時間序列,模型的截距不同,則可以采用在模型中添加虛擬變量的方法估計回歸參數(shù)。一種是隨機效應(yīng)模型(Random Effects Regression Model)。如果固定效應(yīng)模型中的截距項包括了截面隨機誤差項和時間隨機誤差項的平均效應(yīng),并且這兩個隨機誤差項都服從正態(tài)分布,則固定效應(yīng)模型就變成了隨機效應(yīng)模型。

在面板數(shù)據(jù)模型形式的選擇方法上,我們經(jīng)常采用F檢驗決定選用混合模型還是固定效應(yīng)模型,然后用Hausman檢驗確定應(yīng)該建立隨機效應(yīng)模型還是固定效應(yīng)模型。

檢驗完畢后,我們也就知道該選用哪種模型了,然后我們就開始回歸:

在回歸的時候,權(quán)數(shù)可以選擇按截面加權(quán)(cross-section weights)的方式,對于橫截面?zhèn)€數(shù)大于時序個數(shù)的情況更應(yīng)如此,表示允許不同的截面存在異方差現(xiàn)象。估計方法采用PCSE(Panel Corrected Standard Errors,面板校正標準誤)方法。Beck和Katz(1995)引入的PCSE估計方法是面板數(shù)據(jù)模型估計方法的一個創(chuàng)新,可以有效的處理復(fù)雜的面板誤差結(jié)構(gòu),如同步相關(guān),異方差,序列相關(guān)等,在樣本量不夠大時尤為有用。

最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
平臺聲明:文章內(nèi)容(如有圖片或視頻亦包括在內(nèi))由作者上傳并發(fā)布,文章內(nèi)容僅代表作者本人觀點,簡書系信息發(fā)布平臺,僅提供信息存儲服務(wù)。
  • 序言:七十年代末,一起剝皮案震驚了整個濱河市,隨后出現(xiàn)的幾起案子,更是在濱河造成了極大的恐慌,老刑警劉巖,帶你破解...
    沈念sama閱讀 227,882評論 6 531
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件,死亡現(xiàn)場離奇詭異,居然都是意外死亡,警方通過查閱死者的電腦和手機,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 98,208評論 3 414
  • 文/潘曉璐 我一進店門,熙熙樓的掌柜王于貴愁眉苦臉地迎上來,“玉大人,你說我怎么就攤上這事。” “怎么了?”我有些...
    開封第一講書人閱讀 175,746評論 0 373
  • 文/不壞的土叔 我叫張陵,是天一觀的道長。 經(jīng)常有香客問我,道長,這世上最難降的妖魔是什么? 我笑而不...
    開封第一講書人閱讀 62,666評論 1 309
  • 正文 為了忘掉前任,我火速辦了婚禮,結(jié)果婚禮上,老公的妹妹穿的比我還像新娘。我一直安慰自己,他們只是感情好,可當我...
    茶點故事閱讀 71,477評論 6 407
  • 文/花漫 我一把揭開白布。 她就那樣靜靜地躺著,像睡著了一般。 火紅的嫁衣襯著肌膚如雪。 梳的紋絲不亂的頭發(fā)上,一...
    開封第一講書人閱讀 54,960評論 1 321
  • 那天,我揣著相機與錄音,去河邊找鬼。 笑死,一個胖子當著我的面吹牛,可吹牛的內(nèi)容都是我干的。 我是一名探鬼主播,決...
    沈念sama閱讀 43,047評論 3 440
  • 文/蒼蘭香墨 我猛地睜開眼,長吁一口氣:“原來是場噩夢啊……” “哼!你這毒婦竟也來了?” 一聲冷哼從身側(cè)響起,我...
    開封第一講書人閱讀 42,200評論 0 288
  • 序言:老撾萬榮一對情侶失蹤,失蹤者是張志新(化名)和其女友劉穎,沒想到半個月后,有當?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體,經(jīng)...
    沈念sama閱讀 48,726評論 1 333
  • 正文 獨居荒郊野嶺守林人離奇死亡,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點故事閱讀 40,617評論 3 354
  • 正文 我和宋清朗相戀三年,在試婚紗的時候發(fā)現(xiàn)自己被綠了。 大學時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片。...
    茶點故事閱讀 42,807評論 1 369
  • 序言:一個原本活蹦亂跳的男人離奇死亡,死狀恐怖,靈堂內(nèi)的尸體忽然破棺而出,到底是詐尸還是另有隱情,我是刑警寧澤,帶...
    沈念sama閱讀 38,327評論 5 358
  • 正文 年R本政府宣布,位于F島的核電站,受9級特大地震影響,放射性物質(zhì)發(fā)生泄漏。R本人自食惡果不足惜,卻給世界環(huán)境...
    茶點故事閱讀 44,049評論 3 347
  • 文/蒙蒙 一、第九天 我趴在偏房一處隱蔽的房頂上張望。 院中可真熱鬧,春花似錦、人聲如沸。這莊子的主人今日做“春日...
    開封第一講書人閱讀 34,425評論 0 26
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽。三九已至,卻和暖如春,著一層夾襖步出監(jiān)牢的瞬間,已是汗流浹背。 一陣腳步聲響...
    開封第一講書人閱讀 35,674評論 1 281
  • 我被黑心中介騙來泰國打工, 沒想到剛下飛機就差點兒被人妖公主榨干…… 1. 我叫王不留,地道東北人。 一個月前我還...
    沈念sama閱讀 51,432評論 3 390
  • 正文 我出身青樓,卻偏偏與公主長得像,于是被迫代替她去往敵國和親。 傳聞我的和親對象是個殘疾皇子,可洞房花燭夜當晚...
    茶點故事閱讀 47,769評論 2 372

推薦閱讀更多精彩內(nèi)容