第七章序列相關(guān)性1 - 副本

1、計量經(jīng)濟學計量經(jīng)濟學理論理論方法方法EViewsEViews應用應用 郭存芝郭存芝 杜延軍杜延軍 李春吉李春吉 編著編著電子教案第七章第七章 序列相關(guān)性序列相關(guān)性 學習目的學習目的 通過本章的學習，你可以知道什么是序列相關(guān)性，序列通過本章的學習，你可以知道什么是序列相關(guān)性，序列相關(guān)性產(chǎn)生的原因是什么，序列相關(guān)性導致什么樣的后果，相關(guān)性產(chǎn)生的原因是什么，序列相關(guān)性導致什么樣的后果，怎樣檢驗和處理具有序列相關(guān)性的模型。怎樣檢驗和處理具有序列相關(guān)性的模型。 基本要求基本要求1）掌握序列相關(guān)性的概念、序列相關(guān)性的后果和檢驗方法；）掌握序列相關(guān)性的概念、序列相關(guān)性的后果和檢驗方法；2）了解廣義最小二乘法

2、和廣義差分法原理；）了解廣義最小二乘法和廣義差分法原理；3）能運用廣義差分法和廣義最小二乘法估計線性回歸模型。）能運用廣義差分法和廣義最小二乘法估計線性回歸模型。序列相關(guān)性及其產(chǎn)生原因序列相關(guān)性及其產(chǎn)生原因 序列相關(guān)性的影響序列相關(guān)性的影響序列相關(guān)性的檢驗序列相關(guān)性的檢驗序列相關(guān)的補救序列相關(guān)的補救第七章第七章 序列相關(guān)性序列相關(guān)性案例分析案例分析第一節(jié)第一節(jié) 序列相關(guān)性及其產(chǎn)生原因序列相關(guān)性及其產(chǎn)生原因、序列相關(guān)性的含義、序列相關(guān)性的含義對于多元線性回歸模型對于多元線性回歸模型 01122 1,2,iiikkiiYXXXin(7-1)在其他假設(shè)仍然成立的條件下，隨機干擾項序列相關(guān)意味著在其他

3、假設(shè)仍然成立的條件下，隨機干擾項序列相關(guān)意味著如果僅存在如果僅存在則稱為則稱為一階序列相關(guān)或自相關(guān)一階序列相關(guān)或自相關(guān)（簡寫為（簡寫為AR(1)，這是常見的一種序列相關(guān)問題。，這是常見的一種序列相關(guān)問題。1()0,1,2,.,iiEin （7-3）(,)()0CovEijij （7-2）自相關(guān)往往可以寫成如下形式：自相關(guān)往往可以寫成如下形式： （7-4）1 , 11iii 其中其中稱為自協(xié)方差系數(shù)或一階自回歸系數(shù)，稱為自協(xié)方差系數(shù)或一階自回歸系數(shù)，i是滿足以下標準是滿足以下標準OLS假定的隨機干擾項：假定的隨機干擾項：2( )0, ( ), ( ,)0(0)iiii sEVarCovs 由于序

4、列相關(guān)性經(jīng)常出現(xiàn)在以由于序列相關(guān)性經(jīng)常出現(xiàn)在以時間序列時間序列數(shù)據(jù)為樣本的模型中，因此，數(shù)據(jù)為樣本的模型中，因此，本節(jié)下面將代表不同樣本點的下表本節(jié)下面將代表不同樣本點的下表i用用t 表示。表示。二、序列相關(guān)的原因二、序列相關(guān)的原因1 1經(jīng)濟數(shù)據(jù)序列慣性經(jīng)濟數(shù)據(jù)序列慣性2 2模型設(shè)定的偏誤模型設(shè)定的偏誤3 3滯后效應滯后效應4 4蛛網(wǎng)現(xiàn)象蛛網(wǎng)現(xiàn)象5 5數(shù)據(jù)的編造數(shù)據(jù)的編造1 1經(jīng)濟數(shù)據(jù)序列慣性經(jīng)濟數(shù)據(jù)序列慣性 GDP、價格指數(shù)、消費等時間序列數(shù)據(jù)通常表現(xiàn)為周期循環(huán)。當經(jīng)、價格指數(shù)、消費等時間序列數(shù)據(jù)通常表現(xiàn)為周期循環(huán)。當經(jīng)濟衰退的谷底開始復蘇時，大多數(shù)經(jīng)濟序列開始上升，在上升期間，序濟衰退的谷

5、底開始復蘇時，大多數(shù)經(jīng)濟序列開始上升，在上升期間，序列在每一時刻的值都高于前一時刻的值?？磥碛幸环N內(nèi)在的動力驅(qū)使這列在每一時刻的值都高于前一時刻的值?？磥碛幸环N內(nèi)在的動力驅(qū)使這一勢頭繼續(xù)下去，直至某些情況出現(xiàn)（如利率或稅收提高）才把它拖慢一勢頭繼續(xù)下去，直至某些情況出現(xiàn)（如利率或稅收提高）才把它拖慢下來。下來。 因此，在涉及時間序列的回歸中，相繼的觀測值很可能是相互依賴的。因此，在涉及時間序列的回歸中，相繼的觀測值很可能是相互依賴的。比如：比如：2 2模型設(shè)定的偏誤模型設(shè)定的偏誤定義：定義： 指所設(shè)定的模型指所設(shè)定的模型“不正確不正確”，主要表現(xiàn)在模型中丟掉了重要的解釋，主要表現(xiàn)在模型中丟掉了

6、重要的解釋變量或模型函數(shù)形式有偏誤。變量或模型函數(shù)形式有偏誤。 例例1：本來應該估計的模型為本來應該估計的模型為 0111233tttttYXXX（7-5） 但在進行回歸時，卻把模型設(shè)定為如下形式：但在進行回歸時，卻把模型設(shè)定為如下形式：t011t22ttY = + X + X + （7-6） （丟掉了重要的解釋變量）（丟掉了重要的解釋變量）2 2模型設(shè)定的偏誤模型設(shè)定的偏誤 如果如果(7-5)(7-5)式是正確的模型，那做（式是正確的模型，那做（7-67-6）式的回歸就相）式的回歸就相當于令當于令33tttvX于是誤差項于是誤差項v v將表現(xiàn)出一種系統(tǒng)性模式，從而形成了將表現(xiàn)出一種系統(tǒng)性模式

7、，從而形成了自相關(guān)自相關(guān)。例例1：（丟掉了重要的解釋變量）（丟掉了重要的解釋變量）2 2模型設(shè)定的偏誤模型設(shè)定的偏誤定義：定義： 指所設(shè)定的模型指所設(shè)定的模型“不正確不正確”，主要表現(xiàn)在模型中丟掉了重要的解釋，主要表現(xiàn)在模型中丟掉了重要的解釋變量或模型函數(shù)形式有偏誤。變量或模型函數(shù)形式有偏誤。 例例2：（模型函數(shù)形式有偏誤）（模型函數(shù)形式有偏誤）（7-7） 在成本在成本產(chǎn)出研究中，如果真實的邊際成本的模型為：產(chǎn)出研究中，如果真實的邊際成本的模型為： 0122ttttY = + X + X + 其中其中Y代表邊際成本，代表邊際成本，X代表產(chǎn)出。代表產(chǎn)出。 （7-8） 01tttYXv但是如果建模

