朱 旭,龐天曉

(浙江大學(xué) 數(shù)學(xué)科學(xué)學(xué)院,浙江杭州 310058)

§1 引言

許多經(jīng)濟(jì)和金融時(shí)間序列的數(shù)據(jù)都具有單變點(diǎn)或者多變點(diǎn)的特征,參見(jiàn)Bai和Perron[1],Hansen[2],Lee等[3],以及Perron和Zhu[4].此外眾所周知,經(jīng)濟(jì)和金融數(shù)據(jù)往往還具有長(zhǎng)相依的特征.例如日元兌美元的匯率數(shù)據(jù)被廣泛認(rèn)為存在長(zhǎng)相依,參見(jiàn)Horv′ath和Kokoszka[5].因此在處理實(shí)際數(shù)據(jù)時(shí),往往面臨結(jié)構(gòu)變點(diǎn)和長(zhǎng)相依同時(shí)存在的情形.例如Jaruˇskov′a[6]在非平穩(wěn)的長(zhǎng)周期的水文數(shù)據(jù)中發(fā)現(xiàn)了長(zhǎng)相依的證據(jù),同時(shí)認(rèn)為數(shù)據(jù)中存在結(jié)構(gòu)變點(diǎn).大量的研究表明,長(zhǎng)相依的存在會(huì)使傳統(tǒng)的變點(diǎn)理論失效.因此對(duì)長(zhǎng)相依數(shù)據(jù)進(jìn)行變點(diǎn)分析既具有理論意義,也具有實(shí)際意義.

有相當(dāng)多的文章研究了長(zhǎng)相依時(shí)間序列中的結(jié)構(gòu)變點(diǎn).例如Kuan和Hsu[7]研究了均值變點(diǎn)并給出了估計(jì)量的收斂速度.Chang和Perron[8]研究了參數(shù)為d的分整(fractionally integrated)過(guò)程中的斜率和均值變點(diǎn),并給出了估計(jì)量的收斂速度和極限分布.Perron和Qu[9]提出用周期圖來(lái)研究股票指數(shù)回報(bào)的波動(dòng)率,發(fā)現(xiàn)存在長(zhǎng)相依和均值漂移.Shao[10]提出一種檢驗(yàn)流程來(lái)檢驗(yàn)長(zhǎng)相依時(shí)間序列中的均值變點(diǎn).Iacone等作者[11]提出一種檢驗(yàn)流程來(lái)檢驗(yàn)參數(shù)為d的分整時(shí)間序列中的斜率變點(diǎn).

面板數(shù)據(jù)中存在公共變點(diǎn)的現(xiàn)象在實(shí)際生活中非常普遍.例如稅收政策的變化可能會(huì)影響每個(gè)人的收入.另一個(gè)著名的例子是2008年發(fā)生的全球金融危機(jī)幾乎影響了所有國(guó)家的GDP(Gross Domestic Product)數(shù)據(jù)和資本市場(chǎng)表現(xiàn),并帶來(lái)一定時(shí)間內(nèi)持續(xù)的影響.類(lèi)似地,新的科學(xué)技術(shù)的誕生,新藥的研發(fā)生產(chǎn),以及政府層面新的政策措施都會(huì)對(duì)經(jīng)濟(jì)社會(huì)造成或多或少的影響.

