第15章時(shí)間序列回歸

1、1 本章我們討論分析時(shí)間序列數(shù)據(jù)（檢驗(yàn)序列相關(guān)性，估 計(jì)arma模型，使用分布滯后，非平穩(wěn)時(shí)間序列的單位根檢 驗(yàn)）的單方程回歸方法。 2 時(shí)間序列回歸中的一個(gè)普遍現(xiàn)象是：殘差和它自己的滯后值相關(guān)。這種 序列相關(guān)性違背了回歸理論的標(biāo)準(zhǔn)假設(shè)：不同時(shí)點(diǎn)的擾動(dòng)項(xiàng)互不相關(guān)。與序 列相關(guān)相聯(lián)系的主要問(wèn)題有： 在線性估計(jì)中ols不再是有效的； 使用ols公式計(jì)算出的標(biāo)準(zhǔn)差不正確； 如果在方程右邊有滯后因變量，ols估計(jì)是有偏的且不一致。 eviews提供了檢測(cè)序列相關(guān)和估計(jì)方法的工具。但首先必須排除虛假序 列相關(guān)。虛假序列相關(guān)是指模型的序列相關(guān)是由于省略了顯著的解釋變量而 引起的。例如，在生產(chǎn)函數(shù)模型中，如

2、果省略了資本這個(gè)重要的解釋變量， 資本對(duì)產(chǎn)出的影響就被歸入隨機(jī)誤差項(xiàng)。由于資本在時(shí)間上的連續(xù)性，以及 對(duì)產(chǎn)出影響的連續(xù)性，必然導(dǎo)致隨機(jī)誤差項(xiàng)的序列相關(guān)。所以在這種情況下， 要把顯著的變量引入到解釋變量中。 3 如果隨機(jī)過(guò)程 的均值和方 差、自協(xié)方差都不取決于 t，則稱 y t 是協(xié)方差平穩(wěn)的或弱平穩(wěn)的： , 12101 ttt yyyyyyy )( t ye 2 )( t yvar 對(duì)所有的 t 對(duì)所有的 t 對(duì)所有的 t 和 s 注意，如果一個(gè)隨機(jī)過(guò)程是弱平穩(wěn)的，則y t與y t- s之間的協(xié)方差僅取決 于s ，即僅與觀測(cè)值之間的間隔長(zhǎng)度s有關(guān)，而與時(shí)期t 無(wú)關(guān)。一般所說(shuō)的 “平穩(wěn)性”含義就

3、是上述的弱平穩(wěn)定義。給定一個(gè)樣本值為t 的時(shí)間序列可 以看作是隨機(jī)過(guò)程 y t 的一個(gè)實(shí)現(xiàn)，仍記為 。 , 21tt yyyy sstt yye )( 4 一般地，我們考慮如下形式： ttt uxy ttt zu 1 是在t時(shí)刻的解釋變量向量； 是前期已知變量向量； 是參數(shù)向量； 是殘差； 是殘差的擾動(dòng)項(xiàng)； 可能包含 的滯后值或 的滯后值。 是無(wú)條件殘差,它是基于結(jié)構(gòu)成分 的殘差，但它不使用 中包 含的信息。 是一步預(yù)測(cè)誤差，它是因變量真實(shí)值和以解釋變量以及以前預(yù)測(cè)誤 差為基礎(chǔ)的預(yù)測(cè)值之差。 t x 1t z, t u 1t z ),( t x t u t t t u 1t z t 5 最簡(jiǎn)單

4、且最常用的序列相關(guān)模型是一階自回歸ar(1)模型。定義如下： ttt uxy ttt uu 1 參數(shù) 是一階序列相關(guān)系數(shù)，實(shí)際上，ar(1)模型是將以前觀測(cè)值的殘差 包含到現(xiàn)觀測(cè)值的回歸模型中。 更為一般，帶有p階自回歸的模型，ar(p)誤差由下式給出： ttt uxy tptpttt uuuu 2211 ar(p)的自相關(guān)將漸漸衰減至零，同時(shí)高于p階的偏自相關(guān)也是零。 6 在使用估計(jì)方程進(jìn)行統(tǒng)計(jì)推斷（如假設(shè)檢驗(yàn)和預(yù)測(cè)）之前，一般應(yīng)檢驗(yàn) 殘差（序列相關(guān)的證據(jù)），eviews提供了幾種方法來(lái)檢驗(yàn)當(dāng)前序列相關(guān)。 eviews將d-w統(tǒng)計(jì)量視為標(biāo)準(zhǔn)回歸輸出的一部分。 還可估算回歸模型鄰近殘差的線性

5、聯(lián)系。d-w統(tǒng)計(jì)量是在下面定義中檢驗(yàn)原假設(shè): 0 ttt uu 1 如果序列不相關(guān)，d-w值在2附近。如果存在正序列相關(guān)，d-w值將小于 2（最小為0），如果存在負(fù)序列相關(guān)，d-w值將在2 - 4之間。 正序列相關(guān)最為普遍，根據(jù)經(jīng)驗(yàn)，對(duì)于有大于50個(gè)數(shù)據(jù)和較少的解釋變 量，d-w值小于1.5的情況，說(shuō)明存在強(qiáng)正一階序列相關(guān)。參考johnston and dinardo（1997版6.6.1章）關(guān)于d-w檢驗(yàn)和統(tǒng)計(jì)量顯著性的論述。 7 1d-w統(tǒng)計(jì)量的擾動(dòng)項(xiàng)在原假設(shè)下依賴于數(shù)據(jù)矩陣x。 2回歸方程右邊如果存在滯后因變量，d-w檢驗(yàn)不再有效。 3僅僅檢驗(yàn)原假設(shè)（無(wú)序列相關(guān)）與備選假設(shè)（一階序列相關(guān)

6、）。 其他兩種檢驗(yàn)序列相關(guān)方法：q-統(tǒng)計(jì)量和breush-godfrey lm檢驗(yàn)克 服了上述不足，應(yīng)用于大多數(shù)場(chǎng)合。 例子：工作文件15_1eq_cs 8 在方程工具欄選擇view/residual tests/correlogram-q-statistics 。 eviews將顯示殘差的自相關(guān)和偏自相關(guān)函數(shù)以及對(duì)應(yīng)于高階序列相關(guān)的 ljung-box q統(tǒng)計(jì)量。如果殘差不存在序列相關(guān)，在各階滯后的自相關(guān)和 偏自相關(guān)值都接近于零。所有的q-統(tǒng)計(jì)量不顯著，并且有大的p值。 k 階滯后的。計(jì)算 式如下 k j j lb jt r ttq 1 2 2 是 j 階自相關(guān)系數(shù)，t是觀測(cè)值的個(gè)數(shù)。 j

7、r 9 下面是這些檢驗(yàn)程序應(yīng)用的例子，考慮用普通最小二乘估計(jì)的簡(jiǎn)單消費(fèi) 函數(shù)的結(jié)果： 10 瀏覽這些結(jié)果：系數(shù)在統(tǒng)計(jì)上是很顯著的，并且擬合得很好。但是，如果 誤差項(xiàng)是序列相關(guān)的，那么估計(jì)ols標(biāo)準(zhǔn)誤差將是無(wú)效的，并且估計(jì)系數(shù)由于 在方程右端有滯后因變量會(huì)發(fā)生偏倚和不一致。在這種情況下d-w統(tǒng)計(jì)量作為 序列相關(guān)的檢驗(yàn)是不合適的，因?yàn)樵诜匠逃叶舜嬖谥粋€(gè)滯后因變量。選擇 view/residual test/correlogram-q-statistice會(huì)產(chǎn)生如下情況 11 選擇view/residual tests/serial correlation lm test，一般地 對(duì)高階的，含有a

