周丹文

[摘 要]從股票指數(shù)的基礎(chǔ)時(shí)間序列出發(fā)，分析了時(shí)間序列的性質(zhì)以及股票指數(shù)的特性。并以此為基礎(chǔ)，運(yùn)用經(jīng)濟(jì)統(tǒng)計(jì)與計(jì)量經(jīng)濟(jì)學(xué)，選取恒生指數(shù)2009年10月至2019年10月十年間日收盤(pán)價(jià)進(jìn)行分析，同時(shí)利用MATLAB的金融工具分析箱，建立GARCH模型，對(duì)恒生指數(shù)的日收益率進(jìn)行建模與預(yù)測(cè)。結(jié)果表明，GARCH模型能夠較好地?cái)M合樣本數(shù)據(jù)的對(duì)數(shù)收益率;在短期內(nèi)，模型同樣擁有較好的預(yù)測(cè)效果。

[關(guān)鍵詞]股票指數(shù);時(shí)間序列分析;GARCH;MATLAB

Abstract： Based on the basic time series of stock index， this paper analyzes the properties of time series and stock index. On this basis， with the knowledge of economic statistics and econometrics， the daily closing price of Hang Seng Index from October 2009 to October 2019 is selected to analysis. GARCH models are established by using the Financial Toolbox of MATLAB to model and predict the daily yield of Hang Seng Index. The results show that GARCH model can fit the logarithmic rate of return of sample data well， and in the short term， the model also has a good prediction effect.

Key Words： Stock Index; Time Series Analysis; GARCH; MATLAB

一、研究背景

（一）股票指數(shù)

股票指數(shù)是用來(lái)度量股票行情的一種指標(biāo)， 反映了股票市場(chǎng)總體價(jià)格水平及其變動(dòng)趨勢(shì)，一般由證券交易所或其他金融服務(wù)機(jī)構(gòu)編制，用以為股民提供一個(gè)衡量股市價(jià)值變化的參考依據(jù)。股票指數(shù)與其編制范圍內(nèi)的股票的平均價(jià)格水平同向變化，能靈敏反映市場(chǎng)所在國(guó)（或地區(qū)）社會(huì)、政治、經(jīng)濟(jì)變化狀況。

（二）時(shí)間序列

時(shí)間序列是指將同一統(tǒng)計(jì)指標(biāo)的數(shù)值按其發(fā)生的時(shí)間先后順序排列而成的數(shù)列。即對(duì)某一個(gè)或者一組變量X（t）進(jìn)行觀(guān)察測(cè)量，在一系列時(shí)刻t1， t2， …， tn所得的離散序列集合（X（ti）為隨機(jī)變量）為時(shí)間序列。其主要目的是根據(jù)已有的歷史數(shù)據(jù)對(duì)未來(lái)進(jìn)行預(yù)測(cè)。根據(jù)觀(guān)察時(shí)間的不同，時(shí)間序列中的時(shí)間可以是任何時(shí)間形式。股票指數(shù)即為時(shí)間序列數(shù)據(jù)。

（三）時(shí)間序列的平穩(wěn)性

平穩(wěn)性是指時(shí)間序列的行為不隨時(shí)間的改變而變化，可分為嚴(yán)平穩(wěn)與寬平穩(wěn)兩種。嚴(yán)平穩(wěn)表示時(shí)間序列的分布不隨時(shí)間的改變而變化，序列中隨機(jī)變量的性質(zhì)是一樣的;寬平穩(wěn)則沒(méi)有分布不隨時(shí)間變化的性質(zhì)，其強(qiáng)調(diào)的是相關(guān)系數(shù)取決于時(shí)間的間隔，而非時(shí)間的起點(diǎn)。寬平穩(wěn)的特性有兩點(diǎn)，一為均值函數(shù)為常數(shù)函數(shù);二為協(xié)方差函數(shù)僅與時(shí)間差相關(guān)。