結(jié)構(gòu)變點(diǎn)從時(shí)間序列的研究延申到面板數(shù)據(jù)的研究可以追溯到1990年.Joseph和Wolfson是最早的一批研究面板數(shù)據(jù)中的結(jié)構(gòu)變點(diǎn)的學(xué)者[12-13],他們提出了N個(gè)面板上存在獨(dú)立同分布的N個(gè)變點(diǎn)的隨機(jī)變點(diǎn)模型,并給出了變點(diǎn)分布的相合估計(jì)量.不久后隨機(jī)變點(diǎn)模型被Joseph等[14]推廣到了自回歸模型.Bai[15]利用最小二乘的方法估計(jì)了面板數(shù)據(jù)中的公共均值變點(diǎn),同時(shí)證實(shí)了變點(diǎn)估計(jì)量是相合的,這與平穩(wěn)時(shí)間序列中變點(diǎn)估計(jì)量是不相合的結(jié)論截然相反.值得一提的是,以上文章中的面板數(shù)據(jù)都是假設(shè)橫截面獨(dú)立的.此后,面板數(shù)據(jù)中的公共變點(diǎn)研究得到了計(jì)量經(jīng)濟(jì)學(xué)界和統(tǒng)計(jì)學(xué)界的大量關(guān)注,橫截面獨(dú)立的面板數(shù)據(jù)模型也在一些文章中被推廣到了橫截面相依情形.例如,Horv′ath和Huˇskov′a[16]研究了存在橫截面相依的面板數(shù)據(jù)中公共均值變點(diǎn)的檢驗(yàn)問(wèn)題,Kim[17-18]研究了橫截面相依的面板數(shù)據(jù)中存在斜率變點(diǎn)以及斜率與截距變點(diǎn)同時(shí)存在時(shí)的變點(diǎn)估計(jì)問(wèn)題,Li等作者[19]利用自適應(yīng)group fused LASSO懲罰的主成分方法估計(jì)橫截面相依的面板數(shù)據(jù)中的變點(diǎn),Qian和Su[20]提出了橫截面相依面板數(shù)據(jù)中的公共變點(diǎn)的收縮估計(jì)量,Baltagi等作者[21-22]提出了CCE(common correlated effects)估計(jì)量并把它應(yīng)用于橫截面相依面板數(shù)據(jù)中變點(diǎn)的研究,Westerlund[23]采用了基于CCE的方法對(duì)時(shí)間長(zhǎng)度固定的橫截面相依面板數(shù)據(jù)中的變點(diǎn)進(jìn)行研究.值得注意的是,大部分現(xiàn)有的文獻(xiàn)都假設(shè)面板數(shù)據(jù)在時(shí)間維度上是弱相依的,研究長(zhǎng)相依面板數(shù)據(jù)中的結(jié)構(gòu)變點(diǎn)的文章相對(duì)較少.

本文假設(shè)面板數(shù)據(jù)具有橫截面的相依性,同時(shí)每個(gè)時(shí)間序列都是長(zhǎng)相依的(記憶參數(shù)d ∈(0,0.5))且?guī)в行甭首凕c(diǎn).這種模型在宏觀(guān)經(jīng)濟(jì)和金融中具有廣泛的應(yīng)用.用最小二乘的方法估計(jì)斜率變點(diǎn)發(fā)生的時(shí)刻以及發(fā)生時(shí)刻的分?jǐn)?shù),并研究了當(dāng)(T,N)(T表示時(shí)間長(zhǎng)度,N表示個(gè)體數(shù))聯(lián)合趨于無(wú)窮時(shí)估計(jì)量的漸近性質(zhì): 估計(jì)量的相合性,收斂速度以及極限分布.得到了一些有趣的結(jié)論: 大部分情況下斜率變點(diǎn)估計(jì)量的收斂速度隨著記憶參數(shù)d的增大而減緩,但當(dāng)T和N滿(mǎn)足T2d=o(N),且公因子與變點(diǎn)的變化幅度存在交互效應(yīng)時(shí),變點(diǎn)時(shí)刻估計(jì)量的收斂速度僅與T有關(guān).

本文結(jié)構(gòu)如下:§2給出模型的具體形式以及一些假設(shè),§3介紹變點(diǎn)估計(jì)量的漸近性質(zhì),§4給出Monte Carlo模擬結(jié)果,用于評(píng)估變點(diǎn)估計(jì)量在有限樣本情形下的表現(xiàn),§5對(duì)論文進(jìn)行了總結(jié),文章中的理論證明放在§6.

