公司股票對(duì)數(shù)收益率的分析與預(yù)測(cè)
——以T能源公司為例

程瀚哲 首都經(jīng)濟(jì)貿(mào)易大學(xué)

股票是現(xiàn)代金融市場(chǎng)中最常見(jiàn)、最重要的金融工具。股票市場(chǎng)在市場(chǎng)經(jīng)濟(jì)活動(dòng)中起到了至關(guān)重要的作用。股票市場(chǎng)可以讓企業(yè)發(fā)行股票籌集生產(chǎn)資本。企業(yè)可以通過(guò)股市傳播企業(yè)信息以起到宣傳作用、能夠反映企業(yè)自身的真實(shí)運(yùn)營(yíng)狀況及優(yōu)化經(jīng)濟(jì)市場(chǎng)資源配置。對(duì)于投資者來(lái)說(shuō)，股票市場(chǎng)有了更多的投資原則，擴(kuò)大了投資的選擇范圍，股票市場(chǎng)上可以隨時(shí)進(jìn)行股票交易也降低了投資者投資所需要承擔(dān)的風(fēng)險(xiǎn)。對(duì)于國(guó)家來(lái)說(shuō)，股票市場(chǎng)一直都是市場(chǎng)經(jīng)濟(jì)的晴雨表，此外很重要的是能使社會(huì)的閑散資金合理流動(dòng)。而沒(méi)有股票市場(chǎng)，常規(guī)的商品交易很難像股市一樣集成全國(guó)乃至全球范圍的，因?yàn)槠涮摂M資本運(yùn)作的形式而擺脫了地域的限制來(lái)統(tǒng)合整體的經(jīng)濟(jì)活動(dòng)[1]。

股票市場(chǎng)的研究和預(yù)測(cè)難度很大，因?yàn)橛绊懝善笔袌?chǎng)的因素眾多，國(guó)家政策、企業(yè)發(fā)展、投資情緒以及普遍存在大量的隨機(jī)事件導(dǎo)致股票市場(chǎng)變幻莫測(cè)，從因素上對(duì)股市進(jìn)行分析很難有準(zhǔn)確的預(yù)測(cè)。因此時(shí)間序列模型往往是用來(lái)進(jìn)行股票價(jià)格預(yù)測(cè)的普遍選項(xiàng)。時(shí)間序列分析是通過(guò)數(shù)據(jù)本身的變化來(lái)反應(yīng)數(shù)據(jù)走向規(guī)律的一種數(shù)據(jù)分析模型，常用的模型包括自回歸AR模型、移動(dòng)平均MA模型、自回歸移動(dòng)平均ARMA模型、差分自回歸移動(dòng)平均ARIMA模型、自回歸條件異方差A(yù)RCH模型、廣義自回歸條件異方差GARCH模型等。應(yīng)用合適的模型是時(shí)間序列分析的重要步驟之一，將股票數(shù)據(jù)擬合為恰當(dāng)?shù)哪Ｐ湍軌蚪o股票價(jià)格帶來(lái)有價(jià)值的預(yù)測(cè)結(jié)果。

一、數(shù)據(jù)背景

本文所使用的數(shù)據(jù)集是來(lái)源于Kaggle平臺(tái)的公開(kāi)數(shù)據(jù)集：https：//www.kaggle.com/rpaguirre/tesla-stockprice。該數(shù)據(jù)集記錄了T公司從2010年6月29日到2017年3月17日每天的股票價(jià)格數(shù)據(jù)。本文選取其中的收盤(pán)價(jià)序列（Close）進(jìn)行分析，計(jì)算其更具研究意義的對(duì)數(shù)收益率，并對(duì)此對(duì)數(shù)收益率序列進(jìn)行時(shí)間序列建模和預(yù)測(cè)。

二、收盤(pán)價(jià)序列的描述性分析

首先，繪制出每日收盤(pán)價(jià)序列的時(shí)序圖（見(jiàn)圖1）。

從收盤(pán)價(jià)序列時(shí)序圖中可以看出，數(shù)據(jù)總體呈現(xiàn)出上升趨勢(shì)，并且具有隨著股價(jià)增加波動(dòng)幅度更大的跡象，該序列的標(biāo)準(zhǔn)差近似與其水平值成正比。因此，本文根據(jù)此特征，需對(duì)此收盤(pán)價(jià)序列進(jìn)行對(duì)數(shù)變換，以調(diào)整該序列的方差至近似一致。