根據(jù)平穩(wěn)性的定義，可以進(jìn)行平穩(wěn)性檢驗(yàn)。如果序列有明顯的趨勢(shì)或者周期性，那么其就是不平穩(wěn)的，因?yàn)槠椒€(wěn)序列具有常數(shù)均值和方差。

（四）股票指數(shù)的特點(diǎn)

股票指數(shù)的收益率一般隨時(shí)間的變化而變化，因此，它是非平穩(wěn)的。同時(shí)，股票指數(shù)的收益率還存在持續(xù)偏高或偏低的現(xiàn)象，即波動(dòng)聚集性。另一方面，股票指數(shù)的收益率序列一般不符合高斯假定中的“同方差”假定，即其序列一般具有異方差。

二、研究對(duì)象

恒生指數(shù)（HSI），由香港恒生銀行全資附屬的恒生指數(shù)公司編制的指數(shù)。入選樣板股為香港股票市場(chǎng)33家上市公司，是以發(fā)行量為權(quán)數(shù)的加權(quán)平均股價(jià)指數(shù)。為香港股票市場(chǎng)最有影響的指數(shù)。

本文選取恒生指數(shù)2009年10月8日至2019年10月8日十年間的日收益率作為研究對(duì)象，建立模型，并依據(jù)其2019年10月9日至2019年12月17日間的收益率對(duì)模型進(jìn)項(xiàng)檢驗(yàn)。

三、模型選擇

GARCH模型（廣義自回歸條件異方差模型）由Bollerslev等人在于1986年提出，該模型建立在A(yíng)RCH模型的基礎(chǔ)上，并在當(dāng)今演化出了多種形式。

傳統(tǒng)的計(jì)量經(jīng)濟(jì)學(xué)假設(shè)時(shí)間序列的方差是固定的，這種假設(shè)與實(shí)際情況相差較大。對(duì)于股票而言， 其收益的波動(dòng)幅度就是隨時(shí)間而變化的， 并非常數(shù)。這使得傳統(tǒng)的時(shí)間序列分析對(duì)實(shí)際問(wèn)題無(wú)效[1]。1982年，Engle在研究英國(guó)通貨膨脹率的波動(dòng)性時(shí)首次提出了ARCH模型，。

ARCH模型在實(shí)際運(yùn)用時(shí)， 誤差項(xiàng)的條件方差依賴(lài)于之前多期的變化量，存在較難精確估計(jì)的缺陷。GARCH模型則在此基礎(chǔ)上用較為簡(jiǎn)單的低階GARCH模型代替高階ARCH模型， 降低參數(shù)估計(jì)的復(fù)雜性，適用于有波動(dòng)性的樣本的分析與預(yù)測(cè)。GARCH模型可表示為GARCH（p，q），當(dāng)p=0時(shí)，GARCH模型則降至ARCH模型。GARCH模型假定當(dāng)前的條件方差依賴(lài)于其滯后項(xiàng)σt-1與殘差滯后項(xiàng)εt-1;而當(dāng)期方差則取決于常數(shù)項(xiàng)，前一期的殘差與前一期的預(yù)測(cè)方差。

在高斯分布下，GARCH（p，q）的函數(shù)形式如下：

四、模型建立

本次實(shí)驗(yàn)運(yùn)用MATLAB R2018b軟件進(jìn)行分析驗(yàn)證。在2014年之后的版本中，MATLAB中有關(guān)GARCH模型的函數(shù)部分發(fā)生了較大的變化。因此，本文采取了最新的函數(shù)形式進(jìn)行模型建立。

（一）數(shù)據(jù)導(dǎo)入

本文的數(shù)據(jù)來(lái)源為RESSET金融研究數(shù)據(jù)庫(kù)。首先從該數(shù)據(jù)庫(kù)中下載相應(yīng)日期及當(dāng)日恒生指數(shù)的日收盤(pán)價(jià)，并保存在EXCEL中。同時(shí)，為了便于MATLAB的讀取與識(shí)別，將EXCEL中的日期列的格式修改為“日期”。通過(guò)導(dǎo)入功能，將日期與日收盤(pán)價(jià)兩列以“數(shù)值矩陣”的形式導(dǎo)入MATLAB中。