§2 模型與假設(shè)

首先介紹本文接下來(lái)需要用到的一些記號(hào).“?”表示概率測(cè)度弱收斂(參見(jiàn)Billingsley[24]),“”表示依概率收斂,“”表示依分布收斂,“iid”表示獨(dú)立同分布,aTbT表示存在兩個(gè)正常數(shù)c1和c2使得對(duì)于所有足夠大的T,c1≤aT/bT ≤c2成立,其中aT和bT是兩個(gè)關(guān)于T的取值恒為正的函數(shù).設(shè)D[0,1]是[0,1]上的函數(shù)空間,且是存在左極限和右連續(xù)的Skorohod拓?fù)?對(duì)于矩陣A,記‖A‖=[tr(A′A)]1/2,其中tr[·]表示矩陣的跡.此外用M表示正常數(shù),它的取值與T和N無(wú)關(guān),但在不同的地方可取不同的值.除非另有說(shuō)明,本文中的極限均理解為(T,N)→∞,即T和N聯(lián)合趨于無(wú)窮.

本文研究的面板數(shù)據(jù)模型為

注假設(shè)1是為了保證變點(diǎn)的可識(shí)別性,是變點(diǎn)研究領(lǐng)域中的常見(jiàn)假設(shè).假設(shè)2(1)假設(shè)公因子是一個(gè)線(xiàn)性過(guò)程并可以應(yīng)用泛函中心極限定理(參見(jiàn)Phillips和Solo[26]).假設(shè)2(2)假設(shè)誤差ei,t在時(shí)間維度上是長(zhǎng)相依的,但其方差在橫截面維度上可以不同.假設(shè)2(3)是一個(gè)技術(shù)性假設(shè),是為了應(yīng)用面板數(shù)據(jù)的聯(lián)合中心極限定理(參見(jiàn)Phillips和Moon[27]).假設(shè)2(4)假設(shè)公因子和模型誤差是相互獨(dú)立的.假設(shè)3和4描述了模型的總體變點(diǎn)強(qiáng)度,公因子的強(qiáng)度,以及變點(diǎn)的變化幅度與公因子之間的相互影響程度.

§3 極限性質(zhì)

本節(jié)介紹估計(jì)量?T1的極限性質(zhì).后文中,所有的Op(·)將按照其嚴(yán)格的定義理解,即隨機(jī)變量不會(huì)是op(·).

定理3.1若假設(shè)1-4成立,當(dāng)(T,N)→∞,d ∈(0,0.5)時(shí),有下面的結(jié)論.

注以上結(jié)論都表明變點(diǎn)估計(jì)量是相合的.結(jié)論(1)說(shuō)明在不存在橫截面相依時(shí),收斂速度與記憶參數(shù)d,T和N有關(guān),且d越大收斂速度越慢,這與直覺(jué)吻合,因?yàn)閐 ≥0.5時(shí)時(shí)間序列將變成非平穩(wěn).結(jié)論(2)說(shuō)明當(dāng)AHγ/0時(shí),即公因子與變點(diǎn)的變化幅度存在交互效應(yīng),如果個(gè)體數(shù)N相對(duì)于時(shí)間長(zhǎng)度T并沒(méi)有足夠大時(shí),收斂速度依然由d,T和N共同決定,且隨著d的增大而變慢.結(jié)論(3)說(shuō)明當(dāng)T2d=o(N),即個(gè)體數(shù)N相對(duì)于時(shí)間長(zhǎng)度T足夠大時(shí),的收斂速度僅與T有關(guān).

有了收斂速度,接下來(lái)推導(dǎo)估計(jì)量的極限分布.有下面的定理.

定理3.2若假設(shè)1-4成立,當(dāng)(T,N)→∞,d ∈(0,0.5)時(shí),有下面的結(jié)論.

