劉麗萍

（貴州財(cái)經(jīng)大學(xué)數(shù)學(xué)與統(tǒng)計(jì)學(xué)院，貴陽 550025）

引言

協(xié)方差陣在實(shí)際生活中有著廣泛的應(yīng)用，尤其是在資產(chǎn)組合及風(fēng)險(xiǎn)管理中，其扮演著不可替代的角色。協(xié)方差陣的估計(jì)方法大都是基于低頻數(shù)據(jù)的。近年來，高頻數(shù)據(jù)越來越容易獲得，較低頻數(shù)據(jù)而言，其包含了更加豐富的信息。高頻數(shù)據(jù)的信息量更高，但抽樣頻率提高的同時(shí)會(huì)帶來微觀噪聲影響越來越大的問題。由于短期內(nèi)的價(jià)格波動(dòng)等因素，資產(chǎn)收益率可能出現(xiàn)幅度較大的變化，數(shù)據(jù)呈現(xiàn)跳躍狀態(tài)，該跳躍對(duì)高頻數(shù)據(jù)協(xié)方差陣估計(jì)帶來更多難題。信息技術(shù)的發(fā)展促使更多高頻交易場(chǎng)景以及更高頻率的交易場(chǎng)景的出現(xiàn)。人們對(duì)高頻協(xié)方差陣的研究也隨之深入，對(duì)微觀噪聲和跳躍對(duì)其估計(jì)的影響也進(jìn)行了探索。研究指出，微觀噪聲以及跳躍存在時(shí)，積分協(xié)方差陣的一致估計(jì)量不是已實(shí)現(xiàn)協(xié)方差陣，并提出降低微觀噪聲以及跳躍影響的估計(jì)量（比如Zhang（2011）[1]、Griffin和Oomen（2011）[2]、Barndorff-Nielsen和Hansen等（2011）[3]、以及Nole和Voev（2008，2012）[4～5]）。但之前的研究往往分別研究噪聲或者跳躍的影響，劉麗萍（2013）[6]提出了修正的門限預(yù)平均已實(shí)現(xiàn)協(xié)方差陣（MTPCOV），MTPCOV估計(jì)量同時(shí)考慮了噪聲和跳躍的影響，以期彌補(bǔ)之前研究的不足。

頻率高、維度大是當(dāng)今金融機(jī)構(gòu)所研究的數(shù)據(jù)的顯著特征。這種情況下，維數(shù)詛咒導(dǎo)致一些協(xié)方差陣的估計(jì)和預(yù)測(cè)模型的適用性較差。為了克服這些問題，本文將MTPCOV和VAR-LASSO結(jié)合起來，來解決噪聲、跳躍和維數(shù)詛咒對(duì)協(xié)方差陣估計(jì)的影響。本文后續(xù)內(nèi)容如下：第一部分詳細(xì)介紹了MTPCOV估計(jì)量；第二部分將MTPCOV和VAR-LASSO模型相結(jié)合來估計(jì)和預(yù)測(cè)高維高頻的金融協(xié)方差陣；第三部分是實(shí)證分析；最后一部分是本文的結(jié)論。

一、MTPCOV估計(jì)量

修正的門限預(yù)平均已實(shí)現(xiàn)協(xié)方差陣（MTPCOV）是由劉麗萍（2013）[6]提出的，其構(gòu)造形式如下：

在上式（1）中，由馬丹和尹優(yōu)平（2012）[7]的研究可知，MTPCOV的主對(duì)角線元素MTPRV的形式為：

由劉麗萍（2013）[6]的研究可知，MTPCOV的副對(duì)角線元素MTPCV的形式為：

劉麗萍（2013）[6]的研究指出，協(xié)方差陣MTPCOV解決了維數(shù)詛咒和噪聲的影響，是積分協(xié)方差陣的一致估計(jì)量。但是，當(dāng)資產(chǎn)的維度較高時(shí)，該估計(jì)量將面臨著嚴(yán)重的維數(shù)詛咒問題，估計(jì)效果很差。為此，本文考慮將MTPCOV估計(jì)量和VAR-LASSO模型相結(jié)合，來解決維數(shù)詛咒問題。

二、考慮噪聲和跳躍影響的高維高頻協(xié)方差陣的估計(jì)和預(yù)測(cè)

前文提到的MTPCOV估計(jì)量同時(shí)考慮了噪聲和跳躍的影響，但是當(dāng)資產(chǎn)的維度較高時(shí)，該估計(jì)量將面臨著嚴(yán)重的維數(shù)詛咒問題。為此，我們考慮將MTPCOV估計(jì)量和VAR-LASSO模型相結(jié)合，首先采用MTPCOV估計(jì)量估計(jì)出高頻數(shù)據(jù)的協(xié)方差陣，將其記為∑。VAR-LASSO模型是由Callot和Kock等（2017）[8]針對(duì)高維高頻數(shù)據(jù)而提出的模型，該模型認(rèn)為，前p期的協(xié)方差陣的值會(huì)影響到當(dāng)期的值。也就是說，第t期的協(xié)方差陣∑t依賴于∑t-1，…，∑t-p。假定資產(chǎn)的維度為n，則協(xié)方差陣∑t為n×n，共有不同元素n（n+1）/2個(gè)。對(duì)采用MTPCOV估計(jì)量估計(jì)的t時(shí)期的協(xié)方差陣∑t的上三角元素進(jìn)行拉直操作，符號(hào)vech（）表示拉直操作，得到向量yt=vech（∑t）。假定yt服從VAR（p）模型，即：