8、rma誤差項(xiàng)的情況執(zhí)行breush-godfrey lm （lagrange multiplier，拉格朗日乘數(shù)檢驗(yàn)）。在滯后定義對(duì)話 框，輸入要檢驗(yàn)序列的最高階數(shù)。 eviews將給出兩個(gè)統(tǒng)計(jì)量：f統(tǒng)計(jì)量和nr2（觀測(cè)值個(gè)數(shù)乘 以r2），nr2在原假設(shè)下服從 分布。f統(tǒng)計(jì)量分布未知，但常 用來(lái)對(duì)原假設(shè)進(jìn)行非正規(guī)檢驗(yàn)。 2 12 上一例子中相關(guān)圖在滯后值3時(shí)出現(xiàn)峰值。q統(tǒng)計(jì)量在各階滯后值中都 具有顯著性，它顯示的是殘差中的顯著序列相關(guān)。 進(jìn)行序列相關(guān)的lm檢驗(yàn)，選擇view/residual tests/serial correlation lm test，輸入滯后2產(chǎn)生如下結(jié)果： 此檢驗(yàn)拒絕

9、直至2階的無(wú)序列相關(guān)的假設(shè)。q-統(tǒng)計(jì)和lm檢驗(yàn)都表明： 殘差是序列相關(guān)的，并且方程在被用于假設(shè)檢驗(yàn)和預(yù)測(cè)之前應(yīng)該重新定 義。 13 在使用本章描述的工具之前，可以首先檢驗(yàn)?zāi)Ｐ推渌矫娴腻e(cuò)誤。誤差存 在序列相關(guān)是模型定義存在的嚴(yán)重問(wèn)題。特別地，應(yīng)注意使用ols得出的過(guò)分 限制的定義。有時(shí)，在回歸方程中添加不應(yīng)被排除的變量會(huì)消除序列相關(guān)。 在eviews中估計(jì)一個(gè)ar(1)模型，選擇quick/estimate equation打開(kāi)一個(gè)方 程，用列表法輸入方程后，最后將ar(1)項(xiàng)加到列表中。例如：估計(jì)一個(gè)帶有 ar(1)誤差的簡(jiǎn)單消費(fèi)函數(shù) ttt uu 1 應(yīng)定義方程為： cs c gdp cs

10、(-1) ar(1)。例子：工作文件15_1eq_cs_ar1 tttt ucscgdpcccs 1321 cst = -22.35 + 0.0924 * gdpt + 0.874 * cst-1 ut = 0.2789 * ut-1 14 估計(jì)高階ar模型稍稍復(fù)雜些，為估計(jì)ar(k)，應(yīng)輸入模型的定義和所包 括的各階ar值。如果想估計(jì)一個(gè)有1-5階自回歸的模型 tttt uuu 5511 應(yīng)輸入： cs c gdp cs(-1) ar(1) ar(2) ar(3) ar(4) ar(5) 例子：工作文件15_1eq_cs_ar5 可以輸入在模型中想包括的各個(gè)自回歸，eviews在消除序列相關(guān)

11、時(shí)給予 很大靈活性。例如，如果有季度數(shù)據(jù)而且想用一個(gè)單項(xiàng)來(lái)說(shuō)明季節(jié)自回歸， 可以輸入：cs c gdp cs(-1) ar(4)。 tttt ucscgdpcccs 1321 15 15.3.3 存在序列相關(guān)的非線性模型存在序列相關(guān)的非線性模型 eviews可以估計(jì)帶有ar誤差項(xiàng)的非線性回歸模型。例如：估計(jì)如下的帶 有附加ar(2)誤差的非線性方程 t c tt ugdpccs 2 1 tttt ucucu 2413 使用eviews表達(dá)式定義模型，在后面的方括號(hào)內(nèi)描述ar修正項(xiàng)，對(duì)每 一階ar滯后項(xiàng)都應(yīng)包括一個(gè)系數(shù)，每項(xiàng)之間用逗號(hào)隔開(kāi)。 cs=c(1)+gdpc(2)+ar(1)=c(3),

