1、協(xié)整檢驗(yàn),協(xié)整與誤差修正模型,一、長(zhǎng)期均衡與協(xié)整分析二、協(xié)整檢驗(yàn)—EG檢驗(yàn)三、協(xié)整檢驗(yàn)—JJ檢驗(yàn)四、誤差修正模型,一、長(zhǎng)期均衡與協(xié)整分析Equilibrium and Cointegration,1、問題的提出,經(jīng)典回歸模型（classical regression model）是建立在平穩(wěn)數(shù)據(jù)變量基礎(chǔ)上的，對(duì)于非平穩(wěn)變量，不能使用經(jīng)典回歸模型，否則會(huì)出現(xiàn)虛假回歸等諸多問題。由于許多經(jīng)濟(jì)變量是非平穩(wěn)的，這就給經(jīng)典的回歸分析方法

2、帶來了很大限制。但是，如果變量之間有著長(zhǎng)期的穩(wěn)定關(guān)系，即它們之間是協(xié)整的（cointegration)，則是可以使用經(jīng)典回歸模型方法建立回歸模型的。例如，中國(guó)居民人均消費(fèi)水平與人均GDP變量的例子, 從經(jīng)濟(jì)理論上說，人均GDP決定著居民人均消費(fèi)水平，它們之間有著長(zhǎng)期的穩(wěn)定關(guān)系，即它們之間是協(xié)整的。,經(jīng)濟(jì)理論指出，某些經(jīng)濟(jì)變量間確實(shí)存在著長(zhǎng)期均衡關(guān)系，這種均衡關(guān)系意味著經(jīng)濟(jì)系統(tǒng)不存在破壞均衡的內(nèi)在機(jī)制，如果變量在某時(shí)期受到干擾后偏離其

3、長(zhǎng)期均衡點(diǎn)，則均衡機(jī)制將會(huì)在下一期進(jìn)行調(diào)整以使其重新回到均衡狀態(tài)。 假設(shè)X與Y間的長(zhǎng)期“均衡關(guān)系”由式描述,2、長(zhǎng)期均衡,該均衡關(guān)系意味著:給定X的一個(gè)值，Y相應(yīng)的均衡值也隨之確定為??0+?1X。,在t-1期末，存在下述三種情形之一：Y等于它的均衡值：Yt-1= ?0+?1Xt ；Y小于它的均衡值：Yt-1 ?0+?1Xt ；,在時(shí)期t，假設(shè)X有一個(gè)變化量?Xt，如果變量X與Y在時(shí)期t與t-1末期仍滿足它們間的長(zhǎng)期均衡關(guān)

4、系，即上述第一種情況，則Y的相應(yīng)變化量為:,vt=?t-?t-1,如果t-1期末，發(fā)生了上述第二種情況，即Y的值小于其均衡值，則t期末Y的變化往往會(huì)比第一種情形下Y的變化大一些；反之，如果t-1期末Y的值大于其均衡值，則t期末Y的變化往往會(huì)小于第一種情形下的?Yt ?？梢姡绻鸜t=?0+?1Xt+?t正確地提示了X與Y間的長(zhǎng)期穩(wěn)定的“均衡關(guān)系”，則意味著Y對(duì)其均衡點(diǎn)的偏離從本質(zhì)上說是“臨時(shí)性”的。一個(gè)重要的假設(shè)就是:隨機(jī)擾動(dòng)項(xiàng)?

5、t必須是平穩(wěn)序列。如果?t有隨機(jī)性趨勢(shì)（上升或下降），則會(huì)導(dǎo)致Y對(duì)其均衡點(diǎn)的任何偏離都會(huì)被長(zhǎng)期累積下來而不能被消除。,式Y(jié)t=?0+?1Xt+?t中的隨機(jī)擾動(dòng)項(xiàng)也被稱為非均衡誤差（disequilibrium error），它是變量X與Y的一個(gè)線性組合：,如果X與Y間的長(zhǎng)期均衡關(guān)系正確，該式表述的非均衡誤差應(yīng)是一平穩(wěn)時(shí)間序列，并且具有零期望值，即是具有0均值的I(0)序列。,非穩(wěn)定的時(shí)間序列，它們的線性組合也可能成為平穩(wěn)的。稱變量X與Y

