序列相關(guān)性教材

序列相關(guān)性

SerialCorrelation一、序列相關(guān)性概念二、實際經(jīng)濟(jì)問題中的序列相關(guān)性

三、序列相關(guān)性的后果四、序列相關(guān)性的檢驗五、具有序列相關(guān)性模型的估計六、案例序列相關(guān)性

一、序列相關(guān)性概念

如果對于不同的樣本點，隨機(jī)誤差項之間不再是不相關(guān)的，而是存在某種相關(guān)性，則認(rèn)為出現(xiàn)了序列相關(guān)性。

對于模型

Yi=

0+

1X1i+

2X2i+…+

kXki+

i

i=1,2,…,n隨機(jī)項互不相關(guān)的基本假設(shè)表現(xiàn)為

Cov(

i

,

j)=0

i

j,i,j=1,2,…,n或稱為一階列相關(guān)，或自相關(guān)（autocorrelation）其中：

被稱為自協(xié)方差系數(shù)（coefficientofautocovariance）或一階自相關(guān)系數(shù)（first-ordercoefficientofautocorrelation）

i是滿足以下標(biāo)準(zhǔn)的OLS假定的隨機(jī)干擾項：如果僅存在

E(

i

i+1)0

i=1,2,…,n

自相關(guān)往往可寫成如下形式：

i=

i-1+

i-1<

<1

由于序列相關(guān)性經(jīng)常出現(xiàn)在以時間序列為樣本的模型中，因此，本節(jié)將用下標(biāo)t代表i。

（a）正自相關(guān)6（b）負(fù)自相關(guān)7（c）無自相關(guān)8

二、實際經(jīng)濟(jì)問題中的序列相關(guān)性

大多數(shù)經(jīng)濟(jì)時間數(shù)據(jù)都有一個明顯的特點:慣性，表現(xiàn)在時間序列不同時間的前后關(guān)聯(lián)上。由于消費習(xí)慣的影響被包含在隨機(jī)誤差項中，則可能出現(xiàn)序列相關(guān)性（往往是正相關(guān)）。例如，絕對收入假設(shè)下居民總消費函數(shù)模型：

Ct=

0+1Yt+tt=1,2,…,n

1、經(jīng)濟(jì)變量固有的慣性

2、模型設(shè)定的偏誤

所謂模型設(shè)定偏誤（Specificationerror）是指所設(shè)定的模型“不正確”。主要表現(xiàn)在模型中丟掉了重要的解釋變量或模型函數(shù)形式有偏誤。

例如，本來應(yīng)該估計的模型為

Yt=

0+1X1t+2X2t+3X3t+t但在模型設(shè)定中做了下述回歸：

Yt=

0+1X1t+1X2t+vt因此，vt=

3X3t+t，如果X3確實影響Y，則出現(xiàn)序列相關(guān)。

但建模時設(shè)立了如下模型：

Yt=

0+1Xt+vt

因此，由于vt=

2Xt2+t,

，包含了產(chǎn)出的平方對隨機(jī)項的系統(tǒng)性影響，隨機(jī)項也呈現(xiàn)序列相關(guān)性。又如：如果真實的邊際成本回歸模型應(yīng)為：

Yt=

0+1Xt+2Xt2+t其中：Y=邊際成本，X=產(chǎn)出，

3、數(shù)據(jù)的“編造”

例如：季度數(shù)據(jù)來自月度數(shù)據(jù)的簡單平均，這種平均的計算減弱了每月數(shù)據(jù)的波動性，從而使隨機(jī)干擾項出現(xiàn)序列相關(guān)。還有就是兩個時間點之間的“內(nèi)插”技術(shù)往往導(dǎo)致隨機(jī)項的序列相關(guān)性。

在實際經(jīng)濟(jì)問題中，有些數(shù)據(jù)是通過已知數(shù)據(jù)生成的。因此，新生成的數(shù)據(jù)與原數(shù)據(jù)間就有了內(nèi)在的聯(lián)系，表現(xiàn)出序列相關(guān)性。

計量經(jīng)濟(jì)學(xué)模型一旦出現(xiàn)序列相關(guān)性，如果仍采用OLS法估計模型參數(shù)，會產(chǎn)生下列不良后果：

