【正文】 Economic Statistics[J], ， 1992， 467470. 。 15 參考文獻(xiàn) [1]張曉峒：《計量經(jīng)濟(jì)分析》 [M]，經(jīng)濟(jì)科學(xué)出版社， 2020. [2]張曉峒：《 Eviews使用指南與案例》 [M]，機(jī)械工業(yè)出版社， 2020. [3]Dickey， D． A． and W． A． Fuller， Distribution of the Estimators for Autoregressive Time Series With a Unit Root， Journal of the American statistical association[J]， 1979， 427431. [4]Hamilton． J． D， Time Series Analysis[M]， Princeton University Press． 1994. [5]Phillips， P． C B， Time Series Regression with a Unit Root， Econometrica[J]， 1987， 277301. [6]Pierre Perron, The Great Crash, the Oil Price Shock, and the Unit Root Hypothesis， Econometrica[J]， ， 1989， 13611401. [7]Pierre Perron, Testing for a Unit Root in a Time Series with a Changing Mean，Journal of Business amp。 ③ 見 Hamilton(1994)， 第 511 頁。 ① 見 Hamilton(1994)， 第 506 頁。而 在現(xiàn)實中 ， 除差分平穩(wěn)過程、退勢平穩(wěn)過程具有時間趨勢外，分?jǐn)?shù)單整過程、帶結(jié)構(gòu)突變的單位根過程均顯示出一定的時間趨勢。 為退勢平穩(wěn)過程檢驗提供理論支撐和檢驗基礎(chǔ)。 六、結(jié)語 本文給出 退勢平穩(wěn)過程檢驗 的 主要 步驟及 相關(guān) 檢驗 統(tǒng)計量的主要性質(zhì)。 很明顯， 檢驗式 (1)、 (2)、 (3)對應(yīng)的 各非參數(shù)檢驗統(tǒng)計量 具有相同的漸近 分布。 根據(jù)相同的分析方法，可以分別考慮運(yùn)用檢驗式 (2)、 (3)進(jìn)行退勢平穩(wěn)過程檢驗時檢驗統(tǒng)計量 的 性質(zhì)。 因而有 ? ? ??????? 10222 2)1(j j?????， 其中： 2? 、 )1(? 定義如前，j? 是檢驗式 (1)隨機(jī)擾動項的 j 階自協(xié)方差② 。其中 tu? 為檢驗式 (1)的回歸殘差。 14 212101010210210101010102? }))((4)()(12))((12)({)()1()]1(21)([])(21)([12)1)1((211 drrWdrrrWdrrWdrrrWdrrWdrrWWWdrrWdrrWdrrrWWZ P??????????????????? ??? 在實際檢驗時，要計算1??Z，除已知 2Ts 、1??t、1????外，還必須知道總體參數(shù) 2? 及 0? 的 估計值，可以通過以下方法給出它們的一致估計量。故對于檢驗式 (1)，可以建立 形式相同的 非參數(shù)檢驗統(tǒng)計量： TT sTtsZ 111 ?02?2122? ?}{21}{ ??? ?? ??? ?????，其中 2Ts 、1??t、1????分別表示檢驗式 (1)回歸結(jié)果中回歸方差、檢驗統(tǒng)計量1??t及估計量 1?? 的標(biāo)準(zhǔn)差。 2Ts 表示模型(5)的回歸方差、*?Tt?表示模型 (5)中估計量 *?T? 對應(yīng)的 t 值、 T 表示樣本容量、*??T??表示模型 (5)中估計量 *?T? 的標(biāo)準(zhǔn)差。 將殘差 te 代入檢驗式 (1),與前面的分析相同， 可得 模型 ： ttt yty ???? ???? ?1 (4) 其中： 10 ??)1( ????? ??? 、 )1(?1 ??? ?? 、 1??? 。 同樣， 先對 ty 運(yùn)用模型tt ty ??? ??? 10 進(jìn)行 退勢， 得殘差序列 tye tt 10 ?? ?? ??? ，再運(yùn)用不 同 檢驗式對 te 進(jìn)行單位根檢驗。 11 定義： )1(0 ????? ?? ???j j， ,. .. .)2,1,0(,),(02 ??? ? ?? ?? juuE s jssjttj ????。給定樣本容量為 5000， 模擬 30000 次，根據(jù)模擬的漸近 分布確定各分位點的 漸近 臨界值，結(jié)果如表所示。 ③ 根據(jù)推導(dǎo)的漸近 分布可如下模擬漸近臨界值：給定ty 的數(shù)據(jù)生成過程 ttt uyy ?? ?1 、 )1,0(...~ Ndiiut 、 00?y 。 ② 3DF指單位根檢驗式 ttt ytay ??? ???? ?1，對應(yīng)原假設(shè) 1:0 ??H 時 ，檢驗統(tǒng)計量 ??t 的漸近臨界值。 其中： ???????0j jj ?， ),0(...~ 2?? Ndiit ， ??)( 4tE? 。 為此，需要對以上分析框架 做相應(yīng)擴(kuò)展： 放松 檢驗式 隨機(jī)擾動項無自相關(guān)的假設(shè) 。 表 7 不同方法 計算 檢驗統(tǒng)計量 漸近臨界值的比較 計算方法 檢驗式 3DF ② 極限分布 模擬 ③ (1) (2) (3) 檢 驗 水 平 五、檢驗統(tǒng)計量的擴(kuò)展 以上 分析是在 數(shù)據(jù)生成過程： ttt uyy ??? ?10? ， ),(~0 U? 、 00?y 、 ~tu )1,0(... Ndii 下得到的結(jié)果 ， 它 隱含的 假設(shè) 是檢驗式的 隨機(jī)擾動項是無自相關(guān)的 。 所以當(dāng)運(yùn)用 DF單位根檢驗臨界值進(jìn)行退勢平穩(wěn)過程檢驗時，會增大犯第二類錯誤(取偽 )的概率。 (Ⅳ ) 與 DF單位根檢驗的臨界值表 比較， 在 給定顯著性水平、 樣本容量 條件 下 ， 對應(yīng)檢驗式 的 經(jīng)驗 臨界值均小于單位根檢驗臨界值。 (Ⅲ ) 給定顯著性水 平，對于特定的樣本容量而言，檢驗式 (1)的檢驗統(tǒng)計量 經(jīng)驗 臨界值最小，檢驗式(3)的檢驗統(tǒng)計量 經(jīng)驗 臨界值最大。 而對于顯著性水平 、 、 、 ， 隨著樣本容量增大，其 經(jīng)驗臨界值逐漸變小。 有限樣本臨界值是根據(jù) 表 5 中 各 響應(yīng) 面函數(shù) 得 到 。 9 5432100 1 0 0 2 0 0 3 0 0 4 0 0 5 0 0 6 0 0TC V 0 .0 1C V 0 .0 2 5C V 0 .0 5C V 0 .1C V 0 .9C V 0 .9 5C V 0 .9 7 5C V 0 .9 9C V T 543210.0 0 .0 1 .0 2 .0 3 .0 4 .0 5 .0 61 / TC V 0 .0 1C V 0 .0 2 5C V 0 .0 5C V 0 .1C V 0 .9C V 0 .9 5C V 0 .9 7 5C V 0 .9 9C V 1 /T 圖 5 T 與 統(tǒng)計量3??t分位數(shù)的散點圖 圖 6 T/1 與 統(tǒng)計量3??t分位數(shù)的散點圖 表 5 檢驗統(tǒng)計量3??t不同分位點 響應(yīng)面函數(shù)估計 結(jié)果 檢驗統(tǒng)計量 檢驗水平 ?3? 31? 2R se 3??t 注意： 檢驗統(tǒng)計量 臨界值計算公式： 313133 ???. wTVC ??? ?? ?? ，其中 T 為樣本容量。 圖 5是 檢驗統(tǒng)計量3??t臨界值關(guān)于樣本容量 T 的散點圖，圖 6是 檢驗統(tǒng)計量3??t臨界值關(guān)于樣本容量倒數(shù) ( T/1 )的散點圖。 有限樣本臨界值是根據(jù) 表 3 中 各 響應(yīng)面函數(shù) 得 到 。 5432100 10 0 20 0 30 0 40 0 50 0 60 0TC V 0 .0 1C V 0 .0 2 5C V 0 .0 5C V 0 .1C V 0 .9C V 0 .9 5C V 0 .9 7 5C V 0 .9 9C V T 543210.00 .01 .02 .03 .04 .05 .061 /TC V 0 .0 1C V 0 .0 2 5C V 0 .0 5C V 0 .1C V 0 .9C V 0 .9 5C V 0 .9 7 5C V 0 .9 9C V 1 /T 圖 3 T 與 統(tǒng)計量2??t分位數(shù)的散點圖 圖 4 T/1 與 統(tǒng)計量2??t分位數(shù)的散點圖 表 3 檢驗統(tǒng)計量2??t不同分位點響應(yīng)面函數(shù)估計結(jié)果 8 檢驗統(tǒng)計量 檢驗水平 ?2? 21? 2R se 2??t 注意： 檢驗統(tǒng)計量 臨界值計算公式： 212122 ???. wTVC ??? ?? ?? ，其中 T 為樣本容量。 圖 3是 檢驗統(tǒng)計量2??t臨界值關(guān)于樣本容量 T 的散點圖，圖 4是 檢驗統(tǒng)計量2??t臨界值關(guān)于樣本容量倒數(shù) ( T/1 )的散點圖。 有限樣本 臨界值 是根據(jù) 表 1 中 各 響應(yīng) 面函數(shù) 得 到 。 7 注意： 檢驗統(tǒng)計量 臨界值計算公式： 111111 ???. wTVC ??? ?? ?? ，其中 T 為樣本容量。 常用 分位點下檢驗統(tǒng)計量 經(jīng)驗 臨界值 如表 2所示 。 5432100 10 0 20 0 30 0 40 0 50 0 60 0TC V 0 .0 1C V 0 .0 2 5C V 0 .0 5C V 0 .1C V 0 .9C V 0 .9 5C V 0 .9 7 5C V 0 .9 9C V T 543210.0 0 .0 1 .0 2 .0 3 .0 4 .0 5 .0 61 /T0 10 2 5C V 0 .0 5