（二）繪制收益率曲線(xiàn)

1.日收盤(pán)價(jià)

通過(guò)導(dǎo)入的日期與日收盤(pán)價(jià)，運(yùn)用“plot”函數(shù)，可以繪制出恒生指數(shù)的日收盤(pán)價(jià)波動(dòng)曲線(xiàn)。

對(duì)收盤(pán)價(jià)進(jìn)行分析可知，十年間恒生指數(shù)的收盤(pán)價(jià)總體而言呈上升趨勢(shì)，但波動(dòng)較大。在2012年、2015年與2019年，發(fā)生了三次劇烈的下降。其中，2012年是由于歐債危機(jī)惡化、歐美外圍股市下跌等因素的影響;2015是由于股災(zāi);2018年則是因?yàn)槿蚴袌?chǎng)的動(dòng)蕩加劇。由于港元采取釘住美元的匯率政策，香港市場(chǎng)對(duì)于美國(guó)利率的波動(dòng)尤為敏感。恒生指數(shù)也因此較易產(chǎn)生大幅的波動(dòng)。

2.對(duì)數(shù)收益率

恒生指數(shù)的日收盤(pán)價(jià)呈現(xiàn)出劇烈的波動(dòng)趨勢(shì)，因此，考慮在日收盤(pán)價(jià)的基礎(chǔ)上，計(jì)算其對(duì)數(shù)收益率。

對(duì)數(shù)收益率是兩個(gè)時(shí)期資產(chǎn)價(jià)值取對(duì)數(shù)后的差額，即資產(chǎn)多個(gè)時(shí)期的對(duì)數(shù)收益率等于其各時(shí)期對(duì)數(shù)收益率之和。通過(guò)取對(duì)數(shù)，數(shù)據(jù)可以變得更加平穩(wěn)，同時(shí)，也利于直接觀(guān)察收益率。在MATLAB中，可以通過(guò)“price2ret”函數(shù)對(duì)列向量進(jìn)行對(duì)數(shù)收益率的計(jì)算，并通過(guò)“plot”函數(shù)繪制出其圖像。在本例中，將錄有取對(duì)數(shù)后恒生指數(shù)十年間日收益率的列向量命名為“y”。

從圖2可以看出，恒生指數(shù)的對(duì)數(shù)收益率具有波動(dòng)聚集性，集中在（-0.04，0.04）之間，但在波動(dòng)較大的點(diǎn)達(dá)到了+0.06的水平。達(dá)到峰值后又迅速回落。但相較于日收益率，其數(shù)據(jù)已更加平穩(wěn)。

（三）分布檢驗(yàn)

通過(guò)MATLAB中的函數(shù)，可以計(jì)算峰度、偏度等數(shù)值，用以判斷數(shù)據(jù)是否符合正態(tài)分布。在此次實(shí)驗(yàn)中，相關(guān)數(shù)值如下：

通過(guò)計(jì)算樣本峰度、樣本偏差可知，數(shù)據(jù)的偏度為-0.3121，峰度為5.2792，數(shù)據(jù)分布左偏且呈現(xiàn)出厚尾特征。說(shuō)明從整體上來(lái)看，收益率低于其均值的時(shí)候較多，數(shù)據(jù)的波動(dòng)性強(qiáng)。通過(guò)jbtest檢驗(yàn)，亦可知數(shù)據(jù)拒絕了服從正態(tài)分布的假定。

（四）相關(guān)性檢驗(yàn)

1.平穩(wěn)性檢驗(yàn)