二、序列相關(guān)性的后果1、參數(shù)估計量非有效

因為，在有效性證明中利用了

E(NN’)=

2I

即同方差性和互相獨立性條件。而且，在大樣本情況下，參數(shù)估計量雖然具有一致性，但仍然不具有漸近有效性。2、變量的顯著性檢驗失去意義

在變量的顯著性檢驗中，統(tǒng)計量是建立在參數(shù)方差正確估計基礎(chǔ)之上的，這只有當(dāng)隨機(jī)誤差項具有同方差性和互相獨立性時才能成立。

其他檢驗也是如此。3、模型的預(yù)測失效

區(qū)間預(yù)測與參數(shù)估計量的方差有關(guān)，在方差有偏誤的情況下，使得預(yù)測估計不準(zhǔn)確，預(yù)測精度降低。所以，當(dāng)模型出現(xiàn)序列相關(guān)性時，它的預(yù)測功能失效。三、序列相關(guān)性的檢驗

然后，通過分析這些“近似估計量”之間的相關(guān)性，以判斷隨機(jī)誤差項是否具有序列相關(guān)性。

序列相關(guān)性檢驗方法有多種，但基本思路相同：

基本思路:

三、序列相關(guān)性的檢驗1、圖示法2、回歸檢驗法……

如果存在某一種函數(shù)形式，使得方程顯著成立，則說明原模型存在序列相關(guān)性。

回歸檢驗法的優(yōu)點是：（1）能夠確定序列相關(guān)的形式，（2）適用于任何類型序列相關(guān)性問題的檢驗。3、杜賓-瓦森（Durbin-Watson）檢驗法

D-W檢驗是杜賓（J.Durbin）和瓦森(G.S.Watson)于1951年提出的一種檢驗序列自相關(guān)的方法，該方法的假定條件是：（1）解釋變量X非隨機(jī)；（2）隨機(jī)誤差項

i為一階自回歸形式：

i=i-1+i（3）回歸模型中不應(yīng)含有滯后應(yīng)變量作為解釋變量，即不應(yīng)出現(xiàn)下列形式：

Yi=

0+1X1i+kXki+Yi-1+i（4）回歸含有截距項

該統(tǒng)計量的分布與出現(xiàn)在給定樣本中的X值有復(fù)雜的關(guān)系，因此其精確的分布很難得到。但是，他們成功地導(dǎo)出了臨界值的下限dL和上限dU

，且這些上下限只與樣本的容量n和解釋變量的個數(shù)k有關(guān)，而與解釋變量X的取值無關(guān)。

杜賓和瓦森針對原假設(shè)：H0:=0，即不存在一階自回歸，構(gòu)如下造統(tǒng)計量：

D.W.統(tǒng)計量:D.W檢驗步驟:（1）計算DW值（2）給定，由n和k的大小查DW分布表，得臨界值dL和dU（3）比較、判斷

若0<D.W.<dL

存在正自相關(guān)

dL<D.W.<dU

不能確定

dU<D.W.<4－dU

無自相關(guān)

4－dU<D.W.<4－dL

不能確定

4－dL<D.W.<4存在負(fù)自相關(guān)

0dLdU24-dU4-dL

正相關(guān)不能確定無自相關(guān)不能確定負(fù)相關(guān)

當(dāng)D.W.值在2左右時，模型不存在一階自相關(guān)。

證明：展開D.W.統(tǒng)計量：

(*)如果存在完全一階正相關(guān)，即

=1，則D.W.0

完全一階負(fù)相關(guān)，即

=-1,則D.W.4

完全不相關(guān)，即

=0，則D.W.2這里，為一階自回歸模型

i=

i-1+

i的參數(shù)估計。4、拉格朗日乘數(shù)（Lagrangemultiplier）檢驗

拉格朗日乘數(shù)檢驗克服了DW檢驗的缺陷，適合于高階序列相關(guān)以及模型中存在滯后被解釋變量的情形。它是由布勞殊（Breusch）與戈弗雷（Godfrey）于1978年提出的，也被稱為GB檢驗。

對于模型如果懷疑隨機(jī)擾動項存在p階序列相關(guān)：

GB檢驗可用來檢驗如下受約束回歸方程

約束條件為：

H0:

1=

2=…=

