胡俊娟

(浙江科技學院 理學院,杭州 310023)

?

浙江科技學院學報,第27卷第1期,2015年2月

Journal of Zhejiang University of Science and Technology

Vol.27 No.1, Feb. 2015

基于ESTAR模型的單位根檢驗——Wald統(tǒng)計量的研究

胡俊娟

(浙江科技學院 理學院,杭州 310023)

摘要:已有的基于ESTAR模型的單位根檢驗通常假設α=0,但實踐研究表明α顯著。對α≠0的情況進行單位根檢驗,研究該情形下的Wald檢驗統(tǒng)計量。從理論上推導了Wald檢驗統(tǒng)計量的極限分布,并通過Monte Carlo隨機模擬,結(jié)果表明該統(tǒng)計量有效。

關(guān)鍵詞:單位根檢驗;Wald統(tǒng)計量;ESTAR模型

自Ter?svirta[1]提出STAR(smooth transition autoregressive)模型以來,非線性模型一直受到國內(nèi)外學者的關(guān)注。STAR模型很好地描述了模型的狀態(tài)轉(zhuǎn)化和結(jié)構(gòu)變化,因而被廣泛地應用于預測工業(yè)產(chǎn)出、實際匯率、失業(yè)率等主要宏觀時間序列[2-3]。由于模型擬合前需要進行單位根檢驗,眾多學者對指數(shù)STAR(ESTAR)模型進行了研究。Kapetanios等[4]提出了Dickey-Fuller型t統(tǒng)計量，Kruse[5]用修正的t統(tǒng)計量進行了研究，Hanck[6]對帶均值的情況進行了研究。這些檢驗都在假設α=0的情況下進行,但實踐研究表明α顯著[7]。筆者在此基礎(chǔ)上放寬了條件,在α≠0的情況下進行單位根檢驗的研究。

1基于ESTAR模型的單位根檢驗

ESTAR模型是STAR類模型中最常用的模型之一,首先給出ESTAR模型的具體形式:

yt=βyt-1+γyt-1G(yt-1;θ)+εt,t=1,…,T,

(1)

該過程也可以寫成:

Δyt=αyt-1+γyt-1G(yt-1;θ)+εt,

(2)

式(2)中:α=β-1,Δ—差分算子。

1.1　Dickey-Fuller型t統(tǒng)計量

Kapetanios等[4]在假設α=0的條件下,考慮θ是否為零。即對于方程:

Δyt=γyt-1G(yt-1;θ)+εt

給出了原假設H0:θ=0,備擇假設H1:θ>0來進行單位根檢驗?？紤]到參數(shù)γ未知,常用的做法是用一階泰勒展開來逼近非線性函數(shù)G(yt-1;θ),則檢驗方程就變成:

從而對方程檢驗H0:δ=0。Kapetanios等給出了Dickey-Fuller型t統(tǒng)計量:

(3)

1.2　Wald檢驗統(tǒng)計量

考慮到ESTAR模型(式(1))中α不一定為零,則對ESTAR模型進行一階泰勒展開后得到的輔助線性方程是:

(4)

考慮到當α=δ=0時,模型含有單位根,所以對模型檢驗:H0:α=δ=0。

Wald檢驗是計量經(jīng)濟學三大檢驗方法之一,只要求計算無約束模型,在實證分析中得到廣泛應用。在原假設H0:Rβ=0情況下的Wald統(tǒng)計量[8]為:

(5)

2Wald統(tǒng)計量的極限分布

為了給出統(tǒng)計量W的極限分布,先假設序列yt的六階矩存在。有關(guān)W極限分布的定理如下。

定理1考慮模型(2)在H0:α=δ=0的條件下,統(tǒng)計量W有如下的極限分布:

W?A′B-1A,

即

進一步,有:

并且

所以,可以得到:

現(xiàn)對模型(2)進行推廣,考慮誤差項εt序列相關(guān)。根據(jù)ADF檢驗,通過加入p階滯后項來代替誤差項的序列相關(guān),則模型(2)變成

(6)

式(6)中:εt～i.i.d.(0,σ2)。則相應地可以得到輔助方程:

對于統(tǒng)計量W有如下分布:

定理2考慮模型(6)在H0:α=δ=0的條件下,統(tǒng)計量W有與定理1一樣的極限分布。

證明:設T×p的數(shù)據(jù)矩陣Z=(ΔY-1,…,ΔY-p),其中ΔY-i=(Δy-i+1,…,Δy-i+T),設T×T矩陣M=IT-Z(Z′Z)-1Z′,向量ε=ε1,…,εT′,則

和

所以,類似地可以得到W?A′B-1A,

3Monte Carlo模擬

