計(jì)量經(jīng)濟(jì)學(xué)第四章序列相關(guān)性詳解演示文稿

計(jì)量經(jīng)濟(jì)學(xué)第四章序列相關(guān)性詳解演示文稿目前一頁(yè)\總數(shù)四十八頁(yè)\編于十九點(diǎn)優(yōu)選計(jì)量經(jīng)濟(jì)學(xué)第四章序列相關(guān)性目前二頁(yè)\總數(shù)四十八頁(yè)\編于十九點(diǎn)一、序列相關(guān)性的概念如果對(duì)于不同的樣本點(diǎn)，隨機(jī)誤差項(xiàng)之間不再是不相關(guān)的，而是存在某種相關(guān)性，即：Cov(i

,j)≠0ij,i,j=1,2,…,n則認(rèn)為出現(xiàn)了序列相關(guān)性（serialcorrelation）。對(duì)于模型：

Yi=0+1X1i+2X2i+…+kXki+i

i=1,2,…,n目前三頁(yè)\總數(shù)四十八頁(yè)\編于十九點(diǎn)#序列相關(guān)性下的方差－協(xié)方差陣此時(shí)，隨機(jī)誤差項(xiàng)之間的方差－協(xié)方差陣為：目前四頁(yè)\總數(shù)四十八頁(yè)\編于十九點(diǎn)#自相關(guān)（autocorrelation）序列相關(guān)經(jīng)常出現(xiàn)在以時(shí)間序列數(shù)據(jù)為樣本的模型中，此時(shí)，不同樣本點(diǎn)的區(qū)別僅在于時(shí)間的不同這意味著，此時(shí)的序列相關(guān)性表現(xiàn)為不同時(shí)間上的隨機(jī)誤差項(xiàng)存在相關(guān)，這一情形下的序列相關(guān)也通常稱(chēng)之為自相關(guān)為此，本節(jié)將表示不同樣本點(diǎn)的下標(biāo)i

改為t。目前五頁(yè)\總數(shù)四十八頁(yè)\編于十九點(diǎn)如果僅存在：cov(t,t－1)=E(tt－1)0t=2,…,n即：隨機(jī)誤差項(xiàng)只與其前一期值有關(guān)（或者說(shuō)，僅是相鄰的隨機(jī)誤差項(xiàng)之間存在相關(guān)），則稱(chēng)為一階自相關(guān)。一階序列相關(guān)時(shí)，隨機(jī)誤差項(xiàng)可以表示為：t=t-1+t-1<<1

稱(chēng)為一階自回歸形式，記為AR(1)，其中：

：被稱(chēng)為一階自相關(guān)系數(shù)（first-ordercoefficientofautocorrelation）

i：滿足標(biāo)準(zhǔn)的OLS假定的隨機(jī)干擾項(xiàng)#一階自相關(guān)（first-orderautocorrelation）目前六頁(yè)\總數(shù)四十八頁(yè)\編于十九點(diǎn)序列相關(guān)的一般形式可以表示成：稱(chēng)為P階自回歸形式，記為AR(p)，表示模型存在P階自相關(guān)。

t－1、t－2、…、t－p分別表示t的前1期、前2期、…、前p期項(xiàng)，又稱(chēng)為滯后1期、滯后2期、…、滯后p期項(xiàng)。1、2、…，p稱(chēng)為1階、2階、…，p階自相關(guān)系數(shù)。#高階自相關(guān)（high-orderautocorrelation）目前七頁(yè)\總數(shù)四十八頁(yè)\編于十九點(diǎn)二、實(shí)際經(jīng)濟(jì)問(wèn)題中的序列相關(guān)性大多數(shù)經(jīng)濟(jì)時(shí)間數(shù)據(jù)都有一個(gè)明顯的特點(diǎn):慣性，表現(xiàn)在時(shí)間序列不同時(shí)間的前后關(guān)聯(lián)上。由于消費(fèi)習(xí)慣的影響被包含在隨機(jī)誤差項(xiàng)中，則可能出現(xiàn)序列相關(guān)性（往往是正相關(guān)）。例如：絕對(duì)收入假設(shè)下居民總消費(fèi)函數(shù)模型：