12、 ar(2)=c(4) eviews通過(guò)差分來(lái)轉(zhuǎn)換這種非線性模型且使用gauss-newton迭代法來(lái)估 計(jì)轉(zhuǎn)換后的非線性模型。 16 15.3.4 存在序列相關(guān)的二階段回歸模型存在序列相關(guān)的二階段回歸模型 通過(guò)把二階段最小二乘法或二階段非線性最小二乘法和ar項(xiàng)結(jié)合起來(lái)， 對(duì)于在回歸因子和擾動(dòng)項(xiàng)存在相關(guān)性的情況和殘差存在序列相關(guān)一樣估計(jì)模 型。 如果原始回歸模型是線性的，eviews使用marquardt算法來(lái)估計(jì)變形后模 型的參數(shù)。 如果原始回歸模型是非線性的，eviews使用gauss-newton算法來(lái)估計(jì)ar 修正后的模型。 對(duì)于存在序列相關(guān)的情況，可以通過(guò)向方程添加ar項(xiàng)來(lái)調(diào)整tsl

13、s。 eviews會(huì)自動(dòng)將模型轉(zhuǎn)化為非線性最小二乘問(wèn)題，并用工具變量估計(jì)模型。 估計(jì)對(duì)話框中的options 鈕用來(lái)改變非線性工具變量過(guò)程的迭代次數(shù)限制和收 斂標(biāo)準(zhǔn)。 17 例子：例子：15_1eq_cs_tsls_ar 假設(shè)用二階段最小二乘估計(jì)消費(fèi)函數(shù)，考慮存在一階序列相關(guān)。 二階段最小二乘變量列表為： cs c gdp ar (1) 工具變量列表為： c gov log(m1) cs(-1) gdp(-1) 注意因變量的滯后（cs(-1)）和內(nèi)生變量的滯后（gdp(-1)）都包括 在工具變量表中。 類(lèi)似地，考慮消費(fèi)函數(shù)， cs c cs(-1) gdp ar(1) 有效的工具變量表為： c

14、 gov log(m1) cs(-1) cs(-2) gdp(-1) 18 當(dāng)估計(jì)某個(gè)含有ar項(xiàng)的模型時(shí)，在解釋結(jié)果時(shí)一定要小心。在用通常的方 法解釋估計(jì)系數(shù)，系數(shù)標(biāo)準(zhǔn)誤差和t-統(tǒng)計(jì)量時(shí)，涉及殘差的結(jié)果會(huì)不同于ols的 估計(jì)結(jié)果。 要理解這些差別，記住一個(gè)含有ar項(xiàng)的模型有兩種殘差： 第一種是 bxyu ttt 通過(guò)原始變量以及估計(jì)參數(shù) 算出。在用同期信息對(duì)y t值進(jìn)行預(yù)測(cè)時(shí)， 這些殘差是可以觀測(cè)出的誤差，但要忽略滯后殘差中包含的信息。 19 第二種殘差是估計(jì)的 。如名所示，這種殘差代表預(yù) 測(cè)誤差。如果使用前期數(shù)據(jù)殘差和當(dāng)前信息作預(yù)測(cè)，實(shí)際上，通過(guò)利用滯后 殘差的預(yù)測(cè)能力，改善了無(wú)條件預(yù)測(cè)和殘

15、差。 對(duì)于含有ar項(xiàng)的模型，基于殘差的回歸統(tǒng)計(jì)量，如r2 (回歸標(biāo)準(zhǔn)誤差)和 d-w值都是以一期向前預(yù)測(cè)誤差為基礎(chǔ)的。含有ar項(xiàng)的模型獨(dú)有的統(tǒng)計(jì)量是 估計(jì)的ar系數(shù) 。對(duì)于簡(jiǎn)單ar(1)模型， 是無(wú)條件殘差的序列相關(guān)系數(shù)。 對(duì)于平穩(wěn)ar(1)模型， 在-1（極端負(fù)序列相關(guān)）和+1（極端正序列相關(guān)）之 間。 eviews在回歸輸出的底部給出這些根：inverted ar roots。如果存在虛 根，根的模應(yīng)該小于1。 i 20 課本上經(jīng)常描述估計(jì)ar模型的技術(shù)。探討最多的方法，如cochrane- orcutt (科克蘭內(nèi)-奧克特) 、prais-winsten、hatanaka以及hildre

16、th-lu程序 都是使用標(biāo)準(zhǔn)線性回歸進(jìn)行估計(jì)的多步方法。當(dāng)使用滯后因變量作為回 歸自變量或使用高階ar項(xiàng)定義模型時(shí)所有這些方法都有嚴(yán)重的缺點(diǎn)。見(jiàn) davidson& mackinnon (1994, pp.329-341), greene(1997, p.600-607)。 eviews估計(jì)ar模型采用非線性回歸方法。這種方法的優(yōu)點(diǎn)在于：易 被理解，應(yīng)用廣泛，易被擴(kuò)展為非線性定義的模型。注意：非線性最小 二乘估計(jì)漸進(jìn)等于極大似然估計(jì)且漸進(jìn)有效。 21 為估計(jì)ar(1)模型，eviews通過(guò)將線性模型 ttt ttt uu uxy 1 轉(zhuǎn)換成非線性模型。將第二個(gè)方程代入第一個(gè)方程，整理 tttt