8、時設(shè)立了如下回歸模型但是如果建模時設(shè)立了如下回歸模型：2 2模型設(shè)定的偏誤模型設(shè)定的偏誤例例2：（模型函數(shù)形式有偏誤）（模型函數(shù)形式有偏誤） 因此在因此在(7-8)(7-8)中，中， 它包含了產(chǎn)出的平方對隨機它包含了產(chǎn)出的平方對隨機 干擾項的系統(tǒng)性影響，隨機干擾項呈現(xiàn)序列相關(guān)性。干擾項的系統(tǒng)性影響，隨機干擾項呈現(xiàn)序列相關(guān)性。tttXv22 （7-8） 01tttYXv 但是如果建模時設(shè)立了如下回歸模型但是如果建模時設(shè)立了如下回歸模型：3 3滯后效應滯后效應 考慮一個消費支出對收入進行回歸的時間序列模型，人們常常發(fā)考慮一個消費支出對收入進行回歸的時間序列模型，人們常常發(fā)現(xiàn)當期的消費支出除了依賴其

9、他當期收入外，還依賴前期的消費支出，現(xiàn)當期的消費支出除了依賴其他當期收入外，還依賴前期的消費支出，即回歸模型為：即回歸模型為：0121ttttCYC （7-9）其中，其中，C是消費，是消費，Y是收入。是收入。 類似（類似（7-9）式的回歸模型被稱為）式的回歸模型被稱為自回歸模型自回歸模型 由于心理上、技術(shù)上以及制度上的原因，消費者不會輕易改變其消費由于心理上、技術(shù)上以及制度上的原因，消費者不會輕易改變其消費習慣，如果我們忽視（習慣，如果我們忽視（7-9）式中的滯后消費對當前消費的影響，那所帶來）式中的滯后消費對當前消費的影響，那所帶來的誤差項就會體現(xiàn)出一種系統(tǒng)性的模式。的誤差項就會體現(xiàn)出一種系

10、統(tǒng)性的模式。注意：注意：4 4蛛網(wǎng)現(xiàn)象蛛網(wǎng)現(xiàn)象例如：例如：假定某農(nóng)產(chǎn)品的供給模型為：假定某農(nóng)產(chǎn)品的供給模型為：01-1tttSP (7-10)假設(shè)t時期的價格Pt低于t-1時期的價格Pt-1，農(nóng)民就很可能決定在時期t+1生產(chǎn)比t時期更少的東西。顯然在這種情形中，農(nóng)民由于在年度t的過量生產(chǎn)很可能在年度t+1消減他們的產(chǎn)量。諸如此類的現(xiàn)象，就不能期望干擾t t是隨機，從而出現(xiàn)蛛網(wǎng)式的序列相關(guān)。5 5數(shù)據(jù)的編造數(shù)據(jù)的編造新生成的數(shù)據(jù)與原數(shù)據(jù)間就有了內(nèi)在的聯(lián)系，表現(xiàn)出序列相關(guān)性。新生成的數(shù)據(jù)與原數(shù)據(jù)間就有了內(nèi)在的聯(lián)系，表現(xiàn)出序列相關(guān)性。 例如：例如： 季度數(shù)據(jù)來自月度數(shù)據(jù)的簡單平均，這種平均的計算減季

11、度數(shù)據(jù)來自月度數(shù)據(jù)的簡單平均，這種平均的計算減弱了每月數(shù)據(jù)的波動而引進了數(shù)據(jù)中的勻滑性，這種勻滑性弱了每月數(shù)據(jù)的波動而引進了數(shù)據(jù)中的勻滑性，這種勻滑性本身就能使隨機干擾項中出現(xiàn)系統(tǒng)性的因素，從而出現(xiàn)序列本身就能使隨機干擾項中出現(xiàn)系統(tǒng)性的因素，從而出現(xiàn)序列相關(guān)性。相關(guān)性。 利用數(shù)據(jù)的內(nèi)插或外推技術(shù)構(gòu)造的數(shù)據(jù)也會呈現(xiàn)某種系統(tǒng)性的模式。利用數(shù)據(jù)的內(nèi)插或外推技術(shù)構(gòu)造的數(shù)據(jù)也會呈現(xiàn)某種系統(tǒng)性的模式。 一般經(jīng)驗表明，對于采用時間序列數(shù)據(jù)做樣本的計量經(jīng)濟學模型，一般經(jīng)驗表明，對于采用時間序列數(shù)據(jù)做樣本的計量經(jīng)濟學模型，由于在不同樣本點上解釋變量意外的其他因素在時間上的連續(xù)性，由于在不同樣本點上解釋變量意外的

12、其他因素在時間上的連續(xù)性，帶來了他們對被解釋變量的影響的連續(xù)性，所以往往存在序列相關(guān)性。帶來了他們對被解釋變量的影響的連續(xù)性，所以往往存在序列相關(guān)性。第二節(jié)第二節(jié) 序列相關(guān)性序列相關(guān)性的影響的影響1 1參數(shù)估計量非有效參數(shù)估計量非有效2 2隨機誤差項方差估計量是有偏的隨機誤差項方差估計量是有偏的3 3擬合優(yōu)度檢驗擬合優(yōu)度檢驗R2統(tǒng)計量和方程顯著性檢驗統(tǒng)計量和方程顯著性檢驗F F統(tǒng)計量無效統(tǒng)計量無效4 4變量的顯著性檢驗變量的顯著性檢驗t t檢驗統(tǒng)計量和相應的參數(shù)置信區(qū)間估計失去意義檢驗統(tǒng)計量和相應的參數(shù)置信區(qū)間估計失去意義5 5模型的預測失效模型的預測失效(,)()0OLSttjttjCovE

13、 如果我們在干擾中通過假定引進序列相關(guān)，但保留經(jīng)典模型的全部其他假定，對估計量及其方差來說會出現(xiàn)什么情況呢？1 1參數(shù)估計量非有效參數(shù)估計量非有效 根據(jù)根據(jù)OLS估計中關(guān)于參數(shù)估計量的無偏性和有效性的證明過程估計中關(guān)于參數(shù)估計量的無偏性和有效性的證明過程可以看出，當計量經(jīng)濟學模型出現(xiàn)序列相關(guān)性時，其可以看出，當計量經(jīng)濟學模型出現(xiàn)序列相關(guān)性時，其OLS參數(shù)估計參數(shù)估計量仍然具有線性無偏性，但不具有有效性。因為在有效性證明中我量仍然具有線性無偏性，但不具有有效性。因為在有效性證明中我們利用了們利用了2()E I（7-11） 即同方差和相互獨立性條件。而且在大樣本情況下，參數(shù)估計量雖然即同方差和相互