6、是協(xié)整的（cointegrated）。,3、協(xié)整,如果序列{X1t,X2t,…,Xkt}都是d階單整，存在向量?=(?1,?2,…,?k)，使得Zt=?XT ~ I(d-b)， 其中，b>0，X=(X1t,X2t,…,Xkt)T，則認(rèn)為序列{X1t,X2t,…,Xkt}是(d,b)階協(xié)整，記為Xt~CI(d,b)，?為協(xié)整向量（cointegrated vector）。如果兩個(gè)變量都是單整變量，只有當(dāng)它們的單整階數(shù)相同時(shí)，才可

7、能協(xié)整；如果它們的單整階數(shù)不相同，就不可能協(xié)整。,3個(gè)以上的變量，如果具有不同的單整階數(shù)，有可能經(jīng)過線性組合構(gòu)成低階單整變量。,,,（d,d）階協(xié)整是一類非常重要的協(xié)整關(guān)系，它的經(jīng)濟(jì)意義在于：兩個(gè)變量，雖然它們具有各自的長(zhǎng)期波動(dòng)規(guī)律，但是如果它們是（d,d）階協(xié)整的，則它們之間存在著一個(gè)長(zhǎng)期穩(wěn)定的比例關(guān)系。例如，中國(guó)CPC和GDPPC，它們各自都是2階單整，如果它們是(2,2)階協(xié)整，說明它們之間存在著一個(gè)長(zhǎng)期穩(wěn)定的比例關(guān)系，從計(jì)量經(jīng)

8、濟(jì)學(xué)模型的意義上講，建立如下居民人均消費(fèi)函數(shù)模型是合理的。,盡管兩個(gè)時(shí)間序列是非平穩(wěn)的，也可以用經(jīng)典的回歸分析方法建立回歸模型。,從這里，我們已經(jīng)初步認(rèn)識(shí)到：檢驗(yàn)變量之間的協(xié)整關(guān)系，在建立計(jì)量經(jīng)濟(jì)學(xué)模型中是非常重要的。 而且，從變量之間是否具有協(xié)整關(guān)系出發(fā)選擇模型的變量，其數(shù)據(jù)基礎(chǔ)是牢固的，其統(tǒng)計(jì)性質(zhì)是優(yōu)良的。,二、協(xié)整檢驗(yàn)—EG檢驗(yàn),1、兩變量的Engle-Granger檢驗(yàn),為了檢驗(yàn)兩變量Yt,Xt是否為協(xié)整，Engle和G

9、ranger于1987年提出兩步檢驗(yàn)法，也稱為EG檢驗(yàn)。 第一步，用OLS方法估計(jì)方程 Yt=?0+?1Xt+?t并計(jì)算非均衡誤差，得到：,稱為協(xié)整回歸(cointegrating)或靜態(tài)回歸(static regression)。,非均衡誤差的單整性的檢驗(yàn)方法仍然是DF檢驗(yàn)或者ADF檢驗(yàn)。需要注意是，這里的DF或ADF檢驗(yàn)是針對(duì)協(xié)整回歸計(jì)算出的誤差項(xiàng)，而非真正的非均衡誤差。而OLS法采用了殘差最小平方和原理，因此估計(jì)量?