Ct=0+1Yt+tt=1,2,…,n1、經(jīng)濟(jì)變量固有的慣性序列相關(guān)性往往出現(xiàn)在以時(shí)間序列數(shù)據(jù)為樣本的模型中，產(chǎn)生這一問(wèn)題的原因主要來(lái)自三個(gè)方面：目前八頁(yè)\總數(shù)四十八頁(yè)\編于十九點(diǎn)許多經(jīng)濟(jì)行為存在滯后效應(yīng)，即當(dāng)期的經(jīng)濟(jì)行為不僅影響當(dāng)期的有關(guān)結(jié)果，而且也會(huì)對(duì)以后若干期的結(jié)果存在影響，這使得作為結(jié)果變量的經(jīng)濟(jì)變量在不同時(shí)間上呈現(xiàn)出序列相關(guān)性。例如：

固定資產(chǎn)的形成，不僅與當(dāng)期的固定資產(chǎn)投資有關(guān)，也與前期多年的固定資產(chǎn)投資有關(guān)

今年的家庭消費(fèi)水平，不僅與今年的收入有關(guān)，也與前期多年的收入有關(guān)以及前期多年的消費(fèi)支出有關(guān)

企業(yè)當(dāng)期的銷(xiāo)售收入，同樣會(huì)受到前期的商品銷(xiāo)售水平有關(guān)2、經(jīng)濟(jì)行為的滯后性目前九頁(yè)\總數(shù)四十八頁(yè)\編于十九點(diǎn)所謂模型設(shè)定偏誤（Specificationerror）是指所設(shè)定的模型“不正確”。主要表現(xiàn)在模型中丟掉了重要的解釋變量或模型函數(shù)形式有偏誤。例如：本來(lái)應(yīng)該估計(jì)的模型為：

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

Yt=0+1X1t+1X2t+vt因此：vt=3X3t+t，如果X3確實(shí)影響Y，則出現(xiàn)序列相關(guān)。3、模型設(shè)定的偏誤這是橫截面數(shù)據(jù)也可能存在序列相關(guān)性的重要原因目前十頁(yè)\總數(shù)四十八頁(yè)\編于十九點(diǎn)

例如：

季度數(shù)據(jù)來(lái)自月度數(shù)據(jù)的簡(jiǎn)單平均，這種平均的計(jì)算減弱了每月數(shù)據(jù)的波動(dòng)性，從而使隨機(jī)干擾項(xiàng)出現(xiàn)序列相關(guān)。還有就是兩個(gè)時(shí)間點(diǎn)之間的“內(nèi)插”技術(shù)往往導(dǎo)致隨機(jī)項(xiàng)的序列相關(guān)性。在實(shí)際經(jīng)濟(jì)問(wèn)題中，有些數(shù)據(jù)是通過(guò)已知數(shù)據(jù)生成的。因此，新生成的數(shù)據(jù)與原數(shù)據(jù)間就有了內(nèi)在的聯(lián)系，表現(xiàn)出序列相關(guān)性。4、數(shù)據(jù)的處理目前十一頁(yè)\總數(shù)四十八頁(yè)\編于十九點(diǎn)三、序列相關(guān)性的后果1、參數(shù)估計(jì)量仍然無(wú)偏，但非有效因?yàn)椋涸谟行宰C明中利用了：

E(NN’)=2I