14、獨立性條件。而且在大樣本情況下，參數(shù)估計量雖然具有一致性，但仍然不具有漸近有效性。具有一致性，但仍然不具有漸近有效性。為了具體說明這一點，我們回到簡單的一元回歸模型為了具體說明這一點，我們回到簡單的一元回歸模型（7-12）01iiiYX為方便我們不妨假定干擾項為為方便我們不妨假定干擾項為(7-4)所示的一階序列相關(guān)：所示的一階序列相關(guān)：1ttt (7-13)12tttx yx (7-14)對于干擾項為一階序列相關(guān)的一元回歸模型采用對于干擾項為一階序列相關(guān)的一元回歸模型采用OLS估計，如以前估計，如以前一樣，一樣，1 1的的OLS估計量為：估計量為： 但給定干擾項為一階序列相關(guān)時，但給定干擾項為

15、一階序列相關(guān)時，1 的方差估計量現(xiàn)在為：的方差估計量現(xiàn)在為：11()ARVar1式中式中為一階序列相關(guān)時為一階序列相關(guān)時的方差。的方差。212()tVarx （7-16）把該式與沒有干擾項自相關(guān)情形的通常公式把該式與沒有干擾項自相關(guān)情形的通常公式11212212211111222221112()nnntttttnntttARnnnttttttttx xx xx xVarxxxxx （7-15）X相比，可以看出前者等于后者加上另一與自相關(guān)系數(shù)相比，可以看出前者等于后者加上另一與自相關(guān)系數(shù)和各期和各期的樣本協(xié)方差有關(guān)的項。的樣本協(xié)方差有關(guān)的項。 2 2隨機誤差項方差估計量是有偏的隨機誤差項方差估計

16、量是有偏的 在存在干擾項序列相關(guān)的情況下，在存在干擾項序列相關(guān)的情況下，隨機誤差方差的隨機誤差方差的OLS估計量偏離估計量偏離了真實的隨機誤差項的方差了真實的隨機誤差項的方差2。 以一元回歸模型為例，在經(jīng)以一元回歸模型為例，在經(jīng)典假設(shè)情況下，干擾項的典假設(shè)情況下，干擾項的OLS方差估計量方差估計量2212ntten222()E是真實的是真實的的無偏估計，即有的無偏估計，即有。但若隨機誤差項存在一階序列相關(guān)但若隨機誤差項存在一階序列相關(guān) 則可以證明：則可以證明： 222/(1)2()2nrEn式中式中1221nt ttnttx xxr為為X的相繼觀測值之間的樣本相關(guān)系數(shù)。的相繼觀測值之間的樣本相

17、關(guān)系數(shù)。 3 3擬合優(yōu)度檢驗擬合優(yōu)度檢驗R2統(tǒng)計量和方程顯著性檢驗統(tǒng)計量和方程顯著性檢驗F統(tǒng)計量無效統(tǒng)計量無效 由于在序列相關(guān)時由于在序列相關(guān)時OLS對隨機誤差方差估計有偏，結(jié)果基于對隨機誤差方差估計有偏，結(jié)果基于OLS殘差平方和計算出來的擬合優(yōu)度檢驗統(tǒng)計量殘差平方和計算出來的擬合優(yōu)度檢驗統(tǒng)計量R2也失去意義，也失去意義，相應的方程顯著性檢驗統(tǒng)計量相應的方程顯著性檢驗統(tǒng)計量F統(tǒng)計量也無效。統(tǒng)計量也無效。4 4變量的顯著性檢驗變量的顯著性檢驗t t 檢驗統(tǒng)計量和相應的參數(shù)置檢驗統(tǒng)計量和相應的參數(shù)置 信區(qū)間估計失去意義信區(qū)間估計失去意義 用用OLS法估計序列相關(guān)的模型得到的隨機誤差項的方差不僅是法

18、估計序列相關(guān)的模型得到的隨機誤差項的方差不僅是有偏的，而且這一偏誤也將傳遞到用有偏的，而且這一偏誤也將傳遞到用OLS方法得到的參數(shù)估計方法得到的參數(shù)估計量的方差中來，從而使得建立在量的方差中來，從而使得建立在OLS參數(shù)估計量方差基礎(chǔ)上的參數(shù)估計量方差基礎(chǔ)上的變量顯著性檢驗失去意義。變量顯著性檢驗失去意義。212()tVarx沒有被低估，通常沒有被低估，通常OLS參數(shù)估計量的方差式（參數(shù)估計量的方差式（7-16）2即使隨機誤差的方差即使隨機誤差的方差也是存在一階序列相關(guān)時參數(shù)估計量方差的偏誤估計量。也是存在一階序列相關(guān)時參數(shù)估計量方差的偏誤估計量。 以一元回歸模型為例，以一元回歸模型為例，01i

19、iiYX5 5模型的預測失效模型的預測失效 在存在序列相關(guān)時在存在序列相關(guān)時OLS估計的隨機誤差項方差有偏，參數(shù)估計量方估計的隨機誤差項方差有偏，參數(shù)估計量方差非有效，這樣回歸模型的被解釋變量的預測值及預測區(qū)間就不準確，差非有效，這樣回歸模型的被解釋變量的預測值及預測區(qū)間就不準確，預測精度降低。預測精度降低。 被解釋變量預測值區(qū)間與模型參數(shù)和隨機誤差的估計量的方差有關(guān)。被解釋變量預測值區(qū)間與模型參數(shù)和隨機誤差的估計量的方差有關(guān)。所以，當模型出現(xiàn)序列相關(guān)時，它的預測功能失效。所以，當模型出現(xiàn)序列相關(guān)時，它的預測功能失效。tolsttYYe第三節(jié)第三節(jié) 序列相關(guān)性的序列相關(guān)性的檢驗檢驗不同的檢驗方

20、法的共同思路不同的檢驗方法的共同思路: 序列相關(guān)性的檢驗方法有多種，如馮諾曼比檢驗法、回歸檢驗法、序列相關(guān)性的檢驗方法有多種，如馮諾曼比檢驗法、回歸檢驗法、D.W.檢驗法等。檢驗法等。 te 首先采用普通最小二乘法估計模型，首先采用普通最小二乘法估計模型，以得到隨機干擾項的近似估以得到隨機干擾項的近似估計量，我們用計量，我們用表示表示近似估計量：近似估計量： （7-19） 然后通過分析這些近似估計量之間的相關(guān)性以達到判斷隨機干擾然后通過分析這些近似估計量之間的相關(guān)性以達到判斷隨機干擾項是否具有序列相關(guān)性的目的。項是否具有序列相關(guān)性的目的。序列相關(guān)性的檢驗方法序列相關(guān)性的檢驗方法一、圖示法一、圖

21、示法二、回歸檢驗法二、回歸檢驗法三、杜賓三、杜賓沃森檢驗沃森檢驗四、拉格朗日乘子檢驗四、拉格朗日乘子檢驗一、圖示法一、圖示法tttete 由于殘差由于殘差 可以作為隨機誤差可以作為隨機誤差 的估計，因此，如果的估計，因此，如果 存在序列相關(guān)性，存在序列相關(guān)性， 反映出來，因此可以利用反映出來，因此可以利用 的變化來判斷隨機干擾項的序列的變化來判斷隨機干擾項的序列必然會由殘差項必然會由殘差項相關(guān)性，如圖相關(guān)性，如圖71所示。所示。te二二、回歸檢驗法回歸檢驗法，（7-20）（7-21） 以以te21 ,ttee等為解釋變量，等為解釋變量， 為解釋變量，以各種可能的相關(guān)變量，諸如為解釋變量，以各種