10、是向下偏倚的，這樣將導(dǎo)致拒絕零假設(shè)的機(jī)會(huì)比實(shí)際情形大。于是對(duì)et平穩(wěn)性檢驗(yàn)的DF與ADF臨界值應(yīng)該比正常的DF與ADF臨界值還要小。,MacKinnon(1991)通過模擬試驗(yàn)給出了協(xié)整檢驗(yàn)的臨界值。,例 檢驗(yàn)中國(guó)居民人均消費(fèi)水平CPC與人均國(guó)內(nèi)生產(chǎn)總值GDPPC的協(xié)整關(guān)系。,已知CPC與GDPPC都是I(2)序列，已知它們的回歸式,R2=0.9981,對(duì)該式計(jì)算的殘差序列作ADF檢驗(yàn)，適當(dāng)檢驗(yàn)?zāi)Ｐ蜑椋?（-4.47） (3

11、.93) (3.05) LM(1)=0.00 LM(2)=0.00,t=-4.47<-3.75=ADF0.05，拒絕存在單位根的假設(shè)，殘差項(xiàng)是平穩(wěn)的。因此中國(guó)居民人均消費(fèi)水平與人均GDP是(2,2)階協(xié)整的，說明了該兩變量間存在長(zhǎng)期穩(wěn)定的“均衡”關(guān)系。,2、多變量協(xié)整關(guān)系的檢驗(yàn)—擴(kuò)展的E-G檢驗(yàn),多變量協(xié)整關(guān)系的檢驗(yàn)要比雙變量復(fù)雜一些，主要在于協(xié)整變量

12、間可能存在多種穩(wěn)定的線性組合。 假設(shè)有4個(gè)I(1)變量Z、X、Y、W，它們有如下的長(zhǎng)期均衡關(guān)系：,非均衡誤差項(xiàng)?t應(yīng)是I(0)序列：,然而，如果Z與W，X與Y間分別存在長(zhǎng)期均衡關(guān)系：,則非均衡誤差項(xiàng)v1t、v2t一定是穩(wěn)定序列I(0)。于是它們的任意線性組合也是穩(wěn)定的。例如,由于vt象?t一樣，也是Z、X、Y、W四個(gè)變量的線性組合，由此vt 式也成為該四變量的另一穩(wěn)定線性組合。 （1, -?0,-?1,-?2,-?3）是對(duì)應(yīng)于

13、?t 式的協(xié)整向量，（1,-?0-?0,-?1,1,-?1）是對(duì)應(yīng)于vt式的協(xié)整向量。,一定是I(0)序列。,檢驗(yàn)程序：對(duì)于多變量的協(xié)整檢驗(yàn)過程，基本與雙變量情形相同，即需檢驗(yàn)變量是否具有同階單整性，以及是否存在穩(wěn)定的線性組合。在檢驗(yàn)是否存在穩(wěn)定的線性組合時(shí)，需通過設(shè)置一個(gè)變量為被解釋變量，其他變量為解釋變量，進(jìn)行OLS估計(jì)并檢驗(yàn)殘差序列是否平穩(wěn)。如果不平穩(wěn)，則需更換被解釋變量，進(jìn)行同樣的OLS估計(jì)及相應(yīng)的殘差項(xiàng)檢驗(yàn)。當(dāng)所有的變

14、量都被作為被解釋變量檢驗(yàn)之后，仍不能得到平穩(wěn)的殘差項(xiàng)序列，則認(rèn)為這些變量間不存在（d,d）階協(xié)整。,檢驗(yàn)殘差項(xiàng)是否平穩(wěn)的DF與ADF檢驗(yàn)臨界值要比通常的DF與ADF檢驗(yàn)臨界值小，而且該臨界值還受到所檢驗(yàn)的變量個(gè)數(shù)的影響。,MacKinnon(1991)通過模擬試驗(yàn)得到的不同變量協(xié)整檢驗(yàn)的臨界值。,三、協(xié)整檢驗(yàn)—JJ檢驗(yàn),⒈ JJ檢驗(yàn)的原理,Johansen于1988年，以及與Juselius一起于1990年提出了一種用向量自回歸模型進(jìn)行