在MATLAB中，可以通過(guò)”adftest”函數(shù)判斷數(shù)據(jù)是否平穩(wěn)。經(jīng)過(guò)檢驗(yàn)，可以發(fā)現(xiàn)，雖然十年間恒生指數(shù)的日收盤(pán)價(jià)序列并不平穩(wěn)，符合股票指數(shù)的特征;但adftest（y）‘的結(jié)果為1，說(shuō)明其對(duì)數(shù)收益率序列是平穩(wěn)的。因此可以利用其十年間的對(duì)數(shù)收益率建立相關(guān)模型。

2.自相關(guān)與偏自相關(guān)圖像

通過(guò)”autocorr（y）”與”parcorr（y）”函數(shù)，可以分別表示出自相關(guān)與偏自相關(guān)的收斂階數(shù)。

通過(guò)上述自相關(guān)圖與偏自相關(guān)圖像，可知兩者均約在16階收斂，表明日收益率沒(méi)有明顯的日相關(guān)性。

3.Q檢驗(yàn)與ARCH檢驗(yàn)

在本例中，用“y-mean（y）”計(jì)算恒生指數(shù)日對(duì)數(shù)收益率的殘差。MATLAB中分別用于Q檢驗(yàn)與ARCH的函數(shù)如下：

經(jīng)過(guò)Q檢驗(yàn)，可知實(shí)驗(yàn)數(shù)據(jù)在10，15，20滯后階數(shù)內(nèi)，均呈現(xiàn)出h=0。說(shuō)明日收益率殘差與收益波動(dòng)之間不存在自相關(guān)性。ARCH檢驗(yàn)中，可知實(shí)驗(yàn)數(shù)據(jù)在10，15，20滯后階數(shù)內(nèi)，均有h=1，說(shuō)明收益率序列有高階ARCH效應(yīng)即GARCH效應(yīng)。因此考慮建立GARCH模型。

（五）具體模型的選擇

在GARCH模型中，一般有GARCH（1，1），GARCH（1，2），GARCH（2，2）等多種形式。每種形式中，又可分為有平均偏移量和無(wú)平均偏移量?jī)煞N。模型形式的選擇可以通過(guò)” 赤池信息準(zhǔn)則（AIC）”判斷，具體而言，計(jì)算出的“AIC”的數(shù)值的絕對(duì)值越大，則擬合效果越好。

其中，y指恒生指數(shù)的對(duì)數(shù)收益率，aicbic函數(shù)中的”3”，為損失的自由度，（具體指常數(shù)項(xiàng)，GARCH{1}項(xiàng)與ARCH{1}項(xiàng)），會(huì)隨著模型中參數(shù)的改變而改變。

分別建立GARCH（1，1），GARCH（1，2），GARCH（2，2）沒(méi)有平均偏移量的模型，并通過(guò)赤池信息準(zhǔn)則進(jìn)行檢驗(yàn)。經(jīng)過(guò)檢驗(yàn)可知在5%的顯著性水平下，沒(méi)有平均偏移量的GARCH（1，2）模型擬合較好（通過(guò)aicbic（）函數(shù)計(jì)算，可知其絕對(duì)值最大）。再將其與有平均偏移量的GARCH（1，2）模型進(jìn)行比較，可以發(fā)現(xiàn)仍是沒(méi)有平均偏移量的模型擬合較好。各模型的計(jì)算結(jié)果如下：

模型的γ1的系數(shù)接近于1且通過(guò)了顯著性檢驗(yàn)， 說(shuō)明恒生指數(shù)過(guò)去價(jià)格的波動(dòng)與其無(wú)限長(zhǎng)期價(jià)格波動(dòng)的大小都有關(guān)系。條件方差方程中， 系數(shù)α1，α2和γ1都為正且顯著，說(shuō)明過(guò)去的波動(dòng)對(duì)市場(chǎng)未來(lái)波動(dòng)有影響， 且影響是正向的。γ1+α1+α2的值接近于1， 這說(shuō)明股市波動(dòng)對(duì)外部沖擊的反應(yīng)函數(shù)以一個(gè)相對(duì)較慢的速度遞減， 一旦出現(xiàn)大的波動(dòng)在短期內(nèi)很難消除[2]。