即同方差性和互相獨(dú)立性條件。而且：在大樣本情況下，參數(shù)估計(jì)量雖然具有一致性，但仍然不具有漸近有效性。*通常情形下，采用OLS將會(huì)低估參數(shù)估計(jì)量的標(biāo)準(zhǔn)差，也會(huì)低估隨機(jī)誤差項(xiàng)的方差б2目前十二頁(yè)\總數(shù)四十八頁(yè)\編于十九點(diǎn)在變量的顯著性檢驗(yàn)中，統(tǒng)計(jì)量是建立在參數(shù)方差正確估計(jì)基礎(chǔ)之上的，這只有當(dāng)隨機(jī)誤差項(xiàng)具有同方差性和互相獨(dú)立性時(shí)才能成立。2、變量的顯著性檢驗(yàn)失去意義通常情況下，存在序列相關(guān)性時(shí)，參數(shù)估計(jì)值的樣本方差往往會(huì)被低估，此時(shí)變量t檢驗(yàn)和方程F檢驗(yàn)的顯著性容易被夸大！目前十三頁(yè)\總數(shù)四十八頁(yè)\編于十九點(diǎn)參數(shù)估計(jì)值非有效（真實(shí)方差往往被低估），失去最優(yōu)性，樣本估計(jì)式失準(zhǔn)隨機(jī)誤差項(xiàng)的方差一般會(huì)被低估區(qū)間預(yù)測(cè)與參數(shù)估計(jì)量的方差和隨機(jī)誤差項(xiàng)的方差均有關(guān)在方差有偏誤的情況下，使得預(yù)測(cè)估計(jì)不準(zhǔn)確，預(yù)測(cè)可信度降低。所以，當(dāng)模型出現(xiàn)序列相關(guān)性時(shí)，它的預(yù)測(cè)功能失效。3、模型的預(yù)測(cè)失效目前十四頁(yè)\總數(shù)四十八頁(yè)\編于十九點(diǎn)然后，通過(guò)分析這些“近似估計(jì)量”之間的相關(guān)性，以判斷隨機(jī)誤差項(xiàng)是否具有序列相關(guān)性?；舅悸罚核?、序列相關(guān)性的檢驗(yàn)?zāi)壳笆屙?yè)\總數(shù)四十八頁(yè)\編于十九點(diǎn)（一）圖示檢驗(yàn)法目前十六頁(yè)\總數(shù)四十八頁(yè)\編于十九點(diǎn)（二）回歸檢驗(yàn)法……

如果存在某一種函數(shù)形式，使得方程顯著成立，則說(shuō)明原模型存在序列相關(guān)性。優(yōu)點(diǎn)：（1）能夠確定序列相關(guān)的形式;（2）適用于任何類(lèi)型序列相關(guān)性問(wèn)題的檢驗(yàn)。缺點(diǎn)：工作量大，計(jì)算復(fù)雜，檢驗(yàn)繁瑣目前十七頁(yè)\總數(shù)四十八頁(yè)\編于十九點(diǎn)（三）杜賓-瓦森檢驗(yàn)法（DW檢驗(yàn)）D-W檢驗(yàn)是杜賓（J.Durbin）和瓦森(G.S.Watson)于1951年提出的一種檢驗(yàn)序列自相關(guān)的方法該方法只適用于檢驗(yàn)一階自相關(guān)（1）解釋變量X非隨機(jī)；（2）隨機(jī)誤差項(xiàng)t為一階自回歸形式：

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

Yt=0+1X1t+kXkt+Yt-1+t（4）回歸含有截距項(xiàng)假定條件目前十八頁(yè)\總數(shù)四十八頁(yè)\編于十九點(diǎn)該統(tǒng)計(jì)量的分布與出現(xiàn)在給定樣本中的X值有復(fù)雜的關(guān)系，因此其精確的分布很難得到。但是，他們成功地導(dǎo)出了臨界值的下限dL和上限dU，且這些上下限只與樣本的容量n

和解釋變量的個(gè)數(shù)k

有關(guān)，而與解釋變量X的取值無(wú)關(guān)。杜賓和瓦森針對(duì)原假設(shè)：H0:=0，即不存在一階自回歸，構(gòu)造如下統(tǒng)計(jì)量：＃D.W.檢驗(yàn)統(tǒng)計(jì)量目前十九頁(yè)\總數(shù)四十八頁(yè)\編于十九點(diǎn)D.W.檢驗(yàn)步驟（1）提出假設(shè)：H0:ρ=0（不存在一階自相關(guān)）H1:ρ≠0（2）計(jì)算DW值（3）給定，由n和（k＋1）的大小查DW分布表，得臨界值dL和dU（4）比較、判斷0<D.W.<dL

存在正自相關(guān)dL<D.W.<dU

不能確定dU<D.W.<4－dU無(wú)自相關(guān)4－dU<D.W.<4－dL