17、t xxyy )( 11 參數(shù)通過(guò)應(yīng)用marquarat非線性最小二乘法估計(jì)。 對(duì)于非線性定義，eviews將非線性模型 ttt ttt uu uxfy 1 ),( 轉(zhuǎn)換成： ttttt xfxfyy ),(),( 11 使用gauss-newton算法來(lái)估計(jì)參數(shù)。 22 高階ar定義情況也類(lèi)似。例如，在方程中運(yùn)用非線性最小二乘估計(jì)的 非線性ar(3)如下： ),(),()( 11332211 ttttyt xfxfyyyy ttt xfxf ),(),( 3322 23 arma估計(jì)理論都是基于平穩(wěn)時(shí)間序列。如果一個(gè)序列的均值和自協(xié)方 差不依賴于時(shí)間，就說(shuō)它是平穩(wěn)的。 非平穩(wěn)序列的典型例子是

18、隨機(jī)游動(dòng): ttt yy 1 是平穩(wěn)隨機(jī)擾動(dòng)項(xiàng)。序列y的方差隨時(shí)間增長(zhǎng)，若設(shè) ，則 的方差 是 。但是隨機(jī)游動(dòng)是差分平穩(wěn)序列，因?yàn)閥一階差分后平穩(wěn)： t tttt ylyy )1 ( 1 差分平穩(wěn)序列稱為單整，記為i(d)，d為單整階數(shù)。單整階數(shù)是使序列平穩(wěn) 而差分的階數(shù)。對(duì)于上面的隨機(jī)游動(dòng)，有一個(gè)單位根，所以是i(1)，同樣，平 穩(wěn)序列是i(0)。檢查序列平穩(wěn)性的標(biāo)準(zhǔn)方法是單位根檢驗(yàn)。 0 0 y t y 2 t 24 arima（autoregressive integrated moving average, 自回歸單整動(dòng)平均）模 型是簡(jiǎn)單的ar模型的一般化，使用三種工具來(lái)為擾動(dòng)項(xiàng)的序列

19、相關(guān)建模。 上述的ar(1)模型只運(yùn)用了一階ar項(xiàng)，一般地，可以使用高 階ar項(xiàng)。每一個(gè)ar項(xiàng)對(duì)應(yīng)于在無(wú)條件殘差預(yù)測(cè)方法中的滯后值。p 階的自回歸 模型ar(p)有下面的形式： tptpt t t uuuu 22 1 1 每一單整階數(shù)對(duì)應(yīng)于對(duì)序列進(jìn)行差分。一階單整意味著對(duì)原 始序列進(jìn)行一次差分，二階單整對(duì)應(yīng)于進(jìn)行兩次差分，依此類(lèi)推。 25 動(dòng)平均預(yù)測(cè)模型使用預(yù)測(cè)誤差的滯后值來(lái)改善當(dāng)前預(yù)測(cè)。一 階動(dòng)平均利用前期預(yù)測(cè)誤差，二階動(dòng)平均利用前兩期預(yù)測(cè)誤差，以此類(lèi)推。 ma(q)有如下形式： qtqttt u 11 請(qǐng)注意，一些教科書(shū)和軟件包的系數(shù)使用的是與常規(guī)相反的符號(hào)，因此 ma系數(shù)也可能是相反的。

20、 自回歸和動(dòng)平均可以結(jié)合在一起形成arma(p , q)定義為： tptpttt uuuu 2211 qtqtt 2211 盡管計(jì)量經(jīng)濟(jì)學(xué)家常應(yīng)用arma模型于回歸模型殘差分析，它也可直接應(yīng) 用于序列。后種方法提供了一種單變量模型，將條件序列均值設(shè)定為一個(gè)常數(shù)， 將殘差估計(jì)成均值的序列差分。 26 在arima預(yù)測(cè)中，用上述的三種程序塊的結(jié)合來(lái)建立一個(gè)完整的預(yù) 測(cè)模型。建立殘差序列arima模型的是看其自相關(guān)性，可以使用 相關(guān)圖來(lái)達(dá)到這一目的。arima模型建模過(guò)程被稱為識(shí)別（勿與有關(guān)聯(lián) 立方程的著作中使用的同一詞混淆）。殘差當(dāng)前值與過(guò)去值之間的相關(guān)性 為選擇arima形式提供了導(dǎo)向。自相關(guān)很

21、容易解釋每一值是序列當(dāng) 前值和滯后一定區(qū)間的值的相關(guān)系數(shù)。偏自相關(guān)有些復(fù)雜，它們是考慮了 序列所有值滯后后的預(yù)測(cè)能力后，計(jì)算當(dāng)前和滯后序列的相關(guān)性。例如滯 后6階偏自相關(guān)是計(jì)算當(dāng) 已在預(yù)測(cè)模型中時(shí)， 的預(yù)測(cè)能力。 實(shí)際上， 如果懷疑左側(cè)的因變量和其他的預(yù)測(cè)值之間存在一個(gè)分布滯后關(guān)系， 那么可以在執(zhí)行估計(jì)前看它們的交叉相關(guān)。 51tt uu 6t u 6t u t u 27 是決定使用何種arima模型。如果自相關(guān)函數(shù)以幾何速率衰減， 偏自相關(guān)函數(shù)一階滯后后為零，即為一階自回歸模型。同樣，如果自相關(guān)一 階滯后后為零，而偏自相關(guān)以幾何速率衰減，即為一階動(dòng)平均。如果自相關(guān) 有季節(jié)特征，這說(shuō)明存在季節(jié)

22、arma結(jié)構(gòu)。例如：選擇eviews example filesdatahshs序列工作欄中view/correlogram便可以檢驗(yàn)基本住房序列相 關(guān)圖。 28 上圖中存在季節(jié)頻度的波動(dòng)循環(huán)，建議用季節(jié)arma模型擬合hs序列。 arima分析的目的是控制殘差過(guò)程的一種過(guò)度節(jié)約的表示法。也可 僅用ar或ma來(lái)擬和殘差的特性。使用aic 準(zhǔn)則和schwarz準(zhǔn)則來(lái)選擇滯 后階數(shù)。 建立合適的arima模型后，應(yīng)當(dāng)確認(rèn)模型沒(méi)有殘差自相關(guān)。檢查擾 動(dòng)項(xiàng)的自相關(guān)和偏自相關(guān)，還要考慮是否有重要的預(yù)測(cè)能力被忽略。 eviews提供了估計(jì)之后的診斷檢查方案。 29 eviews估計(jì)的是考慮到具有右側(cè)解釋變

23、量的arima(p, d, q)設(shè)定形式，盡 管這樣的模型有時(shí)叫做arimax模型： arima模型是arimax模型的特例，但是我們把這一類(lèi)的模型都叫做arima 模型。 為建立arima模型，首先需要做以下兩點(diǎn)： 1確定單整階數(shù)，差分因變量序列； 2描述結(jié)構(gòu)回歸模型（因變量和解釋變量），用前面介紹的方式加入ar 或ma項(xiàng)。 qtqttptpttt tktkttt uuuu uxxxy 112211 2211 30 d算子被用來(lái)定義序列差分。定義一階差分，僅把序列名寫(xiě)入d后的括 號(hào)。例如，d(gdp)定義gdp的一階差分，或gdp-gdp(-1)。更復(fù)雜的差分形 式可以使用兩個(gè)參數(shù) n，s。

24、d(x, n)定義序列x的n 階差分 xlnxd n )1 (),( l是滯后算子。例如：d(gdp, 2)定義了gdp的2階差分： d(gdp, 2)=gdp - 2*gdp(-1)+gdp(-2) d(x, n, s)定義序列x的n階普通差分，帶有滯后s階的季節(jié)差分： xllsnxd sn )1 ()1 (),( 例如：d(gdp, 0, 4)定義帶有滯后4階季節(jié)差分的零階普通差分，即gdp- gdp(-4)。 如果需要對(duì)數(shù)形式，可以使用dlog算子，它以對(duì)數(shù)值返回差分。例如： dlog(gdp)定義log(gdp)的一階差分，即 log(gdp)-log(gdp(-1)。 31 可以直接

25、在估計(jì)定義式中包含差分算子d。例如：gdpi(1)，即gdp是 一階單整序列。對(duì)gdp估計(jì)arima(1,1,1)模型，可以輸入列表(15_1eq_dy)： d(gdp) c ar(1) ma(1) 使用因變量差分因子d(gdp)定義模型， eviews將提供水平變量gdp的預(yù)測(cè)值。 模型中ar和ma部分應(yīng)使用關(guān)鍵詞ar和ma定義。在上面ar定義中，我們 已見(jiàn)過(guò)這種方法的例子。這對(duì)ma也同樣適用。 32 例如，估計(jì)一個(gè)2階自回歸和1階動(dòng)平均過(guò)程arma(2,1)，應(yīng)將ar(1), ma(1), ar(2)和其它解釋變量一起包含在回歸因子列表中： y c gov ar(1) ar(2) ma(1

26、) 不必連續(xù)使用ar和ma項(xiàng)。例如想用4階季節(jié)自回歸模型來(lái)擬合季節(jié) 變化，可以僅使用ar(4)： y c gov ar(4 ) 也可僅用ma項(xiàng)來(lái)定義純動(dòng)平均模型。如可以表示出殘差的ma(2)模型。 y c gov ma(1) ma(2) 在這種情況下，解釋變量將僅包含一個(gè)c加上ar，ma項(xiàng)，例如： y c ar(1) ar(2) ma(1) ma(2) 這是標(biāo)準(zhǔn)的box-jenkins arma(2, 2)模型。 33 對(duì)于帶有季節(jié)因素的季度數(shù)據(jù)，box and jenkins(1976) 建議使用季節(jié)自 回歸sar和季節(jié)動(dòng)平均sma。sar(p)定義為帶有p階滯后的季節(jié)自回歸項(xiàng)。 估計(jì)中使用

27、的滯后多項(xiàng)式是ar項(xiàng)和sar項(xiàng)定義的結(jié)合。 與此類(lèi)似，sma(q)定義為帶有q階滯后的季節(jié)動(dòng)平均。估計(jì)中使用的滯 后多項(xiàng)式是ma項(xiàng)和sma項(xiàng)定義的結(jié)合。存在sar項(xiàng)則允許建立一個(gè)滯后多項(xiàng) 式。 例如：沒(méi)有季節(jié)項(xiàng)的2階ar過(guò)程 tttt uuu 2211 ntt n xxll , 用滯后算子 ，則上式可表示為： tt ull)1 ( 2 21 可以通過(guò)回歸自變量的ar(1)，ar(2)項(xiàng)來(lái)估計(jì)這個(gè)過(guò)程。 34 對(duì)于季度數(shù)據(jù)，可以加入sar(4)來(lái)表示季節(jié)因素，定義方程： y c x ar(1) ar(2) sar(4) 估計(jì)誤差結(jié)構(gòu)為： tt ulll)1)(1 ( 42 21 等價(jià)于 tttt

28、ttt uuuuuu 625142211 參數(shù) 和季節(jié)因素相聯(lián)系。注意：這是對(duì)系數(shù)有非線性約束的ar(6)模型。 在另一個(gè)例子中，無(wú)季節(jié)性的二階ma過(guò)程如下 , 2211 tttt u 可以通過(guò)包含ma(1)和ma(2)來(lái)估計(jì)二階ma過(guò)程。 35 對(duì)季度數(shù)據(jù)，可以添加sma(4)考慮季節(jié)性。例如定義方程： y c x ma(1 ) ma(2) sma(4) 估計(jì)模型為： tt lllu)1)(1 ( 42 21 等價(jià)于： 625142211 ttttttt u 參數(shù) 和季節(jié)因素相聯(lián)系。這是對(duì)系數(shù)有非線性約束的ma(6)模型。還可以 在方程說(shuō)明中同時(shí)包括sar，sma項(xiàng)。 ttt uxy 21

29、36 含有ar或ma項(xiàng)的模型的估計(jì)輸出和ols模型一樣，只是在底部增加了一 個(gè)ar，ma多項(xiàng)式的根的倒數(shù)。如果我們利用滯后多項(xiàng)式 和 寫(xiě)一般 的arma模型： )(l)(l tt lul)()( 輸出表中報(bào)告的結(jié)果相當(dāng)于下列多項(xiàng)式 0)( 1 x和0)( 1 x 如果 有絕對(duì)值大于1的實(shí)根或一對(duì)復(fù)根的逆在單位圓外（即模大于1）， 這意味著自回歸過(guò)程是發(fā)散的。如果 的根的倒數(shù)在單位圓外，說(shuō)明ma過(guò)程是 不可逆的，應(yīng)使用不同的初值重新估計(jì)模型，直到得到滿足可逆性的動(dòng)平均。 如果估計(jì)的ma模型的根的模接近于1，有可能是對(duì)數(shù)據(jù)差分過(guò)多，這就很難估 計(jì)和預(yù)測(cè) 。如果可能的話，應(yīng)減少差分階數(shù)重新估計(jì)。 的

30、根。這些根（可能是虛根）的模應(yīng)小于1，如果不滿足這個(gè)條件，輸出表中將 顯示警告信息。 37 38 這個(gè)arma估計(jì)輸出例子的結(jié)果對(duì)應(yīng)于如下定義： 或等同于： 注意：ma項(xiàng)的符號(hào)和教科書(shū)中的符號(hào)可能相反。倒根的模接近于1，這對(duì) 于許多宏觀經(jīng)濟(jì)序列是很典型的。 541 57. 091. 063. 0 ttt tt tt llull ur )91. 01)(63. 01 ()93. 01)(86. 01 ( 21. 7 44 ttttt rrrr 541 80. 093. 086. 007. 0 39 如前所述，帶有ar或ma的模型用非線性最小二乘法估計(jì)。非線性估計(jì) 方法對(duì)所有系數(shù)估計(jì)都要求初值。e

31、views自行確定初值。 有時(shí)當(dāng)?shù)畲笾颠_(dá)到時(shí)，方程終止迭代，盡管還未達(dá)到收斂。從前一 步初值重新開(kāi)始方程，使方程從中止處開(kāi)始而不是從開(kāi)始處開(kāi)始。也可以試 試不同的初值來(lái)保證估計(jì)是全部而不是局部平方誤差最小，可以通過(guò)提供初 值加速估計(jì)過(guò)程。 為控制arma估計(jì)初值，在方程定義對(duì)話框單擊options。在eviews提供 的選項(xiàng)中，有幾項(xiàng)設(shè)置初值的選擇。 eviews缺省方法是ols/tsls，這種方法先進(jìn)行沒(méi)有arma項(xiàng)的預(yù)備估計(jì)， 再?gòu)倪@些值開(kāi)始非線性估計(jì)。另一選擇是使用ols或tsls系數(shù)的一部分作為 初值?？梢赃x擇0.3，0.5，0.8或者可以將所有初值設(shè)為零。 用戶確定初值選項(xiàng)是us