15、檢驗(yàn)的方法，通常稱為Johansen檢驗(yàn)，或JJ檢驗(yàn)，是一種進(jìn)行多重I(1)序列協(xié)整檢驗(yàn)的較好方法。,沒有移動(dòng)平均項(xiàng)的向量自回歸模型表示為：,,,,,,差分Yt為M個(gè)I(1)過程構(gòu)成的向量,I(0)過程,I(0)過程,只有產(chǎn)生協(xié)整，才能保證新生誤差是平穩(wěn)過程,將y的協(xié)整問題轉(zhuǎn)變?yōu)橛懻摼仃嚘暗男再|(zhì)問題,,,于是，將yt中的協(xié)整檢驗(yàn)變成對(duì)矩陣Π的分析問題。這就是JJ檢驗(yàn)的基本原理。,兩種檢驗(yàn)方法：特征值軌跡檢驗(yàn)最大特征值檢驗(yàn)☆,⒉ J

16、J檢驗(yàn)的預(yù)備工作,第一步：用OLS分別估計(jì)下式中的每一個(gè)方程，計(jì)算殘差，得到殘差矩陣S0，為一個(gè)(M×T)階矩陣。,,,第一步：用OLS分別估計(jì)下式中的每一個(gè)方程，計(jì)算殘差，得到殘差矩陣S1，也為一個(gè)(M×T)階矩陣。,,,,第三步：構(gòu)造上述殘差矩陣的積矩陣：,,,,,第四步：計(jì)算有序特征值和特征向量。,第五步：設(shè)定似然函數(shù)。,⒊ JJ檢驗(yàn)之一—特征值軌跡檢驗(yàn),服從Johansen分布。被稱為特征值軌跡統(tǒng)計(jì)量。,,…

17、，一直檢驗(yàn)下去，直到出現(xiàn)第一個(gè)不顯著的η(M－r)為止，說明存在r個(gè)協(xié)整向量。這r個(gè)協(xié)整向量就是對(duì)應(yīng)于最大的r個(gè)特征值的經(jīng)過正規(guī)化的特征向量。,⒋ JJ檢驗(yàn)之一—最大特征值檢驗(yàn),該統(tǒng)計(jì)量被稱為最大特征值統(tǒng)計(jì)量。于是該檢驗(yàn)被稱為最大特征值檢驗(yàn)。,,由 Johansen和Juselius于1990年計(jì)算得到 Johansen分布臨界值表。,⒌JJ檢驗(yàn)實(shí)例,GDP、CONSR、CONSP、INV取對(duì)數(shù)后為I(1)序列。即lnGDP、lnCON

18、SR、lnCONSP、lnINV。對(duì)它們之間的協(xié)整關(guān)系進(jìn)行檢驗(yàn)。,兩種方法的結(jié)論是一致的。,如何處理高階單整序列？,從理論上講。JJ 檢驗(yàn)只適用于多個(gè)1階單整序列。多個(gè)同階高階單整序列，差分為1階后再檢驗(yàn)，顯然是可行的。但是意義發(fā)生變化。沒有看到關(guān)于高階多重協(xié)整檢驗(yàn)的文獻(xiàn)，難度太大。能否先檢驗(yàn)，然后建立均衡方程，通過對(duì)誤差項(xiàng)的單位根檢驗(yàn)以判斷發(fā)生何種協(xié)整？未見經(jīng)典。,如何選擇截距和時(shí)間趨勢(shì)項(xiàng)？,分別考慮CE和VAR中是否有截距和

19、時(shí)間趨勢(shì)項(xiàng)作為假設(shè)顯著性檢驗(yàn)重新檢驗(yàn)對(duì)協(xié)整關(guān)系檢驗(yàn)結(jié)果無顯著影響（檢驗(yàn)統(tǒng)計(jì)量發(fā)生變化，但臨界值同時(shí)發(fā)生變化）,如何在多個(gè)協(xié)整關(guān)系中作出選擇？,一般選擇對(duì)應(yīng)于最大特征值的第1個(gè)協(xié)整關(guān)系從應(yīng)用的目的出發(fā)選擇,四、誤差修正模型Error Correction Model, ECM,1、一般差分模型的問題,對(duì)于非穩(wěn)定時(shí)間序列，可通過差分的方法將其化為穩(wěn)定序列，然后才可建立經(jīng)典的回歸分析模型。,模型只表達(dá)了X與Y間的短期關(guān)系，而沒有揭