不能確定4－dL<D.W.<4存在負(fù)自相關(guān)目前二十頁(yè)\總數(shù)四十八頁(yè)\編于十九點(diǎn)#DW檢驗(yàn)的圖示0dLdU24-dU4-dL

正相關(guān)不能確定無(wú)自相關(guān)不能確定負(fù)相關(guān)目前二十一頁(yè)\總數(shù)四十八頁(yè)\編于十九點(diǎn)

證明：展開(kāi)D.W.統(tǒng)計(jì)量：(*)#D.W.檢驗(yàn)統(tǒng)計(jì)量的說(shuō)明DW檢驗(yàn)表明：當(dāng)D.W.值在2左右時(shí)，模型不存在一階自相關(guān)其中：ρ為一階自相關(guān)系數(shù)一階自回歸模型：i=i-1+i

的參數(shù)估計(jì)。目前二十二頁(yè)\總數(shù)四十八頁(yè)\編于十九點(diǎn)由于自相關(guān)系數(shù)的值介于－1和＋1之間，因此：0≤DW≈2(1-ρ)≤4如果存在完全一階正相關(guān)，即=1，則D.W.0

完全一階負(fù)相關(guān)，即=-1,則D.W.4

完全不相關(guān)，即=0，則D.W.20dLdU24-dU4-dL

正相關(guān)不能確定無(wú)自相關(guān)不能確定負(fù)相關(guān)目前二十三頁(yè)\總數(shù)四十八頁(yè)\編于十九點(diǎn)DW檢驗(yàn)是最常用的自相關(guān)性的檢驗(yàn)方法，在報(bào)告回歸分析的結(jié)果時(shí)，一般將DW值連同R2、t值等一起標(biāo)明。但在應(yīng)用DW檢驗(yàn)時(shí)需要注意：1）DW值接近于2時(shí)，只能說(shuō)明模型不存在一階線性自相關(guān)，但并不意味著模型不存在高階自相關(guān)或者非線性相關(guān)2）DW值落入兩個(gè)無(wú)法判斷的區(qū)域時(shí)，需要采用其它檢驗(yàn)方法3）不適用于聯(lián)立方程組模型中各單一方程隨機(jī)誤差項(xiàng)序列相關(guān)的檢驗(yàn)4）DW檢驗(yàn)不適用于模型中含有滯后被解釋變量的情況，即不適用于如下模型Yt=0+1X1t++kXkt+Yt-1+t＃使用D.W.檢驗(yàn)時(shí)需要注意的問(wèn)題目前二十四頁(yè)\總數(shù)四十八頁(yè)\編于十九點(diǎn)針對(duì)滯后變量模型：Yt=0+1X1t++kXkt+Yt-1+t上述模型，Durbin提出Durbin－h(huán)統(tǒng)計(jì)量：#DH統(tǒng)計(jì)量目前二十五頁(yè)\總數(shù)四十八頁(yè)\編于十九點(diǎn)（四）拉格朗日乘數(shù)檢驗(yàn)（LagrangeMultiplier）LM檢驗(yàn)是由布勞殊（Breusch）與戈弗雷（Godfrey）于1978年提出的，也被稱(chēng)為GB檢驗(yàn)。拉格朗日乘數(shù)檢驗(yàn)克服了DW檢驗(yàn)的缺陷，適合于高階序列相關(guān)以及模型中存在滯后被解釋變量的情形。對(duì)于模型如果懷疑隨機(jī)擾動(dòng)項(xiàng)存在p階序列相關(guān)，即隨機(jī)誤差項(xiàng)存在：目前二十六頁(yè)\總數(shù)四十八頁(yè)\編于十九點(diǎn)則構(gòu)造以下輔助回歸模型：在原假設(shè)：H0:1=2=…=p=0（無(wú)序列相關(guān)）成立時(shí)，有：其中：n為輔助回歸樣本容量，R2為輔助回歸的可決系數(shù)：給定，查臨界值2(p)，與LM值比較，如果超出則拒絕H0實(shí)際檢驗(yàn)中，可從1階、2階、…逐次向更高階檢驗(yàn)。目前二十七頁(yè)\總數(shù)四十八頁(yè)\編于十九點(diǎn)檢驗(yàn)時(shí)需要事先確定準(zhǔn)備檢驗(yàn)的階數(shù)P,實(shí)際檢驗(yàn)中，可從1階、2階、…逐次向更高階檢驗(yàn)。檢驗(yàn)結(jié)果顯著時(shí)，可以說(shuō)明存在序列相關(guān)，但是并不一定代表序列相關(guān)的階數(shù)一定能夠達(dá)到所檢驗(yàn)的階數(shù)。低階序列相關(guān)的存在往往會(huì)導(dǎo)致高階序列相關(guān)檢驗(yàn)的顯著性具體階數(shù)的判斷，需要結(jié)合輔助回歸中自相關(guān)系數(shù)的顯著性＃使用GB檢驗(yàn)時(shí)需要注意的問(wèn)題(0.22)(-0.497)(4.541)（-1.842）（0.087）