22、可能的相關(guān)變量，諸如建立各種方程：建立各種方程： 對方程進行估計并進行顯著性檢驗，如果存在某一種函數(shù)形式，使得方對方程進行估計并進行顯著性檢驗，如果存在某一種函數(shù)形式，使得方程顯著成立，則說明原模型存在序列相關(guān)性。程顯著成立，則說明原模型存在序列相關(guān)性。 一旦確定了模型存在序列相關(guān)性，也就同時知道了相關(guān)的形式，而且它一旦確定了模型存在序列相關(guān)性，也就同時知道了相關(guān)的形式，而且它適用于任何類型的序列相關(guān)性問題的檢驗。適用于任何類型的序列相關(guān)性問題的檢驗。優(yōu)點：優(yōu)點：),.,3( ,.,2 21121nteeen) (teettttttt三、三、杜賓杜賓沃森檢驗沃森檢驗 D-W檢驗是杜賓（檢驗是杜

23、賓（J.Durbin）和沃森（）和沃森（G.S.Watson）于）于1951年提出年提出的一種檢驗序列自相關(guān)的方法。雖然該方法很常用，但它有一些的一種檢驗序列自相關(guān)的方法。雖然該方法很常用，但它有一些基本假定基本假定：（1）回歸含有截距項。）回歸含有截距項。（2）解釋變量）解釋變量X是非隨機的，或者在重復抽樣中被固定的。是非隨機的，或者在重復抽樣中被固定的。t1ttt（3）隨機干擾項）隨機干擾項 為一階自回歸形式：為一階自回歸形式：。（4）回歸模型中不應把滯后應變量作為解釋變量之一，即不應出現(xiàn)如下形式模型：）回歸模型中不應把滯后應變量作為解釋變量之一，即不應出現(xiàn)如下形式模型：011221 tt

24、tkktttYXXXY（5）沒有缺失數(shù)據(jù)。）沒有缺失數(shù)據(jù)。 tntntttteeeDW12212)(0:0H杜賓杜賓沃森針對原假設(shè)沃森針對原假設(shè) ，即，即 不存在一階自相關(guān)，構(gòu)造如下統(tǒng)計量：不存在一階自相關(guān)，構(gòu)造如下統(tǒng)計量： （7-22）2UdLd檢驗，它沒有唯一的臨界值可以導出拒絕或檢驗，它沒有唯一的臨界值可以導出拒絕或和下限和下限，且這些上下限只與，且這些上下限只與因此因此D-W檢驗不同于檢驗不同于t、F或或接受原假設(shè)。但他們成功導出了臨界值的上限接受原假設(shè)。但他們成功導出了臨界值的上限樣本容量樣本容量n和解釋變量的個數(shù)有關(guān)，而與解釋變量的取值無關(guān)。和解釋變量的個數(shù)有關(guān)，而與解釋變量的取值

25、無關(guān)。 杜賓杜賓沃森證明該統(tǒng)計量的分布與出現(xiàn)在給定樣本中的沃森證明該統(tǒng)計量的分布與出現(xiàn)在給定樣本中的X值有復雜的關(guān)系，值有復雜的關(guān)系，其準確的抽樣或概率分布很難得到；其準確的抽樣或概率分布很難得到； 又依賴于給定的又依賴于給定的X的值。的值。因為因為D.W.值要從值要從中算出，而中算出，而teteLdUd 因此，在運用因此，在運用D-W檢驗時，只須計算該統(tǒng)計量的值，再根據(jù)樣本容量檢驗時，只須計算該統(tǒng)計量的值，再根據(jù)樣本容量n 和和 ，然后按下列準則，然后按下列準則考察考察 和解釋變量數(shù)目和解釋變量數(shù)目k查查D.W.分布表，得到臨界值分布表，得到臨界值 計算得到的計算得到的D.W.值，以判斷模型

26、的自相關(guān)狀態(tài)：值，以判斷模型的自相關(guān)狀態(tài)： 若若0.LDWd ，則存在正自相關(guān)；，則存在正自相關(guān)； 若若.LUdDWd ，則不確定；，則不確定； 若若.4UUdDWd ，則無自相關(guān)；，則無自相關(guān)； 若若4.4ULdDWd ，則不確定；，則不確定； 若若4.4LdDW ，則存在負自相關(guān)。，則存在負自相關(guān)。也就是說，當也就是說，當D.W.值在值在2附近時，模型不存在一階自相關(guān)。附近時，模型不存在一階自相關(guān)。 例例7-1 7-1 給定一個含有給定一個含有50個觀測值的樣本和個觀測值的樣本和3個解釋變量個解釋變量,如果如果（a）D.W.=1.05，（，（b）D.W.=1.40，（c）D.W.=2.50

27、，（，（d）D.W.=3.97你能對自相關(guān)的問題說些什么？你能對自相關(guān)的問題說些什么？解：解： 根據(jù)根據(jù)D-W檢驗判斷準則可知檢驗判斷準則可知 （b）D.W.=1.40Ld，隨機誤差項存在一階正自相關(guān)；，隨機誤差項存在一階正自相關(guān)； （d）4 Ld =2.58 D.W.=3.97，隨機誤差項存在負一階自相關(guān)。，隨機誤差項存在負一階自相關(guān)。LdUd查查D.W.分布表可知，當樣本數(shù)為分布表可知，當樣本數(shù)為n=50，解釋變量數(shù)，解釋變量數(shù)k=3時，在時，在5%的的為為1.42，為為1.67。 顯著性水平下顯著性水平下D.W.統(tǒng)計量臨界值的下界統(tǒng)計量臨界值的下界（a）D.W.=1.05 D.W.=2.

28、504 =2.33，不能確定隨機誤差項是否，不能確定隨機誤差項是否存在一階自相關(guān)；存在一階自相關(guān)； 在許多情況下，人們發(fā)現(xiàn)上限在許多情況下，人們發(fā)現(xiàn)上限Ud 差不多就是真實的顯著性界限，因而，差不多就是真實的顯著性界限，因而， 如果如果D.W.的估計值落入不能確定的區(qū)域，人們可以使用以下修正的的估計值落入不能確定的區(qū)域，人們可以使用以下修正的D-W 檢驗程序。給定顯著性水平檢驗程序。給定顯著性水平：0:0H1:0H （2）原假設(shè)為）原假設(shè)為 ，備擇假設(shè)為，備擇假設(shè)為0:0H1:0H （1）原假設(shè)為）原假設(shè)為 ，備擇假設(shè)為，備擇假設(shè)為.UDWd 如果有如果有 ，則在顯著性水平，則在顯著性水平上拒

29、絕原假設(shè)上拒絕原假設(shè)H0，接受，接受備擇假設(shè)備擇假設(shè)H1，也就是存在統(tǒng)計上顯著的正相關(guān)。，也就是存在統(tǒng)計上顯著的正相關(guān)。4.UDWd 如果有如果有 ，則在顯著性水平，則在顯著性水平上拒絕原假設(shè)上拒絕原假設(shè)H0，接受，接受備擇假設(shè)備擇假設(shè)H1，也就是存在統(tǒng)計上顯著的負相關(guān)。，也就是存在統(tǒng)計上顯著的負相關(guān)。 在許多情況下，人們發(fā)現(xiàn)上限在許多情況下，人們發(fā)現(xiàn)上限Ud 差不多就是真實的顯著性界限，因而，差不多就是真實的顯著性界限，因而， 如果如果D.W.的估計值落入不能確定的區(qū)域，人們可以使用以下修正的的估計值落入不能確定的區(qū)域，人們可以使用以下修正的D-W 檢驗程序。給定顯著性水平檢驗程序。給定顯著

30、性水平：0:0H1:0H （3）原假設(shè)為）原假設(shè)為 ，備擇假設(shè)為，備擇假設(shè)為.UDWd4.UDWd 如果有如果有 或者或者則在顯著性水平則在顯著性水平上拒絕原假設(shè)上拒絕原假設(shè)H0，接受備擇假設(shè)，接受備擇假設(shè)H1，也就是存在統(tǒng)計上顯著的自相關(guān)。也就是存在統(tǒng)計上顯著的自相關(guān)。四、四、拉格朗日乘子檢驗拉格朗日乘子檢驗 拉格朗日乘子檢驗克服了拉格朗日乘子檢驗克服了D-W檢驗的缺陷，適合于高階序列相關(guān)檢驗的缺陷，適合于高階序列相關(guān)及模型中存在滯后被解釋變量的情形。它是由布勞殊（及模型中存在滯后被解釋變量的情形。它是由布勞殊（Breusch）與）與戈弗雷（戈弗雷（Godfrey）于）于1978年提出的，也

31、稱為年提出的，也稱為GB檢驗檢驗。對于模型對于模型01122tttkkttYXXX（7-24）如果要檢驗隨機誤差項是否存在如果要檢驗隨機誤差項是否存在p階序列相關(guān)：階序列相關(guān)： （7-25）1212.ttttptp 那么檢驗如下受約束回歸方程就是拉格朗日乘子檢驗：那么檢驗如下受約束回歸方程就是拉格朗日乘子檢驗：011221212.tttkkttttptpYXXX （7-26）約束條件為約束條件為 （7-27）012:0pH 如果約束條件為真，則如果約束條件為真，則LM統(tǒng)計量服從大樣本下自由度為統(tǒng)計量服從大樣本下自由度為p的漸近的漸近2分布：分布：22()pLMn pR（7-28）其中其中n p

32、和和2R分別為如下輔助回歸方程的樣本容量和可決系數(shù)：分別為如下輔助回歸方程的樣本容量和可決系數(shù)：tptptktktteeXXe.11110（7-29） (7-29)中的被解釋變量中的被解釋變量是對原模型（是對原模型（7-24）進行）進行OLS回歸后得到的殘差。回歸后得到的殘差。 te p值即滯后的長度無法預先給定，因此實踐操作中可從值即滯后的長度無法預先給定，因此實踐操作中可從1階、階、2階階逐次相更高階檢驗，并用輔助回歸方程（逐次相更高階檢驗，并用輔助回歸方程（7-29）式中各個殘差項前面的）式中各個殘差項前面的參數(shù)的顯著性來幫助判斷序列相關(guān)的階數(shù)。參數(shù)的顯著性來幫助判斷序列相關(guān)的階數(shù)。LM

33、檢驗的一個缺陷檢驗的一個缺陷給定顯著性水平給定顯著性水平，查自由度為，查自由度為p p的分布的相應臨界值，如果計的分布的相應臨界值，如果計算的算的LMLM統(tǒng)計量的值超過該臨界值，則拒絕約束條件為真的原假統(tǒng)計量的值超過該臨界值，則拒絕約束條件為真的原假設(shè)，表明存在隨機誤差項存在直到設(shè)，表明存在隨機誤差項存在直到p p階的序列相關(guān)性。階的序列相關(guān)性。例例7-2 7-2 假定用假定用32個樣本做個樣本做Y對對X(包含截距包含截距)的回歸的回歸而這樣的而這樣的2數(shù)值對應的概率數(shù)值對應的概率p為為0.0003，這是一個很低的概率。，這是一個很低的概率。 因此我們可以拒絕輔助回歸方程中原始回歸殘差序列的全

34、部因此我們可以拒絕輔助回歸方程中原始回歸殘差序列的全部1到到5階滯后階滯后序列系數(shù)均為零的假設(shè)，至少有一個滯后殘差序列的系數(shù)不為零。序列系數(shù)均為零的假設(shè)，至少有一個滯后殘差序列的系數(shù)不為零。 這表明原始回歸的殘差中至少存在這表明原始回歸的殘差中至少存在1到到5階中的某一滯后的自相關(guān)，當然階中的某一滯后的自相關(guān)，當然要確定到底是幾階序列相關(guān)還必須進一步進行要確定到底是幾階序列相關(guān)還必須進一步進行4階、階、3階階等不同階數(shù)的拉格等不同階數(shù)的拉格朗日乘子檢驗。朗日乘子檢驗。如果我們懷疑回歸殘差序列有如果我們懷疑回歸殘差序列有5階滯后相關(guān)，那么輔助回歸方程中我們階滯后相關(guān)，那么輔助回歸方程中我們可以用

35、殘差對可以用殘差對X以及殘差序列的以及殘差序列的1到到5階滯后序列進行回歸，假定從輔助回歸方階滯后序列進行回歸，假定從輔助回歸方程中回歸得到的擬合優(yōu)度程中回歸得到的擬合優(yōu)度R2為為0.8860。 由于原始回歸中有由于原始回歸中有32個樣本，而輔助回歸中用了個樣本，而輔助回歸中用了5個滯后值，這樣輔助個滯后值，這樣輔助2()np R等于等于(32-5)0.886即等于即等于23.382。 回歸方程中僅有回歸方程中僅有27個樣本，因此個樣本，因此第四節(jié)第四節(jié) 序列相關(guān)的補救序列相關(guān)的補救 由于序列相關(guān)出現(xiàn)時由于序列相關(guān)出現(xiàn)時OLS估計量是非有效的，因此如果回歸模型被證明估計量是非有效的，因此如果回

36、歸模型被證明存在序列相關(guān)性，則應該發(fā)展新的方法來估計模型。類似于處理異方差的情存在序列相關(guān)性，則應該發(fā)展新的方法來估計模型。類似于處理異方差的情況，在大樣本下我們也可以用與自相關(guān)相一致的況，在大樣本下我們也可以用與自相關(guān)相一致的OLS回歸殘差的方差協(xié)方差回歸殘差的方差協(xié)方差矩陣來處理隨機誤差項的自相關(guān)情況，這樣矩陣來處理隨機誤差項的自相關(guān)情況，這樣OLS估計也仍然是有效的，只是估計也仍然是有效的，只是我們需要報告相應的自相關(guān)穩(wěn)健標準差和相應的統(tǒng)計量，其處理方法完全類我們需要報告相應的自相關(guān)穩(wěn)健標準差和相應的統(tǒng)計量，其處理方法完全類似于異方差穩(wěn)健推斷，這里我們不再對自相關(guān)穩(wěn)健推斷詳細論述，我們詳

37、細介似于異方差穩(wěn)健推斷，這里我們不再對自相關(guān)穩(wěn)健推斷詳細論述，我們詳細介紹一般情況下處理序列相關(guān)最常用的紹一般情況下處理序列相關(guān)最常用的廣義最小二乘法廣義最小二乘法（GLS）和）和廣義差分法廣義差分法。一、廣義最小二乘法一、廣義最小二乘法定義定義: : 最具有普遍意義的最小二乘法最具有普遍意義的最小二乘法.普通最小二乘法普通最小二乘法和和加權(quán)最小二乘法加權(quán)最小二乘法是它的特例。是它的特例。 一般情況下，對于模型一般情況下，對于模型（7-30）Y = X + 如果存在序列相關(guān)性，同時存在異方差，即有如果存在序列相關(guān)性，同時存在異方差，即有21121222122212( )()nnnnnCovE

38、顯然，顯然， 是一對稱矩陣，因此存在一可逆矩陣，使得是一對稱矩陣，因此存在一可逆矩陣，使得DD 用用1D左乘（左乘（7-30）式兩邊，得到一個新的模型）式兩邊，得到一個新的模型111D YD XD（7-31）即即 YX該模型具有同方差性和隨機干擾項相互獨立性。因為該模型具有同方差性和隨機干擾項相互獨立性。因為1111)()()(DEDDDEEIDDDDDD2121121則則111111111 ()() ()X XX YX DD XX DD YXXXY這就是原模型（這就是原模型（7-30）式的廣義最小二乘估計量，它是無偏有效的估計量。）式的廣義最小二乘估計量，它是無偏有效的估計量。于是，可以用普

39、通最小二乘法估計模型（于是，可以用普通最小二乘法估計模型（7-31）式，記參數(shù)估計量為）式，記參數(shù)估計量為 ，2由上面的推導過程可知，只要知道隨機干擾項的方差由上面的推導過程可知，只要知道隨機干擾項的方差-協(xié)方差矩陣協(xié)方差矩陣，就可以采用廣義最小二乘法得到參數(shù)的最佳線性無偏估計量。就可以采用廣義最小二乘法得到參數(shù)的最佳線性無偏估計量。j(1)/2n n 然而若只有然而若只有n個樣本點，要對包括各個個樣本點，要對包括各個 在內(nèi)的在內(nèi)的進行估計是困難的，在實踐操作中，往往通過廣義差分法來實現(xiàn)廣義最小二乘估計。進行估計是困難的，在實踐操作中，往往通過廣義差分法來實現(xiàn)廣義最小二乘估計。+k+1個未知參

40、數(shù)個未知參數(shù)二、廣義差分法二、廣義差分法 廣義差分法需要對隨機干擾項自相關(guān)系數(shù)事先給出必要的假設(shè)，廣義差分法需要對隨機干擾項自相關(guān)系數(shù)事先給出必要的假設(shè)，可區(qū)分為兩種情形：自相關(guān)系數(shù)已知和未知。可區(qū)分為兩種情形：自相關(guān)系數(shù)已知和未知。1）自相關(guān)系數(shù)已知時）自相關(guān)系數(shù)已知時t 由于干擾項由于干擾項 是不可觀測的，關(guān)于序列相關(guān)的性質(zhì)往往是一種猜測是不可觀測的，關(guān)于序列相關(guān)的性質(zhì)往往是一種猜測 遵循形如遵循形如(7-4)式那樣的一階自回歸方式，式那樣的一階自回歸方式，t 或?qū)嶋H體驗。實踐中，常假定或?qū)嶋H體驗。實踐中，常假定1ttt （7-32）即：即：（7-32）式中自回歸系數(shù)和隨機干擾項滿足）式中