32、er supplied。在這個(gè)選項(xiàng)下，eviews使用c系數(shù)向 量中的值。為設(shè)置初值，雙擊圖標(biāo)，為c系數(shù)向量開(kāi)一窗口，進(jìn)行編輯。 40 為適當(dāng)?shù)卦O(shè)置初值，需對(duì)eviews如何為arma設(shè)置系數(shù)多些了解。eviews 使用c系數(shù)向量。它按下列規(guī)則為變量安排系數(shù)： 1. 變量系數(shù)，以輸入為序。 2. 定義的ar項(xiàng)，以輸入為序。 3sar，ma，sma系數(shù)（以輸入為序） 這樣，下面兩種定義將有同樣規(guī)格的系數(shù) y c x ma(2) ma(1) sma(4) ar(1) y sma(4 ) c ar(1) ma(2) x ma(1) 也可使用程序指令安排c向量值 param c(1) 50 c(2 )

33、 0.8 c(3) 0.2 c(4) 0.6 c(5) 0.1 c(6) 0.5 初值：常數(shù)是50， x系數(shù)的初值是0.8， ar(1)、ma(2)、ma(1)、sma(4) 系數(shù) 的初值分別是0.2 ， 0.6，0.1，0.5。 估計(jì)后，可在方程表達(dá)式representation選項(xiàng)見(jiàn)到系數(shù)安排。也可以從估計(jì) 方程中填寫(xiě)c向量，選擇pros/update/ coefs from equations。 41 對(duì)于arma模型估計(jì)，要考慮一些問(wèn)題。首先，ma模型很難估計(jì)。 特別的，應(yīng)避免高階ma，除非模型非常需要，因?yàn)樗鼈兛赡芤鸸烙?jì)困 難。例如，相關(guān)圖上滯后57有一個(gè)大波峰并不要求必須在模型中

