• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    偏方差波動(dòng)率預(yù)測(cè)模型

    2023-12-21 12:02:58陳梓榮周瑤
    上海管理科學(xué) 2023年6期

    陳梓榮 周瑤

    摘 要:論文研究了閾值數(shù)量和大小對(duì)已實(shí)現(xiàn)波動(dòng)率預(yù)測(cè)的影響,并提出了偏方差已實(shí)現(xiàn)波動(dòng)率預(yù)測(cè)模型——HAR-PV(G),該模型進(jìn)一步提高了已實(shí)現(xiàn)波動(dòng)率的預(yù)測(cè)效果??紤]到不同大小收益對(duì)已實(shí)現(xiàn)波動(dòng)率的影響具有非對(duì)稱性,以HAR-RS模型為基礎(chǔ),選取不同數(shù)量和不同大小的閾值組合對(duì)日內(nèi)收益進(jìn)行分割,并計(jì)算對(duì)應(yīng)的偏方差,從而構(gòu)建HAR-PV(G)模型。論文以滬深300指數(shù)為研究對(duì)象,比較了不同HAR-PV(G)模型的樣本外預(yù)測(cè)能力。樣本外分析表明,閾值數(shù)量為3的平分偏方差模型具有比傳統(tǒng)HAR、HAR-RS以及其他閾值組合的偏方差模型更好的預(yù)測(cè)能力。全樣本的參數(shù)分析也顯示閾值數(shù)量為3的平分偏方差模型對(duì)數(shù)據(jù)的擬合效果更出眾。

    關(guān)鍵詞:偏方差;已實(shí)現(xiàn)波動(dòng)率;日內(nèi)收益;高頻數(shù)據(jù)

    中圖分類號(hào):F 830.9

    文獻(xiàn)標(biāo)志碼:A

    Partial Variance Volatility Forecasting Model

    CHEN Zirong ZHOU Yao

    (Antai College of Economic and Management, Shanghai Jiao Tong University, Shanghai 200030, China)

    Abstract:This paper proposes a partial variance prediction model, HAR-PV(G), which can test the effects of different quantities and different sizes of thresholds on the prediction of realized volatility, aiming at improving the predicting effect of realized volatility. Due to the asymmetric effects of different intraday returns on the realized volatility, this paper selects threshold combinations of different quantities and sizes to segment intraday returns based on HAR-RS model. Then, this paper calculates the corresponding partial variance to construct HAR-PV(G) model. This paper uses 5-minute high-frequency trading data of CSI 300 index to compare in-sample and out-of-sample performances of HAR-PV(G) model. The results show an equal division of partial variance model with 3 thresholds could achieve a better out-of-sample performance compared with traditional HAR, HAR-RS and other HAR-PV(G) models with different threshold combinations. Meanwhile, this equal division of partial variance model with 3 thresholds also has an excellent fitting effect in the in-sample analysis.

    Key words:partial variance; realized volatility; intraday return; high-frequency data

    0 引言

    金融市場(chǎng)中,資產(chǎn)具有波動(dòng)性是一種被普遍認(rèn)可的金融事實(shí)。用于衡量資產(chǎn)風(fēng)險(xiǎn)情況的波動(dòng)率在投資組合構(gòu)建、資產(chǎn)定價(jià)以及風(fēng)險(xiǎn)管理等金融領(lǐng)域扮演著極其重要的角色,因此對(duì)波動(dòng)率的統(tǒng)計(jì)分析和建模預(yù)測(cè)一直都是金融計(jì)量領(lǐng)域關(guān)注的熱點(diǎn)問(wèn)題。隨著我國(guó)衍生品市場(chǎng)快速發(fā)展,市場(chǎng)對(duì)衍生品定價(jià)的準(zhǔn)確性要求強(qiáng)調(diào)了波動(dòng)率預(yù)測(cè)研究的重要性。此外,宏觀政策層面的波動(dòng)必然帶來(lái)微觀個(gè)體的變動(dòng),進(jìn)而引起金融市場(chǎng)的波動(dòng);微觀市場(chǎng)層面的波動(dòng)與金融環(huán)境和實(shí)體經(jīng)濟(jì)的穩(wěn)定密切相關(guān),正確預(yù)測(cè)波動(dòng)率對(duì)市場(chǎng)從業(yè)人員和監(jiān)管者都有著不可忽視的作用。

    波動(dòng)率通常具有較為復(fù)雜的特性,Bollerslev和Taylor提出的廣義自回歸條件異方差(GARCH) 模型和隨機(jī)波動(dòng)率(SV) 模型是傳統(tǒng)的兩類刻畫(huà)波動(dòng)率的模型。隨著數(shù)據(jù)存儲(chǔ)技術(shù)的進(jìn)步,日內(nèi)高頻數(shù)據(jù)的獲取成本也大幅降低,因此包含豐富信息的日內(nèi)高頻數(shù)據(jù)為波動(dòng)率的估計(jì)和預(yù)測(cè)提供了新的手段。Andersen等首次提出將已實(shí)現(xiàn)波動(dòng)率(Realized Volatility, RV)作為日內(nèi)波動(dòng)率的估計(jì)量,相較于傳統(tǒng)的GARCH和SV等模型,已實(shí)現(xiàn)波動(dòng)率不僅不需要復(fù)雜的模型計(jì)算和數(shù)理推導(dǎo),而且能夠充分利用日內(nèi)高頻交易數(shù)據(jù)所包含的信息,更精確地刻畫(huà)金融市場(chǎng)波動(dòng)率,因此已實(shí)現(xiàn)波動(dòng)率一經(jīng)提出便得到了廣泛的研究和應(yīng)用。

    在已實(shí)現(xiàn)波動(dòng)率預(yù)測(cè)的研究中,最具代表性的是Corsi提出的HAR模型。Corsi以“異質(zhì)市場(chǎng)假說(shuō)”為理論基礎(chǔ),將波動(dòng)率劃分為短期波動(dòng)、中期波動(dòng)和長(zhǎng)期波動(dòng),提出了異質(zhì)自回歸已實(shí)現(xiàn)波動(dòng)率預(yù)測(cè)模型(Heterogeneous Autoregressive,HAR)。HAR模型不僅降低了傳統(tǒng)模型的估計(jì)復(fù)雜度,也具有很強(qiáng)的擴(kuò)展性和明確的經(jīng)濟(jì)含義,因此逐漸成為已實(shí)現(xiàn)波動(dòng)率預(yù)測(cè)的基準(zhǔn)模型。Andersen等將已實(shí)現(xiàn)波動(dòng)率分解成連續(xù)波動(dòng)(Continuous Volatility, CV)和跳躍波動(dòng)(Jump Volatility,JV),并分別考慮這兩種波動(dòng)對(duì)未來(lái)已實(shí)現(xiàn)波動(dòng)率的影響,提出了HAR-RV-J和HAR-CJ模型。作者通過(guò)實(shí)證發(fā)現(xiàn)分解波動(dòng)率可以顯著提高已實(shí)現(xiàn)波動(dòng)率的預(yù)測(cè)能力。Chen和Ghysels發(fā)現(xiàn)好消息和壞消息對(duì)波動(dòng)率的影響是非對(duì)稱的,因此把已實(shí)現(xiàn)波動(dòng)率分解為好消息驅(qū)動(dòng)的波動(dòng)率成分和壞消息驅(qū)動(dòng)的波動(dòng)率成分可以提高已實(shí)現(xiàn)波動(dòng)率的預(yù)測(cè)能力。遵循這一思路,Patton和Sheppard將日內(nèi)收益按照正負(fù)進(jìn)行劃分,分別平方求和得到已實(shí)現(xiàn)半變方差(Realized Semi-Variance,RS)和符號(hào)跳躍變差(Signed Jump Variation, SJV),在此基礎(chǔ)上提出了HAR-RS、HAR-RV-SJV和HAR-RV-SJVd等模型,實(shí)證發(fā)現(xiàn)這三種模型的預(yù)測(cè)能力都優(yōu)于傳統(tǒng)的HAR模型。

    國(guó)內(nèi)對(duì)已實(shí)現(xiàn)波動(dòng)率預(yù)測(cè)的研究起步較晚。文鳳華、劉曉群、唐海如等在HAR模型的基礎(chǔ)上考慮市場(chǎng)波動(dòng)的杠桿效應(yīng)和量?jī)r(jià)關(guān)系,提出了LHAR-RV-V模型,并使用滬深300指數(shù)的高頻數(shù)據(jù)進(jìn)行實(shí)證分析,證明該模型可以提高傳統(tǒng)HAR模型的預(yù)測(cè)能力。馬鋒、魏宇、黃登仕等將符號(hào)跳躍變差應(yīng)用到中國(guó)市場(chǎng),比較了HAR、HAR-RV-J、HAR-RV-CJ和HAR-RV-TCJ等已實(shí)現(xiàn)波動(dòng)率預(yù)測(cè)模型在加入符號(hào)變差后的預(yù)測(cè)能力。陳聲利、關(guān)濤和李一軍將百度指數(shù)引入HAR模型中,實(shí)證發(fā)現(xiàn)新模型可以提高對(duì)滬深300股指期貨波動(dòng)率的預(yù)測(cè)能力。李俊儒、汪壽陽(yáng)和魏云捷考慮波動(dòng)率 的測(cè)量誤差,發(fā)現(xiàn)加入波動(dòng)率測(cè)量誤差后的HAR模型預(yù)測(cè)能力的持續(xù)性有所提高,且其預(yù)測(cè)效果優(yōu)于傳統(tǒng)的HAR模型。龔旭、曹杰、文鳳華等在4個(gè)經(jīng)典的HAR族模型的基礎(chǔ)上,探究了杠桿效應(yīng)和結(jié)構(gòu)突變因素對(duì)已實(shí)現(xiàn)波動(dòng)率的預(yù)測(cè)能力。

    盡管上述HAR族模型對(duì)已實(shí)現(xiàn)波動(dòng)率已有較強(qiáng)的預(yù)測(cè)能力,但是考慮到波動(dòng)率在金融市場(chǎng)的重要作用,進(jìn)一步探究如何提高已實(shí)現(xiàn)波動(dòng)率的預(yù)測(cè)能力仍有必要。因此,本文基于Patton和Sheppard提出的HAR-RS模型,提出偏方差已實(shí)現(xiàn)波動(dòng)率預(yù)測(cè)模型——HAR-PV(G),并研究了不同閾值組合對(duì)已實(shí)現(xiàn)波動(dòng)率模型預(yù)測(cè)能力的影響。

    在HAR-RS模型中,作者對(duì)日內(nèi)收益以0作為閾值進(jìn)行分割,將已實(shí)現(xiàn)波動(dòng)率分為上行波動(dòng)率和下行波動(dòng)率并加入HAR模型,最終發(fā)現(xiàn)HAR-RS模型比傳統(tǒng)的HAR模型有更強(qiáng)的預(yù)測(cè)能力。基于HAR-RS的研究,本文進(jìn)一步探究以下兩個(gè)具體問(wèn)題:第一,在中國(guó)金融市場(chǎng),選取非0的閾值對(duì)日內(nèi)收益進(jìn)行分割,然后計(jì)算對(duì)應(yīng)的波動(dòng)率成分且加入HAR模型中,新的模型是否比HAR-RS模型有更好的預(yù)測(cè)能力;第二,如果閾值數(shù)量大于1,即將日內(nèi)收益率分成三段或者更多,然后計(jì)算對(duì)應(yīng)的波動(dòng)率成分并加入HAR模型中,新的模型對(duì)已實(shí)現(xiàn)波動(dòng)率的預(yù)測(cè)能力是否會(huì)得到進(jìn)一步提升。

    1 計(jì)量模型

    1.1 已實(shí)現(xiàn)波動(dòng)率預(yù)測(cè)模型

    假設(shè)金融資產(chǎn)在交易日t的對(duì)數(shù)價(jià)格p=log (P),服從以下擴(kuò)散模型:

    p=∫ μ ds+∫ σ dW+J.??? (1)

    其中:μ是有限方差的連續(xù)跳躍過(guò)程;σ代表瞬時(shí)波動(dòng)率,且是一個(gè)嚴(yán)格為正的càdlaàg過(guò)程;J是跳躍過(guò)程。如果用Δp=pp代表s時(shí)點(diǎn)的跳躍大小,那么p的二次變差為[p,p]=∫ σ ds+∑(Δp) 。其中:∫ σ ds為積分波動(dòng)率,又叫作連續(xù)波動(dòng)率;∑(Δp)為跳躍波動(dòng)。

    假設(shè)在交易日t,擁有間隔相等的M+1個(gè)觀察值p,p,…,p。令r代表交易日t內(nèi)第i時(shí)段的對(duì)數(shù)收益率,r=pp,i=1,…,M。Andersen等提出可以使用日內(nèi)收益率的平方和作為二次變差的估計(jì)量,該估計(jì)量又被叫作已實(shí)現(xiàn)波動(dòng)率,其計(jì)算方法如下:

    在眾多已實(shí)現(xiàn)波動(dòng)率預(yù)測(cè)模型中,Corsi等提出的HAR模型因能夠刻畫(huà)已實(shí)現(xiàn)波動(dòng)率的長(zhǎng)尾效應(yīng)獲得最為廣泛的應(yīng)用。HAR模型的形式非常簡(jiǎn)潔,只需使用上一天、上一周以及上一個(gè)月的平均已實(shí)現(xiàn)波動(dòng)率作為預(yù)測(cè)下一期已實(shí)現(xiàn)波動(dòng)率的預(yù)測(cè)變量:

    其中,RV是t-4到t日的平均已實(shí)現(xiàn)波動(dòng)率,RV是t-19到t日的平均已實(shí)現(xiàn)波動(dòng)率。

    Patton和Sheppard考慮了波動(dòng)率的杠桿效應(yīng),提出以0為閾值對(duì)日內(nèi)收益進(jìn)行分割,并將已實(shí)現(xiàn)波動(dòng)率分解成上行已實(shí)現(xiàn)半方差RS=∑ r I{r>0}和下行已實(shí)現(xiàn)半方差RS= ∑ r I{r<0},并在此基礎(chǔ)之上提出了HAR-RS模型:

    在Patton和Sheppard、馬鋒、魏宇、黃登仕等的研究中,下行已實(shí)現(xiàn)波動(dòng)率和上行已實(shí)現(xiàn)波動(dòng)率對(duì)未來(lái)已實(shí)現(xiàn)波動(dòng)率的影響具有顯著的非對(duì)稱性,下行已實(shí)現(xiàn)波動(dòng)率的影響更為顯著,且HAR-RS模型擁有比HAR模型更好的預(yù)測(cè)能力。

    1.2 偏方差波動(dòng)率預(yù)測(cè)模型

    HAR-RS的研究表明正負(fù)收益對(duì)未來(lái)已實(shí)現(xiàn)波動(dòng)率的影響是不同的。既然不同大小的收益對(duì)未來(lái)已實(shí)現(xiàn)波動(dòng)率的影響是不同的,那么是否可以通過(guò)對(duì)日內(nèi)收益進(jìn)行更細(xì)致的劃分來(lái)獲得更好的預(yù)測(cè)效果將是本文試圖回答的問(wèn)題。

    在Patton和Sheppard的研究中,周(RV)和月(RV)已實(shí)現(xiàn)波動(dòng)率的分解對(duì)預(yù)測(cè)效果的影響較小,因此本文的研究集中在日已實(shí)現(xiàn)波動(dòng)率(RV)的分解對(duì)波動(dòng)率預(yù)測(cè)的影響上。

    從公式(5)中可以發(fā)現(xiàn),偏方差預(yù)測(cè)模型中閾值的選擇面臨兩個(gè)問(wèn)題:閾值數(shù)量和閾值大小。為了探究如何選取合適的閾值組合進(jìn)行預(yù)測(cè),本文擬從實(shí)證的角度來(lái)研究閾值的大小和數(shù)量對(duì)波動(dòng)率預(yù)測(cè)的影響。本文將通過(guò)比較不同閾值組合的樣本外表現(xiàn),最終對(duì)閾值組合的選擇提出建議。

    2 數(shù)據(jù)

    本文選取滬深300股票指數(shù)高頻交易數(shù)據(jù)作為研究對(duì)象,數(shù)據(jù)來(lái)源于微盛投資數(shù)據(jù)庫(kù)。為了盡可能多地運(yùn)用高頻數(shù)據(jù)包含的信息,同時(shí)避免因?yàn)槭袌?chǎng)結(jié)構(gòu)噪聲帶來(lái)的影響,本文采用5分鐘頻率的日內(nèi)收益數(shù)據(jù)。數(shù)據(jù)樣本時(shí)間跨度為2011年1月4日至2021年5月17日(數(shù)據(jù)庫(kù)時(shí)間截至2021年5月17日)。

    從圖1的已實(shí)現(xiàn)波動(dòng)率時(shí)序圖來(lái)看,波動(dòng)率在2015年出現(xiàn)了一個(gè)高峰時(shí)期,其正好對(duì)應(yīng)了2015年的股災(zāi)時(shí)期,符合金融事實(shí)。此外,從表1的描述性統(tǒng)計(jì)來(lái)看,已實(shí)現(xiàn)波動(dòng)率的偏度和峰度都很高,說(shuō)明已實(shí)現(xiàn)波動(dòng)率存在右偏和尖峰厚尾等特征。

    3 實(shí)證分析

    為了進(jìn)一步深入研究閾值組合對(duì)已實(shí)現(xiàn)波動(dòng)率預(yù)測(cè)的影響,本節(jié)將從以下四個(gè)角度具體展開(kāi)研究:

    1)閾值數(shù)量固定時(shí),不同閾值組合下的偏方差模型的預(yù)測(cè)能力比較。

    2)放開(kāi)閾值數(shù)量固定的限制,不同閾值數(shù)量的偏方差模型的預(yù)測(cè)能力比較。

    3)放開(kāi)閾值數(shù)量和組合的限制,允許閾值的數(shù)量和大小在每一期都發(fā)生變化,并和其他偏方差預(yù)測(cè)模型的預(yù)測(cè)能力進(jìn)行比較。

    4)通過(guò)全樣本的參數(shù)估計(jì)分析不同大小收益對(duì)未來(lái)已實(shí)現(xiàn)波動(dòng)率的影響。

    前三個(gè)分析使用樣本外表現(xiàn)來(lái)比較不同模型的預(yù)測(cè)能力,第四個(gè)分析使用樣本內(nèi)分析探究不同大小收益對(duì)未來(lái)波動(dòng)率的影響。

    在實(shí)際選取閾值時(shí),考慮到每天的已實(shí)現(xiàn)波動(dòng)率的變化是不同的且波動(dòng)較大,本文并不直接選取固定的數(shù)值作為閾值,而是選取一些特定的分位點(diǎn)作為閾值:

    其中,Q(z;q)是z序列的q分位點(diǎn)。同時(shí),本文將用HAR-PV()代表當(dāng)閾值數(shù)量為G且分位點(diǎn)將日內(nèi)收益率進(jìn)行平分的偏方差預(yù)測(cè)模型。

    本文主要通過(guò)樣本外分析對(duì)不同偏方差模型的預(yù)測(cè)能力進(jìn)行比較。樣本外預(yù)測(cè)采用一步向前滾動(dòng)窗估計(jì),估計(jì)樣本窗口設(shè)定為1000個(gè)交易日,即預(yù)測(cè)從第1001個(gè)交易日開(kāi)始。由于文章篇幅的限制,本文的實(shí)證分析主要集中于預(yù)測(cè)下一天的已實(shí)現(xiàn)波動(dòng)率RV。

    對(duì)于樣本外預(yù)測(cè)能力的評(píng)判標(biāo)準(zhǔn),現(xiàn)有研究并未對(duì)損失函數(shù)的選擇達(dá)成共識(shí)。在傳統(tǒng)的HAR模型研究中,不少文獻(xiàn)同時(shí)采用多個(gè)損失函數(shù)來(lái)比較不同模型的樣本外表現(xiàn)。本文遵循Audrino和Hu的做法,選取兩個(gè)常用的損失函數(shù):

    其中,RV是已實(shí)現(xiàn)波動(dòng)率的預(yù)測(cè)值,RV是已實(shí)現(xiàn)波動(dòng)率的真實(shí)值,n代表預(yù)測(cè)樣本數(shù)量。預(yù)測(cè)模型的損失函數(shù)值越小,表明模型的預(yù)測(cè)能力越強(qiáng)。

    然而,通過(guò)比較不同損失函數(shù)值的大小可以直接比較不同模型預(yù)測(cè)的精度,但是無(wú)法判斷一個(gè)模型的預(yù)測(cè)能力是否顯著優(yōu)于另一個(gè)或另一組模型,因此本文還將采用模型可信集(MCS)檢驗(yàn)方法來(lái)對(duì)比不同模型的預(yù)測(cè)能力。模型可信集檢驗(yàn)是由Hansen等基于Hansen改進(jìn)之后提出的用于對(duì)不同模型的預(yù)測(cè)能力優(yōu)劣進(jìn)行排序的檢驗(yàn)方法。MCS檢驗(yàn)統(tǒng)計(jì)量的構(gòu)造和比較過(guò)程非常復(fù)雜,由于篇幅原因在此不過(guò)多展開(kāi),詳細(xì)步驟可參考Hansen等。

    在后續(xù)的結(jié)果里,本文將重點(diǎn)報(bào)告MCS檢驗(yàn)中的P值,根據(jù)Hansen等論文里的解釋,該P(yáng)值可以簡(jiǎn)單理解為在所有進(jìn)行檢測(cè)的模型中,該模型擁有最好預(yù)測(cè)能力的概率。本文的MCS檢驗(yàn)采用較為常用的范圍統(tǒng)計(jì)量(Range)作為檢驗(yàn)統(tǒng)計(jì)量,窗口長(zhǎng)度為25,重復(fù)抽樣次數(shù)為5000,臨界值的選擇為0.1(置信度為90%)。實(shí)證中如果模型的Range統(tǒng)計(jì)量對(duì)應(yīng)的檢驗(yàn)P值大于0.1,說(shuō)明該模型具有較好的預(yù)測(cè)能力;P值越大,表明該模型具有較好預(yù)測(cè)能力的概率越高。

    3.1 固定數(shù)量不同閾值組合的預(yù)測(cè)效果比較

    本節(jié)先從比較簡(jiǎn)單的情形入手,比較閾值數(shù)量固定的閾值組合對(duì)下一期已實(shí)現(xiàn)波動(dòng)率的預(yù)測(cè)效果。在下文的分析中,閾值分位點(diǎn)的選擇區(qū)間為0.05至0.95,步長(zhǎng)為0.05。

    首先,在給定一個(gè)閾值的情況下,圖2繪制了該情形下的樣本外表現(xiàn)。其中,橫軸代表了不同的閾值分位點(diǎn),縱軸代表了對(duì)應(yīng)模型的樣本外平均MSE的值。兩條虛線HAR-RV和HAR-RS分別代表了傳統(tǒng)的HAR-RV模型和HAR-RS模型的樣本外平均MSE的值。從圖2中可以發(fā)現(xiàn),HAR-PV(1)模型在閾值分位點(diǎn)處于0.5附近時(shí)取到最小值,即是說(shuō),已實(shí)現(xiàn)波動(dòng)率處在中值水平時(shí)取得最小值,此時(shí)已實(shí)現(xiàn)波動(dòng)率的中值非常接近于0,因此閾值分位點(diǎn)為0.5的偏方差模型非常接近于HAR-RS模型。該結(jié)果表明當(dāng)閾值數(shù)量只有1個(gè)的時(shí)候,在樣本外表現(xiàn)上擊敗HAR-RS模型是非常困難的。

    此外,從圖2所顯示的兩側(cè)尾部的表現(xiàn)來(lái)看,當(dāng)閾值分位點(diǎn)非常偏左時(shí),預(yù)測(cè)的樣本外效果表現(xiàn)并不理想,不僅弱于傳統(tǒng)的HAR-RV模型,也遠(yuǎn)比不上在右側(cè)極端分位點(diǎn)上進(jìn)行分割后的表現(xiàn),該結(jié)果類似于Bollerslev等對(duì)美國(guó)S&P500的研究結(jié)果,但并不符合通常的認(rèn)知。在傳統(tǒng)的金融研究里,一般認(rèn)為左側(cè)的極端收益對(duì)未來(lái)波動(dòng)率的影響相較于右側(cè)收益更大,因此單獨(dú)考慮左側(cè)極端風(fēng)險(xiǎn)通常能夠提高對(duì)未來(lái)波動(dòng)率的預(yù)測(cè)能力?;诖耍疚膶⒃?.4節(jié)中進(jìn)一步深入探討不同大小的收益對(duì)未來(lái)波動(dòng)率的影響。

    接著,本文繼續(xù)考慮兩個(gè)閾值的情況,這里閾值分位點(diǎn)同樣選擇從0.05至0.95,步長(zhǎng)為0.05。這里有兩個(gè)閾值,因此共有171個(gè)不同的閾值組合。

    圖3刻畫(huà)了給定兩個(gè)閾值下的樣本外表現(xiàn)。兩個(gè)橫軸分別代表了兩個(gè)閾值的選擇,縱軸代表了不同模型樣本外的平均MSE值。從圖3中可以發(fā)現(xiàn),圖形中存在兩個(gè)相對(duì)低點(diǎn),而低點(diǎn)的位置非常接近于[0.35,0.65]的閾值組合,說(shuō)明 HAR-PV(2)模型在閾值組合為[0.35,0.65]附近的樣本外表現(xiàn)相對(duì)較好。

    綜合圖3和圖4可知,在閾值數(shù)量固定為1或者2的兩種情況下,具有最好的樣本外表現(xiàn)的模型都接近于對(duì)日內(nèi)收益進(jìn)行平分后的偏方差模型,即HAR-PV()模型。因此,本文進(jìn)一步推測(cè)當(dāng)閾值數(shù)量固定時(shí),平分日內(nèi)收益率的偏方差模型具有相對(duì)較好的預(yù)測(cè)能力。

    為了進(jìn)一步驗(yàn)證這一推測(cè),下文繼續(xù)考慮給定閾值數(shù)量更多的情況。這一部分一共對(duì)5種不同閾值數(shù)量的情況進(jìn)行考慮,即閾值數(shù)量為1、2、3、4、5的情況。當(dāng)閾值數(shù)量大于2時(shí),其樣本外的平均MSE的值的分布難以通過(guò)繪圖的方式表現(xiàn)出來(lái),為了讓閾值組合的結(jié)果可視化,在每個(gè)閾值數(shù)量下選取了三組有代表性的閾值組合,表2給出了具體的閾值選擇的情況。從表2中不難發(fā)現(xiàn),組合1是對(duì)日內(nèi)收益平分的閾值組合;當(dāng)閾值數(shù)量大于等于2時(shí),組合2側(cè)重考慮兩側(cè)尾部風(fēng)險(xiǎn);組合3側(cè)重考慮左側(cè)尾部風(fēng)險(xiǎn)。不對(duì)右側(cè)尾部風(fēng)險(xiǎn)單獨(dú)考慮,主要是因?yàn)樵趥鹘y(tǒng)的金融風(fēng)險(xiǎn)研究中,左側(cè)風(fēng)險(xiǎn)對(duì)波動(dòng)率預(yù)測(cè)、資產(chǎn)配置等有更為顯著的影響。本節(jié)將分別計(jì)算每一種閾值分位點(diǎn)組合的樣本外平均MSE和QLIKE的值,并對(duì)每一個(gè)閾值數(shù)量確定下的三種組合進(jìn)行MCS檢驗(yàn),進(jìn)而判斷哪一種閾值組合是較好的預(yù)測(cè)模型。

    表3報(bào)告了不同組合下的樣本外平均MSE和QLIKE的值,其中MSE列的值乘以了10。MCS列報(bào)告了置信水平為90%的MCS檢驗(yàn)的P值,N代表了閾值的數(shù)量。從表3的結(jié)果可以看出,在絕大多數(shù)情況下,無(wú)論損失函數(shù)選擇MSE 還是QLIKE,平分閾值都有著最小的平均損失值。MCS的結(jié)果也更進(jìn)一步驗(yàn)證了這一點(diǎn),除去閾值數(shù)量是4個(gè),且損失函數(shù)采用QLIKE的情況,其他情況下,平分閾值模型的MCS檢驗(yàn)下的P值都是1.000,該結(jié)果說(shuō)明平分閾值模型是三種模型中預(yù)測(cè)能力相對(duì)較好的模型。

    3.2 不同數(shù)量閾值的預(yù)測(cè)效果比較

    上節(jié)比較了在閾值數(shù)量固定時(shí),不同閾值組合的樣本外表現(xiàn)。本節(jié)將比較不同閾值數(shù)量的樣本外表現(xiàn)。為了方便比較,本節(jié)直接應(yīng)用3.1中觀察到的結(jié)果,直接對(duì)比不同數(shù)量閾值下的平分閾值組合的偏方差模型的樣本外表現(xiàn)。同樣閾值數(shù)量的選擇為1到5,每種閾值數(shù)量下都采用平分日內(nèi)收益率的方法來(lái)計(jì)算對(duì)應(yīng)的偏方差波動(dòng)率,即HAR-PV()模型。

    圖4繪制了不同閾值數(shù)量下的樣本外表現(xiàn)。其中,橫軸代表了閾值的數(shù)量,縱軸代表了預(yù)測(cè)模型的樣本外的平均MSE的值,兩條虛線HAR-RV和HAR-RS,分別代表了傳統(tǒng)HAR-RV和HAR-RS模型的樣本外平均MSE的值。從圖4中可以清晰地看到,樣本外表現(xiàn)和閾值數(shù)量呈現(xiàn)一個(gè)正“U”形,樣本外表現(xiàn)在閾值數(shù)量取到3的時(shí)候取到最小值,MSE為3.9E-08。與傳統(tǒng)的HAR-RV模型相比,所有對(duì)日內(nèi)波動(dòng)率進(jìn)行平均分割的偏方差模型的樣本外表現(xiàn)都遠(yuǎn)好于傳統(tǒng)的HAR-RV模型。由此可見(jiàn),對(duì)日內(nèi)收益率的分割是必要的,這也足以說(shuō)明偏方差模型可以顯著提高對(duì)已實(shí)現(xiàn)波動(dòng)率的預(yù)測(cè)能力。

    為了進(jìn)一步驗(yàn)證上述結(jié)論,表4統(tǒng)計(jì)了不同閾值數(shù)量下不同模型的樣本外表現(xiàn)。從表中不難發(fā)現(xiàn),不管損失函數(shù)采用MSE還是QLIKE,HAR-PV(3)的樣本外損失函數(shù)值都是最好的之一,這也與圖5一致。此外,在MCS檢驗(yàn)中,HAR-PV(3)的MCS檢驗(yàn)P值都是1.000,說(shuō)明HAR-PV(3)模型是擁有最好預(yù)測(cè)能力的模型集合中的一個(gè)。綜上所述,在不固定閾值數(shù)量的不同組合中,三個(gè)閾值數(shù)量的偏方差模型具有較好的波動(dòng)率預(yù)測(cè)效果。

    3.3 閾值動(dòng)態(tài)模型

    從上述兩節(jié)的分析中可知,閾值的數(shù)量和大小都會(huì)對(duì)已實(shí)現(xiàn)波動(dòng)率的預(yù)測(cè)產(chǎn)生影響。本節(jié)將繼續(xù)探究一個(gè)問(wèn)題:如果允許閾值的數(shù)量和大小都隨時(shí)間的變化而變化,或者稱為閾值動(dòng)態(tài)模型,該模型是否具有更好的波動(dòng)率預(yù)測(cè)能力。本文將用HAR-PV()代表閾值數(shù)量和組合均可隨時(shí)間變化而變化的閾值動(dòng)態(tài)模型。

    本節(jié)通過(guò)交叉驗(yàn)證的方法來(lái)確定每一期的閾值組合;具體方法操作如下:

    首先,選取1000個(gè)交易日的數(shù)據(jù)作為訓(xùn)練集。具體地,假設(shè)在交易日t,選取t-1000到t日的交易數(shù)據(jù)作為訓(xùn)練集。

    其次,從訓(xùn)練集中隨機(jī)選取600個(gè)連續(xù)交易日。

    接著,將隨機(jī)選取的600個(gè)連續(xù)交易日分為400和200兩個(gè)區(qū)間。

    然后,對(duì)于每一種閾值組合,利用上文提到的滾動(dòng)預(yù)測(cè)的方差預(yù)測(cè)200個(gè)交易日的已實(shí)現(xiàn)波動(dòng)率。

    隨后,根據(jù)預(yù)測(cè)的已實(shí)現(xiàn)波動(dòng)率計(jì)算每一種閾值組合的樣本外MSE值,作為評(píng)估閾值組合的標(biāo)準(zhǔn)。

    重復(fù)步驟2至步驟5五次。

    最后,選取5次平均MSE值最小的閾值組合對(duì)應(yīng)的偏方差模型作為預(yù)測(cè)t+1日的偏方差模型。

    為了節(jié)省計(jì)算能力且不失一般性,本節(jié)的閾值分位點(diǎn)選取為0.1至0.9,步長(zhǎng)為0.1。

    表5報(bào)告了HAR-PV()閾值動(dòng)態(tài)模型的樣本外表現(xiàn)。這里將HAR-PV()模型和其他HAR-PV()模型進(jìn)行了比較,便于發(fā)現(xiàn)最優(yōu)的波動(dòng)率預(yù)測(cè)模型。首先,將HAR-PV()模型和傳統(tǒng)的HAR模型的樣本外表現(xiàn)進(jìn)行比較,HAR-PV()的MSE值為4.061,低于HAR模型的4.693,但其QLIKE值為0.201,略高于HAR的0.195,因此,結(jié)合兩種損失函數(shù)的結(jié)果,并不能得出HAR-PV()完全優(yōu)于傳統(tǒng)的HAR模型。接著,將HAR-PV()模型和前兩節(jié)得到論證的日內(nèi)收益率平均模型HAR-PV()進(jìn)行比較,類似地,并不能說(shuō)明HAR-PV()比HAR-PV(3)有更好的樣本外表現(xiàn)。此外,雖然HAR-PV()的MCS檢驗(yàn)P值均為1.000,但是從平均損失值的角度來(lái)看,并不能認(rèn)為HAR-PV()是比HAR-PV(3)更好的預(yù)測(cè)模型。

    綜合3.1、3.2和3.3的分析,可以得到以下三個(gè)結(jié)論:第一,在固定閾值數(shù)量的情況下,平分日內(nèi)收益的偏方差模型往往具有更好的預(yù)測(cè)能力。第二,偏方差模型的預(yù)測(cè)能力和閾值數(shù)量呈現(xiàn)一個(gè)倒“U”形的關(guān)系,模型的預(yù)測(cè)能力在閾值數(shù)量為3附近達(dá)到最優(yōu)。第三,即便放開(kāi)閾值數(shù)量和組合的限制,閾值動(dòng)態(tài)模型的預(yù)測(cè)能力并未優(yōu)于HAR-PV(3)。綜上所述,閾值數(shù)量為3的平分偏方差模型可以擊敗傳統(tǒng)的HAR-RS模型和其他偏方差模型,具有更好的預(yù)測(cè)能力。

    3.4 參數(shù)分析

    正如前三節(jié)所分析的,日內(nèi)收益率的狀況會(huì)對(duì)波動(dòng)率產(chǎn)生影響,為了進(jìn)一步理解不同大小的收益率對(duì)下一期已實(shí)現(xiàn)波動(dòng)率的影響,本節(jié)通過(guò)分析模型的回歸系數(shù)研究不同的偏方差對(duì)下一期已實(shí)現(xiàn)波動(dòng)率的具體影響。

    同樣,本節(jié)采用HAR-PV()模型作為研究對(duì)象。因?yàn)椴煌钠讲钕嗖钶^大,本文對(duì)所有的回歸變量進(jìn)行歸一處理,從而可以更好地在不同變量的參數(shù)之間進(jìn)行比較。與上文一致,本節(jié)依舊只考慮預(yù)測(cè)下一天的已實(shí)現(xiàn)波動(dòng)率。

    表6報(bào)告了全樣本內(nèi)不同HAR-PV()模型的估計(jì)結(jié)果,括號(hào)中是對(duì)應(yīng)參數(shù)的Newey-West調(diào)整后的t檢驗(yàn)值(West),表中的PV(i)對(duì)應(yīng)了不同模型下的第i段偏方差的回歸系數(shù)。因?yàn)榉治霾簧婕俺?shù)項(xiàng),因此表中忽略了常數(shù)項(xiàng)。

    首先,關(guān)注HAR-PV(1)模型的結(jié)果。在HAR-PV(1)模型中,負(fù)的偏方差的系數(shù)為0.519,而對(duì)應(yīng)的正偏方差只有-0.154,且負(fù)偏方差的t值達(dá)到了4.332,遠(yuǎn)高于正偏方差的2.190??梢?jiàn)在預(yù)測(cè)下一天的已實(shí)現(xiàn)波動(dòng)率時(shí),負(fù)日內(nèi)收益包含了比正日內(nèi)收益更多的信息,且其對(duì)下一期已實(shí)現(xiàn)波動(dòng)率的影響幾乎是正偏方差的四倍。該結(jié)果進(jìn)一步驗(yàn)證了龔旭、曹杰、文鳳華等的研究結(jié)果,即中國(guó)股票市場(chǎng)的波動(dòng)率具有杠桿效應(yīng)。

    HAR-PV(2)模型類似于Andersen等提出的HAR-J模型。HAR-J模型將已實(shí)現(xiàn)波動(dòng)率分解為連續(xù)波動(dòng)和跳躍波動(dòng),不同于HAR-J模型中考慮極端收益的影響,HAR-PV(2)模型重點(diǎn)考慮兩側(cè)較大的尾部收益的影響。HAR-PV(2)模型的回歸結(jié)果顯示,三段偏方差都是顯著的。右側(cè)尾部收益的影響和左側(cè)尾部收益的影響相反,左側(cè)尾部收益對(duì)下一期已實(shí)現(xiàn)波動(dòng)率的影響為正,而右側(cè)影響為負(fù)。

    從HAR-PV(3)至HAR-PV(5)的參數(shù)估計(jì)中可以總結(jié)出兩個(gè)主要結(jié)論。第一,中位偏方差在預(yù)測(cè)波動(dòng)率中占據(jù)主導(dǎo)地位。從估計(jì)參數(shù)的符號(hào)和顯著性情況可以看到,在零附近的偏方差的系數(shù)始終是正向且顯著的,該結(jié)果說(shuō)明常規(guī)收益對(duì)下一期已實(shí)現(xiàn)波動(dòng)率的影響是要大于極端收益的。第二,左側(cè)極端收益對(duì)下一期已實(shí)現(xiàn)波動(dòng)率的影響通常高于右側(cè)極端收益,但右側(cè)較大收益對(duì)下一期已實(shí)現(xiàn)波動(dòng)率的影響通常大于左側(cè)較大收益。從HAR-PV(5)的回歸結(jié)果可以看出,在兩側(cè)極端的偏方差中,PV(1)的系數(shù)為0.159,t值為1.822,而對(duì)應(yīng)的PV(6)的系數(shù)僅為-0.091,t值為-1.399??梢?jiàn)左側(cè)的極端風(fēng)險(xiǎn)對(duì)預(yù)測(cè)下一期已實(shí)現(xiàn)波動(dòng)率的作用更大。這也符合金融市場(chǎng)的風(fēng)險(xiǎn)認(rèn)知,在傳統(tǒng)的金融研究里,左側(cè)風(fēng)險(xiǎn)確實(shí)是更為重要的一種。但是,這種情況在考慮兩側(cè)較大的收益時(shí)發(fā)生了反轉(zhuǎn)。同樣在PV(5)模型的回歸結(jié)果中,PV(5)的系數(shù)為-0.306,t值為-2.824,而PV(2)的系數(shù)僅為0.216,t值為1.547。 可見(jiàn)如果考慮較大的收益對(duì)下一期已實(shí)現(xiàn)波動(dòng)率的影響,右側(cè)收益的影響要大于左側(cè)收益。

    最后,從回歸結(jié)果調(diào)整后的R看,HAR-PV(4)模型對(duì)數(shù)據(jù)的擬合效果最好,其調(diào)整后的R值為0.651。此外,與傳統(tǒng)HAR模型相比,對(duì)已實(shí)現(xiàn)波動(dòng)率進(jìn)行分割后對(duì)數(shù)據(jù)的擬合效果都更好(高于0.536),該結(jié)果進(jìn)一步驗(yàn)證了對(duì)已實(shí)現(xiàn)波動(dòng)率進(jìn)行分割的必要性。

    4 結(jié)語(yǔ)

    本文研究了在HAR-RS模型中,是否可以通過(guò)改變閾值的選取提高已實(shí)現(xiàn)波動(dòng)率的預(yù)測(cè)能力的問(wèn)題。本文以HAR-RS模型為基礎(chǔ),構(gòu)造了HAR-PV(G)模型,并使用滬深300股票指數(shù)5分鐘高頻交易數(shù)據(jù)進(jìn)行實(shí)證分析,結(jié)果表明恰當(dāng)選取合適的閾值確實(shí)能夠提高對(duì)已實(shí)現(xiàn)波動(dòng)率的預(yù)測(cè)能力。

    本文對(duì)不同模型的樣本外表現(xiàn)進(jìn)行比較,發(fā)現(xiàn)閾值數(shù)量固定時(shí),平分閾值組合的偏方差模型樣本外表現(xiàn)最優(yōu)。在閾值數(shù)量可以變動(dòng)的情況下,偏方差模型的樣本外表現(xiàn)和閾值數(shù)量的關(guān)系呈現(xiàn)一個(gè)“U”形,閾值數(shù)量為3的偏方差模型的樣本外表現(xiàn)最好。隨后,本文進(jìn)一步考慮了時(shí)變閾值模型即閾值動(dòng)態(tài)模型,利用交叉檢驗(yàn)的方法在每一期選出最優(yōu)的閾值組合用以預(yù)測(cè)下一期已實(shí)現(xiàn)波動(dòng)率,并將最終的預(yù)測(cè)效果和平分閾值模型進(jìn)行比較后發(fā)現(xiàn),閾值動(dòng)態(tài)模型并不能顯著擊敗閾值數(shù)量為3的平分偏方差模型。

    為了進(jìn)一步探究不同大小的收益對(duì)下一期已實(shí)現(xiàn)波動(dòng)率的影響,本文通過(guò)平分模型的樣本內(nèi)估計(jì)發(fā)現(xiàn)負(fù)收益對(duì)下一期已實(shí)現(xiàn)波動(dòng)率的影響更大且更顯著。除此之外,0附近的收益在波動(dòng)率預(yù)測(cè)中占據(jù)主導(dǎo)地位。同時(shí),在極端的收益中,左側(cè)極端收益對(duì)預(yù)測(cè)的影響強(qiáng)于右側(cè)的極端收益;但在較大的收益中,右側(cè)較大的收益的影響卻強(qiáng)于左側(cè)較大的收益。

    雖然本文取得了一些創(chuàng)新型的結(jié)果,但還存在一些問(wèn)題需要進(jìn)一步拓展研究。具體表現(xiàn)在三個(gè)方面。第一,本文的研究主要局限在股票市場(chǎng),文中的結(jié)論能否應(yīng)用于其他金融市場(chǎng),比如期貨、大宗商品等,有待進(jìn)一步驗(yàn)證。第二,本文的研究主要局限在預(yù)測(cè)下一天的已實(shí)現(xiàn)波動(dòng)率,在預(yù)測(cè)下一周甚至下個(gè)月的已實(shí)現(xiàn)波動(dòng)率時(shí),是否發(fā)生閾值選擇的變化有待考證。第三,本文僅研究了波動(dòng)率預(yù)測(cè)的問(wèn)題,后續(xù)可以將偏方差模型運(yùn)用到資產(chǎn)配置、期權(quán)定價(jià)和風(fēng)險(xiǎn)管理等實(shí)際問(wèn)題中。

    參考文獻(xiàn):

    [1]BOLLERSLEV T. Generalized autoregressive conditional heteroskedasticity[J]. Journal of Econometrics, 1986, 31(3):307-27.

    [2]TAYLOR S J. Modeling stochastic volatility:a review and comparative study[J]. Mathematical Finance, 1994, 4(2):183-204.

    [3]ANDERSEN T G, BOLLERSLEV T, DIEBOLD F X, et al. Modeling and forecasting realized volatility[J]. Econometrica, 2003, 71(2):579-625.

    [4]CORSI F. A Simple approximate long-memory model of realized volatility[J]. Journal of Financial Econometrics, 2009, 7(2):174-196.

    [5]ANDERSEN T G, BOLLERSLEV T, DIEBOLD F X. Roughing it up:including jump components in the measurement, modeling, and forecasting of return volatility[J]. The Review of Economics and Statistics, 2007, 89(4):701-720.

    [6]CHEN X, GHYSELS E. News—good or bad—and its impact on volatility predictions over multiple horizons[J]. The Review of Financial Studies, 2010, 24(1):46-81.

    [7]PATTON A J, SHEPPARD K. Good volatility, bad volatility:signed jumps and the persistence of volatility[J]. The Review of Economics and Statistics, 2015, 97(3):683-697.

    [8]文鳳華, 劉曉群, 唐海如,等.基于LHAR-RV-V模型的中國(guó)股市波動(dòng)性研究 [J]. 管理科學(xué)學(xué)報(bào), 2012(15):59-67.

    [9]馬鋒, 魏宇, 黃登仕,等.基于跳躍和符號(hào)跳躍變差的HAR-RV預(yù)測(cè)模型及其MCS檢驗(yàn) [J]. 系統(tǒng)管理學(xué)報(bào), 2015, 24:700-710.

    [10]陳聲利, 關(guān)濤,李一軍.基于跳躍、好壞波動(dòng)率與百度指數(shù)的股指期貨波動(dòng)率預(yù)測(cè) [J]. 系統(tǒng)工程理論與實(shí)踐, 2018(38):299-316.

    [11]李俊儒, 汪壽陽(yáng),魏云捷.基于波動(dòng)率測(cè)量誤差的波動(dòng)率預(yù)測(cè)模型 [J]. 系統(tǒng)工程理論與實(shí)踐, 2018(38):1905-1918.

    [12]龔旭, 曹杰, 文鳳華,等.基于杠桿效應(yīng)和結(jié)構(gòu)突變的HAR族模型及其對(duì)股市波動(dòng)率的預(yù)測(cè)研究 [J]. 系統(tǒng)工程理論與實(shí)踐, 2020(40):1113-1133.

    [13]CORSI F, PIRINO D, RENò R. Threshold bipower variation and the impact of jumps on volatility forecasting [J]. Journal of Econometrics, 2010, 159(2):276-288.

    [14]AUDRINO F, HU Y. Volatility forecasting:downside risk, jumps and leverage effect[J]. Econometrics, 2016, 4(1):8.

    [15]HANSEN P R, LUNDE A, NASON J M. The model confidence set[J]. Econometrica, 2011, 79(2):453-497.

    [16]HANSEN P R. A test for superior predictive ability[J]. Journal of Business & Economic Statistics, 2005, 23(4):365-380.

    [17]BOLLERSLEV T, MEDEIROS M C, PATTON A J, et al. From zero to hero:realized partial (co)variances[J]. Journal of Econometrics, 2021.

    [18]WEST N. A simple, positive semi-definite, heteroskedasticity and autocorrelation consistent covariance matrix[J]. Econometrica, 1987, 55(3):703-708.

    国产熟女欧美一区二区| 熟妇人妻久久中文字幕3abv| 成人午夜高清在线视频| 1000部很黄的大片| 免费黄频网站在线观看国产| 欧美性感艳星| 免费看a级黄色片| 人妻一区二区av| 午夜激情欧美在线| 嫩草影院精品99| 国产精品综合久久久久久久免费| 亚洲三级黄色毛片| 国产亚洲av嫩草精品影院| 国产高潮美女av| 高清在线视频一区二区三区| 日韩av在线大香蕉| 中文欧美无线码| 亚洲成色77777| 永久网站在线| av在线蜜桃| 亚洲av日韩在线播放| 夜夜看夜夜爽夜夜摸| 日韩国内少妇激情av| 日韩电影二区| 内地一区二区视频在线| 日本欧美国产在线视频| 久久韩国三级中文字幕| 淫秽高清视频在线观看| 黑人高潮一二区| 亚洲精品日韩av片在线观看| 日本-黄色视频高清免费观看| 在线播放无遮挡| 亚洲国产成人一精品久久久| 亚洲最大成人手机在线| 日本爱情动作片www.在线观看| 美女黄网站色视频| 国产午夜精品久久久久久一区二区三区| 国内精品一区二区在线观看| 人人妻人人看人人澡| 国产综合懂色| 久久久久久久国产电影| 亚洲一级一片aⅴ在线观看| 白带黄色成豆腐渣| 直男gayav资源| 久久久久久久大尺度免费视频| 国产成人aa在线观看| 美女内射精品一级片tv| 嫩草影院入口| 欧美另类一区| 久久99热这里只有精品18| 国内少妇人妻偷人精品xxx网站| 国产精品一区二区性色av| 夜夜看夜夜爽夜夜摸| 国产免费福利视频在线观看| 成人亚洲精品一区在线观看 | 日韩一区二区视频免费看| 亚洲图色成人| 日本一二三区视频观看| 欧美zozozo另类| 乱人视频在线观看| 五月玫瑰六月丁香| 男人舔女人下体高潮全视频| 久久精品人妻少妇| 大陆偷拍与自拍| 国产毛片a区久久久久| 日本三级黄在线观看| 国产精品一二三区在线看| 亚洲av免费高清在线观看| 毛片一级片免费看久久久久| 偷拍熟女少妇极品色| 亚洲自偷自拍三级| 免费播放大片免费观看视频在线观看| 日韩av在线免费看完整版不卡| 九九在线视频观看精品| 九草在线视频观看| 人人妻人人看人人澡| 最近最新中文字幕大全电影3| 久久久久久久大尺度免费视频| 亚洲一区高清亚洲精品| 插阴视频在线观看视频| 91久久精品国产一区二区成人| 亚洲四区av| 久久精品夜色国产| 性色avwww在线观看| 丝袜美腿在线中文| 校园人妻丝袜中文字幕| 国产精品爽爽va在线观看网站| 亚洲国产欧美在线一区| 亚洲av.av天堂| 国产精品av视频在线免费观看| 精品久久久久久久久亚洲| av在线观看视频网站免费| 人妻夜夜爽99麻豆av| 精品国产露脸久久av麻豆 | 中文天堂在线官网| 亚洲欧美一区二区三区黑人 | 好男人在线观看高清免费视频| 日韩欧美精品免费久久| 一二三四中文在线观看免费高清| 观看美女的网站| 偷拍熟女少妇极品色| 综合色av麻豆| 久久久精品94久久精品| 日韩av免费高清视频| 五月玫瑰六月丁香| av又黄又爽大尺度在线免费看| 免费av毛片视频| 久久久成人免费电影| 精品久久久久久电影网| 热99在线观看视频| 日韩不卡一区二区三区视频在线| 亚洲精品日韩av片在线观看| 亚洲自拍偷在线| 97人妻精品一区二区三区麻豆| 国产成年人精品一区二区| 国产极品天堂在线| 99视频精品全部免费 在线| 赤兔流量卡办理| 国产久久久一区二区三区| 国产黄色视频一区二区在线观看| 欧美成人一区二区免费高清观看| 听说在线观看完整版免费高清| 免费电影在线观看免费观看| 中文欧美无线码| 欧美xxⅹ黑人| .国产精品久久| 中国美白少妇内射xxxbb| 99热这里只有是精品50| 欧美成人午夜免费资源| 成人性生交大片免费视频hd| 国产精品久久久久久av不卡| 久久久精品欧美日韩精品| 亚洲精品乱久久久久久| 国产精品无大码| 国产精品综合久久久久久久免费| 国产真实伦视频高清在线观看| 国产精品国产三级国产av玫瑰| 免费av毛片视频| 我的女老师完整版在线观看| 午夜激情久久久久久久| 激情五月婷婷亚洲| 国产黄片视频在线免费观看| 一个人观看的视频www高清免费观看| 午夜福利在线观看吧| 在线 av 中文字幕| 久久人人爽人人片av| 极品教师在线视频| 尤物成人国产欧美一区二区三区| 欧美zozozo另类| 精品久久久久久久久av| 成人漫画全彩无遮挡| 久久6这里有精品| 少妇猛男粗大的猛烈进出视频 | 亚洲av男天堂| 五月玫瑰六月丁香| 欧美一区二区亚洲| 欧美性感艳星| 18禁裸乳无遮挡免费网站照片| 成人av在线播放网站| 女人被狂操c到高潮| 青春草国产在线视频| 欧美成人精品欧美一级黄| 成人亚洲精品av一区二区| 欧美一级a爱片免费观看看| 欧美变态另类bdsm刘玥| 少妇的逼水好多| 亚洲三级黄色毛片| 亚洲色图av天堂| 亚洲精品日韩av片在线观看| 国产精品日韩av在线免费观看| 欧美成人一区二区免费高清观看| 成人二区视频| 国产淫语在线视频| 成人亚洲欧美一区二区av| 日本与韩国留学比较| 日韩精品有码人妻一区| 美女xxoo啪啪120秒动态图| 麻豆精品久久久久久蜜桃| 成年女人在线观看亚洲视频 | 51国产日韩欧美| 免费黄色在线免费观看| 免费看av在线观看网站| 九九爱精品视频在线观看| 欧美成人精品欧美一级黄| 内地一区二区视频在线| 亚洲人成网站高清观看| 天堂网av新在线| 国产精品人妻久久久久久| 久久精品国产亚洲av涩爱| 黄色日韩在线| 成年女人看的毛片在线观看| 永久免费av网站大全| 国产黄a三级三级三级人| 久久久精品欧美日韩精品| 亚洲欧美日韩卡通动漫| 美女cb高潮喷水在线观看| 最近2019中文字幕mv第一页| 日本av手机在线免费观看| 丝袜美腿在线中文| 日本欧美国产在线视频| 综合色丁香网| 成人亚洲精品一区在线观看 | 久久鲁丝午夜福利片| 免费av不卡在线播放| 亚洲电影在线观看av| 亚洲色图av天堂| 最近中文字幕2019免费版| 亚洲精品aⅴ在线观看| 国产黄频视频在线观看| 看非洲黑人一级黄片| 国内揄拍国产精品人妻在线| 成人国产麻豆网| 街头女战士在线观看网站| 永久免费av网站大全| 国产一区二区三区av在线| 色综合亚洲欧美另类图片| 高清av免费在线| 久久综合国产亚洲精品| 成人无遮挡网站| 97精品久久久久久久久久精品| 亚洲图色成人| 亚洲av中文字字幕乱码综合| 免费大片黄手机在线观看| 日本黄大片高清| 亚洲熟女精品中文字幕| 日韩欧美三级三区| 亚洲av一区综合| 日本一二三区视频观看| 成年免费大片在线观看| videos熟女内射| 中国美白少妇内射xxxbb| 国产av码专区亚洲av| 国产av在哪里看| 最近最新中文字幕免费大全7| 亚洲最大成人av| 好男人视频免费观看在线| 日韩 亚洲 欧美在线| 亚洲天堂国产精品一区在线| 日韩在线高清观看一区二区三区| av天堂中文字幕网| 亚洲成人中文字幕在线播放| 亚洲国产欧美在线一区| 久久久色成人| 国产亚洲av嫩草精品影院| 欧美xxxx性猛交bbbb| 久久草成人影院| 极品少妇高潮喷水抽搐| 国产久久久一区二区三区| 欧美zozozo另类| 国产精品人妻久久久影院| 日产精品乱码卡一卡2卡三| 99久久九九国产精品国产免费| 国产色爽女视频免费观看| 最近2019中文字幕mv第一页| 最近中文字幕高清免费大全6| 色吧在线观看| 精品人妻偷拍中文字幕| 久久久久国产网址| 又黄又爽又刺激的免费视频.| 欧美xxxx性猛交bbbb| 久久99热这里只频精品6学生| 美女国产视频在线观看| 成人毛片60女人毛片免费| 80岁老熟妇乱子伦牲交| 一级av片app| 国语对白做爰xxxⅹ性视频网站| 国产成人福利小说| 成人漫画全彩无遮挡| 一夜夜www| 99热全是精品| 18+在线观看网站| 成人毛片60女人毛片免费| 精品久久久久久久久av| 简卡轻食公司| 国产久久久一区二区三区| 人妻一区二区av| 中文天堂在线官网| 亚洲国产精品专区欧美| 少妇裸体淫交视频免费看高清| 国产老妇伦熟女老妇高清| 丝瓜视频免费看黄片| 日韩成人伦理影院| 人妻夜夜爽99麻豆av| av免费在线看不卡| 日本av手机在线免费观看| 国产黄色视频一区二区在线观看| 天堂影院成人在线观看| 免费大片黄手机在线观看| 亚洲精品久久久久久婷婷小说| 精品午夜福利在线看| 天天躁日日操中文字幕| 久久精品国产亚洲网站| 免费观看精品视频网站| 国产老妇女一区| 久久精品国产自在天天线| 舔av片在线| 简卡轻食公司| 亚洲av成人精品一区久久| av福利片在线观看| 国产真实伦视频高清在线观看| 激情五月婷婷亚洲| 亚洲精品日本国产第一区| 99久久人妻综合| 日韩国内少妇激情av| 婷婷色麻豆天堂久久| 久久久精品94久久精品| 禁无遮挡网站| 亚洲精品一区蜜桃| 国产熟女欧美一区二区| 亚洲伊人久久精品综合| 免费看美女性在线毛片视频| 超碰97精品在线观看| 大香蕉久久网| 免费av毛片视频| 在线 av 中文字幕| 成年女人在线观看亚洲视频 | 欧美激情国产日韩精品一区| 午夜福利高清视频| 国产精品人妻久久久影院| 91在线精品国自产拍蜜月| 99热全是精品| 80岁老熟妇乱子伦牲交| 精品熟女少妇av免费看| 热99在线观看视频| 91在线精品国自产拍蜜月| 99热这里只有是精品在线观看| 午夜免费观看性视频| 麻豆久久精品国产亚洲av| 亚洲熟妇中文字幕五十中出| 男人和女人高潮做爰伦理| 亚洲国产精品国产精品| 亚洲色图av天堂| 日韩中字成人| 我要看日韩黄色一级片| 在线 av 中文字幕| 亚洲最大成人手机在线| 成人av在线播放网站| 最近2019中文字幕mv第一页| 建设人人有责人人尽责人人享有的 | 免费黄色在线免费观看| 99热全是精品| 国产v大片淫在线免费观看| videos熟女内射| 欧美97在线视频| 色5月婷婷丁香| 国产毛片a区久久久久| 美女脱内裤让男人舔精品视频| 亚洲欧美一区二区三区黑人 | 七月丁香在线播放| 18禁裸乳无遮挡免费网站照片| 又爽又黄a免费视频| 亚洲国产高清在线一区二区三| 午夜亚洲福利在线播放| 亚洲精品日本国产第一区| 免费看不卡的av| 男女啪啪激烈高潮av片| 国产一区二区亚洲精品在线观看| 水蜜桃什么品种好| 九草在线视频观看| 日韩精品青青久久久久久| www.色视频.com| 国精品久久久久久国模美| 国产精品一二三区在线看| 国产激情偷乱视频一区二区| 国产精品久久久久久久电影| 男人舔女人下体高潮全视频| 日韩成人av中文字幕在线观看| 国产乱来视频区| 人人妻人人看人人澡| 纵有疾风起免费观看全集完整版 | 一级a做视频免费观看| 亚洲综合精品二区| 色尼玛亚洲综合影院| 色网站视频免费| 国产黄a三级三级三级人| 久久99热6这里只有精品| 女人久久www免费人成看片| 91av网一区二区| 毛片一级片免费看久久久久| 国产综合懂色| 亚洲精品aⅴ在线观看| 国产亚洲最大av| 超碰97精品在线观看| 女人久久www免费人成看片| 国产成人91sexporn| 午夜精品在线福利| kizo精华| 精品一区二区免费观看| 日韩 亚洲 欧美在线| 91狼人影院| 在线天堂最新版资源| 高清av免费在线| 国产精品一二三区在线看| 免费大片黄手机在线观看| 国产精品久久久久久精品电影| 男人和女人高潮做爰伦理| 日韩欧美一区视频在线观看 | 精品熟女少妇av免费看| 男女边吃奶边做爰视频| 亚洲一区高清亚洲精品| 国产白丝娇喘喷水9色精品| 国产 亚洲一区二区三区 | 好男人视频免费观看在线| 国产成人a区在线观看| 亚洲av中文av极速乱| 国产黄片视频在线免费观看| 国产在线一区二区三区精| 国产免费福利视频在线观看| a级毛片免费高清观看在线播放| 成人鲁丝片一二三区免费| 亚洲成色77777| 天堂√8在线中文| 欧美日韩亚洲高清精品| 国产有黄有色有爽视频| 成人特级av手机在线观看| 蜜桃亚洲精品一区二区三区| 人体艺术视频欧美日本| or卡值多少钱| 日日啪夜夜爽| 日韩av不卡免费在线播放| 色播亚洲综合网| av专区在线播放| 亚洲成人中文字幕在线播放| av国产久精品久网站免费入址| 免费观看精品视频网站| 免费看a级黄色片| av卡一久久| 国产大屁股一区二区在线视频| 日本黄大片高清| 美女xxoo啪啪120秒动态图| 国产久久久一区二区三区| 国产精品一区二区性色av| 中国国产av一级| 日日啪夜夜爽| 男人舔奶头视频| 国产成人a∨麻豆精品| 啦啦啦韩国在线观看视频| 黄片无遮挡物在线观看| 久久韩国三级中文字幕| 久久精品熟女亚洲av麻豆精品 | 色哟哟·www| 亚洲国产成人一精品久久久| 国内精品美女久久久久久| 青青草视频在线视频观看| 天堂影院成人在线观看| 亚洲精品乱码久久久久久按摩| 日韩av在线大香蕉| 久久精品熟女亚洲av麻豆精品 | 18禁在线播放成人免费| 高清在线视频一区二区三区| 久久久亚洲精品成人影院| 亚洲婷婷狠狠爱综合网| 80岁老熟妇乱子伦牲交| 久久久久精品久久久久真实原创| 成年人午夜在线观看视频 | 国产精品国产三级专区第一集| 性插视频无遮挡在线免费观看| 成人欧美大片| 日韩精品有码人妻一区| 水蜜桃什么品种好| 三级经典国产精品| av一本久久久久| 精品亚洲乱码少妇综合久久| 国产乱人视频| 欧美成人a在线观看| 少妇的逼水好多| 亚洲av男天堂| 美女高潮的动态| 国产成人精品久久久久久| 日本免费a在线| 精品人妻一区二区三区麻豆| 最近视频中文字幕2019在线8| 精品久久久精品久久久| 中文字幕久久专区| 日产精品乱码卡一卡2卡三| 亚洲国产最新在线播放| 中国美白少妇内射xxxbb| av国产免费在线观看| 2021少妇久久久久久久久久久| 国产精品99久久久久久久久| 亚洲怡红院男人天堂| 99热6这里只有精品| 国产精品1区2区在线观看.| 国产高清国产精品国产三级 | 国产精品福利在线免费观看| 18禁裸乳无遮挡免费网站照片| 久久久久免费精品人妻一区二区| 97超视频在线观看视频| 欧美丝袜亚洲另类| 网址你懂的国产日韩在线| 一边亲一边摸免费视频| 亚洲一区高清亚洲精品| av在线蜜桃| 久久99热这里只频精品6学生| 天堂网av新在线| 少妇人妻一区二区三区视频| 亚洲婷婷狠狠爱综合网| 少妇的逼水好多| 国产 亚洲一区二区三区 | 亚洲欧美一区二区三区国产| 少妇熟女aⅴ在线视频| 午夜激情久久久久久久| 嘟嘟电影网在线观看| 纵有疾风起免费观看全集完整版 | 亚洲精品第二区| .国产精品久久| or卡值多少钱| 精品人妻一区二区三区麻豆| eeuss影院久久| 免费看av在线观看网站| 成人毛片a级毛片在线播放| 一夜夜www| 少妇的逼水好多| 国产欧美另类精品又又久久亚洲欧美| 男人狂女人下面高潮的视频| 免费观看av网站的网址| 一区二区三区四区激情视频| 少妇的逼水好多| 日韩三级伦理在线观看| 18禁动态无遮挡网站| 亚洲乱码一区二区免费版| 国产美女午夜福利| 一级毛片黄色毛片免费观看视频| 久久6这里有精品| 国产伦理片在线播放av一区| 国产精品国产三级国产av玫瑰| 欧美丝袜亚洲另类| 肉色欧美久久久久久久蜜桃 | 最近最新中文字幕免费大全7| 嘟嘟电影网在线观看| 亚洲精品aⅴ在线观看| 久久韩国三级中文字幕| 精品久久久久久久人妻蜜臀av| 在线免费观看的www视频| 中文字幕制服av| 91午夜精品亚洲一区二区三区| av国产久精品久网站免费入址| 内射极品少妇av片p| 丰满乱子伦码专区| 在线观看人妻少妇| av免费观看日本| 熟女人妻精品中文字幕| 伦理电影大哥的女人| 国产精品日韩av在线免费观看| 成年女人在线观看亚洲视频 | 久久国产乱子免费精品| 成人性生交大片免费视频hd| 亚洲成色77777| 国产综合精华液| 99久久精品一区二区三区| 亚洲国产精品专区欧美| 亚洲不卡免费看| 欧美日韩在线观看h| 国产成人aa在线观看| 一级毛片黄色毛片免费观看视频| 久久97久久精品| 国产伦在线观看视频一区| 国产国拍精品亚洲av在线观看| 真实男女啪啪啪动态图| 少妇丰满av| 成人亚洲精品av一区二区| 免费高清在线观看视频在线观看| 五月天丁香电影| 国产精品.久久久| 我要看日韩黄色一级片| 国产精品人妻久久久久久| 成人特级av手机在线观看| 亚洲丝袜综合中文字幕| av免费在线看不卡| 内地一区二区视频在线| 亚洲成人精品中文字幕电影| 精品人妻熟女av久视频| 大话2 男鬼变身卡| 男女视频在线观看网站免费| 国产伦在线观看视频一区| 欧美另类一区| 欧美97在线视频| 国产成人91sexporn| 亚洲久久久久久中文字幕| 亚洲成人精品中文字幕电影| 国产精品美女特级片免费视频播放器| 天堂中文最新版在线下载 | 日本午夜av视频| 国产精品国产三级国产专区5o| 男的添女的下面高潮视频| 国产亚洲av嫩草精品影院| 夫妻性生交免费视频一级片| 欧美成人a在线观看| 黄色一级大片看看| 嫩草影院新地址| 成人亚洲精品av一区二区| 国产黄频视频在线观看| 午夜激情福利司机影院| 亚洲激情五月婷婷啪啪| 97精品久久久久久久久久精品| 纵有疾风起免费观看全集完整版 | 中文字幕制服av| 哪个播放器可以免费观看大片| 非洲黑人性xxxx精品又粗又长| 2022亚洲国产成人精品| 狂野欧美白嫩少妇大欣赏| 日本免费在线观看一区| 精品久久久久久久末码| 熟女电影av网| 高清欧美精品videossex| 免费观看在线日韩| 99九九线精品视频在线观看视频| 国产在线一区二区三区精| 日本wwww免费看| 九九爱精品视频在线观看| 日韩伦理黄色片| 51国产日韩欧美| 色综合色国产| 日本猛色少妇xxxxx猛交久久| 天堂av国产一区二区熟女人妻|