式（4）中，Φi（i=1，…，p）是k×k的系數(shù)矩陣（k=n（n+1）/2）；yt-i=vech（∑t-i）。式（4）可以簡(jiǎn)化為如下形式：

式（4）中，yi=（yT，i，…，y1，i）′是觀測(cè)值i（i=1，…，k）的向量；?i=（?T，i，…，?1，i）′是由對(duì)應(yīng)的殘差構(gòu)成的向量；X=（ZT，…，Z1）′，是kp維的系數(shù)向量。根據(jù)式（5）可知，我們所建立的每一個(gè)方程，都需要估計(jì)個(gè)個(gè)參數(shù)，并且待估參數(shù)的個(gè)數(shù)會(huì)隨著資產(chǎn)產(chǎn)維度的增加而急劇地增加，資產(chǎn)維度的增高會(huì)加劇式（5）的估計(jì)困難。所以，Callot和Kock等又進(jìn)一步在VAR（p）模型中引入LASSO方法，將不顯著變量的回歸系數(shù)壓縮至0。其思想可以由式（6）表示，通過最小化式（6）來估計(jì)最優(yōu)的回歸系數(shù)。

估計(jì)出VAR-LASSO模型的回歸系數(shù)后，便可以根據(jù)其進(jìn)一步做預(yù)測(cè)，假定我們做的是滾動(dòng)的一步向前預(yù)測(cè)。則：

三、實(shí)證研究

在該部分的實(shí)證研究中，我們選擇的是2011年1月4日至2014年9月30日的上證180指數(shù)成份股的高頻數(shù)據(jù)，剔除交易缺失的數(shù)據(jù)，共有交易的天數(shù)為906天。將906天劃分為估計(jì)和預(yù)測(cè)兩個(gè)部分，估計(jì)區(qū)間的長(zhǎng)度為806天，預(yù)測(cè)區(qū)間的長(zhǎng)度為100天。采用滾動(dòng)預(yù)測(cè)法來預(yù)測(cè)后100天的協(xié)方差陣。首先采用MTPCOV和Andersen（2003）[9]提出的RCOV估計(jì)量來估計(jì)前806天的協(xié)方差陣，接著采用VAR-LASSO模型和劉麗萍等（2015）[10]提出的CF-ARMA模型來預(yù)測(cè)第807天的協(xié)方差陣；然后采用相同的協(xié)方差陣估計(jì)量來估計(jì)2—807天的協(xié)方差陣，從而預(yù)測(cè)出第808天的協(xié)方差陣；依此類推，直到選用101—905天的數(shù)據(jù)來估計(jì)協(xié)方差陣，最后預(yù)測(cè)出第906天的協(xié)方差陣。我們將預(yù)測(cè)出的100天的協(xié)方差陣應(yīng)用在投資組合中，來驗(yàn)證它的應(yīng)用效果。這里，我們采用的投資組合是等比例風(fēng)險(xiǎn)投資組合。下表給出了我們采用不同的模型預(yù)測(cè)得到的資產(chǎn)組合的收益、波動(dòng)和夏普比例值。

組合收益、波動(dòng)及Sharpe比率值

從上表可以得到三條結(jié)論：首先，無論采用的預(yù)測(cè)模型是CF-ARMA模型還是VAR-LASSO模型，由MTPCOV估計(jì)量構(gòu)造的投資組合的收益明顯高于RCOV，并且MTPCOV的組合波動(dòng)更小，Sharpe比率值更高，這是因?yàn)镸TPCOV估計(jì)量克服了噪聲和跳躍的影響。其次，預(yù)測(cè)模型VAR-LASSO明顯要優(yōu)于CF-ARMA。最后，由VAR-LASSO模型預(yù)測(cè)的MTPCOV估計(jì)量構(gòu)造的投資組合最優(yōu)，其對(duì)應(yīng)的Sharpe比率值最高，因?yàn)閂AR-LASSO模型和MTOCOV相結(jié)合來估計(jì)和預(yù)測(cè)金融數(shù)據(jù)的協(xié)方差陣時(shí)，不僅解決了噪聲和跳躍的問題，還解決了維數(shù)災(zāi)難問題。

下頁圖是VAR-LASSO和CF-ARMA模型預(yù)測(cè)的資產(chǎn)組合的邊界圖，下頁圖的結(jié)果基本與上表相符，進(jìn)一步說明考慮了噪聲和跳躍的MTPCOV估計(jì)量明顯優(yōu)于RCOV估計(jì)量，采用VAR-LASSO模型預(yù)測(cè)的MTPCOV估計(jì)量應(yīng)用在投資組合中有較好的表現(xiàn)。

組合的有效邊界圖

結(jié)語

高維高頻協(xié)方差陣的估計(jì)和預(yù)測(cè)是統(tǒng)計(jì)領(lǐng)域研究的熱點(diǎn)問題。噪聲、跳躍以及維數(shù)詛咒問題為協(xié)方差陣的估計(jì)帶來了眾多的挑戰(zhàn)。本文將MTPCOV和VAR-LASSO模型相結(jié)合來來估計(jì)和預(yù)測(cè)高維高頻協(xié)方差陣，在考慮了噪聲和跳躍對(duì)協(xié)方差陣估計(jì)影響的同時(shí)，也解決了維數(shù)詛咒問題。將VAR-LASSO模型預(yù)測(cè)的MTPCOV估計(jì)量應(yīng)用在投資組合中，發(fā)現(xiàn)與其他估計(jì)量相比，由其構(gòu)造的投資組合有更好的表現(xiàn)，收益更高，風(fēng)險(xiǎn)更小，其對(duì)應(yīng)的Sharpe比率值更高。