34、包括 ma(57)，除非知道每57期都有特別事情發(fā)生。相關(guān)圖中的突起很可能是序 列中的一個(gè)或多個(gè)奇異值的結(jié)果。模型中含有許多ma項(xiàng)，將會(huì)喪失自由 度，并且可能犧牲估計(jì)的穩(wěn)定性和可靠性。 如果ma過(guò)程的根的模接近于1，可能會(huì)遇到估計(jì)困難。eviews會(huì)報(bào)告 在迭代到最大次數(shù)時(shí)，不能收斂或不能提高平方和，這說(shuō)明可能對(duì)數(shù)據(jù)差 分過(guò)多，應(yīng)檢查序列相關(guān)圖來(lái)看是否可以減少差分階數(shù)來(lái)重新估計(jì)。 42 對(duì)arima進(jìn)行二階段最小二乘法或工具變量法沒(méi)有什么特殊困難。 eviews將估計(jì)帶有自回歸項(xiàng)的非線性最小二乘模型。 eviews目前不估計(jì)有ma誤差的非線性模型。然而，可以運(yùn)用狀態(tài)空 間方法來(lái)定義估計(jì)這些模型

35、。 eviews不會(huì)自動(dòng)估計(jì)帶有arma誤差項(xiàng)的加權(quán)模型如果對(duì)一加權(quán) 模型加入ar項(xiàng)，加權(quán)序列會(huì)被忽略。 43 如果arma模型定義正確，模型殘差將為白噪聲。這意味著殘差中應(yīng)不存 在序列相關(guān)。d-w統(tǒng)計(jì)量是當(dāng)方程右邊沒(méi)有滯后變量時(shí)對(duì)一階序列相關(guān)的檢驗(yàn)。 如上所述，對(duì)殘差中序列相關(guān)更多的檢驗(yàn)可以如：view/residual tests/ corre- logram-q-statistic和view/residual tests/serial correlation lm test。 44 在經(jīng)濟(jì)分析中人們發(fā)現(xiàn)，一些經(jīng)濟(jì)變量，它們的數(shù)值是由自身的滯后量 或者其他變量的滯后量所決定的，表現(xiàn)在計(jì)量經(jīng)