20、示它們間的長(zhǎng)期關(guān)系。關(guān)于變量水平值的重要信息將被忽略。,誤差項(xiàng)?t不存在序列相關(guān)， ?t是一個(gè)一階移動(dòng)平均時(shí)間序列，因而是序列相關(guān)的。,2、誤差修正模型,是一種具有特定形式的計(jì)量經(jīng)濟(jì)學(xué)模型，它的主要形式是由Davidson、 Hendry、Srba和Yeo于1978年提出的，稱為DHSY模型。,,,由于現(xiàn)實(shí)經(jīng)濟(jì)中很少處在均衡點(diǎn)上，假設(shè)具有（1, 1）階分布滯后形式,Y的變化決定于X的變化以及前一時(shí)期的非均衡程度。一階誤差修正模型(fi

21、rst-order error correction model)的形式：,若(t-1)時(shí)刻Y大于其長(zhǎng)期均衡解?0+?1X，ecm為正，則(-?ecm)為負(fù)，使得?Yt減少；若(t-1)時(shí)刻Y小于其長(zhǎng)期均衡解?0+?1X ，ecm為負(fù)，則(-?ecm)為正，使得?Yt增大。體現(xiàn)了長(zhǎng)期非均衡誤差對(duì)短期變化的控制。,復(fù)雜的ECM形式，例如：,,,誤差修正模型的優(yōu)點(diǎn)：如： a）一階差分項(xiàng)的使用消除了變量可能存在的趨勢(shì)因素，從而

22、避免了虛假回歸問題； b）一階差分項(xiàng)的使用也消除模型可能存在的多重共線性問題； c）誤差修正項(xiàng)的引入保證了變量水平值的信息沒有被忽視； d）由于誤差修正項(xiàng)本身的平穩(wěn)性，使得該模型可以用經(jīng)典的回歸方法進(jìn)行估計(jì)，尤其是模型中差分項(xiàng)可以使用通常的t檢驗(yàn)與F檢驗(yàn)來進(jìn)行選取；等等。,3、誤差修正模型的建立,Granger 表述定理（Granger representaion theorem） Engl

23、e 與 Granger 1987年提出 如果變量X與Y是協(xié)整的，則它們間的短期非均衡關(guān)系總能由一個(gè)誤差修正模型表述。,模型中沒有明確指出Y與X的滯后項(xiàng)數(shù)，可以是多階滯后；由于一階差分項(xiàng)是I(0)變量，因此模型中允許采用X的非滯后差分項(xiàng)?Xt 。,建立誤差修正模型，需要： 首先對(duì)變量進(jìn)行協(xié)整分析，以發(fā)現(xiàn)變量之間的協(xié)整關(guān)系，即長(zhǎng)期均衡關(guān)系，并以這種關(guān)系構(gòu)成誤差修正項(xiàng)。 然后建立短期模型，將誤差修正項(xiàng)看作一個(gè)解釋變量，連

24、同其它反映短期波動(dòng)的解釋變量一起，建立短期模型，即誤差修正模型。,Engle-Granger兩步法 第一步，進(jìn)行協(xié)整回歸（OLS法），檢驗(yàn)變量間的協(xié)整關(guān)系，估計(jì)協(xié)整向量（長(zhǎng)期均衡關(guān)系參數(shù)）； 第二步，若協(xié)整性存在，則以第一步求到的殘差作為非均衡誤差項(xiàng)加入到誤差修正模型中，并用OLS法估計(jì)相應(yīng)參數(shù)。 需要注意的是：在進(jìn)行變量間的協(xié)整檢驗(yàn)時(shí)，如有必要可在協(xié)整回歸式中加入趨勢(shì)項(xiàng)，這時(shí)，對(duì)殘差項(xiàng)的穩(wěn)定性檢驗(yàn)就無須再設(shè)趨勢(shì)