41、自回歸系數(shù)和隨機干擾項滿足(7-4)的假定。若假定的假定。若假定(7-32)是是為已知時，序列相關(guān)問題就可以圓滿解決。為已知時，序列相關(guān)問題就可以圓滿解決。 真實的，當自相關(guān)系數(shù)真實的，當自相關(guān)系數(shù)為說明這一點，考慮以下多元回歸模型為例：為說明這一點，考慮以下多元回歸模型為例：101.ttktktYXX（7-33）如果（如果（7-33）在時刻）在時刻t成立，則在時刻成立，則在時刻t-1也成立，因此有：也成立，因此有：1101111.tktkttYXX （7-34）用用乘（乘（7-34）兩邊，得到：）兩邊，得到：1011111.ttkkttYXX （7-35）011ttkkttYXX （7-37

42、）其中，其中，0011(1), , , 1,2,tttitititYYYXXXik, 由于由于t 滿足全部滿足全部OLS假定，故可以直接對方程（假定，故可以直接對方程（7-37）進行）進行OLS回歸回歸得到具有得到具有BLUE性質(zhì)的估計量。性質(zhì)的估計量。將（將（7-36）式簡寫為）式簡寫為用（用（7-33）減去（）減去（7-35）得到）得到1111011111011(1). (1)(). ()ttttktktkkttttktkttkYYXXXXXXXX （7-36）更一般地如果多元回歸模型更一般地如果多元回歸模型01122tttkkttYXXX （7-38）中的隨機干擾項存在中的隨機干擾項存在