36、濟(jì)模型中，解釋變量中經(jīng)常 包含某些滯后變量。以投資函數(shù)為例，分析中國(guó)的投資問(wèn)題發(fā)現(xiàn)，當(dāng)年的投 資額除了取決于當(dāng)年的收入（即國(guó)內(nèi)生產(chǎn)總值）外，由于投資的連續(xù)性，它 還受到前1 個(gè)、2個(gè)、3個(gè)時(shí)期投資額的影響。已經(jīng)開(kāi)工的項(xiàng)目總是要繼續(xù)下 去的，而每個(gè)時(shí)期的投資額又取決于每個(gè)時(shí)期的收入，所以可以建立如下關(guān) 于投資的計(jì)量經(jīng)濟(jì)方程 其中i 表示投資額，y 表示國(guó)內(nèi)生產(chǎn)總值。 ttttt yyyi 22110 45 對(duì)于有限滯后長(zhǎng)度的情形，分布滯后模型的一般形式如下 tktktttt xxxy 110 系數(shù) 描述x對(duì)y作用的滯后。在模型中解釋變量與隨機(jī)誤差項(xiàng)不相關(guān)的 情況下，可以直接使用ols估計(jì)參數(shù)。但

37、是，一個(gè)顯然的問(wèn)題是解釋變量之 間，即x的當(dāng)前和滯后值之間具有高度共線性，而共線性問(wèn)題的一個(gè)直接后 果是參數(shù)估計(jì)量失去意義，不能揭示x的各個(gè)滯后量對(duì)因變量的影響，所以 必須尋求另外的估計(jì)方法。 可以使用多項(xiàng)式分布滯后（polynomial distributed lags , pdls）來(lái)減少 要估計(jì)的參數(shù)個(gè)數(shù)，以此來(lái)平滑滯后系數(shù)。平滑就是要求系數(shù)服從一個(gè)相對(duì) 低階的多項(xiàng)式。p 階pdls模型限制 系數(shù)服從如下形式的p階多項(xiàng)式 p pj cjcjcj)()()( 1 2 321 j = 0 , 1 , 2 , , k (15.2) (15.1) 46 是事先定義常數(shù)： c 是偶數(shù) 是奇數(shù)（ p

38、k pk c 2/ )( 2/ ) 1( pdls有時(shí)被稱為almon分布滯后模型。常數(shù) 僅用來(lái)避免共線性引起的數(shù) 值問(wèn)題，不影響 的估計(jì)。這種定義允許僅使用參數(shù) p 來(lái)估計(jì)一個(gè)x 的 k 階 滯后的模型（如果 p k，將顯示“近似奇異“錯(cuò)誤信息）。 定義一個(gè)pdl模型，eviews用(15.2)式代入到(15.1)式，將產(chǎn)生如下形式方 程 c tppt zzzy 112211 其中 kt p t p t p p kttt kttt xckxcxcz xckxcxcz xxxz )()1 ()( )()1 ( 11 12 11 (15.3) 47 一旦從(15.3)式估計(jì)出 ，利用(15.2)

39、式就可得到 的各系數(shù)。這一過(guò)程 很明了，因?yàn)槭?的 線性變換。定義一個(gè)pdls要有三個(gè)元素：滯后長(zhǎng)度k， 多項(xiàng)式階數(shù)（多項(xiàng)式最高次冪數(shù)）p和附加的約束條件。 一個(gè)近端約束限制x對(duì)y一期超前作用為零： 0)1()1()1( 1 2 3211 p p ccc 一個(gè)遠(yuǎn)端約束限制x對(duì)y的作用在大于定義滯后的數(shù)目衰減： 0)1()1()1( 1 2 3211 p pk ckckck 如果限制滯后算子的近端或遠(yuǎn)端，參數(shù)個(gè)數(shù)將減少一個(gè)來(lái)解釋這種約束。 如果對(duì)近端和遠(yuǎn)端都約束，參數(shù)個(gè)數(shù)將減少二個(gè)。 eviews缺省不加任何約束。 48 通過(guò)pdl項(xiàng)定義一個(gè)多項(xiàng)式分布滯后，信息在隨后的括號(hào)內(nèi)，按下列規(guī)則 用逗號(hào)

40、隔開(kāi)： 1、序列名 2、滯后長(zhǎng)度（序列滯后數(shù)） 3、多項(xiàng)式階數(shù) 4、一個(gè)數(shù)字限制碼來(lái)約束滯后多項(xiàng)式： 1 = 限制滯后近端為零 2 = 限制遠(yuǎn)端 3 = 兩者都限制 如果不限制滯后多項(xiàng)式，可以省略限制碼。方程中可以包含多個(gè)pdl項(xiàng)。 例如： sales c pdl(y , 8 , 3 )是用常數(shù)，解釋變量y的當(dāng)前和8階分布滯后來(lái) 擬合因變量sales，這里解釋變量y的滯后系數(shù)服從沒(méi)有約束的3階多項(xiàng)式。 49 類(lèi)似地， y c pdl(x , 12 , 4 , 2) 用常數(shù)，解釋變量x的當(dāng)前和12階分布滯后擬合因變量y，這里解釋變量x的 系數(shù)服從帶有遠(yuǎn)端約束的4階多項(xiàng)式。 pdl也可用于二階段最

41、小二乘法tsls。如果pdl序列是外生變量，應(yīng)當(dāng)在工 具表中也包括序列的pdl項(xiàng)。為此目的，可以定義pdl(*)作為一個(gè)工具變量，則 所有的pdl變量都將被作為工具變量使用。例如：如果定義tsls方程為 sales c inc pdl(y(-1) , 12 , 4) 使用工具變量 z z(-1) pdl(*) 則y的分布滯后和z，z(-1)都被用作工具變量。 50 三、例子三、例子 投資inv關(guān)于關(guān)于gdp的 分布滯后模型的結(jié)果如下 51 逐個(gè)觀察，inv滯后的系數(shù)統(tǒng)計(jì)上都不顯著。但總體上講回歸具有一個(gè)合 理的r2， (盡管dw統(tǒng)計(jì)量很低)。這是回歸自變量中多重共線的典型現(xiàn)象， 建議擬合一個(gè)多

42、項(xiàng)式分布滯后模型。估計(jì)一個(gè)無(wú)限制的5階多項(xiàng)式滯后模型， 輸入變量列表：log(inv) c pdl(log(gdp), 12, 5)，窗口中顯示的多項(xiàng)式估計(jì)系 數(shù)，pdl01, pdl02, pdl03 , 分別對(duì)應(yīng)方程(15.3)中z1, z 2 , 的系數(shù)。 52 方程（15.1）中的系數(shù) 在表格底部顯示。 j 表格底部的滯后值是分布滯后的估計(jì)系數(shù)值，并且在平穩(wěn)的假設(shè)下有 gdp對(duì)inv的長(zhǎng)期影響的解釋。 53 eviews提供兩種單位根檢驗(yàn)：dickey-fuller(df)、增廣df(augmented df)檢驗(yàn)和phillips - perron（pp）檢驗(yàn)。本節(jié)提供這兩種檢驗(yàn)的一

43、些理論背 景。 為說(shuō)明adf檢驗(yàn)的使用，先考慮一個(gè)ar(1)過(guò)程 ttt yy 1 是參數(shù)， 假設(shè)為白噪聲。如果-1 1，y是平穩(wěn)序列。如果 =1， y是 非平穩(wěn)序列（帶漂移的隨機(jī)游動(dòng)）。如果這一過(guò)程在某一點(diǎn)開(kāi)始， y的方差隨 時(shí)間增長(zhǎng)趨于無(wú)窮。如果 的絕對(duì)值大于1，序列發(fā)散。因此，一個(gè)序列是否平 穩(wěn)，可以檢驗(yàn) 是否嚴(yán)格小于1。df和pp都用單位根作為原假設(shè)， 。 因?yàn)榘l(fā)散序列沒(méi)有經(jīng)濟(jì)學(xué)含義，所以備選假設(shè)為單邊假設(shè) 。 , t 1: 0 h 1: 1 h （15.4） 54 從方程兩邊同時(shí)減去 1t y ttt yy 1 其中 1 所以原假設(shè)和備選假設(shè)可改為 0: 0: 1 0 h h （15

