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

    基于插值修正的上證50ETF期權(quán)隱含高階矩研究

    2024-12-16 00:00:00王欣王俊葉五一
    中國證券期貨 2024年6期

    摘要:期權(quán)隱含高階矩是衡量金融資產(chǎn)波動性和非對稱性的關(guān)鍵指標(biāo)。然而,真實市場中期權(quán)的執(zhí)行價格具有離散性和有界性,直接應(yīng)用理論方法基于市場離散期權(quán)數(shù)據(jù)計算的隱含高階矩存在近似誤差,影響指標(biāo)可靠性。本文利用上證50ETF期權(quán)數(shù)據(jù),采用無模型方法計算隱含波動率和隱含偏度指標(biāo),并通過“內(nèi)插-外推”方法進(jìn)行修正。研究結(jié)果表明:①未修正的隱含高階矩通常被高估;②隱含波動率指數(shù)與上證50ETF收益之間存在顯著的負(fù)相關(guān)性,這支持了市場存在杠桿效應(yīng)和波動率反饋效應(yīng)的理論;③隱含偏度指數(shù)與上證50ETF收益之間存在顯著的正相關(guān)性,這可以理解為投資者更偏好右偏資產(chǎn),而右偏資產(chǎn)價值傾向于被高估,從而導(dǎo)致負(fù)的收益;④在大多數(shù)時期,已實現(xiàn)波動率比隱含波動率低,這表明大多數(shù)投資者大部分時間傾向于規(guī)避風(fēng)險,并且愿意為對沖風(fēng)險支付費用。

    關(guān)鍵詞:期權(quán)隱含高階矩;上證50ETF;已實現(xiàn)高階矩;風(fēng)險溢價

    一、引言

    長期以來,期權(quán)價格中蘊含的信息一直是學(xué)術(shù)界和業(yè)界關(guān)注的焦點。早期研究表明,通過Black-Scholes(BS)公式可以獲得隱含波動率(Black和Scholes,1973)。隨著期權(quán)定價模型與技術(shù)手段的不斷發(fā)展,研究者發(fā)現(xiàn)期權(quán)面板數(shù)據(jù)中還包含了更為復(fù)雜的信息,如隱含偏度。這些期權(quán)隱含信息被廣泛認(rèn)為是對真實測度下歷史信息的重要補充,因為期權(quán)價格反映了標(biāo)的資產(chǎn)在風(fēng)險中性測度下的預(yù)期走勢。本文主要探討期權(quán)的隱含信息,包括隱含波動率和隱含偏度。這兩者分別代表了風(fēng)險中性測度下標(biāo)的資產(chǎn)收益率的預(yù)期波動性和偏度。

    提取期權(quán)價格中的隱含信息對學(xué)術(shù)界和業(yè)界都具有重要意義。在實際應(yīng)用中,隱含信息可作為市場信號使用。例如,芝加哥期權(quán)交易所(CBOE)的市場波動率指數(shù)(VIX),通常被稱為“恐慌指數(shù)”,用于預(yù)測未來市場波動率;中國波動率指數(shù)(iVIX)用于預(yù)警市場潛在風(fēng)險,并作為投資者調(diào)整倉位的重要信號。隱含偏度可以反映股票的預(yù)期未來走勢,隱含偏度值越小預(yù)示著未來資產(chǎn)收益分布越左偏,即資產(chǎn)收益出現(xiàn)下跌的概率會增加。除此之外,這些信息可以與真實測度下的信息結(jié)合,用于構(gòu)建風(fēng)險溢價指標(biāo),以反映市場預(yù)期和實際情況的偏差,從而為研究市場收益提供支持。因此,充分利用期權(quán)價格中的隱含信息可以為風(fēng)險管理和金融資產(chǎn)定價提供有效指導(dǎo)。

    一些重要文獻(xiàn)研究了期權(quán)隱含高階矩。例如,Christoffersen等(2013)總結(jié)了期權(quán)隱含信息的相關(guān)研究,介紹了兩種構(gòu)建期權(quán)隱含高階矩的方法:參數(shù)方法(Heston,1993;Christoffersen等,2009)和無模型方法(Carr和Madan,1998)。相比之下,無模型方法無須依賴特定的模型設(shè)定,因此被認(rèn)為更具靈活性。在實際市場中,芝加哥期權(quán)交易所就是采用無模型方法來計算VIX指數(shù)。

    目前,關(guān)于期權(quán)隱含信息的研究已經(jīng)取得了諸多成果。Jiang和Tian(2005)發(fā)現(xiàn),無模型方法計算的隱含波動率可以更有效地預(yù)測未來實際波動率。鄭振龍等(2016)利用S&P500指數(shù)期權(quán)的方差和偏度對股票橫截面序列進(jìn)行了回歸分析。叢明舒(2018)發(fā)現(xiàn),中國市場和美國市場的期權(quán)隱含信息對股票收益率的影響存在顯著差異。

    盡管無模型方法具有多種優(yōu)勢,但仍存在一些局限性,如它依賴期權(quán)執(zhí)行價格連續(xù)性的假設(shè)。然而,在實際市場中,我們只能觀察到離散的期權(quán)執(zhí)行價格,因此這一假設(shè)并不成立。此外,期權(quán)執(zhí)行價格的范圍是有界的,這與計算公式中的無界假設(shè)不符。為了盡可能解決這些局限性,本文在無模型方法的基礎(chǔ)上引入了“內(nèi)插-外推”方法,以修正隱含高階矩的估算值,并將修正后的隱含高階矩信息應(yīng)用于中國市場,進(jìn)一步研究修正后的隱含高階矩對50ETF收益率的影響,以及隱含高階矩與已實現(xiàn)高階矩之間的關(guān)系。此外,本文基于風(fēng)險中性測度和真實測度構(gòu)建了方差風(fēng)險溢價指標(biāo),并研究了該指標(biāo)對50ETF收益率的影響。

    本文剩余的結(jié)構(gòu)安排如下:第二部分詳細(xì)介紹了隱含高階矩指標(biāo)的構(gòu)建方法。第三部分分析了計算隱含高階矩時近似誤差的來源。第四部分提出了近似誤差的修正方法。第五部分進(jìn)行實證分析。第六部分總結(jié)本文的研究結(jié)論,并討論其在實踐中的潛在應(yīng)用和意義。

    二、隱含高階矩構(gòu)建

    從期權(quán)面板數(shù)據(jù)中計算隱含高階矩的方法主要分為有模型和無模型兩種。由于前者容易受模型不確定性的影響,故本文主要借鑒Bakshi等(2003)關(guān)于無模型方法的思想來推導(dǎo)期權(quán)隱含高階矩。設(shè)定R(t,τ)=lnS(t+τ)-lnS(t)為t時刻τ期的收益率,S(t)為t時刻資產(chǎn)的價格,則在t時刻τ期的風(fēng)險中性收益率的隱含波動率和隱含偏度可以分別定義為:

    Vol(t,τ)={Et[R(t,τ)-Et(R(t,τ))]2}1/2(1)

    Skew(t,τ)=Et[R(t,τ)-Et(R(t,τ))]3{Et[R(t,τ)-Et(R(t,τ))]2}3/2(2)

    這里,Et(·)為風(fēng)險中性下的條件期望,即Et(·)=E(·|S(t))。

    根據(jù)CBOE白皮書、Bakshi等(2003)以及Bali和Murray(2013)的研究,可以通過虛值期權(quán)數(shù)據(jù)推導(dǎo)出對應(yīng)資產(chǎn)的隱含高階矩:

    Vol(t,τ)=(erτV(t,τ)-μ(t,τ)2)1/2(3)

    Skew(t,τ)=erτW(t,τ)-3μ(t,τ)erτV(t,τ)+2μ(t,τ)3[erτV(t,τ)-μ(t,τ)2]3/2(4)

    其中,μ(t,τ)=erτ-1-erτ2V(t,τ)-erτ6W(t,τ)-erτ24X(t,τ),r為無風(fēng)險利率,τ為期權(quán)到期期限。隱含高階矩中對應(yīng)的函數(shù)V、W和X的具體形式如下:

    V(t,τ)=∫∞S(t)21-lnKS(t)K2C(t,τ;K)dK+∫S(t)021+lnS(t)KK2P(t,τ;K)dK(5)

    W(t,τ)=∫∞S(t)6lnKS(t)-3lnKS(t)2K2C(t,τ;K)dK-∫S(t)06lnS(t)K+3lnS(t)K2K2P(t,τ;K)dK(6)

    X(t,τ)=∫∞S(t)12lnKS(t)2+4lnKS(t)3K2

    C(t,τ;K)dK+∫S(t)012lnS(t)K2+4lnS(t)K3K2P(t,τ;K)dK(7)

    其中,C(t,τ;K)和P(t,τ;K)分別是t時刻開始t+τ時刻到期且執(zhí)行價格為K的看漲和看跌期權(quán)價格。

    為了計算V、W和X,需要使用連續(xù)執(zhí)行價格下的看漲和看跌期權(quán)的價格。然而,由于期權(quán)面板數(shù)據(jù)是離散的,只能利用離散的期權(quán)數(shù)據(jù)計算V、W和X。

    V︿(t,τ)=∑nci=121-lnKciS(t)*Kci2C(t,τ;Kci)△Kci+∑npi=121+lnS(t)*KpiKpi2P(t,τ;Kpi)△Kpi(8)

    W︿(t,τ)=∑nci=16lnKciS(t)*-3lnKciS(t)*2Kci2C(t,τ;Kci)△Kci-∑npi=16lnS(t)*Kpi+3lnS(t)*Kpi2Kpi2P(t,τ;Kpi)△Kpi(9)

    X︿(t,τ)=∑nci=112lnKciS(t)*2+4lnKciS(t)*3Kci2C(t,τ;Kci)△Kci+∑npi=112lnS(t)*Kpi2+4lnS(t)*Kpi3Kpi2P(t,τ;Kpi)△Kpi(10)

    這里,S(t)*=S(t)-PVDivs為除息后的股價,PVDivs為股息現(xiàn)值,Kci和Kpi分別為第i個虛值看漲期權(quán)和看跌期權(quán)的行權(quán)價,C(t,τ;Kci)和P(t,τ;Kpi)分別為第i個看漲和看跌期權(quán)的價格,nc和np分別為看漲和看跌期權(quán)的數(shù)量,對于2≤i≤nc,△Kci=Kci-1-Kci,△Kc1=S*t-Kc1,對于2≤i≤np,△Kpi=Kpi-1-Kpi,△Kp1=S*t-Kp1。這樣通過離散方法可以計算出隱含波動率和隱含偏度。

    為了估計無模型隱含高階矩,式(8)~式(10)中給出的離散近似V︿、W︿和X︿會產(chǎn)生多種近似誤差,如截斷誤差和離散化誤差。這些誤差可能對估計隱含高階矩產(chǎn)生重要影響。接下來,將詳細(xì)介紹每種近似誤差。

    三、近似誤差來源

    在本節(jié)中,筆者將探討在利用離散期權(quán)數(shù)據(jù)計算隱含高階矩時產(chǎn)生的兩種近似誤差,即截斷誤差和離散化誤差。

    (一)截斷誤差

    截斷誤差是由于執(zhí)行價格K范圍有界導(dǎo)致的。正如V、W和X的定義所示,無模型隱含高階矩的計算需要0到無窮的執(zhí)行價格,然而,當(dāng)使用離散化的數(shù)據(jù)進(jìn)行計算時,只能利用有限范圍的執(zhí)行價格數(shù)據(jù)。記KL和KU分別是給定期限的最低和最高執(zhí)行價格。由于無限范圍的執(zhí)行價格被有限范圍的執(zhí)行價格[KL,KU]所取代,因此可以得到以下近似式(先不考慮給定范圍內(nèi)的離散情況):

    V(t,τ)≈∫KUS(t)21-lnKS(t)K2C(t,τ;K)dK+∫S(t)KL21+lnS(t)KK2P(t,τ;K)dK(11)

    W(t,τ)≈∫KUS(t)6lnKS(t)-3lnKS(t)2K2C(t,τ;K)dK-∫S(t)KL6lnS(t)K+3lnS(t)K2K2P(t,τ;K)dK(12)

    X(t,τ)≈∫KUS(t)12lnKS(t)2+4lnKS(t)3K2C(t,τ;K)dK+∫S(t)KL12lnS(t)K2+4lnS(t)K3K2

    P(t,τ;K)dK(13)

    因此,V、W和X的截斷誤差由以下表達(dá)式給出:

    Vtrunc=-∫∞KU21-lnKS(t)K2C(t,τ;K)dK+∫KL021+lnS(t)KK2P(t,τ;K)dK(14)

    Wtrunc=-∫∞KU6lnKS(t)-3lnKS(t)2K2C(t,τ;K)dK-∫KL06lnS(t)K+3lnS(t)K2K2

    P(t,τ;K)dK(15)

    Xtrunc=-∫∞KU12lnKS(t)2+4lnKS(t)3K2

    C(t,τ;K)dK+∫KL012lnS(t)K(t)2+4lnS(t)K3K2

    P(t,τ;K)dK(16)

    負(fù)號表明截斷誤差會導(dǎo)致計算的高階矩產(chǎn)生負(fù)的偏差。

    (二)離散化誤差

    離散化誤差是由于數(shù)值積分導(dǎo)致的,因為CBOE在計算中采用式(8)、式(9)、式(10)。盡管從理論上看,數(shù)值積分誤差可以通過使用足夠精細(xì)的執(zhí)行價格網(wǎng)絡(luò)來最小化,但實際期權(quán)面板數(shù)據(jù)中可獲取的離散數(shù)據(jù)始終是有限的,且行權(quán)價之間的離散程度往往很大,因此會產(chǎn)生較為明顯的離散化誤差。離散化誤差可以表示為:

    Vdisc=V︿(t,τ)-∫KUS(t)21-lnKS(t)K2C(t,τ;K)dK+∫S(t)KL21+lnS(t)KK2P(t,τ;K)dK(17)

    Wdisc=W︿(t,τ)-∫KUS(t)6lnKS(t)-3lnKS(t)2K2C(t,τ;K)dK-∫S(t)KL6lnS(t)K+3lnS(t)K2K2

    P(t,τ;K)dK(18)

    Xdisc=X︿(t,τ)-∫KUS(t)12lnKS(t)2+4lnKS(t)3K2C(t,τ;K)dK+∫S(t)KL12lnS(t)K2+4lnS(t)K3K2

    P(t,τ;K)dK(19)

    四、近似誤差的處理

    近似誤差主要由于期權(quán)數(shù)據(jù)在不同執(zhí)行價格下的稀疏性引起。為了減少這種誤差,需要獲取更多執(zhí)行價格對應(yīng)的期權(quán)價格數(shù)據(jù)。解決這一問題的一種方法是通過“平滑”技術(shù)填補稀疏的期權(quán)數(shù)據(jù)。通常而言,“平滑”技術(shù)應(yīng)用于期權(quán)合約的波動率,而非直接應(yīng)用于期權(quán)價格。具體步驟如下:首先,利用BS公式計算已有期權(quán)價格對應(yīng)的波動率;其次,使用“平滑”技術(shù)構(gòu)建期權(quán)合約對應(yīng)的波動率曲線,以獲取所需執(zhí)行價格對應(yīng)的波動率;最后,再次利用BS公式計算相應(yīng)執(zhí)行價格的看漲和看跌期權(quán)價格。本文采用了一種有效的平滑方法來提取無模型波動率,即采用“內(nèi)插-外推”方法構(gòu)建波動率曲面。

    (一)已有執(zhí)行價格間的內(nèi)插法

    假定對于給定的到期日T有N個執(zhí)行價格,分別為Ki,i=1,2,…,N,其對應(yīng)的看漲看跌期權(quán)的市場價格分別是CM(Ki,T)和PM(Ki,T)。不失一般性,有0<KL=K1<K2<…<KN=KU<∞。這里僅使用OTM期權(quán)來建立波動率函數(shù)。通過BS模型,結(jié)合N個執(zhí)行價格以及對應(yīng)的期權(quán)價格,可以反推出這N個執(zhí)行價格對應(yīng)的波動率(σ(Ki,T),i=1,2,…,N)。

    構(gòu)建執(zhí)行價格區(qū)間[KL,KU]上的波動率函數(shù)是關(guān)鍵。已有研究中,參數(shù)法和非參數(shù)法都被廣泛應(yīng)用(Bahra,2007;Bliss和Panigirtzoglou,2002;Anagnou等,2002)。然而,參數(shù)方法對模型形式的限制,無法靈活適應(yīng)觀測數(shù)據(jù)。因此,本文采用非參數(shù)方法來構(gòu)建波動率函數(shù),以便更好地處理和利用不同執(zhí)行價格下的期權(quán)價格數(shù)據(jù)。

    采用“內(nèi)插法”在N個執(zhí)行價格之間建立光滑函數(shù),可以很好地擬合這N個已知的波動率。筆者希望在執(zhí)行價格區(qū)間[KL,KU]找到一個可微函數(shù)f(K),有

    f(Ki)=σ(Ki,T),i=1,2,…,N(20)

    這里參考Jiang和Tian(2005),使用三次樣條來擬合已知的波動率。為了使擬合結(jié)果更好,從K2到KN-1所有的內(nèi)部執(zhí)行價格均被用作節(jié)點來擬合三次樣條。接著,對于新的內(nèi)插點便可以通過f(K)獲得其對應(yīng)合約的波動率。

    (二)已有執(zhí)行價格外的外推法

    由于在已列出的執(zhí)行價格區(qū)間之外,沒有已知的期權(quán)價格,而范圍更廣的執(zhí)行價格對于計算隱含高階矩有較大的好處,因此需要根據(jù)已構(gòu)建的波動率函數(shù)對區(qū)間外的波動率進(jìn)行外推。具體來說,考慮兩種外推方法。

    1平推法

    平推法是指對于執(zhí)行價格區(qū)間外的波動率,用最接近的已知波動率值來代替,即對于區(qū)間(0,KL]的波動率都用σ(KL)來代替,對于區(qū)間[KU,∞)的波動率都用σ(KU)來代替,其核心思想如圖1上圖所示,該方法常被用于波動率估計的修正。然而,“平推法”存在兩個缺點:一是為了滿足無套利條件,波動率函數(shù)需要保持連續(xù)性,但是“平推法”可能導(dǎo)致波動率函數(shù)在KL和KU處出現(xiàn)拐點,從而使風(fēng)險中性概率密度函數(shù)變?yōu)樨?fù)數(shù);二是根據(jù)波動率“微笑”或者波動率“偏斜”的觀察,隨著執(zhí)行價格遠(yuǎn)離當(dāng)前股價,波動率往往會上升。因此,對于執(zhí)行價格區(qū)間外的波動率使用常數(shù)波動率代替,可能會導(dǎo)致波動率被低估。

    2線性外推法

    線性外推法指的是在最低KL和最高KU的執(zhí)行價格點上施加平滑條件,并保持一個線性的外推結(jié)構(gòu)。具體來說,通過調(diào)整處于執(zhí)行價格兩邊的線性外推斜率,使其與最低KL或最高KU執(zhí)行價格內(nèi)部端的斜率一致,如圖1下圖所示。這種方法確保了在整個執(zhí)行價格范圍內(nèi)都有平滑的波動率曲線,并且能夠很好地擬合已知的波動率數(shù)據(jù)。“線性外推法”構(gòu)建的波動率曲線更符合通常觀察到的波動率“微笑”或者波動率“偏斜”。

    由此可見,一旦構(gòu)建好波動率函數(shù),適當(dāng)范圍內(nèi)任意執(zhí)行價格處的期權(quán)價格都可以通過使用BS模型根據(jù)相應(yīng)擬合的波動率計算得到。這樣一來,可以利用更多的期權(quán)數(shù)據(jù)對數(shù)值積分V、W和X進(jìn)行離散化計算。

    五、實證分析

    (一)數(shù)據(jù)描述

    本文選取2021年1月4日至2021年12月31日每天9:30到11:30以及13:00到15:00的50ETF指數(shù)和50ETF期權(quán)的1分鐘級數(shù)據(jù),共計243天(交易日),每天包含241個1分鐘數(shù)據(jù)。剔除所有不滿足無套利條件的期權(quán)樣本,選擇期權(quán)到期期限為10~60天,并且每天至少存在兩支認(rèn)購和認(rèn)沽期權(quán)的樣本。當(dāng)數(shù)據(jù)處理完后,可以使用前面章節(jié)的方法計算出風(fēng)險中性高階矩,這里所計算的均為期限τ為30天的風(fēng)險中性矩。

    (二)修正方案間的比較分析

    采用CBOE的指標(biāo)構(gòu)建方法,分別構(gòu)建波動率指數(shù)和偏度指數(shù):

    iVIXt,τ=100*Volt,τ(21)

    iSKEWt,τ=100-10*Skewt,τ(22)

    本節(jié)對比分析了不同的“內(nèi)插-外推”方法對50ETF波動率指數(shù)(iVIX)和偏度指數(shù)(iSKEW)的影響。鑒于截斷誤差和離散化誤差,對不同方法進(jìn)行比較分析,以評估其對iVIX和iSKEW的影響。對“內(nèi)插-外推”方法進(jìn)行兩種不同的內(nèi)插點數(shù)選擇,分別為在已有執(zhí)行價格區(qū)間內(nèi)選擇100個和500個等間距內(nèi)插點。在內(nèi)插計算完成后,還采用了線性外推和平值外推兩種方法。對于線性外推,當(dāng)內(nèi)插點數(shù)為100時,考慮2個、10個和20個等間距外推點;當(dāng)內(nèi)插點數(shù)為500時,則考慮10個、50個和100個等間距外推點,以確保外推的區(qū)間長度一致。對于平值外推,外推點個數(shù)的選擇與線性外推方法一致。

    表1展示了不同插值方法計算的波動率指數(shù)和偏度指數(shù)的描述性統(tǒng)計。筆者觀察到修正后的iVIX指數(shù)相對于未修正的iVIX指數(shù),在均值和最小值上幾乎都顯示出下降趨勢。具體而言,平值外推方法得到的iVIX指數(shù)相對于線性外推方法更小,且內(nèi)插點數(shù)為500時的iVIX指數(shù)比內(nèi)插點數(shù)為100時的iVIX指數(shù)要小。對于iSKEW的計算,平值外推對其影響不大。然而,當(dāng)線性外推點數(shù)過多時會導(dǎo)致iSKEW值顯著增加。

    圖2呈現(xiàn)了不同插值方法修正后的iVIX和iSKEW指數(shù)的時間序列。可以發(fā)現(xiàn),未插值情況下的波動率和偏度指數(shù)一般位于所有插值計算結(jié)果之上,顯示出插值后的指數(shù)普遍偏低。此外,線性外推點數(shù)過多可能會產(chǎn)生異常值。為了更好地顯示不同方法之間的差異,圖3展示了不同插值方法修正后的iVIX和iSKEW指數(shù)的平均日內(nèi)走勢??梢园l(fā)現(xiàn),對于平均日內(nèi)iVIX指數(shù),經(jīng)過平值外推和線性外推修正后的波動率指數(shù)均低于未插值計算結(jié)果,且內(nèi)插點數(shù)對其影響較小,而外推點數(shù)對其影響較大。對于平均日內(nèi)iSKEW指數(shù),平值外推修正的iSKEW指數(shù)低于未插值的計算結(jié)果。然而,線性外推修正的iSKEW指數(shù)在外推點數(shù)過多時會顯著高于未插值的計算結(jié)果。因此,不建議使用過多的線性外推點數(shù)。

    綜上所述,未插值的隱含高階矩往往存在高估的情況,這可能導(dǎo)致市場投資者在期權(quán)交易、金融資產(chǎn)定價、資產(chǎn)收益分析以及金融風(fēng)險傳染度量等方面做出錯誤的決策。并且內(nèi)插點數(shù)越多,適度的線性外推結(jié)果越穩(wěn)健。由于理論上采用插值計算的隱含高階矩更為準(zhǔn)確,且實證結(jié)果顯示插值和不插值結(jié)果差異比較大?;谶@些考量,選擇了內(nèi)插點數(shù)為500、外推點數(shù)為50的方案進(jìn)行后續(xù)數(shù)據(jù)分析和實證研究。

    (三)變量的描述性統(tǒng)計

    圖4分別展示了修正后的iVIX指數(shù)和iSKEW指數(shù)的時間序列。以下是對圖4結(jié)果的分析:①從iSKEW指數(shù)圖中可以觀察到,2021年中國市場的iSKEW指數(shù)有部分?jǐn)?shù)值小于100,處于90~100,顯示出市場預(yù)期收益率出現(xiàn)右偏的現(xiàn)象。這與Chang等(2013)在美國市場計算的隱含偏度指數(shù)始終

    大于100有所不同。②iSKEW指數(shù)通常在達(dá)到115之后會下跌,而不是持續(xù)上升,這是因為市場上存在套利機(jī)會的投資者,他們的介入會使iSKEW指數(shù)回歸其均值附近,從而避免持續(xù)的左偏現(xiàn)象。這些分析結(jié)果揭示了市場波動率和偏度指數(shù)在不同時間段的動態(tài)變化,對理解市場行為和制定投資策略具有重要意義。

    圖5顯示了以下內(nèi)容:左側(cè)為iVIX指數(shù)和iSKEW指數(shù)本身的自相關(guān)圖;右側(cè)則展示了對iVIX指數(shù)和iSKEW指數(shù)進(jìn)行差分后的自相關(guān)圖。在自相關(guān)圖中,陰影部分表示95%的置信區(qū)間。觀察圖可知,iVIX指數(shù)和iSKEW指數(shù)本身具有顯著的自相關(guān)性,然而,經(jīng)過一階差分處理后,大部分iVIX指數(shù)和iSKEW指數(shù)的自相關(guān)系數(shù)均落在95%的置信區(qū)間內(nèi)。這表明對隱含高階矩進(jìn)行差分處理后的時間序列更加平穩(wěn)。

    (四)回歸分析

    1收益率回歸分析

    Jondeau等(2019)在其研究中使用美國數(shù)據(jù)發(fā)現(xiàn),偏度可以影響個股的收益率。Harvey和Siddique(2000,2000)的研究結(jié)果顯示,偏度不僅影響個股的收益率,還對整個股票市場的收益率產(chǎn)生影響。那么,在中國市場,隱含高階矩是否會對50ETF的收益率產(chǎn)生影響呢?

    為了研究隱含高階矩對上證50ETF指數(shù)收益的影響,本節(jié)使用了差分后的iVIX指數(shù)、iSKEW指數(shù)及其組合對上證50ETF指數(shù)的收益率進(jìn)行回歸分析。研究數(shù)據(jù)覆蓋了2021年1月4日至2021年12月31日的243個交易日的50ETF和期權(quán)數(shù)據(jù)。我們建立了以下三個模型:

    Model1:Rt=α+β1△iVIXt+t(23)

    Model2:Rt=α+β2△iSKEWt+t(24)

    Model3:Rt=α+β1△iVIXt+β2△iSKEWt+t(25)

    這里的Rt是第t個時刻50ETF的對數(shù)價格收益率,△iVIXt=iVIXt-iVIXt-1,△iSKEWt=iSKEWt-iSKEWt-1。表2展示了隱含高階矩與市場收益率回歸的結(jié)果。

    從表2的回歸結(jié)果可以得出:①收益率與△iVIX之間表現(xiàn)出顯著的負(fù)相關(guān)性,這種關(guān)系主要是由于杠桿效應(yīng)和波動率反饋效應(yīng)所導(dǎo)致的。②△iSKEW對收益率具有顯著的預(yù)測能力,這可以理解為投資者更偏好右偏的資產(chǎn),而右偏的資產(chǎn)往往被高估,從而導(dǎo)致負(fù)的收益。此外,當(dāng)iSKEW>100時,表明資產(chǎn)收益存在左偏的可能性,進(jìn)一步指示資產(chǎn)收益下跌的風(fēng)險增加。由于風(fēng)險厭惡的投資者會尋求更多的風(fēng)險補償來規(guī)避尾部風(fēng)險的影響,因此更高的△iSKEW通常伴隨著更高的資產(chǎn)收益。③當(dāng)△iVIX和△iSKEW同時加入收益率的回歸模型后,它們對收益率的影響依然顯著。這表明,從期權(quán)數(shù)據(jù)中提取的隱含高階矩均包含了影響50ETF收益率的重要信息。

    2已實現(xiàn)高階矩之間的回歸分析

    通過期權(quán)數(shù)據(jù)計算的隱含波動率指數(shù)和隱含偏度指數(shù)間接反映了投資者對市場風(fēng)險的度量,而已實現(xiàn)波動率(RV)和已實現(xiàn)偏度(RSkew)是直接通過50ETF收益率數(shù)據(jù)計算的,能夠直接反映資產(chǎn)的收益情況,因此它們之間存在一定的關(guān)聯(lián)。在本節(jié)中,建立已實現(xiàn)高階矩與隱含高階矩之間的回歸模型,并分析它們之間的關(guān)系。首先,根據(jù)Andersen和Bollerslev(1998)及后續(xù)文獻(xiàn)提出的方法,構(gòu)建已實現(xiàn)高階矩:

    RVt=∑N~i=0R2t-i(26)

    RSkewt=N~∑N~i=1R3t-i(RVt)3/2(27)

    RSKEWt=100-10×RSkewt(28)

    這里,Rt-i表示第t-i時刻的對數(shù)收益率。使用滑動步長為1分鐘,窗寬大小N~為240分鐘的數(shù)據(jù)集來計算高頻RV和RSkewt。接下來,建立已實現(xiàn)高階矩與隱含高階矩之間的回歸關(guān)系模型:

    RV2t=α+βiVIX2t+t(29)

    RSKEWt=α+βiSKEWt+t(30)

    使用修正后的隱含波動率指數(shù)和隱含偏度指數(shù),得到的回歸結(jié)果如下:

    RV2t=-00294+15606iVIX2t(-969840)(2187779)(31)

    RSKEWt=1570276-05636iSKEWt

    (595898)(-215865)(32)

    其中括號內(nèi)為各個系數(shù)的t-統(tǒng)計量。從隱含高階矩對已實現(xiàn)高階矩的回歸結(jié)果可以看出:①已實現(xiàn)波動率和iVIX指數(shù)之間存在顯著的正相關(guān)關(guān)系;②已實現(xiàn)偏度指數(shù)與iSKEW指數(shù)之間存在顯著的負(fù)相關(guān)關(guān)系。由此可以說明50ETF期權(quán)隱含信息對50ETF收益率的方差和偏度具有一定的解釋作用。

    3方差風(fēng)險溢價回歸分析

    已實現(xiàn)波動率反映了投資者對市場真實風(fēng)險的度量,而隱含波動率則反映了對未來波動率風(fēng)險的預(yù)期。因此,方差風(fēng)險溢價(VRP)可以定義為在風(fēng)險中性測度下的方差減去真實測度下的已實現(xiàn)方差,即VRPt=iVIX2t-RV2t。VRP在一定程度上反映了投資者的風(fēng)險規(guī)避程度。圖6展示了高頻已實現(xiàn)波動率(RV)和隱含波動率(iVIX)的時間序列,可以觀察到在許多時段,iVIX比RV要大。

    圖7呈現(xiàn)了方差風(fēng)險溢價??梢钥闯?,VRP大多時候都是大于0的,這說明中國投資者更傾向于規(guī)避風(fēng)險。當(dāng)未來市場走向不確定增加時,投資者規(guī)避風(fēng)險的意愿會增強(qiáng),導(dǎo)致VRP處于一個較高的水平。然而,當(dāng)市場出現(xiàn)高波動(如2021年2月以及2021年7月)時,VRP會出現(xiàn)很大的負(fù)值。這表明在這些時期,市場相對較為動蕩,導(dǎo)致真實市場中由50ETF價格直接計算的已實現(xiàn)波動率要高于投資者的預(yù)期波動率。

    Bollerslev等(2009)使用美國數(shù)據(jù)發(fā)現(xiàn)VRP在對季度收益率的預(yù)測上表現(xiàn)最為強(qiáng)烈;Qiao等(2024)則利用9個新興市場的指數(shù)構(gòu)建VRP,發(fā)現(xiàn)VRP可以預(yù)測股指收益、貨幣收益以及資本流動,并且在短期內(nèi)預(yù)測效果尤其顯著。本節(jié)利用插值的隱含波動率指數(shù)計算VRP,并研究其對上證50ETF收益率的影響。建立以下回歸模型:

    Rt=α+βVRPt+t(33)

    其中,Rt表示第t個時刻50ETF的對數(shù)價格收益率,VRPt表示第t個時刻的方差風(fēng)險溢價。得到的回歸結(jié)果如下:

    Rt=-00744+57866VRPt

    (-03819)(05615)(34)

    括號內(nèi)為參數(shù)估計的t-統(tǒng)計量?;貧w結(jié)果顯示,截距項和斜率項在95%置信水平下均不顯著,表明VRP對上證50ETF的收益率沒有顯著的影響,這與田鳳平和楊科(2021)的研究結(jié)果一致。

    六、結(jié)論

    本文利用上證50ETF及其期權(quán)數(shù)據(jù),采用“內(nèi)插-外推”方法,計算了50ETF指數(shù)的風(fēng)險中性隱含高階矩,并將其與基于CBOE白皮書未插值計算的隱含高階矩進(jìn)行了比較。結(jié)果顯示,未插值處理的隱含波動率指數(shù)和隱含偏度指數(shù)都存在一定程度的被高估。隱含波動率的被高估通常意味著期權(quán)價格可能也被高估,這可能誤導(dǎo)依賴隱含波動率指數(shù)進(jìn)行期權(quán)交易的投資者,從而導(dǎo)致?lián)p失。鑒于國內(nèi)期權(quán)市場近年來的快速發(fā)展,精確構(gòu)建隱含波動率和偏度指標(biāo)對維護(hù)投資者信心和衍生品市場的穩(wěn)定發(fā)展至關(guān)重要。因此,本文建議采用基于“內(nèi)插-外推”方法修正的隱含波動率和隱含偏度作為新的研究指標(biāo)。

    同時,本文還研究了插值后的隱含高階矩信息對50ETF收益率以及已實現(xiàn)高階矩的影響,并進(jìn)一步基于隱含高階矩構(gòu)建了方差風(fēng)險溢價,分析其對50ETF收益率的影響。研究發(fā)現(xiàn),隱含高階矩與50ETF收益率之間存在顯著的相關(guān)性。一方面,隱含波動率指數(shù)和收益率之間的負(fù)相關(guān)性為中國市場可能存在杠桿效應(yīng)和波動率反饋效應(yīng)提供了一定的實證依據(jù)。另一方面,隱含偏度指數(shù)與收益率之間的正相關(guān)性在一定程度上說明了當(dāng)市場投資者更偏好右偏資產(chǎn)時,由于這類資產(chǎn)被高估,可能導(dǎo)致負(fù)收益。此外,本文還發(fā)現(xiàn)已實現(xiàn)高階矩與隱含高階矩之間存在顯著的相關(guān)性,這表明50ETF期權(quán)隱含信息對50ETF收益率的方差和偏度具有一定的解釋作用。對于方差風(fēng)險溢價,筆者發(fā)現(xiàn)市場平靜期間,隱含波動率通常高于已實現(xiàn)波動率,表明中國投資者大部分時間在規(guī)避風(fēng)險,并且愿意為對沖風(fēng)險支付費用。同時,研究發(fā)現(xiàn)中國市場的方差風(fēng)險溢價對收益率的影響并不顯著。

    07d8b73b6b3a126084e933814ebd3169

    綜上所述,投資者在制定投資策略時應(yīng)充分考慮期權(quán)隱含的高階矩風(fēng)險信息,并合理配置資產(chǎn)。同時,監(jiān)管部門也應(yīng)關(guān)注期權(quán)隱含的高階矩風(fēng)險,以便更好地管理衍生品市場的風(fēng)險。

    參考文獻(xiàn)

    [1]ANAGNOUI,BEDENDOM,HODGESSD,etalTherelationbetweenimpliedandrealisedprobabilitydensityfunctions[J]SSRNElectronicJournal,2002

    [2]ANDERSENTG,BOLLERSLEVTAnsweringtheskeptics:Yes,standardvolatilitymodelsdoprovideaccurateforecasts[J]InternationalEconomicReview,1998,39(4):885-905

    [3]BAHRABImpliedrisk-neutralprobabilitydensityfunctionsfromoptionprices:Acentralbankperspective[J]ForecastingVolatilityintheFinancialMarkets(ThirdEdition),2007:201-226

    [4]BAKSHIG,KAPADIAN,MADANDStockreturncharacteristics,skewlaws,andthedifferentialpricingofindividualequityoptions[J]TheReviewofFinancialStudies,2003,16(1):101-143

    [5]BALITG,MURRAYSDoesrisk-neutralskewnesspredictthecross-sectionofequityoptionportfolioreturns?[J]JournalofFinancialandQuantitativeAnalysis,2013,48(4):1145-1171

    [6]BLISSRR,PANIGIRTZOGLOUNTestingthestabilityofimpliedprobabilitydensityfunctions[J]JournalofBankingandFinance,2002,26(2-3):381-422

    [7]BLACKF,SCHOLESMThepricingofoptionsandcorporateliabilities[J]JournalofPoliticalEconomy,1973,81(3):637-654

    [8]BOLLERSLEVT,TAUCHENG,ZHOUHExpectedstockreturnsandvarianceriskpremia[J]TheReviewofFinancialStudies,2009,22(11):4463-4492

    [9]CHANGBY,CHRISTOFFERSENP,JACOBSKMarketskewnessriskandthecrosssectionofstockreturns[J]JournalofFinancialEconomics,2013,107(1):46-68

    [10]HARVEYCR,SIDDIQUEAConditionalskewnessinassetpricingtests[J]TheJournalofFinance,2000,55(3):1263-1295

    [11]HARVEYCR,SIDDIQUEATime-varyingconditionalskewnessandthemarketriskpremium[J]ResearchinBankingandFinance,2000,1(1):27-60

    [12]JIANGGJ,TIANYSExtractingmodel-freevolatilityfromoptionprices:AnexaminationoftheVIXindex[J]TheJournalofDerivatives,2007,14(3):35-60

    [13]JIANGGJ,TIANYSThemodel-freeimpliedvolatilityanditsinformationcontent[J]TheReviewofFinancialStudies,2005,18(4):1305-1342

    [14]JONDEAUE,ZHANGQ,ZHUXAverageskewnessmatters[J]JournalofFinancialEconomics,2019,134(1):29-47

    [15]CARRP,MADANDTowardsatheoryofvolatilitytrading[J]Volatility:Newestimationtechniquesforpricingderivatives,1998(29):417-427

    [16]CARRP,WULVarianceriskpremiums[J]TheReviewofFinancialStudies,2009,22(3):1311-1341

    [17]CHRISTOFFERSENP,JACOBSK,CHANGBYForecastingwithoptionimpliedinformation[J]HandbookofEconomicForecasting,2013,2(A):581-656

    [18]CHRISTOFFERSENP,HESTONS,JACOBSKTheshapeandtermstructureoftheindexoptionsmirk:Whymultifactorstochasticvolatilitymodelsworksowell[J]ManagementScience,2009,55(12):1914-1932

    19]QIAOF,XUL,ZHANGX,etalVarianceriskpremiumsinemergingmarkets[J]JournalofBanking&Finance,2024(167):107259

    [20]HESTONSLAclosed-formsolutionforoptionswithstochasticvolatilitywithapplicationstobondandcurrencyoptions[J]TheReviewofFinancialStudies,1993,6(2):327-343.

    [21]叢明舒中國場內(nèi)期權(quán)市場研究——基于中美關(guān)于期權(quán)隱含方差的差異[J]金融研究,2018(12):189-206

    [22]李志勇,余湄,汪壽陽方差風(fēng)險溢價和收益率預(yù)測:來自上證50ETF期權(quán)市場的證據(jù)[J]系統(tǒng)工程理論與實踐,2022,42(2):306-319

    [23]田鳳平,楊科我國商品期貨的無模型隱含方差與方差風(fēng)險溢價研究基于豆粕和白糖的分析[J]系統(tǒng)工程理論與實踐,2021,41(8):2015-2029

    [24]王琳玉,倪中新,郭婧上證50ETF隱含高階矩風(fēng)險對股票收益的預(yù)測研究[J]統(tǒng)計研究,2020,37(12):75-90

    [25]鄭振龍,孫清泉,吳強(qiáng)方差和偏度的風(fēng)險價格[J]管理科學(xué)報,2016,19(12):110-123

    日本撒尿小便嘘嘘汇集6| 国产av一区二区精品久久| 亚洲欧美日韩无卡精品| АⅤ资源中文在线天堂| 麻豆国产av国片精品| tocl精华| 日日夜夜操网爽| 精品日产1卡2卡| 91大片在线观看| 久久久久国产精品人妻aⅴ院| 国产精品久久久久久久电影 | e午夜精品久久久久久久| 无人区码免费观看不卡| 国产精品免费一区二区三区在线| 黄色 视频免费看| 婷婷精品国产亚洲av| 国产亚洲av高清不卡| 一边摸一边做爽爽视频免费| 黄色毛片三级朝国网站| 亚洲在线自拍视频| 精品日产1卡2卡| 观看免费一级毛片| 久久久久久国产a免费观看| 悠悠久久av| 亚洲国产精品sss在线观看| 床上黄色一级片| www.www免费av| 操出白浆在线播放| 91九色精品人成在线观看| 亚洲av成人不卡在线观看播放网| 国内毛片毛片毛片毛片毛片| 免费看日本二区| 9191精品国产免费久久| 无遮挡黄片免费观看| 最近视频中文字幕2019在线8| 国产熟女午夜一区二区三区| 国产av麻豆久久久久久久| 精品一区二区三区视频在线观看免费| 久久香蕉国产精品| 国产成人av教育| 可以在线观看毛片的网站| 欧美日韩黄片免| 88av欧美| 99久久国产精品久久久| 女人爽到高潮嗷嗷叫在线视频| 国产精品久久久久久亚洲av鲁大| 午夜亚洲福利在线播放| www.精华液| 老熟妇乱子伦视频在线观看| 日韩欧美一区二区三区在线观看| 午夜福利视频1000在线观看| 久久 成人 亚洲| 999久久久国产精品视频| 日韩三级视频一区二区三区| 狂野欧美白嫩少妇大欣赏| 99精品在免费线老司机午夜| 99国产极品粉嫩在线观看| 国产精品爽爽va在线观看网站| 18禁黄网站禁片午夜丰满| 亚洲精品av麻豆狂野| 一区二区三区高清视频在线| 久久久水蜜桃国产精品网| 亚洲成av人片免费观看| 国产在线精品亚洲第一网站| 精品一区二区三区视频在线观看免费| 欧美日韩中文字幕国产精品一区二区三区| 欧美3d第一页| 久热爱精品视频在线9| 他把我摸到了高潮在线观看| 1024手机看黄色片| 黑人操中国人逼视频| 国产三级中文精品| 亚洲无线在线观看| 人妻久久中文字幕网| 12—13女人毛片做爰片一| 1024香蕉在线观看| 免费在线观看黄色视频的| 真人做人爱边吃奶动态| 日韩欧美精品v在线| 极品教师在线免费播放| 精品午夜福利视频在线观看一区| 三级男女做爰猛烈吃奶摸视频| 免费av毛片视频| 在线观看午夜福利视频| 一级毛片精品| 日本免费一区二区三区高清不卡| av在线天堂中文字幕| 久久久久性生活片| 免费无遮挡裸体视频| 日韩欧美一区二区三区在线观看| 亚洲成av人片免费观看| 香蕉久久夜色| 国产欧美日韩一区二区三| 亚洲精华国产精华精| 欧洲精品卡2卡3卡4卡5卡区| 999久久久精品免费观看国产| 亚洲18禁久久av| 久久久久国内视频| 成人欧美大片| 精品人妻1区二区| 听说在线观看完整版免费高清| 国产97色在线日韩免费| 在线观看免费视频日本深夜| 久久精品影院6| 三级毛片av免费| 99国产精品99久久久久| 国产精品久久久av美女十八| 1024视频免费在线观看| 黄色a级毛片大全视频| 每晚都被弄得嗷嗷叫到高潮| 国产一区二区三区在线臀色熟女| 国产69精品久久久久777片 | 波多野结衣高清无吗| 国产一区二区三区在线臀色熟女| 精品午夜福利视频在线观看一区| 亚洲精品中文字幕在线视频| 变态另类丝袜制服| 久久久久性生活片| 777久久人妻少妇嫩草av网站| 免费av毛片视频| 亚洲成av人片免费观看| netflix在线观看网站| 欧美日本亚洲视频在线播放| a在线观看视频网站| 中文字幕最新亚洲高清| 日本免费一区二区三区高清不卡| 一个人免费在线观看电影 | 国产亚洲精品一区二区www| 日本五十路高清| 免费观看人在逋| 国产精品亚洲av一区麻豆| 婷婷亚洲欧美| 天堂√8在线中文| 国产麻豆成人av免费视频| 欧美高清成人免费视频www| 欧美一区二区精品小视频在线| 日韩三级视频一区二区三区| 亚洲欧美精品综合久久99| 国内精品一区二区在线观看| 欧美在线黄色| 久久九九热精品免费| 国产99久久九九免费精品| 男人舔奶头视频| 制服丝袜大香蕉在线| 精品久久蜜臀av无| 久久婷婷成人综合色麻豆| 90打野战视频偷拍视频| 性欧美人与动物交配| 特级一级黄色大片| 亚洲av美国av| 精品一区二区三区四区五区乱码| 757午夜福利合集在线观看| 国产精品av视频在线免费观看| 色精品久久人妻99蜜桃| 三级毛片av免费| 欧美 亚洲 国产 日韩一| 国产精品av视频在线免费观看| 久久九九热精品免费| 日日摸夜夜添夜夜添小说| 2021天堂中文幕一二区在线观| 免费观看人在逋| 亚洲国产中文字幕在线视频| 日本黄色视频三级网站网址| 国产av不卡久久| 日本一区二区免费在线视频| 国产精品国产高清国产av| 男女做爰动态图高潮gif福利片| 午夜免费激情av| 亚洲 欧美一区二区三区| 久久久久亚洲av毛片大全| а√天堂www在线а√下载| 亚洲精品国产一区二区精华液| 一级黄色大片毛片| 女生性感内裤真人,穿戴方法视频| 国产成人一区二区三区免费视频网站| 亚洲午夜精品一区,二区,三区| 天堂动漫精品| 色综合欧美亚洲国产小说| 国产精品1区2区在线观看.| 99久久精品热视频| 国内精品久久久久久久电影| 国产精品免费视频内射| 久久久久久久午夜电影| 亚洲av电影不卡..在线观看| 三级国产精品欧美在线观看 | e午夜精品久久久久久久| 老司机在亚洲福利影院| 亚洲av电影在线进入| 欧美日韩亚洲国产一区二区在线观看| 在线观看免费视频日本深夜| 叶爱在线成人免费视频播放| 男女那种视频在线观看| 日韩欧美精品v在线| 最近最新免费中文字幕在线| 在线国产一区二区在线| 亚洲自拍偷在线| 99国产综合亚洲精品| 99热这里只有精品一区 | 97碰自拍视频| 岛国在线免费视频观看| 九色国产91popny在线| 又粗又爽又猛毛片免费看| av在线天堂中文字幕| 亚洲国产精品久久男人天堂| 日韩精品青青久久久久久| 大型av网站在线播放| 免费在线观看日本一区| 午夜福利欧美成人| 亚洲人成电影免费在线| 在线观看免费视频日本深夜| 性色av乱码一区二区三区2| 伦理电影免费视频| 精品电影一区二区在线| 免费看十八禁软件| 国产精品乱码一区二三区的特点| 日韩欧美三级三区| 国产成人一区二区三区免费视频网站| 精品国产亚洲在线| 色播亚洲综合网| 成人一区二区视频在线观看| 少妇人妻一区二区三区视频| 一边摸一边抽搐一进一小说| 9191精品国产免费久久| 亚洲国产欧洲综合997久久,| 欧美一级毛片孕妇| 伊人久久大香线蕉亚洲五| 男女视频在线观看网站免费 | 久久人妻福利社区极品人妻图片| 亚洲欧美日韩高清在线视频| 在线永久观看黄色视频| aaaaa片日本免费| 国产黄片美女视频| 精品国产乱码久久久久久男人| 国产片内射在线| 又粗又爽又猛毛片免费看| 99久久无色码亚洲精品果冻| 日韩国内少妇激情av| 99热只有精品国产| 国产激情久久老熟女| 日韩免费av在线播放| 老鸭窝网址在线观看| 中文字幕av在线有码专区| 久久99热这里只有精品18| 国产成人aa在线观看| 国产乱人伦免费视频| 九九热线精品视视频播放| 可以在线观看的亚洲视频| 久久香蕉激情| 美女午夜性视频免费| АⅤ资源中文在线天堂| 非洲黑人性xxxx精品又粗又长| 真人做人爱边吃奶动态| 男人舔女人的私密视频| 村上凉子中文字幕在线| 88av欧美| 精品久久久久久久末码| 欧美一级毛片孕妇| 亚洲欧美精品综合一区二区三区| 亚洲美女视频黄频| 精品久久久久久久末码| 精品第一国产精品| 麻豆久久精品国产亚洲av| 日韩三级视频一区二区三区| 很黄的视频免费| 岛国在线观看网站| 国产欧美日韩一区二区三| 一本一本综合久久| 可以在线观看毛片的网站| 亚洲精品一区av在线观看| 两个人视频免费观看高清| 成人三级做爰电影| av福利片在线观看| 亚洲av美国av| 久久久久九九精品影院| 国产精品电影一区二区三区| 久久精品夜夜夜夜夜久久蜜豆 | 不卡av一区二区三区| www.999成人在线观看| 欧美一级毛片孕妇| 亚洲美女黄片视频| 精品无人区乱码1区二区| 夜夜看夜夜爽夜夜摸| 欧美成人性av电影在线观看| 国产亚洲精品一区二区www| 亚洲精品美女久久久久99蜜臀| 成人三级做爰电影| 国产成人系列免费观看| 国产区一区二久久| 97碰自拍视频| 操出白浆在线播放| 日韩免费av在线播放| 国模一区二区三区四区视频 | 男人舔奶头视频| 国产精华一区二区三区| av福利片在线观看| 亚洲av熟女| 精品久久久久久久人妻蜜臀av| 在线免费观看的www视频| 亚洲av电影不卡..在线观看| 国产精品久久久av美女十八| 国产成人欧美在线观看| 日韩精品青青久久久久久| 十八禁人妻一区二区| 国产又黄又爽又无遮挡在线| 亚洲成人久久性| 久久久久九九精品影院| 又紧又爽又黄一区二区| 亚洲五月天丁香| 老司机深夜福利视频在线观看| 黄色视频不卡| 99精品在免费线老司机午夜| 久久久水蜜桃国产精品网| 中文字幕最新亚洲高清| 久久这里只有精品中国| 精品久久久久久成人av| 国产在线观看jvid| 最新在线观看一区二区三区| 精品电影一区二区在线| 久99久视频精品免费| 麻豆成人av在线观看| 国产成人精品无人区| 9191精品国产免费久久| 亚洲中文av在线| 久久国产精品人妻蜜桃| 午夜a级毛片| 性欧美人与动物交配| 精品人妻1区二区| 免费在线观看黄色视频的| 熟女电影av网| 亚洲av熟女| 亚洲专区字幕在线| 在线观看美女被高潮喷水网站 | 草草在线视频免费看| 国产区一区二久久| 色在线成人网| 国产v大片淫在线免费观看| 国产精品影院久久| x7x7x7水蜜桃| 99国产极品粉嫩在线观看| 免费在线观看黄色视频的| 亚洲中文字幕一区二区三区有码在线看 | 中国美女看黄片| 国产精品一区二区三区四区免费观看 | 男人舔奶头视频| 成人av在线播放网站| 亚洲中文字幕一区二区三区有码在线看 | 午夜精品久久久久久毛片777| 亚洲熟女毛片儿| 久久精品91无色码中文字幕| 777久久人妻少妇嫩草av网站| 人成视频在线观看免费观看| 国语自产精品视频在线第100页| 欧美成人性av电影在线观看| 老司机在亚洲福利影院| av超薄肉色丝袜交足视频| 国产精品国产高清国产av| 欧美日韩乱码在线| 欧美精品啪啪一区二区三区| 母亲3免费完整高清在线观看| 又黄又粗又硬又大视频| 一个人免费在线观看电影 | 日韩 欧美 亚洲 中文字幕| 国产亚洲欧美98| 国产av在哪里看| 97超级碰碰碰精品色视频在线观看| 成人国产一区最新在线观看| 免费在线观看成人毛片| 人人妻人人澡欧美一区二区| av片东京热男人的天堂| 国产野战对白在线观看| 欧美中文日本在线观看视频| 97人妻精品一区二区三区麻豆| 国产视频内射| 啦啦啦免费观看视频1| 精品久久蜜臀av无| 在线观看一区二区三区| 日韩欧美在线二视频| av欧美777| 99热这里只有精品一区 | 制服人妻中文乱码| 久久精品91无色码中文字幕| 成人一区二区视频在线观看| 久久国产乱子伦精品免费另类| 丝袜人妻中文字幕| 麻豆一二三区av精品| 中文字幕人成人乱码亚洲影| 免费在线观看黄色视频的| 欧美 亚洲 国产 日韩一| 国产av一区二区精品久久| 99在线视频只有这里精品首页| 国产激情偷乱视频一区二区| 午夜免费成人在线视频| 亚洲七黄色美女视频| 欧美久久黑人一区二区| 女人高潮潮喷娇喘18禁视频| 日本 欧美在线| 在线观看舔阴道视频| 亚洲激情在线av| 午夜老司机福利片| av福利片在线观看| 亚洲av日韩精品久久久久久密| 国产成人啪精品午夜网站| 天天躁夜夜躁狠狠躁躁| cao死你这个sao货| 日日夜夜操网爽| avwww免费| 午夜精品久久久久久毛片777| 国产亚洲精品久久久久久毛片| 俺也久久电影网| 亚洲 欧美 日韩 在线 免费| 国产精品美女特级片免费视频播放器 | 老司机深夜福利视频在线观看| 国产黄色小视频在线观看| 特级一级黄色大片| 欧美成人性av电影在线观看| 精品不卡国产一区二区三区| 国产av又大| 身体一侧抽搐| 亚洲成人久久爱视频| 亚洲专区中文字幕在线| 亚洲欧美激情综合另类| 国产不卡一卡二| 法律面前人人平等表现在哪些方面| 一区福利在线观看| 伦理电影免费视频| 十八禁网站免费在线| 日日爽夜夜爽网站| 欧美日韩乱码在线| 国产精品99久久99久久久不卡| 一边摸一边抽搐一进一小说| 一区二区三区高清视频在线| 美女 人体艺术 gogo| av天堂在线播放| 在线看三级毛片| 人人妻人人澡欧美一区二区| 午夜福利成人在线免费观看| 少妇被粗大的猛进出69影院| 欧美久久黑人一区二区| 国产三级黄色录像| 好看av亚洲va欧美ⅴa在| 国产成人精品无人区| av福利片在线观看| 99精品久久久久人妻精品| 亚洲av成人不卡在线观看播放网| 美女 人体艺术 gogo| 国产精品久久久av美女十八| 丰满人妻熟妇乱又伦精品不卡| 美女高潮喷水抽搐中文字幕| 高清在线国产一区| 精品午夜福利视频在线观看一区| 久久亚洲精品不卡| 日本一本二区三区精品| 国产一区二区三区在线臀色熟女| 亚洲专区字幕在线| 久久天躁狠狠躁夜夜2o2o| 亚洲电影在线观看av| 久久久国产成人免费| 亚洲国产高清在线一区二区三| 日韩av在线大香蕉| 一边摸一边做爽爽视频免费| 午夜精品久久久久久毛片777| 一区二区三区激情视频| 国产精品久久久久久人妻精品电影| 在线观看66精品国产| 欧美色视频一区免费| 久久久久亚洲av毛片大全| 中出人妻视频一区二区| 亚洲男人天堂网一区| 日本五十路高清| 少妇的丰满在线观看| 久久精品aⅴ一区二区三区四区| 黄色a级毛片大全视频| 舔av片在线| 欧美激情久久久久久爽电影| 午夜福利视频1000在线观看| 午夜免费成人在线视频| 亚洲男人的天堂狠狠| 日本精品一区二区三区蜜桃| 国产成人精品久久二区二区免费| 亚洲av熟女| 久久久久久免费高清国产稀缺| 热99re8久久精品国产| 精华霜和精华液先用哪个| 级片在线观看| 久久久久久大精品| 国产精品一区二区三区四区免费观看 | 在线免费观看的www视频| 国产精品野战在线观看| 亚洲精品中文字幕一二三四区| 亚洲熟女毛片儿| 亚洲人与动物交配视频| 日韩精品免费视频一区二区三区| 九色国产91popny在线| 亚洲av电影在线进入| 日韩精品中文字幕看吧| 欧美成狂野欧美在线观看| av中文乱码字幕在线| 成人三级做爰电影| 天堂√8在线中文| 国产亚洲精品一区二区www| 亚洲国产欧美一区二区综合| 亚洲专区中文字幕在线| 午夜福利成人在线免费观看| 久久精品影院6| 国产伦一二天堂av在线观看| 黄片大片在线免费观看| 99国产极品粉嫩在线观看| 亚洲精品国产一区二区精华液| 亚洲精品美女久久久久99蜜臀| 搡老熟女国产l中国老女人| 1024视频免费在线观看| 老熟妇仑乱视频hdxx| 少妇的丰满在线观看| 久久久久久久久久黄片| 91九色精品人成在线观看| 久久久久国内视频| 亚洲av成人不卡在线观看播放网| 国产激情久久老熟女| 19禁男女啪啪无遮挡网站| 超碰成人久久| 两性夫妻黄色片| 2021天堂中文幕一二区在线观| 熟女少妇亚洲综合色aaa.| 亚洲国产欧洲综合997久久,| 亚洲va日本ⅴa欧美va伊人久久| 一边摸一边抽搐一进一小说| ponron亚洲| 欧美成人免费av一区二区三区| 国产精品一区二区三区四区久久| 香蕉久久夜色| 亚洲精品久久国产高清桃花| 精品电影一区二区在线| 欧美国产日韩亚洲一区| 亚洲精品国产精品久久久不卡| 香蕉久久夜色| 男人的好看免费观看在线视频 | 天堂影院成人在线观看| 久久久久久国产a免费观看| 欧美激情久久久久久爽电影| 日韩有码中文字幕| 一a级毛片在线观看| 欧美久久黑人一区二区| 欧美日韩亚洲综合一区二区三区_| 亚洲午夜精品一区,二区,三区| 欧美成狂野欧美在线观看| 亚洲成人精品中文字幕电影| 天堂影院成人在线观看| 91成年电影在线观看| 老熟妇乱子伦视频在线观看| 免费在线观看黄色视频的| 亚洲美女视频黄频| 日日摸夜夜添夜夜添小说| 亚洲精品中文字幕一二三四区| 悠悠久久av| 久久久久国产一级毛片高清牌| 免费在线观看完整版高清| 麻豆一二三区av精品| 久久精品91无色码中文字幕| 天天躁夜夜躁狠狠躁躁| 校园春色视频在线观看| 国产精品一区二区三区四区久久| 亚洲国产精品合色在线| 在线观看免费午夜福利视频| 18禁黄网站禁片免费观看直播| 亚洲人成网站高清观看| 叶爱在线成人免费视频播放| 欧美三级亚洲精品| 国产区一区二久久| 国产成年人精品一区二区| 免费在线观看视频国产中文字幕亚洲| 在线观看免费日韩欧美大片| av中文乱码字幕在线| 中文在线观看免费www的网站 | 成人国产综合亚洲| 亚洲一卡2卡3卡4卡5卡精品中文| 男人舔女人下体高潮全视频| 在线国产一区二区在线| 午夜福利在线观看吧| 日韩精品青青久久久久久| 最新美女视频免费是黄的| 极品教师在线免费播放| 一边摸一边做爽爽视频免费| 欧美色欧美亚洲另类二区| 天堂动漫精品| 亚洲欧美日韩东京热| 麻豆国产av国片精品| 草草在线视频免费看| 一区二区三区国产精品乱码| 18美女黄网站色大片免费观看| 久久久国产欧美日韩av| 人妻丰满熟妇av一区二区三区| 51午夜福利影视在线观看| 日韩免费av在线播放| 高清在线国产一区| 亚洲精品粉嫩美女一区| 女同久久另类99精品国产91| xxxwww97欧美| 男女下面进入的视频免费午夜| 国产精品电影一区二区三区| 日本 欧美在线| 日韩精品青青久久久久久| 两个人看的免费小视频| 午夜亚洲福利在线播放| 男人舔女人的私密视频| 中文字幕最新亚洲高清| 久久精品国产亚洲av高清一级| 国产爱豆传媒在线观看 | 90打野战视频偷拍视频| 亚洲熟女毛片儿| 性欧美人与动物交配| 久久香蕉国产精品| 在线观看免费视频日本深夜| 精品午夜福利视频在线观看一区| 国产亚洲欧美在线一区二区| 麻豆成人av在线观看|