43、p階序列相關(guān)：階序列相關(guān)：1122.ttttptp （7-39）那么可以將原模型（那么可以將原模型（7-38）式變換為）式變換為110111111111.(1.) (.). (.) ttptppttptpkktktpktptYYYXXXXXX 1, 2 , . , tppn （7-40）（7-40）式即為多元回歸形式的廣義差分模型，該模型不存在序列相關(guān)性。）式即為多元回歸形式的廣義差分模型，該模型不存在序列相關(guān)性。 采用采用OLS法估計該模型得到的參數(shù)估計量即為原模型參數(shù)的無偏有效法估計該模型得到的參數(shù)估計量即為原模型參數(shù)的無偏有效估計量，這樣處理序列相關(guān)的方法就是估計量，這樣處理序列相關(guān)的方

44、法就是廣義差分法廣義差分法。 廣義差分法就是前面我們討論過的廣義最小二乘法（廣義差分法就是前面我們討論過的廣義最小二乘法（GLS），但應注），但應注意滯后的觀測值被排除了。意滯后的觀測值被排除了。 為看清這一點，我們?nèi)匀豢紤]前面的一階序列相關(guān)的情況為看清這一點，我們?nèi)匀豢紤]前面的一階序列相關(guān)的情況我們用矩陣形式把上述估計過程重寫一遍。對于一階序列相關(guān)的隨機誤差項我們用矩陣形式把上述估計過程重寫一遍。對于一階序列相關(guān)的隨機誤差項 1ttt 我們可以證明該隨機干我們可以證明該隨機干擾項的方差和協(xié)方差分別為擾項的方差和協(xié)方差分別為2222221()11(,)1tsstt sVarCov 用矩陣表示為