44、.5） 單位根檢驗(yàn)可以看作對(duì) 進(jìn)行t 檢驗(yàn)。對(duì)于（15.5）可使用常規(guī)的ols方法 得到 的估計(jì)量和t統(tǒng)計(jì)量(即系數(shù)的估計(jì)量除以其標(biāo)準(zhǔn)偏差)，但是這時(shí)t統(tǒng)計(jì)量 不再漸近正態(tài)，甚至不對(duì)稱。dickey和fuller(1979)通過(guò)蒙特卡洛方法研究發(fā)現(xiàn)， 檢驗(yàn) 的統(tǒng)計(jì)量的臨界值依賴于回歸的形式(是否引進(jìn)了常數(shù)項(xiàng)和趨勢(shì)項(xiàng))， 并且和樣本長(zhǎng)度n有關(guān)，一般臨界值隨著n的增大而減小。t 統(tǒng)計(jì)量在原假設(shè)下不 具有一般 t 分布。dickey和fuller對(duì)選擇的樣本序列模擬了臨界值，近年來(lái) mackinnon又進(jìn)行了更大規(guī)模的模擬，允許d-f計(jì)算的臨界值適用于任何樣本和 任何數(shù)目的右邊變量。eviews為單

45、位根檢驗(yàn)提供了這些臨界值。 0 （15.6） 55 上面描述的簡(jiǎn)單單位根檢驗(yàn)稱為df檢驗(yàn)，該檢驗(yàn)只有當(dāng)序列為ar(1)時(shí)才 有效。如果序列存在高階滯后相關(guān)，這就違背了擾動(dòng)項(xiàng)是白噪聲的假設(shè)。adf 和pp檢驗(yàn)使用不同的方法來(lái)控制高階序列相關(guān)。 adf方法通過(guò)在回歸方程右邊加入因變量y的滯后差分項(xiàng)來(lái)控制高階相關(guān) ptptttt yyyy 1111 擴(kuò)展定義將檢驗(yàn) 0: 0: 1 0 h h eviews將df，adf檢驗(yàn)都看成為adf檢驗(yàn)。adf檢驗(yàn)考慮如下三種回歸 形式： t p i ititt yyy 1 1 t p i ititt ytayay 1 210 t p i ititt yyay

46、1 10 56 通過(guò)在模型中增加的滯后項(xiàng)來(lái)消除殘差的序列相關(guān)性 。eviews在檢驗(yàn)回 歸時(shí),詢問(wèn)是否包含有其它外生變量，即上面三種形式：包含截距，即常數(shù)項(xiàng)， 包含截距和線性趨勢(shì)項(xiàng)，或兩者都不包含。 如果序列含有趨勢(shì)（確定的或隨機(jī)的），序列回歸中應(yīng)既有常數(shù)又有趨勢(shì)。 如果序列沒(méi)有表現(xiàn)任何趨勢(shì)且有非零均值，回歸中應(yīng)僅有常數(shù)。如果序列在零 均值波動(dòng)，檢驗(yàn)回歸中應(yīng)既不含有常數(shù)又不含有趨勢(shì)。 phillips和perron（1988）提出一種非參數(shù)方法來(lái)控制序列中高階序列相關(guān)。 對(duì)ar(1)的pp檢驗(yàn)為： ttt yy 1 57 adf檢驗(yàn)通過(guò)在方程右邊添加滯后差分項(xiàng)來(lái)修正高階序列相關(guān)。pp檢驗(yàn)參 數(shù)的t 統(tǒng)計(jì)量來(lái)修正ar(1)的序列相