注以上的三個(gè)結(jié)論表明: 當(dāng)變點(diǎn)分?jǐn)?shù)的真值λ0越接近中間值1/2,或者斜率變點(diǎn)的信號(hào)強(qiáng)度越大(即Aγγ越大),則漸近分布的方差越小.

§4 數(shù)值模擬

在這一節(jié),將通過(guò)Monte Carlo模擬來(lái)評(píng)估估計(jì)量在有限樣本情形下的表現(xiàn).對(duì)所有實(shí)驗(yàn),重復(fù)次數(shù)都為2000次.數(shù)據(jù)由以下過(guò)程產(chǎn)生.設(shè)公因子Ft是一維的且滿(mǎn)足Ft=0.2Ft?1+wt,F0=0,wt獨(dú)立同分布于N(0,1);一維因子載荷hi有兩種情況: 若不存在橫截面相依,則對(duì)所有的i取hi=0,若公因子與變點(diǎn)的變化幅度存在交互效應(yīng),即AHγ/0,設(shè)hi由U(0,1)隨機(jī)生成;參考McLeod和Hipel[28]以及Hosking[29]的方法,隨機(jī)生成N個(gè)長(zhǎng)相依時(shí)間序列{ei,t,t=1···,T},i=1,···,N,使得(1?L)dei,t=εi,t,εi,t=σiηi,t,其中σi由U(1,1.5) 隨機(jī)生成,{ηi,t,t ≥1,i ≥1}獨(dú)立同分布于N(0,1),記憶參數(shù)根據(jù)具體情況從d ∈{0.05,0.1,0.25,0.4,0.45}中選取.為方便起見(jiàn),設(shè)定==0,因?yàn)檫@些參數(shù)的大小并不會(huì)影響Monte Carlo模擬的結(jié)果,影響模擬結(jié)果的是斜率的變化幅度.設(shè)由U(0,0.3)隨機(jī)生成.此外,設(shè)定變點(diǎn)時(shí)刻的真值=0.5T,以及π=0.05.接下來(lái),用直方圖的形式來(lái)評(píng)估在定理3.1中的表現(xiàn).

圖1展示了(T,N)∈ {(100,10),(100,40),(200,10),(200,20)},d ∈ {0.1,0.25,0.4},所有的hi=0時(shí),的直方圖.橫向比較第一行可以看出估計(jì)量的估計(jì)誤差隨記憶參數(shù)d的增大而增加;依次縱向比較可以看出固定d時(shí)估計(jì)誤差隨T和N的增大而減小,這與定理3.1(1)的結(jié)論相符.圖2展示了(T,N)∈{(200,10),(200,20),(400,10),(400,20)},d ∈{0.25,0.4,0.45},公因子與變點(diǎn)的變化幅度存在交互效應(yīng)時(shí),的直方圖,這些T,N,d的設(shè)定是為了保證N=o(T2d)或T2d N能夠被滿(mǎn)足.同樣通過(guò)橫向和縱向的比較,可以看出這與定理3.1(2)中的結(jié)論相符.圖3則展示了(T,N)∈{(40,100),(40,200),(100,100),(100,200)},d ∈{0.05,0.1,0.25},公因子與變點(diǎn)的變化幅度存在交互效應(yīng)時(shí),的直方圖,這些T,N,d的設(shè)定是為了保證T2d=o(N)能夠被滿(mǎn)足.橫向比較第一行可以看出變點(diǎn)估計(jì)量的估計(jì)誤差與d無(wú)關(guān),縱向比較第一列可以看出估計(jì)誤差與N亦無(wú)關(guān),縱向比較第二和三列則可看出變點(diǎn)估計(jì)量的估計(jì)誤差與T有關(guān),T越大估計(jì)誤差越小.以上模擬結(jié)果與定理3.1(3)中的結(jié)論相符.

圖1 當(dāng)所有的hi=0,且N=o(T3?2d)時(shí),的直方圖

圖2 當(dāng)AHγ0,且N=o(T2d)或T2d N時(shí),的直方圖