45、用矩陣表示為1222212( )()1.1.1.1nnnnCovE 根據(jù)線性代數(shù)易知根據(jù)線性代數(shù)易知2212210.0001.00001.000.000. 10000.1000.01 從而有從而有2-1100.00010.00001.000.000.100000.10000.01D用用1D左乘矩陣形式的多元回歸模型左乘矩陣形式的多元回歸模型 Y = X + ,得到得到 111D YD XD （7-41）1111011111011(1). (1)(). ()ttttktktktttttktkttkYYXXXXXXXX 然后展開（然后展開（7-41）式中所有矩陣乘積，去掉展開式的第一行就得到（）式

46、中所有矩陣乘積，去掉展開式的第一行就得到（7-36）一樣的結(jié)果。一樣的結(jié)果。111D YD XD （7-41） 類似地對具有類似地對具有p階序列相關(guān)的多元回歸模型的廣義差分法估計也等同于廣義階序列相關(guān)的多元回歸模型的廣義差分法估計也等同于廣義最小二乘估計，但我們損失了前面最小二乘估計，但我們損失了前面p個樣本觀測值，這一點可以從廣義差分模個樣本觀測值，這一點可以從廣義差分模型（型（7-40）式看出來。在樣本規(guī)模較大而誤差序列相關(guān)階數(shù)較小時，廣義差分）式看出來。在樣本規(guī)模較大而誤差序列相關(guān)階數(shù)較小時，廣義差分法與廣義最小二乘法的估計結(jié)果很接近。但在小樣本或誤差呈現(xiàn)較大的高階序法與廣義最小二乘法的

47、估計結(jié)果很接近。但在小樣本或誤差呈現(xiàn)較大的高階序列相關(guān)時，觀測值的損失可能會對估計結(jié)果有影響。因此在廣義差分變換中，列相關(guān)時，觀測值的損失可能會對估計結(jié)果有影響。因此在廣義差分變換中，有時需彌補這一損失。有時需彌補這一損失。2211,1,(1,2,., )tjtjtYY XXjk這樣廣義差分法的估計結(jié)果就完全等同于廣義最小二乘估計量。這樣廣義差分法的估計結(jié)果就完全等同于廣義最小二乘估計量。例如，在一階序列相關(guān)情況下，對損失的第一次觀測值可進行如下的例如，在一階序列相關(guān)情況下，對損失的第一次觀測值可進行如下的 普萊斯普萊斯-溫斯特（溫斯特（Prais-Winsten）變換）變換：2 2）自相關(guān)系

48、數(shù)未知時的處理）自相關(guān)系數(shù)未知時的處理 盡管廣義差分回歸直接明了，但通常情況下我們并不知道總體模型中隨盡管廣義差分回歸直接明了，但通常情況下我們并不知道總體模型中隨機干擾項的真實自回歸系數(shù)機干擾項的真實自回歸系數(shù)是多少，故廣義差分法一般難以實現(xiàn)。是多少，故廣義差分法一般難以實現(xiàn)。（1 1）一次差分法）一次差分法（2 2）根據(jù)）根據(jù)D.W.D.W.統(tǒng)計量來估計統(tǒng)計量來估計（3 3）科克倫）科克倫- -奧科特（奧科特（Cochrane-OrcuttCochrane-Orcutt）迭代法）迭代法（4 4）杜賓兩步法）杜賓兩步法因此我們需要另想辦法來處理序列相關(guān)問題，我們介紹幾種常用的方法。因此我們需

49、要另想辦法來處理序列相關(guān)問題，我們介紹幾種常用的方法。（1 1）一次差分法）一次差分法 因為自回歸系數(shù)因為自回歸系數(shù)介于（介于（-1，1）之間，我們考慮極端的序列相關(guān)情況，）之間，我們考慮極端的序列相關(guān)情況，即完全的正相關(guān)或負相關(guān)，此時即完全的正相關(guān)或負相關(guān)，此時等于等于1或或 1。 考慮簡單的一元回歸模型：考慮簡單的一元回歸模型： （7-42）01tttYX假定該模型中隨機干擾項為完全一階正相關(guān)，即有假定該模型中隨機干擾項為完全一階正相關(guān)，即有1ttt（7-43）對（對（7-42）進行一次差分得到）進行一次差分得到111111()()()tttttttttYYXXXX即即 （7-44）1tt

50、tYX10（7-44）的差分回歸方程沒有截距，隨機干擾項沒有序列自相關(guān)，因此可以）的差分回歸方程沒有截距，隨機干擾項沒有序列自相關(guān)，因此可以 對它采取過原點對它采取過原點OLS回歸得到回歸得到的的BLUE估計量，注意此時原模型中的截距估計量，注意此時原模型中的截距就不能估計出來了。就不能估計出來了。 如果原模型為包含時間趨勢的模型：如果原模型為包含時間趨勢的模型： 0112tttYXt（7-45） 那么對它進行一次差分后得到那么對它進行一次差分后得到（7-46）12tttYX 該差分模型中含有一截距，因此含有截距的一次差分模型意味著在原模型該差分模型中含有一截距，因此含有截距的一次差分模型意味

51、著在原模型中存在一線性時間趨勢項，而且一次差分模型中的截距就是原模型中時間趨勢中存在一線性時間趨勢項，而且一次差分模型中的截距就是原模型中時間趨勢項的系數(shù)。如果項的系數(shù)。如果2是正的話，這表明原模型中是正的話，這表明原模型中Y除了受除了受X的影響外還有一上升的影響外還有一上升的趨勢。的趨勢。 如果原模型中隨機干擾項是完全一階負相關(guān)的，那么一次差分處理的方法就是相反了。思考思考: : 析析: : 要注意它是以假定要注意它是以假定=1為前提的，如果隨機干擾項不是完全一階為前提的，如果隨機干擾項不是完全一階正相關(guān)，就不能進行這樣的一次差分變換。正相關(guān)，就不能進行這樣的一次差分變換。 怎樣知道假定怎樣

52、知道假定=1是否合理呢？是否合理呢？用貝倫布魯特用貝倫布魯特-韋布（韋布（Belenblutt-Webbtest）統(tǒng)計量來檢驗。）統(tǒng)計量來檢驗。對對進行進行為檢驗假設(shè)為檢驗假設(shè)=1，貝倫布魯特，貝倫布魯特韋布推出如下韋布推出如下g檢驗統(tǒng)計量：檢驗統(tǒng)計量：ntnttteg1222tYjtX各個解釋變量各個解釋變量X的一階差分的一階差分OLS回歸得到的殘差（注意無截距項）。回歸得到的殘差（注意無截距項）。 其中其中是原始模型的是原始模型的OLS殘差，而殘差，而t是被解釋變量是被解釋變量Y的一階差分的一階差分te例例7-3 7-3 假定用假定用32個樣本做個樣本做Y對對X的的OLS回歸得到的殘差平方