通過(guò)MATLAB中的 econometricModeler，可以查看GARCH（1，2）對(duì)應(yīng)的函數(shù)形式，以及其條件方差和標(biāo)準(zhǔn)化殘差。

通過(guò)上述圖像可以看出，標(biāo)準(zhǔn)化殘差變得較為平穩(wěn)，波動(dòng)聚集現(xiàn)象已不再明顯。

從而寫(xiě)出對(duì)應(yīng)的方程：

五、模型預(yù)測(cè)

運(yùn)用MATLAB中的forecast（）函數(shù)可以進(jìn)行模型預(yù)測(cè)。在將上述GARCH（1，2）模型命名為EstMdl2的情況下，具體操作如下：

在執(zhí)行上述操作后，即得到51天內(nèi)觀(guān)測(cè)值的預(yù)測(cè)結(jié)果。再對(duì)foreY取對(duì)數(shù)收益率，以10天為一組，將計(jì)算出的結(jié)果與恒生指數(shù)真實(shí)的收益率進(jìn)行比較，結(jié)果如下：

綜合預(yù)測(cè)與實(shí)際的結(jié)果，可以發(fā)現(xiàn)，總體而言，隨著時(shí)間的推移，預(yù)測(cè)值偏離真實(shí)值的范圍越大，預(yù)測(cè)的準(zhǔn)確率越低。因此，在短期內(nèi)，GARCH模型有較好的擬合效果。

六、總結(jié)

本文選用GARCH模型分析股票指數(shù)，模擬了股票指數(shù)波動(dòng)性變化。結(jié)果說(shuō)明GARCH模型對(duì)于股票指數(shù)的日收益率有較好的擬合效果，同時(shí)短期預(yù)測(cè)有較高的準(zhǔn)確度。股票指數(shù)屬于時(shí)間序列的一種，但其不同于傳統(tǒng)的數(shù)據(jù)，呈現(xiàn)出非平穩(wěn)、波動(dòng)聚集與異方差的特征。因此，對(duì)于它的建模也應(yīng)當(dāng)采用非傳統(tǒng)的方法。GARCH模型用方差預(yù)測(cè)時(shí)間序列，對(duì)股票指數(shù)有較好的擬合與預(yù)測(cè)結(jié)果。因此，對(duì)于投資者的決策能起到一定的指導(dǎo)作用。

[參考文獻(xiàn)]

[1]俞瑩，易榮華.來(lái)自新興市場(chǎng)的公司上市對(duì)目的地市場(chǎng)的影響——以中國(guó)香港市場(chǎng)為例[J].中國(guó)管理科學(xué)，2020，28（2）：37-47.

[2]吳小丹.中美貿(mào)易摩擦以來(lái)對(duì)我國(guó)大豆行業(yè)影響及應(yīng)對(duì)策略[J].對(duì)外經(jīng)貿(mào)，2019（11）：6-11.

[3]徐德貴，李慧，王月月.中證500股指期貨對(duì)股票現(xiàn)貨市場(chǎng)波動(dòng)性影響研究[J].經(jīng)濟(jì)研究導(dǎo)刊，2019（20）：136-140.

[4]黃凝.中國(guó)股票市場(chǎng)指數(shù)波動(dòng)特征研究[J].現(xiàn)代營(yíng)銷(xiāo)（信息版），2019（7）：108-109.

[5]張露. 融資融券對(duì)A股波動(dòng)性的影響研究[D].西北大學(xué)，2019.

[6]宋邵陽(yáng). A股中小板塊股票beta系數(shù)均值回歸特性研究[D].云南大學(xué)，2019.

[7]馬勇，何順.貨幣周期與資產(chǎn)定價(jià)：基于中國(guó)的實(shí)證研究[J].金融評(píng)論，2019，11（2）：1-24+123.

（責(zé)任編輯：顧曉濱 馬 琳）