然后，繪制出對(duì)數(shù)收盤(pán)價(jià)序列的時(shí)序圖及其延遲100階的樣本自相關(guān)圖（見(jiàn)圖2）。

圖2 對(duì)數(shù)收盤(pán)價(jià)序列時(shí)序圖

經(jīng)觀察，對(duì)數(shù)變換后的收盤(pán)價(jià)序列在方差一致性方面確實(shí)有了顯著的改善。此外，對(duì)數(shù)收盤(pán)價(jià)序列仍具有明顯的上升趨勢(shì)，且不具有短期相關(guān)性。因此，本文認(rèn)為此序列非平穩(wěn)，對(duì)于其是否為差分非平穩(wěn)，還需做單位根檢驗(yàn)來(lái)進(jìn)行進(jìn)一步的探究。

圖3：對(duì)數(shù)收盤(pán)價(jià)序列ACF圖

三、模型識(shí)別

（一）ADF檢驗(yàn)

對(duì)上述對(duì)數(shù)收盤(pán)價(jià)序列做一階差分，記差分后的數(shù)據(jù)為對(duì)數(shù)收益率，接下來(lái)，根據(jù)對(duì)數(shù)收益率序列的樣本偏自相關(guān)函數(shù)圖和R語(yǔ)言中的ar函數(shù)選擇ADF檢驗(yàn)時(shí)自回歸模型的階數(shù)，對(duì)數(shù)收益率序列擬合自回歸模型時(shí)AR項(xiàng)階數(shù)選擇為0，于是ADF檢驗(yàn)就退化為DF檢驗(yàn)，對(duì)數(shù)收盤(pán)價(jià)序列DF檢驗(yàn)結(jié)果如下：

表1 ar函數(shù)對(duì)對(duì)數(shù)收益率序列的估計(jì)結(jié)果

由于p值大于0.05，于是在0.05的顯著性水平下不拒絕原假設(shè)，即認(rèn)為對(duì)數(shù)收盤(pán)價(jià)序列單位根非平穩(wěn)。

進(jìn)而需要檢驗(yàn)對(duì)數(shù)收益率序列是否單位根非平穩(wěn)（步驟同上述DF檢驗(yàn)）：

表2 ar函數(shù)對(duì)對(duì)數(shù)收益率一次差分序列的估計(jì)結(jié)果

由于p值小于0.05，于是在0.05的顯著性水平下拒絕原假設(shè)，即認(rèn)為對(duì)數(shù)收益率序列平穩(wěn)，因此不需要對(duì)對(duì)數(shù)收盤(pán)價(jià)序列做二次差分。

綜上，對(duì)數(shù)收盤(pán)價(jià)序列非平穩(wěn)。若想實(shí)現(xiàn)平穩(wěn)，則應(yīng)對(duì)此序列進(jìn)行一次差分處理，形成對(duì)數(shù)收益率序列。

（二）對(duì)數(shù)收益率序列的描述性分析

首先，計(jì)算對(duì)數(shù)收益率序列的描述統(tǒng)計(jì)量。

表3 對(duì)數(shù)收益率序列的描述統(tǒng)計(jì)量

由樣本均值的t檢驗(yàn)結(jié)果知：對(duì)數(shù)收益率序列的總體均值并沒(méi)有顯著地不同于0；正態(tài)性檢驗(yàn)的結(jié)果顯示此序列并不服從于正態(tài)分布；樣本偏度為正，說(shuō)明此序列右偏。

繪制對(duì)數(shù)收益率序列的樣本自相關(guān)函數(shù)圖，結(jié)合對(duì)數(shù)收益率序列的樣本ACF圖和PACF圖，以及其總體均值不顯著的特征，初步認(rèn)定此序列為白噪聲序列，接下來(lái)需進(jìn)一步做Box檢驗(yàn)，檢驗(yàn)結(jié)果p為0.4198。

由于p值顯著大于顯著性水平0.05，所以該序列不能拒絕白噪聲的原假設(shè)。

（三）ARCH效應(yīng)檢驗(yàn)