53、和回歸得到的殘差平方和RSS1=204.2934，再做再做Y對對X的的OLS回歸（注意在此回歸中沒有截距）得到殘差平方和回歸（注意在此回歸中沒有截距）得到殘差平方和RSS2=28.1938。 g=28.1938/204.2934=0.1377 查查D.W.分布表發(fā)現(xiàn)分布表發(fā)現(xiàn)5%的顯著性水平下的顯著性水平下31個樣本和個樣本和1個解釋變量的個解釋變量的D.W.值下界為值下界為1.363，上界為，上界為1.496。 因此因此 這樣計算的這樣計算的g的數(shù)值小于的數(shù)值小于D.W.統(tǒng)計量的下界，我們不能拒絕統(tǒng)計量的下界，我們不能拒絕 基于這一結(jié)果，對原模型進行一次差分后再用基于這一結(jié)果，對原模型進行一

54、次差分后再用OLS估計是合理的。估計是合理的。=1的原假設(shè)。的原假設(shè)。（2 2）根據(jù)）根據(jù)D.W.D.W.統(tǒng)計量來估計統(tǒng)計量來估計根據(jù)該式我們可以得到根據(jù)該式我們可以得到的計算表達式：的計算表達式：12d （7-48）這是從所估計的這是從所估計的D.W.統(tǒng)計量獲得統(tǒng)計量獲得的一個估計值的簡易方法。的一個估計值的簡易方法。 回想我們前面的回想我們前面的D.W.統(tǒng)計量統(tǒng)計量1212121212ntntttteeedDW12d （7-48）由（由（7-48）可見，僅當）可見，僅當 d 等于或接近于等于或接近于0時，一次差分法中假定時，一次差分法中假定1才是對的才是對的01 此外當此外當d=2時時，d

55、=4時，時，因此因此D.W.統(tǒng)計量為我們提供了一個估計統(tǒng)計量為我們提供了一個估計的現(xiàn)成方法。的現(xiàn)成方法。 但要注意的是，（但要注意的是，（7-48）僅提供了一個估計）僅提供了一個估計 的近似式，在小樣本下的近似式，在小樣本下未必可靠，僅在大樣本下才具有最優(yōu)性質(zhì)。未必可靠，僅在大樣本下才具有最優(yōu)性質(zhì)。一旦從（一旦從（7-48）估計出）估計出，我們就可以對原模型進行廣義差分變換，我們就可以對原模型進行廣義差分變換，然后對廣義差分后的模型進行然后對廣義差分后的模型進行OLS估計。估計。 同樣需要注意的是，由于廣義差分法中用的是真實的同樣需要注意的是，由于廣義差分法中用的是真實的，而我們是用，而我們是

56、用來代替真實的來代替真實的，因此就會出現(xiàn)一個問題：，因此就會出現(xiàn)一個問題： 估計的估計的這樣估計的回歸系數(shù)是否有經(jīng)典回歸模型中所說的最優(yōu)性質(zhì)呢？這樣估計的回歸系數(shù)是否有經(jīng)典回歸模型中所說的最優(yōu)性質(zhì)呢？ 當用一個估計的量去代替真值時，當用一個估計的量去代替真值時，OLS估計得到的回歸系數(shù)僅是漸近有估計得到的回歸系數(shù)僅是漸近有效的，就是說僅在大樣本情況下才是最優(yōu)的，而且通常的假設(shè)檢驗統(tǒng)計量也效的，就是說僅在大樣本情況下才是最優(yōu)的，而且通常的假設(shè)檢驗統(tǒng)計量也僅是漸近有效的僅是漸近有效的。 一個一般性的原則：一個一般性的原則：（3 3）科克倫）科克倫- -奧科特（奧科特（Cochrane-Orcutt

57、Cochrane-Orcutt）迭代法）迭代法利用估計的殘差去獲得關(guān)于未知的利用估計的殘差去獲得關(guān)于未知的的信息。的信息。 考慮一元回歸模型：考慮一元回歸模型： 01tttYX （7-49）假定隨機干擾項為一階自相關(guān)，即假定隨機干擾項為一階自相關(guān)，即（7-50）1ttt 按如下步驟來估計自回歸系數(shù)按如下步驟來估計自回歸系數(shù) （7-51） 3. 用（用（7-51）回歸得到的）回歸得到的，對（，對（7-49）做廣義差分方程：）做廣義差分方程：01tttYX （7-52）010011(1),對此式進行對此式進行OLS回歸即可得到回歸即可得到和和的估計值，然后注意到的估計值，然后注意到01,就可以得到

58、原模型（就可以得到原模型（7-49）中系數(shù)）中系數(shù)的估計值。的估計值。 2利用回歸殘差利用回歸殘差做如下做如下OLS回歸：回歸： te 1對（對（7-49）進行）進行OLS回歸得到回歸殘差回歸得到回歸殘差tetttee1 。 5. 現(xiàn)在估計回歸方程：現(xiàn)在估計回歸方程： （ 7-53）的第二次估計值。的第二次估計值。 這樣得到這樣得到 按如下步驟來估計自回歸系數(shù)按如下步驟來估計自回歸系數(shù)如果如果的第二次估計值仍然不能夠令人滿意，我們可以進行的第二次估計值仍然不能夠令人滿意，我們可以進行第四次第四次估計，一直到估計，一直到的估計值達到令人滿意的精度為止。的估計值達到令人滿意的精度為止。 的第三次、

59、的第三次、01, 4. 將第三步得到的將第三步得到的 的估計值重新代入原模型（的估計值重新代入原模型（7-49）并計算得到新）并計算得到新 的殘差的殘差tetttee1（4 4）杜賓兩步法）杜賓兩步法以上面的一元回歸模型為例以上面的一元回歸模型為例 （7-54) 0111(1)()tttttYXXYv 把廣義差分方程把廣義差分方程 改寫為：改寫為：01tttYX以下兩步以下兩步程序程序來估計來估計：1tY 1對（對（7-54）進行）進行OLS回歸，并把對回歸，并把對 的回歸系數(shù)的估計值的回歸系數(shù)的估計值 看作對看作對 的一個估計。雖然這個估計值有偏誤，但它卻是的一個估計。雖然這個估計值有偏誤，

60、但它卻是 的一個一致估計。的一個一致估計。 2求得求得 后，把它代入差分方程（后，把它代入差分方程（7-52），即代入下面的方程），即代入下面的方程10111(1)()()ttttttYYXX該方程改寫為該方程改寫為 （7-55）01tttYX 對（對（7-55）進行）進行OLS回歸得到參數(shù)的估計值。回歸得到參數(shù)的估計值。 由此可見，杜賓兩步法的第一步是要得到由此可見，杜賓兩步法的第一步是要得到的一個估計值，的一個估計值，第二步是要得到回歸的參數(shù)值。第二步是要得到回歸的參數(shù)值。還有一些其他的估計還有一些其他的估計 的方法，這里不再一一介紹。的方法，這里不再一一介紹。其他方法的基本上都是兩步法：