p=0約束條件H0為真時，大樣本下其中，n為樣本容量，R2為如下輔助回歸的可決系數(shù)：

給定

，查臨界值

2(p)，與LM值比較，做出判斷，實際檢驗中，可從1階、2階、…逐次向更高階檢驗。

如果模型被檢驗證明存在序列相關(guān)性，則需要發(fā)展新的方法估計模型。

最常用的方法是廣義最小二乘法（GLS:Generalizedleastsquares）和廣義差分法(GeneralizedDifference)。四、序列相關(guān)的補(bǔ)救

1、廣義最小二乘法

對于模型

Y=X

+

如果存在序列相關(guān)，同時存在異方差，即有

是一對稱正定矩陣，存在一可逆矩陣D，使得

=DD’變換原模型：

D-1Y=D-1X

+D-1

即

Y*=X*

+

*(*)(*)式的OLS估計：

這就是原模型的廣義最小二乘估計量(GLSestimators)，是無偏的、有效的估計量。

該模型具有同方差性和隨機(jī)誤差項互相獨立性：

如何得到矩陣

？

對

的形式進(jìn)行特殊設(shè)定后，才可得到其估計值。

如設(shè)定隨機(jī)擾動項為一階序列相關(guān)形式

i=

i-1+

i則2、廣義差分法

廣義差分法是將原模型變換為滿足OLS法的差分模型，再進(jìn)行OLS估計。如果原模型存在可以將原模型變換為:

該模型為廣義差分模型，不存在序列相關(guān)問題?？蛇M(jìn)行OLS估計。

注意：

廣義差分法就是上述廣義最小二乘法，但是卻損失了部分樣本觀測值。

如：一階序列相關(guān)的情況下,廣義差分是估計這相當(dāng)于去掉第一行后左乘原模型Y=X

+

。即運用了GLS法，但第一次觀測值被排除了。

3、隨機(jī)誤差項相關(guān)系數(shù)的估計

應(yīng)用廣義最小二乘法或廣義差分法，必須已知隨機(jī)誤差項的相關(guān)系數(shù)

1,

2,…,

L

。實際上，人們并不知道它們的具體數(shù)值，所以必須首先對它們進(jìn)行估計。

常用的估計方法有：

科克倫-奧科特（Cochrane-Orcutt）迭代法。杜賓（durbin）兩步法

（1）科克倫-奧科特迭代法。

以一元線性模型為例：首先，采用OLS法估計原模型

Yi=

0+

1Xi+

i得到的

的“近似估計值”，并以之作為觀測值使用OLS法估計下式

i=

1

i-1+

2

i-2+

L

i-L+

i求出

i新的“近擬估計值”，

并以之作為樣本觀測值，再次估計

i=

1

i-1+

2

i-2+

L

i-L+

i

類似地，可進(jìn)行第三次、第四次迭代。

關(guān)于迭代的次數(shù)，可根據(jù)具體的問題來定。一般是事先給出一個精度，當(dāng)相鄰兩次

1,

2,,

L的估計值之差小于這一精度時，迭代終止。實踐中，有時只要迭代兩次，就可得到較滿意的結(jié)果。兩次迭代過程也被稱為科克倫-奧科特兩步法。（2）杜賓（durbin）兩步法

該方法仍是先估計

1,

2,,

l，再對差分模型進(jìn)行估計

第一步，變換差分模型為下列形式進(jìn)行OLS估計，得各Yj（j=i-1,i-2,…,i-l)前的系數(shù)

1,

2,,

l的估計值五、案例：中國商品進(jìn)口模型

經(jīng)濟(jì)理論指出，商品進(jìn)口主要由進(jìn)口國的經(jīng)濟(jì)發(fā)展水平，以及商品進(jìn)口價格指數(shù)與國內(nèi)價格指數(shù)對比因素決定的。由于無法取得中國商品進(jìn)口價格指數(shù)，我們主要研究中國商品進(jìn)口與國內(nèi)生產(chǎn)總值的關(guān)系。（下表）。

1.通過OLS法建立如下中國商品進(jìn)口方程：

（2.32）（20.12）

2.進(jìn)行序列相關(guān)性檢驗。

DW檢驗

取

=5%，由于n=24，k=2(包含常數(shù)項)，查表得：

dl=