ARCH模型是自回歸條件異方差模型，解決了傳統(tǒng)計(jì)量經(jīng)濟(jì)學(xué)時(shí)間序列變量的第二個(gè)假設(shè)方差恒定所引起的問(wèn)題，應(yīng)用于波動(dòng)性有關(guān)的廣泛的研究領(lǐng)域，可以用來(lái)預(yù)測(cè)股票的波動(dòng)率，來(lái)對(duì)風(fēng)險(xiǎn)進(jìn)行控制，等等。不同于大部分預(yù)測(cè)被解釋變量期望著的模型，ARCH模型能預(yù)測(cè)被解釋變量的方差，準(zhǔn)確地模擬時(shí)間序列變量波動(dòng)性的變化。ARCH模型主要解釋序列中明顯的變化是否具有規(guī)律，說(shuō)明這種變化依存的是內(nèi)在傳導(dǎo)而不是外生性結(jié)構(gòu)變化，能夠提高預(yù)測(cè)精度[2]。

繪制對(duì)數(shù)收益率絕對(duì)值序列的樣本自相關(guān)圖，并對(duì)其進(jìn)行白噪聲檢驗(yàn)。檢驗(yàn)結(jié)論為對(duì)數(shù)收益率絕對(duì)值序列存在序列自相關(guān)。白噪聲檢驗(yàn)結(jié)果也表明絕對(duì)值變換后的序列為非白噪聲序列。

繪制對(duì)數(shù)收益率平方序列的樣本自相關(guān)圖，并對(duì)其進(jìn)行白噪聲檢驗(yàn)。

對(duì)數(shù)收益率平方序列存在序列自相關(guān)，白噪聲檢驗(yàn)結(jié)果也表明平方變換后的序列為非白噪聲序列。

非線性變換后不是獨(dú)立的，說(shuō)明原序列也是非獨(dú)立的，因此，對(duì)數(shù)收益率序列具有某種自相關(guān)性，而且這種自相關(guān)性會(huì)造成序列的條件異方差性?；诖擞^點(diǎn)，需對(duì)對(duì)數(shù)收益率序列進(jìn)一步做ARCH效應(yīng)檢驗(yàn)。

對(duì)此對(duì)數(shù)收益率序列做McLeod-Li檢驗(yàn)，McLeod-Li檢驗(yàn)在5%的顯著性水平下都顯著，給數(shù)據(jù)具有ARCH特征提供了強(qiáng)有力的證據(jù)。

對(duì)對(duì)數(shù)收益率序列做拉格朗日乘子檢驗(yàn)，LM檢驗(yàn)在5%的顯著性水平下都顯著，給數(shù)據(jù)具有ARCH特征提供了強(qiáng)有力的證據(jù)。

（四）GARCH模型定階

Bollerslve在Engle的基礎(chǔ)上，借助自回歸移動(dòng)平均ARMA模型的建立，建立了廣義自回歸條件異方差GARCH模型。他彌補(bǔ)了ARCH（q）模型中階數(shù)q過(guò)大時(shí)需要過(guò)多的參數(shù)的缺陷。這是一種允許異方差中同時(shí)存在自回歸項(xiàng)和滑動(dòng)平均項(xiàng)的模型，GARCH（p，q）模型。顯然，如果p=0時(shí)，GARCH（0，q）模型就是ARCH（q）模型，可以將GARCH模型看成是ARCH模型的推廣，或者將ARCH模型看作是GARCH模型的特殊情形。

由于對(duì)數(shù)收益率序列ARCH效應(yīng)的存在，所以應(yīng)對(duì)此序列擬合ARCH模型或GARCH模型，首先需確定模型的階數(shù)。

根據(jù)ARCH模型定階應(yīng)用平方序列的樣本偏自相關(guān)函數(shù)圖和GARCH模型定階應(yīng)用平方序列的樣本EACF圖所示結(jié)果，建議設(shè)定ARCH（2）模型或GARCH（1，1）模型。

四、參數(shù)估計(jì)

（一）ARCH（2）模型的估計(jì)（應(yīng)用正態(tài)分布的新息）

對(duì)數(shù)收益率應(yīng)用正態(tài)分布的新息建立ARCH（2）模型，記作m1，對(duì)m1模型進(jìn)行參數(shù)估計(jì)。

此模型的各參數(shù)估計(jì)值都顯著，模型表達(dá)式為：

檢驗(yàn)結(jié)果顯示：均值方程和波動(dòng)率方程都充分?jǐn)M合了數(shù)據(jù)，模型AIC為-4.05562，但標(biāo)準(zhǔn)化殘差并不滿足正態(tài)性假定。