25、項(xiàng)。 另外，第二步中變量差分滯后項(xiàng)的多少，可以殘差項(xiàng)序列是否存在自相關(guān)性來判斷，如果存在自相關(guān)，則應(yīng)加入變量差分的滯后項(xiàng)。,經(jīng)濟(jì)理論指出，居民消費(fèi)支出是其實(shí)際收入的函數(shù)。 以中國(guó)國(guó)民核算中的居民消費(fèi)支出經(jīng)過居民消費(fèi)價(jià)格指數(shù)縮減得到中國(guó)居民實(shí)際消費(fèi)支出時(shí)間序列（C）； 以支出法GDP對(duì)居民消費(fèi)價(jià)格指數(shù)縮減近似地代表國(guó)民收入時(shí)間序列(GDP)。 時(shí)間段為1978~2000（表9.3.3）,例 中國(guó)

26、居民消費(fèi)的誤差修正模型 ☆,（1）對(duì)數(shù)據(jù)lnC與lnGDP進(jìn)行單整檢驗(yàn),容易驗(yàn)證lnC與lnGDP是一階單整的，它們適合的檢驗(yàn)?zāi)Ｐ腿缦拢?(3.81)（-4.01） （2.66） （2.26） （2.54） LM(1)=0.38 LM(2)=0.67 LM(3)=2.34 LM(4)=2.46,首先，

27、建立lnC與lnGDP的回歸模型,（2）檢驗(yàn)lnC與lnGDP的協(xié)整性，并建立長(zhǎng)期均衡關(guān)系,（0.30） (57.48) R2=0.994 DW=0.744,發(fā)現(xiàn)有殘關(guān)項(xiàng)有較強(qiáng)的一階自相關(guān)性?？紤]加入適當(dāng)?shù)臏箜?xiàng)，得lnC與lnGDP的分布滯后模型,(1.63) (6.62） （4.92） （-2.17）

28、 R2=0.994 DW=1.92 LM(1)=0.00 LM(2)=2.31,自相關(guān)性消除，因此可初步認(rèn)為是lnC與lnGDP的長(zhǎng)期穩(wěn)定關(guān)系。,殘差項(xiàng)的穩(wěn)定性檢驗(yàn)：,（-4.32） R2=0.994 DW=2.01 LM(1)=0.04 LM(2)=1.34,t=-4.32<-3.64=ADF0.05 說明lnC與lnGDP是（1，1）階協(xié)

29、整的，下式即為它們長(zhǎng)期穩(wěn)定的均衡關(guān)系:,以穩(wěn)定的時(shí)間序列,（3）建立誤差修正模型,做為誤差修正項(xiàng)，可建立如下,誤差修正模型:,（6.96） (2.96) (-1.91) (-3.15) R2=0.994 DW=2.06 LM(1)=0.70 LM(2)=2.04,由式,可得lnC關(guān)于lnGDP的長(zhǎng)

30、期彈性： (0.698-0.361)/(1-0.622)=0.892；由（**）式可得lnC關(guān)于lnGDP的短期彈性：0.686,(**),用打開誤差修正項(xiàng)括號(hào)的方法直接估計(jì)誤差修正模型，適當(dāng)估計(jì)式為:,（1.63）(6.62） (-2.99) (2.88) R2=0.791 =0.0064

31、 DW=1.93 LM(2)=2.31 LM(3)=2.78,寫成誤差修正模型的形式如下,由上式知，lnC關(guān)于lnGDP的短期彈性為0.698，長(zhǎng)期彈性為0.892。 可見兩種方法的結(jié)果非常接近。,（4）預(yù)測(cè),由式,給出1998年關(guān)于長(zhǎng)期均衡點(diǎn)的偏差：,=ln(18230)-0.152-0.698ln(39008)-0.662ln(17072) +0.361ln(36684)= 0.0125,由式,預(yù)測(cè)1999年

32、的短期波動(dòng)： ?lnC99=0.686(ln(41400)-ln(39008))+0.784(ln(18230)-ln(17072)) -0.484(ln(39008)-ln(36684))-1.163×0.0125= 0.048,于是,按照式,預(yù)測(cè)的結(jié)果為: ?lnC99=0.698（ln(41400)-ln(39008))-0.378(ln(18230)-0.4