圖3 當(dāng)AHγ/0,且T2d=o(N)時(shí),的直方圖

最后比較估計(jì)量的有限樣本分布與定理3.2的極限分布.圖4和圖5展示T1/2?dN1/2()在不存在橫截面相依和公因子與變點(diǎn)的變化幅度存在交互效應(yīng)這兩種情況下,估計(jì)量的有限樣本分布與極限分布.可以看出d越小,有限樣本分布與極限分布的逼近效果越好.其次,對(duì)于相同的N和d,不存在橫截面相依的情況下有限樣本分布更接近極限分布(即使T更小),這說(shuō)明公因子的存在會(huì)影響估計(jì)量的有限樣本分布.圖6展示了當(dāng)公因子與變點(diǎn)的變化幅度存在交互效應(yīng),且N/T2d收斂到不為零的常數(shù)時(shí),T1/2?dN1/2(?)的有限樣本分布與極限分布.可以看出,兩者的吻合程度較高.最后,圖7展示了當(dāng)公因子與變點(diǎn)的變化幅度存在交互效應(yīng)且T2d=o(N) 時(shí),T1/2()的有限樣本分布與極限分布.可以看出,有限樣本分布與極限分布的吻合度與N和d的關(guān)聯(lián)較小,但當(dāng)T變大時(shí),有限樣本分布更接近極限分布.這些模擬結(jié)果與定理3.2(3)中的結(jié)論相符.

圖4 當(dāng)所有的hi=0,且N=o(T3?2d)時(shí), T1/2?dN1/2( ?)的有限樣本分布和極限分布

圖5 當(dāng)AHγ 0,且N=o(T2d)時(shí), T1/2?dN1/2( ?)的有限樣本分布和極限分布

圖6 當(dāng)AHγ 0,且N/T2d收斂到不為零的常數(shù)時(shí), T1/2?dN1/2( ?)的有限樣本分布和極限分布

圖7 當(dāng)AHγ0,且T2d=o(N)時(shí), T1/2?dN1/2( ?)的有限樣本分布和極限分布

§5 結(jié)論

本文研究了面板數(shù)據(jù)同時(shí)存在結(jié)構(gòu)變點(diǎn),時(shí)間序列長(zhǎng)相依以及橫截面相依的變點(diǎn)估計(jì)問(wèn)題.用最小二乘的方法估計(jì)了變點(diǎn)發(fā)生的時(shí)刻,并討論了估計(jì)量的相合性,收斂速度以及極限分布.本文的結(jié)論表明在不存在橫截面相依時(shí),變點(diǎn)時(shí)刻的估計(jì)量的收斂速度與記憶參數(shù)d,時(shí)間序列長(zhǎng)度T以及個(gè)體數(shù)N有關(guān),且d越大收斂速度越慢;當(dāng)公因子與變點(diǎn)的變化幅度存在交互效應(yīng),且N相對(duì)于T并沒(méi)有足夠大時(shí),收斂速度依然由d,T和N共同決定,并隨著d的增大而變慢;然而當(dāng)T2d=o(N),即N相對(duì)于T足夠大時(shí),變點(diǎn)時(shí)刻估計(jì)量的收斂速度將僅與T有關(guān).最后通過(guò)Monte Carlo模擬比較了變點(diǎn)估計(jì)量的有限樣本分布和極限分布,并印證了本文的理論成果.

§6 證 明

由式(9)可知,若對(duì)所有i都有hi=0,則SXU僅剩其中的第二項(xiàng);若不然,還需要比較式(11)和式(14)階的大小,即比較T2d和N的相對(duì)大小.若N=o(T2d),則SXU的階將由式(9)中的第二項(xiàng)決定;若T2d=o(N),則SXU的階將由式(9)中的第一項(xiàng)決定;若T2d和N同階,則SXU的階將由式(9)整體決定.因此在λ ∈(0,1)上一致地有