（二）ARCH（2）模型的估計(jì)（應(yīng)用學(xué)生t分布的新息）

對(duì)數(shù)收益率應(yīng)用學(xué)生t分布的新息建立ARCH（2）模型，記作m2，對(duì)m2模型進(jìn)行參數(shù)估計(jì)。

此模型的各參數(shù)估計(jì)值都顯著，模型表達(dá)式為：

檢驗(yàn)結(jié)果顯示：均值方程和波動(dòng)率方程都充分?jǐn)M合了數(shù)據(jù)，模型AIC為-4.232647。

（三）GARCH（1，1）模型的估計(jì)（應(yīng)用正態(tài)分布的新息）

對(duì)數(shù)收益率應(yīng)用正態(tài)分布的新息建立GARCH（1，1）模型，記作m3，對(duì)m3模型進(jìn)行參數(shù)估計(jì)。

此模型的各參數(shù)估計(jì)值都顯著，模型表達(dá)式為：

檢驗(yàn)結(jié)果顯示：均值方程和波動(dòng)率方程都充分?jǐn)M合了數(shù)據(jù)，模型AIC為-4.061315，但標(biāo)準(zhǔn)化殘差并不滿足正態(tài)性假定。

（四）GARCH（1，1）模型的估計(jì)（應(yīng)用學(xué)生t分布的新息）

對(duì)數(shù)收益率應(yīng)用正態(tài)分布的新息建立GARCH（1，1）模型，記作m4，對(duì)m4模型進(jìn)行參數(shù)估計(jì)。

此模型的各參數(shù)估計(jì)值除omega外都顯著，但考慮到條件方差的長(zhǎng)期水平不能為0，于是參數(shù)omega不予剔除。模型表達(dá)式為：

檢驗(yàn)結(jié)果顯示：均值方程和波動(dòng)率方程都充分?jǐn)M合了數(shù)據(jù)，模型AIC為-4.254761。

（五）模型比較

匯總上述通過(guò)檢驗(yàn)的模型的AIC值，并進(jìn)行樣本內(nèi)比較：

表4 模型AIC值匯總

m4模型的AIC值最小，因此，認(rèn)為m4模型為相對(duì)最優(yōu)模型。

五、模型診斷

畫(huà)出對(duì)數(shù)收益率m4模型的波動(dòng)率圖像（見(jiàn)圖4）。

圖4 m4模型的波動(dòng)率圖像

畫(huà)出對(duì)數(shù)收益率m4模型的標(biāo)準(zhǔn)化殘差時(shí)序圖（見(jiàn)圖5）。

圖5 m4模型的標(biāo)準(zhǔn)化殘差時(shí)序圖

由上圖可見(jiàn)，標(biāo)準(zhǔn)化殘差序列基本在水平線0的上下波動(dòng)，且波動(dòng)幅度近似一致。

畫(huà)出對(duì)數(shù)收益率的m4模型的標(biāo)準(zhǔn)化殘差樣本自相關(guān)圖（見(jiàn)圖6）。

圖6 m4模型的標(biāo)準(zhǔn)化殘差樣本ACF圖

畫(huà)出對(duì)數(shù)收益率m4模型的標(biāo)準(zhǔn)化殘差平方樣本自相關(guān)圖（見(jiàn)圖7）。

圖7 m4模型的標(biāo)準(zhǔn)化殘差平方樣本ACF圖

圖6不表明標(biāo)準(zhǔn)化殘差過(guò)程中有顯著的序列相關(guān)性，圖7不表明標(biāo)準(zhǔn)化殘差過(guò)程中有顯著的條件異方差性，因此，m4模型通過(guò)了診斷，可以用作預(yù)測(cè)。

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

應(yīng)用m4模型對(duì)T公司的對(duì)數(shù)收益率序列進(jìn)行了10期外推預(yù)測(cè)，預(yù)測(cè)結(jié)果如下：

表5 對(duì)數(shù)收益率的m4模型的預(yù)測(cè)結(jié)果

保留原對(duì)數(shù)收益率序列的后50條數(shù)據(jù)，并結(jié)合上述m4模型的預(yù)測(cè)結(jié)果，繪制出如下預(yù)測(cè)圖像（見(jiàn)圖8）。

圖8 對(duì)數(shù)收益率的m4模型的預(yù)測(cè)圖像