R2=0.6615目前二十八頁(yè)\總數(shù)四十八頁(yè)\編于十九點(diǎn)如果模型被檢驗(yàn)證明存在序列相關(guān)性，則首先需要分析其原因，對(duì)癥下藥：如果產(chǎn)生序列相關(guān)的原因是變量選擇失準(zhǔn)（如遺漏了重要的解釋變量等），則應(yīng)調(diào)整變量；如果是模型設(shè)定不當(dāng)，應(yīng)當(dāng)調(diào)整模型形式。——虛假的序列相關(guān)問(wèn)題如果原因在于客觀經(jīng)濟(jì)現(xiàn)象的自身特點(diǎn)，如經(jīng)濟(jì)變量的慣性作用等，則需要發(fā)展新的估計(jì)方法最常用的方法是廣義最小二乘法（GLS:Generalizedleastsquares）和廣義差分法(GD，GeneralizedDifference)。五、序列相關(guān)性的補(bǔ)救目前二十九頁(yè)\總數(shù)四十八頁(yè)\編于十九點(diǎn)（一）廣義最小二乘法對(duì)于模型：Y=X+

（X為設(shè)計(jì)矩陣，Y、β、μ為列向量）如果存在序列相關(guān)，同時(shí)存在異方差，即有：是一對(duì)稱(chēng)正定矩陣，存在一可逆矩陣D，使得：

=DD’廣義最小二乘法（GLS)是最具有普遍意義的最小二乘法，普通最小二乘法（OLS）和加權(quán)最小二乘法（WLS）是其特例目前三十頁(yè)\總數(shù)四十八頁(yè)\編于十九點(diǎn)變換原模型（D-1左乘）：D-1Y=D-1X+D-1

即：Y*=X*+*

(*)(*)式的OLS估計(jì)：此即原模型的廣義最小二乘估計(jì)量(GLSE)，是無(wú)偏的、有效的估計(jì)量。(*)模型具有同方差性和無(wú)序列相關(guān)性，因?yàn)椋耗壳叭豁?yè)\總數(shù)四十八頁(yè)\編于十九點(diǎn)＃如何得到矩陣？——近似估計(jì)矩陣是原模型隨機(jī)誤差項(xiàng)的方差－協(xié)方差陣。獲得的一種方法是采用隨機(jī)誤差項(xiàng)的近似估計(jì)量構(gòu)造目前三十二頁(yè)\總數(shù)四十八頁(yè)\編于十九點(diǎn)獲取的更精確的方法是根據(jù)原模型序列相關(guān)的具體形式進(jìn)行估計(jì)常見(jiàn)的是假設(shè)隨機(jī)誤差項(xiàng)具有一階序列相關(guān)性，即：