上述分析表明,在單位根檢驗下,參數(shù)約束的Wald統(tǒng)計量的極限分布都不是χ2分布。本研究采用蒙特卡羅模擬方法模擬出上述Wald統(tǒng)計量的臨界值。數(shù)據(jù)生成過程為:Δyt=εt,y0=0,εt～i.i.d.(0,1)?？疾旖y(tǒng)計量W的極限分布,樣本容量T取1 000,針對樣本容量模擬20 000次,從而得到檢驗統(tǒng)計量W常用檢驗水平0.01、0.05、0.10的漸進臨界值(表1)。

表1　統(tǒng)計量W的臨界值

對于小樣本而言,分別針對樣本容量T為50、100、200，考察檢驗統(tǒng)計量的檢驗水平。數(shù)據(jù)生成過程為:Δyt=εt,εt=ρεt-1+ut,ut～i.i.d.(0,1)。為了考察更多情況,ρ取0~0.4的均勻分布。通過表2可以看出,統(tǒng)計量W的實際檢驗水平與名義檢驗水平一致。

表2　統(tǒng)計量W的實際檢驗水平

考慮數(shù)據(jù)生成過程為ESTAR模型,現(xiàn)對式(3)的統(tǒng)計量KSS和統(tǒng)計量W的勢進行對比。不妨假設數(shù)據(jù)生成過程為:

式中,α設為0.1,γ為-1.5、-1、-0.5,θ為0.001~0.01的均勻分布。給定名義水平為0.05,從表3可以看出,統(tǒng)計量W較統(tǒng)計量KSS有更高的勢。

表3　統(tǒng)計量W和統(tǒng)計量KSS的勢

4結(jié)語

本研究主要針對ESTAR模型提出用于單位根檢驗的Wald統(tǒng)計量,通過研究得出以下結(jié)論:無論是否有p階滯后項,含有單位根約束條件的Wald統(tǒng)計量的極限分布均不服從χ2分布,而是服從以Wiener過程泛函表示的分布;蒙特卡羅模擬結(jié)果顯示,有限樣本下Wald統(tǒng)計量W的勢比KSS的勢更高,故認為統(tǒng)計量W更有效。

參考文獻:

[1]Ter?svirta T. Specification, estimation and evaluation of smooth transition autoregressive models[J]. Journal of the American Statistical Association,1994,89:208-218.

[2]Lundbergh S, Ter?svirta T. Forecasting with smooth transition autoregressive models[M]. [S.L.]: Clements M P,2001.

[3]Van Dijk D, Ter?svirta T, Franses P H. Smooth transition autoregressive models: A survey of recent developments[J]. Econometric Reviews,2002,21:1-47.

[4]Kapetanios G, Shin Y, Snell A. Testing for a unit root in the nonlinear STAR framework[J]. Journal of Econometrics,2003,112:359-379.

[5]Kruse R. A new unit root test against ESTAR based on a class of modified statistics[J]. Statistical Papers,2011,52(1):71-85.

[6]Hanck C. On the asymptotic distribution of a unit root test against ESTAR alternatives[J]. Statistics and Probability Letters,2012,82:360-364.

[7]Ubilava D, Helmers C G. Forecasting ENSO with a smooth transition autoregressive model[J].Environmental Modelling and Software,2013,40:181-190.

[8]張凌翔,張曉峒.單位根檢驗中的Wald統(tǒng)計量研究[J].數(shù)量經(jīng)濟技術(shù)經(jīng)濟研究,2009(7):146-158.

[9]聶巧平,張曉峒.ADF單位根檢驗中聯(lián)合檢驗F統(tǒng)計量研究[J].統(tǒng)計研究,2007,24(2):73-80.

Unit root test based on ESTAR process

—Study on Wald statistics

HU Junjuan

(School of Sciences, Zhejiang University of Science and Technology, Hangzhou 310023, China)

Abstract:The alternative hypothesis of exponential smooth transition autoregressive (ESTAR) nonlinearity usually assumes thatα=0, while the Existing tests of the unit root hypothesis is against that. However, empirical work indicates that the estimated value of α is always greater than zero and significant. Hence, the paper relaxes the restriction in the test regression and investigate the Wald-type test for a unit root process against ESTAR process. The asymptotic distributions of the test statistic are derived. Results show that the Wald statistic is effective via Monte Carlo simulation.

Key words:unit root test; Wald-type statistic; exponential smooth transition autoregressive

中圖分類號:O212;F222

文獻標志碼:A

文章編號:1671-8798(2015)01-0006-05

作者簡介:胡俊娟(1979—),女,浙江省蘭溪人,講師,博士研究生,主要從事時間序列分析研究。

收稿日期:2015-01-21

doi:10.3969/j.issn.1671-8798.2015.01.002