i=i-1+i(－1<<1)此時(shí)，可以證明：＃如何得到矩陣？——精確估計(jì)目前三十三頁(yè)\總數(shù)四十八頁(yè)\編于十九點(diǎn)證明：由：i=i-1+i(－1<<1)有：即：由：有：目前三十四頁(yè)\總數(shù)四十八頁(yè)\編于十九點(diǎn)目前三十五頁(yè)\總數(shù)四十八頁(yè)\編于十九點(diǎn)廣義差分法是利用廣義差分變換將原模型變換為滿足基本假設(shè)的差分模型，再進(jìn)行OLS估計(jì)。是一類(lèi)克服序列相關(guān)性的有效方法，被廣泛采用。對(duì)于模型:將模型滯后一期，有：同理，模型滯后p期的形式為：（二）廣義差分法目前三十六頁(yè)\總數(shù)四十八頁(yè)\編于十九點(diǎn)如果模型存在:對(duì)模型施行廣義P階差分變換，有：該模型為原模型的廣義差分模型，不存在序列相關(guān)問(wèn)題，可進(jìn)行OLS估計(jì)，從而獲得原模型的最佳估計(jì)量，即：目前三十七頁(yè)\總數(shù)四十八頁(yè)\編于十九點(diǎn)廣義差分法實(shí)質(zhì)上與廣義最小二乘法是一致的，只是GD法中損失了部分樣本觀測(cè)值。這相當(dāng)于GLS中的D-1去掉第一行后左乘原模型：Y=X+

＃GD和GLS的關(guān)系如：一階序列相關(guān)的情況下,廣義差分是估計(jì)即運(yùn)用了GLS法，但第一次觀測(cè)值被排除了。則GD與GLS完全等價(jià)。（普萊斯－溫斯特變換）目前三十八頁(yè)\總數(shù)四十八頁(yè)\編于十九點(diǎn)（三）隨機(jī)誤差項(xiàng)的自相關(guān)系數(shù)ρ的估計(jì)應(yīng)用廣義最小二乘法或廣義差分法，必須已知隨機(jī)誤差項(xiàng)的自相關(guān)系數(shù)1,

2,…,

L。實(shí)際上，人們并不知道它們的具體數(shù)值，所以必須首先對(duì)它們進(jìn)行估計(jì)。常用的估計(jì)方法有：（1）利用DW統(tǒng)計(jì)量進(jìn)行近似估計(jì)（2）科克倫-奧科特（Cochrane-Orcutt）迭代法。（3）杜賓（durbin）兩步法目前三十九頁(yè)\總數(shù)四十八頁(yè)\編于十九點(diǎn)（1）利用DW統(tǒng)計(jì)量進(jìn)行估計(jì)對(duì)于一階自相關(guān)：由：有：僅適用于一階自相關(guān)情形，用于構(gòu)建一階差分模型。所估計(jì)的為一階自相關(guān)系數(shù)目前四十頁(yè)\總數(shù)四十八頁(yè)\編于十九點(diǎn)（2）科克倫-奧科特迭代法以一元線性模型為例。1）首先采用OLS法估計(jì)原模型：Yt=0+1Xt+t得到的i的“近似估計(jì)值”et(1)，2）以et(1)作為觀測(cè)值使用OLS法估計(jì)下式：

t=1t-1+2t-2+pt-p+t

得到：作為隨機(jī)誤差項(xiàng)的相關(guān)系數(shù)的第一次估計(jì)值目前四十一頁(yè)\總數(shù)四十八頁(yè)\編于十九點(diǎn)求出i新的“近似估計(jì)值”et(2)3）利用構(gòu)造廣義差分模型：進(jìn)行OLS估計(jì)，得到回歸系數(shù)的估計(jì)4）將代回原模型：Yt=0+1Xt+t6）重復(fù)上述步驟，可得相關(guān)系數(shù)的多次迭代值。注：具體迭代次數(shù)，可根據(jù)具體問(wèn)題來(lái)定。一般可事先規(guī)定一個(gè)精度δ，當(dāng)時(shí)，迭代終止。實(shí)際中一般只需要迭代兩次即可。因此上述方法又稱(chēng)為科克倫－奧科特兩步法。5）重復(fù)步驟2），得到相關(guān)系數(shù)的第二次估計(jì)值：目前四十二頁(yè)\總數(shù)四十八頁(yè)\編于十九點(diǎn)#科克倫-奧科特迭代法圖示目前四十三頁(yè)\總數(shù)四十八頁(yè)\編于十九點(diǎn)第一步：變換差分模型為下列形式進(jìn)行OLS估計(jì)，得各Yj（j=t-1,