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

    參數(shù)法、半?yún)?shù)法的動態(tài)VaR模型風險度量

    2016-12-20 12:30:47黃旭東林雪勤
    統(tǒng)計與決策 2016年23期
    關(guān)鍵詞:置信水平度量收益

    張 瓊,黃旭東,林雪勤

    (1.安徽師范大學a.經(jīng)濟管理學院;b.數(shù)學計算機科學學院,安徽蕪湖241000;2.國泰安教育技術(shù)股份有限公司,合肥230088)

    參數(shù)法、半?yún)?shù)法的動態(tài)VaR模型風險度量

    張 瓊1a,黃旭東1b,林雪勤2

    (1.安徽師范大學a.經(jīng)濟管理學院;b.數(shù)學計算機科學學院,安徽蕪湖241000;2.國泰安教育技術(shù)股份有限公司,合肥230088)

    文章采用參數(shù)法和半?yún)?shù)法,分別考慮標準化收益在GED、SGT、GPD分布下以及FSH方法下的GARCH模型、EGARCH模型和PGARCH模型的風險測度的準確性,據(jù)此組建了12種風險測度的動態(tài)VaR模型,并采用道瓊斯股票市場指數(shù)和上證指數(shù)進行實證分析。對收益率進行基本統(tǒng)計分析發(fā)現(xiàn)兩個股票市場的收益率都不服從一般的正態(tài)分布。運用后驗測試的方法,對所有模型的樣本外預(yù)測動態(tài)VaR值采用LR、LRuc和DQ三種方法綜合檢驗,并由損失函數(shù)值可以看出:GARCH模型的風險度量能力最弱,在置信水平99.5%下,EGARCH模型最準確,在置信水平95%下,PGARCH模型最準確;GED分布描述市場的準確程度相對最弱,在較高的置信水平下,半?yún)?shù)模型能更好地度量市場的風險,在較低的置信水平下,參數(shù)模型能更好地度量股票市場的風險。

    動態(tài)VaR模型;風險測度;損失函數(shù)

    0 引言

    做好風險管理,最重要的一個環(huán)節(jié)是進行風險度量,金融風險度量的主要手段之一是在險價值(Value at Risk, VaR)。VaR可以看作是建立在過去和現(xiàn)在信息上的未來投資組合的分位數(shù),由于投資組合的收益分布隨時間在變化,如何準確的估計VaR,在實踐中是一個很大的挑戰(zhàn)。估算VaR方法主要有參數(shù)法、半?yún)?shù)法和非參數(shù)法。

    在參數(shù)法中,GARCH模型是應(yīng)用最為廣泛的時間序列波動率模型,然而在GARCH模型應(yīng)用中,基礎(chǔ)的研究通常假設(shè)資產(chǎn)收益服從正態(tài)分布,而大量的實證分析表明資產(chǎn)收益表現(xiàn)出非對稱性、厚尾、尖峰等特征,因此假設(shè)服從正態(tài)分布來估計投資組合的動態(tài)VaR,對于風險管理者來說,不能產(chǎn)生好的預(yù)測結(jié)果,常常嚴重低估風險。如徐煒和黃炎龍(2008)比較研究GARCH族的11種模型分別在正態(tài)分布和Skewed-t分布下度量動態(tài)VaR值的精確程度,結(jié)果表明Skewed-t分布較好地擬合了金融資產(chǎn)的厚尾特性,Cheng和Hung(2011)在GARCH模型中比較了正態(tài)分布、t分布、SSD和GED下的動態(tài)VaR,結(jié)論是假設(shè)服從SSD和GED下的動態(tài)VaR得出的結(jié)論最好。Polanski和Stoja (2010)比較了正態(tài)分布、t分布、SGT和EGB2分布,發(fā)現(xiàn)僅后兩種分布提供了精確的動態(tài)VaR估計。

    半?yún)?shù)方法主要是將非參數(shù)和參數(shù)法結(jié)合。在非參數(shù)法中,最典型的就是歷史模擬法,Trenca(2009)等結(jié)論表明單純使用歷史模擬法對于估計VaR不準確,歷史模擬法常常會低估風險,并且對于一些重大事件的靈敏性比較低。Barone-Adesi等(1999)提出的一種半?yún)?shù)方法:FSH (Filtered Historical Simulation)方法,即將歷史模擬法和條件波動率模型相結(jié)合。Marimoutou等(2009)認為FSH繼承了歷史模擬法優(yōu)勢的同時提高了模型的靈敏度。Zikovic等(2009)的研究表明FSH方法的預(yù)測效果比歷史模擬法優(yōu)越很多。應(yīng)用極值理論方法(EVT)考慮尾部分布所具有的特征和市場極端變動有密切關(guān)系是另一種半?yún)?shù)方法,極值理論方法不對金融收益整體分布進行建模,僅考慮分布的尾部特質(zhì)。高瑩等(2008)將基于極值理論的GARCH-EVT模型與普通的GARCH-NORMAL模型進行比較,結(jié)果表明GARCH-EVT模型優(yōu)于GARCH-NORMAL模型。Bekiros和Georgoutosos(2005)指出在較高置信水平下,EVT理論體系對于極端事件損失的預(yù)測最準確。

    近年來,已有學者在不同置信水平下研究不同模型的準確性,如葉青(2000)引入基于GARCH模型的方差協(xié)方差法和David Li提出的半?yún)?shù)法對中國的股票市場進行了模擬。Chen等(2013)比較了不同置信水平下不同方法計算的VaR值的準確性。然而對于GARCH族模型下,利用參數(shù)法和半?yún)?shù)法在不同置信水平下對動態(tài)VaR模型的準確性進行比較的文獻相對較少,且用兩種不同市場的數(shù)據(jù)同時實證分析,得到的結(jié)果顯然更具說服力。

    因此本文基于三種有代表性的GARCH模型(GARCH、EGARCH、PGARCH)的參數(shù)與半?yún)?shù)動態(tài)VaR模型,考慮標準化收益在GED、GPD、SGT分布下以及FSH方法下,探究在不同置信水平下,動態(tài)VaR模型預(yù)測的靈敏度,從而找到最優(yōu)預(yù)測模型,為金融市場的風險測度提供合適的度量工具。

    1 動態(tài)VaR模型的定義及其估計

    1.1 動態(tài)VaR模型的定義

    VaR是指在給定的概率α下,在某一段時期內(nèi),預(yù)測投資組合的最大損失值,也是投資組合收益分布的特定分位數(shù)的估計,數(shù)學表達式為:

    其中F(·)是投資資產(chǎn)組合收益的分布,?R是某一段時期內(nèi)投資組合值。VaR按是否考慮波動率的時變性,分為靜態(tài)VaR和動態(tài)VaR,靜態(tài)VaR是沒有考慮波動率的時變性,只是由前一時間段的整體波動率來預(yù)測下一刻的波動率。而動態(tài)VaR考慮了波動率的時變和股票收益序列的波動集群性,下一刻的波動率由前一時刻的波動率來預(yù)測,有更好地準確性和及時性。

    動態(tài)VaR模型的表達式為:

    故在給定置信水平1-α情況下,對于式(2),主要由兩部分組成,一是波動率(即條件均值和條件方差);二是θ,即標準化收益的分位數(shù)。

    動態(tài)VaR模型中對于波動率GARCH模型和θ值的確定都需要知道標準化收益zt的值或者zt服從的分布,本文的研究假定zt服從GED、SGT、GPD分布,將資產(chǎn)收益數(shù)據(jù)進行中心化得到條件標準化收益zt,利用最大似然法估計出zt在不同分布下GARCH、EGARCH和PGARCH三個模型中的未知參數(shù),然后又采用FSH方法對三個模型的未知參數(shù)作出估計。

    1.2 動態(tài)VaR模型波動率的估計

    GARCH模型是反映金融市場時變特征最常用的波動率模型,它能有效捕捉股市的叢集性效應(yīng)、非對稱特征等。

    GARCH族模型有多種不同類型的表達形式,其中最具有代表性的是GARCH(1,1)、EGARCH和PGARCH三個模型。

    Bollerslev在1986年提出的線性對稱GARCH(1,1)模型,具體表達式為:

    Nelson在1991年提出指數(shù)形式的GARCH模型(EGARCH模型),模型的表達式為:

    模型中條件方差采用對數(shù)形式,即使參數(shù)為負數(shù),條件方差也是正的,這就規(guī)避了參數(shù)都為正的條件。若γ≠0時,說明信息作用非對稱。

    PGARCH模型是對GARCH模型作了進一步的擴展,對標準差的冪次進行模擬,這個冪并不需要事先給定,而是通過模型自身來決定。它的形式為:則模型是對稱的;反之不對稱。E(zt)=0,VaR(zt)=1,限制條件為 β0>0,0≤β1<1,0≤β2<1。

    這里的μt和σt是建立在Ωt-1={rt-1,rt-2,…}上的條件均值和條件標準差

    1.3 動態(tài)VaR模型θ值的估計方法

    估計標準化收益zt的分位數(shù)θ,需要先知zt的整個分布或部分分布。

    (1)參數(shù)法

    在傳統(tǒng)的研究中,常常假設(shè)zt服從標準正態(tài)分布,但是在實踐中我們發(fā)現(xiàn)即使對于標準化收益的分布仍然有尖峰、厚尾等特征,需要對zt的分布作出新的假設(shè)。

    GED(Generalised Error Distribution)分布是Newey在1991年提出,其密度函數(shù)的表達式為:

    SGT(the Skewed Generalized T Distribution)分布是由Theodossiou在1998年提出的,標準化收益zt的SGT分布的密度函數(shù)表達為:

    其中sign是符號函數(shù);λ是偏度參數(shù),且|λ|<1;正數(shù)k和n是尖峰參數(shù),且B(·)是 Beta函數(shù)。

    (2)半?yún)?shù)法

    Barone-Adesi等(1999)對于標準化收益分布的刻畫提出了一種半?yún)?shù)模型,即FSH(Filtered Historical Simulation)模型。此模型是非參數(shù)和參數(shù)方法的一種結(jié)合,利用歷史模擬法來模擬標準化收益的分布,然后將歷史模擬法得到的數(shù)據(jù)特征(如捕捉金融數(shù)據(jù)的厚尾、偏態(tài)等非正常特征)和GARCH模型相結(jié)合,從而估計出標準化收益的分布。

    另一種半?yún)?shù)法是依據(jù)條件極值理論(Extreme Value Theory,EVT)進行估計。首先,定義超過閾值u分布的條件概率函數(shù)為:

    Balkerma等(1974)指出GPD(Generalised Pareto Distribution)是Fu(y)的極限分布,并且論證了當閾值u足夠大時,F(xiàn)u(y)=Gξ,σ(u)(y)。

    這里的Gξ,σ(u)(y)表示服從GDP分布,具體表達式為:

    2 動態(tài)VaR模型風險測度準確性的度量

    2.1 動態(tài)VaR模型的可行性檢驗

    對于動態(tài)VaR模型預(yù)測能力的可行性檢驗,主要運用回測檢驗的方法,本文采用非條件覆蓋似然比率檢驗(the LR of Unconditional Coverage testing)、條件覆蓋似然率檢驗(the LR of conditional Coverage test)和條件分位數(shù)回歸檢驗(the Dynamic Quantile test)這三個檢驗來綜合判斷,一個準確的模型均應(yīng)該通過這些檢驗。

    對于異常事件通常有如下描述,將預(yù)測第t+1天的動態(tài) VaRt+1(α)與當日投資組合 rt+1進行比較,若rt+1<VaRt+1(α),則稱異常事件出現(xiàn)。則定義異常指標變量如下:

    非條件覆蓋似然比率檢驗(簡記LRuc)是Kupiec(1995)建立的,這個檢驗是假設(shè)小于動態(tài)VaR(α)的觀測值的實際個數(shù)和期望值是相等的,即原假設(shè)與備擇假設(shè)分別是:似然檢驗統(tǒng)計量為:

    條件覆蓋似然率檢驗(簡記LRcc)是Christoffersen (1998)提出的,指出動態(tài)VaR估計可以看作是金融收益分布在較低尾部的區(qū)間預(yù)測。此檢驗是一種聯(lián)合檢驗,即①是否②It+1是否獨立同分布。LRcc檢驗的原假

    其中:

    nij是It=i時It+1=j出現(xiàn)的個數(shù),t=1,2,…,n-1。相應(yīng)地依據(jù)此檢驗,可以判定當產(chǎn)生過多或者過少的聚類異常事件時,拒絕動態(tài)VaR模型。

    自信息技術(shù)出現(xiàn)以后,憑借著出眾的性能與能力,信息技術(shù)成功的被應(yīng)用到人類社會發(fā)展的方方面面。信息技術(shù)不僅改變了人們對生活工作的看法,同時也提高了人們的生活學習質(zhì)量。尤其是在教育系統(tǒng)當中,信息技術(shù)的出現(xiàn)更是有著豐富教育資源,實現(xiàn)教育公平,不受時間地域限制的線上教育等功能和優(yōu)勢。成人教育是我國教育系統(tǒng)中一項非常重要的教育,和全日制教育有著不同的教育形式和教學形式。是幫助成人取得專業(yè)資格、提高技術(shù)、豐富知識、增長能力的重要途徑。因很多成人缺少足夠的時間參與課堂學習,因此信息技術(shù)在成人教育中有著重大潛力,這就要求成人教育教師必須要掌握信息教育素質(zhì)。

    條件分位數(shù)回歸檢驗(簡記為DQ檢驗)由Engle和Manganelli(2004)提出,基本原理為當動態(tài)VaR已估計,檢驗異常事件指標和屬于Ωt-1的所有變量是否相關(guān)。首先定義一個序列變量指標函數(shù):

    其中I(·)是示性函數(shù),1-α是給定的置信水平。Engle和Manganelli(2004)提出一個優(yōu)的模型應(yīng)該產(chǎn)生一系列無偏性、無相關(guān)性Hit,并且Hit的期望為0,因此DQ檢驗同時檢驗兩個方面內(nèi)容:①E(Hitt)=0;②Hitt是互不相關(guān)的。這兩個檢驗內(nèi)容通過下列回歸進行檢驗:

    其中X=(1,Hitt-1,Hitt-2,…,Hitt-p,X1,…,Xq),B=(β0,β1,…,βp+q)T,Xj是屬于Ωt-1的解釋變量。檢驗的原假設(shè)H0∶B=0。在原假設(shè)基礎(chǔ)上的條件分位數(shù)檢驗統(tǒng)計量為:

    2.2 動態(tài)VaR模型風險估計準確程度的度量

    對于風險管理者來說,關(guān)注的不僅僅是異常事件是否發(fā)生,更關(guān)注當異常事件發(fā)生后損失的大小,因為損失的大小和金融機構(gòu)儲備金多少有直接關(guān)系。Lopez(1998)提出利用損失函數(shù)作為評估異常事件發(fā)生后產(chǎn)生的后果,主要考察當異常事件發(fā)生時,觀察值rt與預(yù)測的動態(tài)VaRt之間的距離,包含了異常事件發(fā)生的次數(shù)和損失大小。表達式為:

    對于動態(tài)VaR模型的優(yōu)劣,可通過lf的大小表現(xiàn)出來,所以當一個模型的lf值小于另一個,這個模型便優(yōu)于后者。

    3 實證分析

    實證分析的數(shù)據(jù)選取了中國股票市場的上證指數(shù)和美國股票市場的道瓊斯指數(shù)。

    道瓊斯指數(shù)選用了1991年1月2日到2013年12月30日的收盤指數(shù)(共5294個樣本)作為實證分析的數(shù)據(jù),中國的上證指數(shù)選用了2012年2月20日到2014年11月7日的每日收盤指數(shù)(共660個樣本)作為實證分析的數(shù)據(jù),所有的數(shù)據(jù)來源于國泰安數(shù)據(jù)庫。

    3.2 股指對數(shù)收益率的統(tǒng)計特征

    表1 道瓊斯指數(shù)(1991/1/2-2013/12/30)描述性統(tǒng)計

    由表1可以看出,道瓊斯股票指數(shù)的每日對數(shù)收益指數(shù)JB=16977.66>5.3995,JB對應(yīng)的概率P=0.00<0.05,說明金融收益序列不服從正態(tài)分布。且偏度S=-0.153461<0,則說明序列分布左尾比右尾密集,峰度K=11.37894>3,,則說明此序列的分布比正態(tài)分布的峰部更尖。

    作出上證綜合指數(shù)的對數(shù)收益率基本統(tǒng)計分析的時間序列圖、分布圖、核密度圖和Q-Q圖,如圖1所示。

    圖1上證指數(shù)對數(shù)收益率的基本統(tǒng)計分析圖

    圖1 直觀的展現(xiàn)了我國股票市場對數(shù)收益率的一些特性。由分布圖和核密度估計圖可以觀察到我國對數(shù)收益率左偏并且有很高的峰度;Q-Q圖上的點在兩端處明顯偏離所示直線,可見上證綜合指數(shù)收益率不服從正態(tài)分布。

    對兩個不同的股票市場的收益率進行基本統(tǒng)計分析發(fā)現(xiàn),兩個市場的收益率都不服從一般的正態(tài)分布。

    3.3 模型的參數(shù)和動態(tài)VaR值估計

    對于道瓊斯指數(shù),本文選用1991—2000年的數(shù)據(jù)作為前段數(shù)據(jù),估計出條件均值-方差模型以及標準化收益的分布函數(shù)的未知參數(shù)。然后用2001年的數(shù)據(jù)來評估模型的好壞,即假設(shè)2001年的條件均值-方差滿足前10年估計出來模型,并且分布函數(shù)在2001年不變。利用2000/12/ 29(2000年交易最后一天數(shù)據(jù))估計2001/1/2(即2001年交易第一天)條件均值-方差,再利用2001/1/2真實收益數(shù)據(jù)和條件均值-方差,預(yù)測出2001/1/3的條件均值-方差,以此類推,得到2001年每天條件均值-方差。接著利用前十年已估計概率分布得到不同模型對應(yīng)的分位數(shù)θ,對伴隨概率α取不同的值,在本文中伴隨概率取0.5%、1%、2.5%和5%,由式(2)算得不同模型基于滾動樣本上2001年每天的VaR估計值,根據(jù)式(12)異常事件的定義,統(tǒng)計出相應(yīng)模型中異常事件的個數(shù)。再接著剔除1991年數(shù)據(jù),加入2001年實際收益數(shù)據(jù),同上面的做法求得2002年的動態(tài)VaR值和對應(yīng)的異常值個數(shù),同理一直滾動求解出后面的年份的動態(tài)VaR值和對應(yīng)的異常值個數(shù),直到2013年,運用每10年的滾動樣本,分別產(chǎn)生了12個不同模型后13年評估的VaR值和異常值個數(shù),得到3269個動態(tài)VaR值的評估數(shù)據(jù)。

    對于上證指數(shù),本文選用了2012/2/20—2013/3/14的數(shù)據(jù)作為觀測樣本,以前220個交易日數(shù)據(jù)作為前段數(shù)據(jù)估計出條件均值-方差模型,以及標準化收益的分布函數(shù)的未知參數(shù)。然后用接下來40個交易日的數(shù)據(jù)進行檢驗。同道瓊斯指數(shù)的估計方法,對12種模型分別產(chǎn)生了440個的評估數(shù)據(jù)。

    以上證指數(shù)為例,圖2給出了2013/1/11—2014/7/15的上證指數(shù)的對數(shù)收益率,并給出了在置信水平99%和95%下,GARCH模型下參數(shù)(SGT)和半?yún)?shù)法(GDP)估計出的每日動態(tài)VaR值。

    圖2 不同置信水平下GARCH模型的每日動態(tài)VaR值

    從圖2可以看出,在不同的置信水平下,基于GARCH模型下參數(shù)(SGT)和半?yún)?shù)法(GDP)估計出的VaR值都刻畫出了金融行業(yè)的市場風險,隨當前波動性產(chǎn)生了很明顯的變化。

    3.4 模型的檢驗和準確程度的度量

    在不同的置信水平下,采用LR、LRuc和DQ三種檢驗方法綜合判斷各個動態(tài)VaR模型的準確性是否達到既定的水平,并根據(jù)式(19)求得各個模型的最大損失函數(shù)值lf,評估出這些動態(tài)VaR模型準確程度的優(yōu)劣。

    在不同置信水平下,12種模型對于道瓊斯指數(shù)和上證指數(shù)的檢驗結(jié)果在表2至表4中,各表中第一列是波動率模型,第二列是標準化收益的分布,Lruc表示對應(yīng)模型給定樣本的Lruc檢驗值,Lrind表示對應(yīng)模型給定樣本的Lrind檢驗值,DQ表示對應(yīng)模型給定樣本的DQ檢驗值,lf表示給定樣本各個模型的最大損失函數(shù)值,排序指的是最大損失函數(shù)值對應(yīng)的排列順序,沒有全部通過三種檢驗的模型不參與排序。

    表2給出了置信水平99.5%下的動態(tài)VaR模型檢驗。當伴隨概率α=0.5%時,Lruc檢驗、Lrind檢驗和DQ檢驗對應(yīng)的拒絕零假設(shè)的邊界值分別為

    表2 置信水平99.5%下的動態(tài)VaR模型檢驗

    由表2可知,在置信水平為99.5%下,12種模型均通過第一階段的三種檢驗,即12種模型預(yù)測的準確性都可以接受。由損失函數(shù)lf的估計值可知:在相同分布下,預(yù)測效果的優(yōu)劣順序是基于EGARCH、PGARCH、GARCH下的動態(tài)VaR模型;在每類波動率模型中,F(xiàn)SH-VaR、GPD-VaR模型優(yōu)于SGT-VaR、GED-VaR模型;12種動態(tài)VaR模型中預(yù)測效果最差的是GARCH-GED模型。

    表3給出了置信水平為97.5%下的動態(tài)VaR模型檢驗。當伴隨概率α=2.5%時,Lruc檢驗、Lrind檢驗和DQ檢驗對應(yīng)的拒絕零假設(shè)的邊界值分別為

    表3 置信水平97.5%下的動態(tài)VaR模型檢驗

    由表3可知,在置信水平為97.5%時,12種模型都通過了準確性檢驗,具有可行性。根據(jù)損失函數(shù)lf的值可以看出,在相同分布下,動態(tài)VaR模型中風險度量的準確程度最差的是基于GARCH下的動態(tài)VaR模型;在EGARCH和PGARCH模型中,SGT-VaR更能準確的描繪市場,優(yōu)于半?yún)?shù)法的FSH-VaR、GPD-VaR模型,但并不顯著,預(yù)測能力最差的仍是GED-VaR模型。

    表4給出了置信水平為95%下的動態(tài)VaR模型檢驗,當伴隨概率α=5%時,Lruc檢驗、Lrind檢驗和DQ檢驗對應(yīng)的拒絕零假設(shè)的邊界值分別為

    表4 置信水平95%下的樣本外動態(tài)VaR模型檢驗

    由表4可知,在置信水平為95%下,對于道瓊斯指數(shù),組合的12種模型都通過了第一階段準確性檢驗,模型具有可行性;對于上證指數(shù),GARCH-SGT-VaR沒有通過檢驗。根據(jù)損失函數(shù)lf的值,在相同分布下,不同GARCH模型的動態(tài)VaR模型優(yōu)劣順序是基于PGARCH、EGARCH、GARCH下的動態(tài)VaR模型;對于EGARCH和PGARCH模型而言,預(yù)測準確程度最差的GED-VaR模型,參數(shù)法的SGT-VaR模型優(yōu)于半?yún)?shù)法的FSH-VaR、GPD-VaR模型;12種模型中風險度量準確程度最高的模型是PGARCH-SGT-VaR模型,最差的是GARCH-GED-VaR模型。

    4 結(jié)論

    本文對GARCH族模型中三個經(jīng)典的模型:GARCH模型、EGARCH模型和PGARCH模型的風險收益率取了三種不同分布,并運用了FSH方法組合了12種動態(tài)VaR模型,并采用道瓊斯股票指數(shù)和上證指數(shù)檢驗了12種模型在不同置信水平下的準確性和風險度量能力的優(yōu)劣,由實證分析得出以下結(jié)論:

    第一,道瓊斯股票市場和上證股票市場的收益率序列都具有尖峰厚尾的特點,不服從一般的正態(tài)分布。因此,在實際操作中,以正態(tài)分布計算VaR值的機構(gòu)或個人需提高風險意識,增強抵抗風險的能力。

    第二,在相同分布假設(shè)下,GARCH模型的風險度量能力最弱,EGARCH模型和PGARCH模型能更好地刻畫市場的波動情況。但在不同的置信水平下,EGARCH模型和PGARCH模型的度量能力也有差異,在置信水平為99.5%時,EGARCH模型的風險度量能力最好,在置信水平為95%,PGARCH模型的風險度量能力最好。

    第三,對于GARCH模型、EGARCH模型和PGARCH模型,在采用參數(shù)法和非參數(shù)法分別對市場風險進行估計時,假定市場收益率服從GED分布的模型市場風險度量能力最弱,這表明在這些方法中GED分布不能準確的描述市場,SGT-VaR、FSH-VaR和GPD-VaR能更好地擬合市場。在置信水平99.5%下,半?yún)?shù)模型(FSH-VaR,GPD-VaR)比參數(shù)模型(SGT-VaR)風險度量能力好;在置信水平97.5%下,參數(shù)模型(SGT-VaR)比半?yún)?shù)模型(FSH-VaR,GPD-VaR)風險度量能力略好,但并不顯著;在置信水平95%下,參數(shù)模型(SGT-VaR)風險度量能力優(yōu)于半?yún)?shù)模型(FSH-VaR,GPD-VaR)??梢?,在參數(shù)模型和非參數(shù)模型都能較好的擬合市場時,置信水平越高,半?yún)?shù)模型(FSH-VaR,GPD-VaR)的風險度量能力越強,置信水平較低時,參數(shù)模型(SGT-VaR)的風險度量的準確程度優(yōu)于半?yún)?shù)模型。

    [1]Bali T,Theodossiou P.A Conditional-SGT-VaR Approach With Al?ternative GARCH Models[J].Annals of Research,2007,(151).

    [2]Balkema A,de Haan L.Residual Life at Great Age[J].Annals of Prob?ability,1974,(2).

    [3]Barone Adesi G,Giannopoulos K,Vosper L.VaR Without Correla?tions for Nonlinear Portfolios[J].Journal of Futures Markets,1999, (19).

    [4]Bekiros S,Georgoutsos D.Estimation of Value at Risk by Extreme Value and Conventional Methods:A Comparative Evaluation of Their Predictive Per-Formance[J].Journal of International Financial Mar?kets,Institutions&Money,2005,15(3).

    [5]Chen Q,Chen R.Method of Value-at-Risk and Empirical Research for Shanghai Stock Market[J].Procedia Computer Science,2013,(17).

    [6]Cheng W,Hung J.Skewness and Leptokurtosis in GARCH-typed VaR Estimation of Petroleum and Metal Asset Returns[J].Journal of Empirical Finance,2011,(18).

    [7]Christoffersen P.Evaluating Interval Forecasting[J].International Eco? nomic Review,1998,(38).

    [8]Engle R F,Manganelli S.CAViaR:Conditional Autoregressive Value at Risk by Regression Quantiles[J].Journal of Business and Economic Statistics,2004,(22).

    [9]Fisher R,Tippett L.Limiting Forms of the Frequency Distribution of The Largest or Smallest Member of a Sample[J].Proceedings of the Cambridge Philosophical Society,1928,(24).

    [10]Kupiec P.Techniques for Verifying the Accuracy of Risk Measure?ment Models[J].Journal of Derivatives,1995,(2).

    [11]Marimoutou V,Raggad B,Trabelsi A.Extreme Value Theory and Val?ue at Risk:Application to Oil Market[J].Energy Economics,2009, (31).

    [12]McNeil A,Frey R.Estimation of Tail-Related Risk Measures for Heteroscedasticity Financial Time Series:An Extreme Value Ap?proach[J].Journal of Empirical Finance,2000,(7).

    [13]Polanski A,Stoja E.Incorporating Higher Moments Into Val?ue-at-risk Forecasting[J].Journal of Forecasting,2010,(29).

    [14]Trenca I.The Use in Banks of VaR Method in Market Risk Manage?ment[J]Economic Sciences Series 2009.

    [15]Zikovic S,Aktan B.Global Financial Crisis and Var Performance in Emerging Markets:A Case of EU Candidate States-Turkey and Croatia[J].Journal of Economics and Business,2009,(27).

    [16]高瑩,周鑫,金秀.GARCH模型在動態(tài)VaR中的應(yīng)用[J].東北大學學報(自然科學報),2008,29(4).

    [17]徐煒,黃炎龍.GARCH模型和VaR的度量研究[J].數(shù)量經(jīng)濟技術(shù)經(jīng)濟研究,2008,(1).

    [18]朱慧明等.基于序數(shù)空間的行業(yè)間戰(zhàn)略相對風險度量[J].統(tǒng)計與決策,2012,(13).

    (責任編輯/易永生)

    F830

    A

    1002-6487(2016)23-0015-06

    全國統(tǒng)計科學研究重點項目(2015LZ54);安徽省高校自然科學重點基金項目(KJ2016A278);安徽省軟科學項目(1502052039);安徽師范大學哲學社會科學繁榮發(fā)展計劃首批重大項目(FRZD201302);2014年安徽師范大學特色優(yōu)勢研究領(lǐng)域建設(shè)項目

    張 瓊(1976—),女,安徽樅陽人,博士,副教授,研究方向:經(jīng)濟統(tǒng)計。

    猜你喜歡
    置信水平度量收益
    有趣的度量
    模糊度量空間的強嵌入
    螃蟹爬上“網(wǎng)” 收益落進兜
    產(chǎn)品控制與市場風險之間的相互作用研究
    迷向表示分為6個不可約直和的旗流形上不變愛因斯坦度量
    單因子方差分析法在卷煙均勻性檢驗中的研究與應(yīng)用
    用VaR方法分析中國A股市場的風險
    2015年理財“6宗最”誰能給你穩(wěn)穩(wěn)的收益
    金色年華(2016年1期)2016-02-28 01:38:19
    地質(zhì)異常的奇異性度量與隱伏源致礦異常識別
    東芝驚爆會計丑聞 憑空捏造1518億日元收益
    IT時代周刊(2015年8期)2015-11-11 05:50:38
    在线观看66精品国产| 搡女人真爽免费视频火全软件 | 久久九九热精品免费| 狂野欧美白嫩少妇大欣赏| av黄色大香蕉| 岛国在线观看网站| 亚洲国产精品999在线| 国产伦精品一区二区三区视频9 | 搡老熟女国产l中国老女人| 国产一级毛片七仙女欲春2| 国产高清激情床上av| 亚洲美女视频黄频| 69人妻影院| av片东京热男人的天堂| 日韩欧美在线乱码| 免费看a级黄色片| 女生性感内裤真人,穿戴方法视频| 国产黄色小视频在线观看| 一级作爱视频免费观看| 国产探花在线观看一区二区| 成人国产一区最新在线观看| 欧美性猛交╳xxx乱大交人| 国产淫片久久久久久久久 | 日韩欧美在线乱码| 午夜影院日韩av| 欧美一区二区国产精品久久精品| 在线观看免费午夜福利视频| 成人一区二区视频在线观看| 成年女人看的毛片在线观看| 免费高清视频大片| 成人性生交大片免费视频hd| 小蜜桃在线观看免费完整版高清| 国产av在哪里看| 九九热线精品视视频播放| 精品熟女少妇八av免费久了| 色综合婷婷激情| 一进一出抽搐gif免费好疼| 久久久久久久久久黄片| 日韩亚洲欧美综合| 在线国产一区二区在线| av天堂在线播放| 淫秽高清视频在线观看| 欧美中文日本在线观看视频| 亚洲va日本ⅴa欧美va伊人久久| 色在线成人网| 色播亚洲综合网| 欧美极品一区二区三区四区| 又紧又爽又黄一区二区| 波多野结衣高清作品| 桃红色精品国产亚洲av| 成人18禁在线播放| 亚洲人成网站在线播放欧美日韩| 丁香欧美五月| 最近最新免费中文字幕在线| 校园春色视频在线观看| 久久香蕉精品热| 欧美午夜高清在线| a在线观看视频网站| 午夜久久久久精精品| av天堂中文字幕网| 搡女人真爽免费视频火全软件 | 久久精品亚洲精品国产色婷小说| 久久精品夜夜夜夜夜久久蜜豆| 成人午夜高清在线视频| 天堂影院成人在线观看| 国产99白浆流出| 88av欧美| 国产中年淑女户外野战色| 亚洲国产欧美网| 久久久久亚洲av毛片大全| 90打野战视频偷拍视频| 成年女人毛片免费观看观看9| 又黄又粗又硬又大视频| 国内揄拍国产精品人妻在线| 久久性视频一级片| 国产亚洲精品av在线| av中文乱码字幕在线| 欧美性感艳星| svipshipincom国产片| 婷婷精品国产亚洲av| 国产精品久久久人人做人人爽| 一区二区三区免费毛片| 免费人成在线观看视频色| 国产精品综合久久久久久久免费| 国内毛片毛片毛片毛片毛片| 成人18禁在线播放| 91久久精品电影网| 免费看美女性在线毛片视频| 欧洲精品卡2卡3卡4卡5卡区| 一区二区三区高清视频在线| 在线观看美女被高潮喷水网站 | 精品一区二区三区av网在线观看| 国产精品久久久久久亚洲av鲁大| 狠狠狠狠99中文字幕| 欧美日韩亚洲国产一区二区在线观看| 欧美三级亚洲精品| 国产一区二区激情短视频| 欧美日韩国产亚洲二区| 亚洲精品日韩av片在线观看 | 国产亚洲精品综合一区在线观看| 欧美bdsm另类| 成人无遮挡网站| 国产aⅴ精品一区二区三区波| 欧美三级亚洲精品| 国产三级黄色录像| 首页视频小说图片口味搜索| 日本成人三级电影网站| 午夜免费成人在线视频| 97碰自拍视频| 亚洲乱码一区二区免费版| 国产男靠女视频免费网站| 亚洲片人在线观看| 深夜精品福利| 亚洲无线在线观看| 免费看美女性在线毛片视频| 久久天躁狠狠躁夜夜2o2o| 亚洲最大成人手机在线| 国产av麻豆久久久久久久| 国产三级中文精品| 一个人观看的视频www高清免费观看| 欧美中文综合在线视频| 人妻夜夜爽99麻豆av| 久久久久久久午夜电影| 国产精品久久电影中文字幕| 在线a可以看的网站| 久久久久久久久大av| 淫秽高清视频在线观看| 国产亚洲av嫩草精品影院| 热99re8久久精品国产| 亚洲国产色片| 成熟少妇高潮喷水视频| 亚洲人成网站高清观看| 91在线观看av| 色哟哟哟哟哟哟| 国产伦精品一区二区三区视频9 | 又爽又黄无遮挡网站| 精品福利观看| 国产视频一区二区在线看| 国产成人系列免费观看| 成人特级av手机在线观看| 国产亚洲精品久久久com| 国产精品永久免费网站| 美女被艹到高潮喷水动态| 1000部很黄的大片| 俄罗斯特黄特色一大片| 国产 一区 欧美 日韩| 亚洲成人久久性| 中文字幕人妻丝袜一区二区| 亚洲成a人片在线一区二区| 欧美激情久久久久久爽电影| 美女高潮喷水抽搐中文字幕| www国产在线视频色| 精品不卡国产一区二区三区| 90打野战视频偷拍视频| 色尼玛亚洲综合影院| 午夜两性在线视频| 欧美日韩乱码在线| 亚洲国产精品久久男人天堂| 午夜两性在线视频| 少妇裸体淫交视频免费看高清| 亚洲人成网站在线播| 中文亚洲av片在线观看爽| 欧美一区二区国产精品久久精品| 亚洲男人的天堂狠狠| 久久精品91蜜桃| 国产精品日韩av在线免费观看| 三级毛片av免费| 日本熟妇午夜| 国产欧美日韩精品亚洲av| 国产视频内射| 亚洲电影在线观看av| 久久久久久大精品| 制服人妻中文乱码| 国产真实伦视频高清在线观看 | а√天堂www在线а√下载| 国产高清videossex| 国产爱豆传媒在线观看| 精品久久久久久久末码| 少妇的逼水好多| 一a级毛片在线观看| 亚洲av熟女| 午夜两性在线视频| 黄色日韩在线| 好男人电影高清在线观看| 久久精品国产清高在天天线| 天天躁日日操中文字幕| 99久久成人亚洲精品观看| 伊人久久精品亚洲午夜| 国产野战对白在线观看| 国产精品自产拍在线观看55亚洲| 操出白浆在线播放| 国产探花在线观看一区二区| 午夜两性在线视频| 国产精品日韩av在线免费观看| 日韩有码中文字幕| 黄色丝袜av网址大全| 好看av亚洲va欧美ⅴa在| 精品久久久久久,| 99国产综合亚洲精品| 18禁黄网站禁片午夜丰满| 99热这里只有精品一区| 在线十欧美十亚洲十日本专区| 色在线成人网| 亚洲av电影在线进入| 午夜激情欧美在线| 亚洲不卡免费看| 免费看光身美女| 成人国产一区最新在线观看| 久久久精品欧美日韩精品| 在线观看日韩欧美| 一本一本综合久久| 麻豆一二三区av精品| 悠悠久久av| 免费电影在线观看免费观看| 国产精品久久久久久人妻精品电影| 欧美成人免费av一区二区三区| 国产探花在线观看一区二区| 动漫黄色视频在线观看| 美女黄网站色视频| 一边摸一边抽搐一进一小说| 亚洲精品456在线播放app | 成人av一区二区三区在线看| 久久这里只有精品中国| 观看美女的网站| 一进一出抽搐动态| 国产成人av激情在线播放| 黄色日韩在线| 国产精品av视频在线免费观看| 超碰av人人做人人爽久久 | 亚洲最大成人手机在线| 精品无人区乱码1区二区| 香蕉av资源在线| 51午夜福利影视在线观看| 精品一区二区三区视频在线 | 中文亚洲av片在线观看爽| 最好的美女福利视频网| 亚洲av第一区精品v没综合| 久久久久精品国产欧美久久久| 欧美最新免费一区二区三区 | 成人精品一区二区免费| 12—13女人毛片做爰片一| 国产不卡一卡二| 高清日韩中文字幕在线| 亚洲一区二区三区色噜噜| 国产三级在线视频| 久久久久久国产a免费观看| 性色av乱码一区二区三区2| 欧美国产日韩亚洲一区| 亚洲色图av天堂| 久久精品国产亚洲av香蕉五月| 亚洲精品影视一区二区三区av| svipshipincom国产片| 国产成+人综合+亚洲专区| 成人三级黄色视频| 亚洲av电影不卡..在线观看| 国内精品久久久久精免费| 搞女人的毛片| 热99在线观看视频| 日韩高清综合在线| 男女那种视频在线观看| 亚洲电影在线观看av| 热99在线观看视频| 久久性视频一级片| 亚洲av成人精品一区久久| 动漫黄色视频在线观看| 久久久久久久午夜电影| 欧美极品一区二区三区四区| 成人国产综合亚洲| 国产一区二区在线观看日韩 | 欧美乱妇无乱码| 欧美中文综合在线视频| 国产精品一区二区三区四区免费观看 | 国产久久久一区二区三区| 可以在线观看毛片的网站| 成年女人看的毛片在线观看| 国产一区二区在线观看日韩 | 亚洲精品乱码久久久v下载方式 | 久久久久九九精品影院| 国内精品久久久久精免费| 久久精品人妻少妇| 在线观看美女被高潮喷水网站 | 18禁黄网站禁片午夜丰满| 国产亚洲精品综合一区在线观看| 3wmmmm亚洲av在线观看| 欧美xxxx黑人xx丫x性爽| 欧美一级a爱片免费观看看| 伊人久久大香线蕉亚洲五| 少妇裸体淫交视频免费看高清| 亚洲男人的天堂狠狠| 国产精品99久久99久久久不卡| 一个人免费在线观看的高清视频| 日日摸夜夜添夜夜添小说| 精品国产美女av久久久久小说| 久久久国产成人精品二区| 亚洲专区中文字幕在线| 好男人在线观看高清免费视频| 亚洲色图av天堂| 天美传媒精品一区二区| 国产成+人综合+亚洲专区| 床上黄色一级片| 色尼玛亚洲综合影院| 亚洲精品成人久久久久久| 国产日本99.免费观看| 一进一出抽搐gif免费好疼| 嫁个100分男人电影在线观看| 色老头精品视频在线观看| 欧美一区二区亚洲| 在线观看舔阴道视频| 美女高潮喷水抽搐中文字幕| 亚洲国产精品久久男人天堂| 国产精品av视频在线免费观看| 国产v大片淫在线免费观看| 国内少妇人妻偷人精品xxx网站| 欧美在线一区亚洲| 国产毛片a区久久久久| 听说在线观看完整版免费高清| 日日摸夜夜添夜夜添小说| 日本精品一区二区三区蜜桃| 国产色爽女视频免费观看| 日本 av在线| 少妇的丰满在线观看| 午夜福利在线在线| 97超级碰碰碰精品色视频在线观看| 悠悠久久av| 给我免费播放毛片高清在线观看| 999久久久精品免费观看国产| 久久久久国内视频| 久久国产精品影院| 在线观看美女被高潮喷水网站 | 国产日本99.免费观看| 波野结衣二区三区在线 | 在线免费十八禁| 九草在线视频观看| xxx大片免费视频| 亚洲美女搞黄在线观看| 在线观看一区二区三区| 久热久热在线精品观看| 久久精品国产自在天天线| 欧美3d第一页| 中文欧美无线码| 国产视频内射| 婷婷色麻豆天堂久久| 亚洲人与动物交配视频| 久久久午夜欧美精品| 亚洲欧美一区二区三区国产| 神马国产精品三级电影在线观看| 最近的中文字幕免费完整| 人妻夜夜爽99麻豆av| 午夜福利视频1000在线观看| 免费观看在线日韩| 99久久精品国产国产毛片| 高清午夜精品一区二区三区| 美女cb高潮喷水在线观看| 熟女电影av网| 欧美xxxx黑人xx丫x性爽| 国产一区二区三区综合在线观看 | 久久精品综合一区二区三区| 欧美变态另类bdsm刘玥| 乱人视频在线观看| 日韩成人伦理影院| 国产黄片视频在线免费观看| 国产成人精品福利久久| 91久久精品国产一区二区成人| 内射极品少妇av片p| 一个人免费在线观看电影| 久久亚洲国产成人精品v| 中国国产av一级| 亚洲国产日韩欧美精品在线观看| 国产精品精品国产色婷婷| 亚洲精品一二三| 听说在线观看完整版免费高清| 亚洲欧美成人精品一区二区| 深爱激情五月婷婷| 国产一区二区三区综合在线观看 | 偷拍熟女少妇极品色| 91精品一卡2卡3卡4卡| 亚洲av电影在线观看一区二区三区 | 国产伦一二天堂av在线观看| 免费观看性生交大片5| 久久鲁丝午夜福利片| 免费观看性生交大片5| 久久鲁丝午夜福利片| 69av精品久久久久久| 欧美激情久久久久久爽电影| 欧美xxxx黑人xx丫x性爽| 精品国产一区二区三区久久久樱花 | 22中文网久久字幕| 日韩一区二区三区影片| 精品人妻视频免费看| 男人狂女人下面高潮的视频| 热99在线观看视频| 久久久a久久爽久久v久久| 97超视频在线观看视频| 可以在线观看毛片的网站| 欧美另类一区| 国产精品99久久久久久久久| 国产av国产精品国产| 国产精品一区二区三区四区免费观看| 成人无遮挡网站| 肉色欧美久久久久久久蜜桃 | 啦啦啦韩国在线观看视频| 亚洲内射少妇av| 小蜜桃在线观看免费完整版高清| 亚洲av成人精品一区久久| 免费av毛片视频| 日韩欧美国产在线观看| 一级片'在线观看视频| 国产探花在线观看一区二区| 嫩草影院精品99| 久久久久久久久中文| 18禁在线播放成人免费| 久久鲁丝午夜福利片| 午夜福利视频精品| 99视频精品全部免费 在线| 亚洲四区av| 日韩欧美 国产精品| 欧美xxxx性猛交bbbb| 国产精品国产三级专区第一集| 女人被狂操c到高潮| 日韩伦理黄色片| av黄色大香蕉| 国产成人福利小说| 国产伦在线观看视频一区| 国精品久久久久久国模美| 夫妻性生交免费视频一级片| 精品欧美国产一区二区三| 成人毛片a级毛片在线播放| 99热全是精品| 男人舔女人下体高潮全视频| 国产亚洲av片在线观看秒播厂 | 成人亚洲精品av一区二区| 看非洲黑人一级黄片| 高清在线视频一区二区三区| 少妇的逼好多水| 日韩在线高清观看一区二区三区| 日本爱情动作片www.在线观看| 一个人看的www免费观看视频| 插阴视频在线观看视频| 天天躁夜夜躁狠狠久久av| 国产精品一二三区在线看| 国产 一区精品| av线在线观看网站| 高清av免费在线| 国产亚洲精品av在线| 国精品久久久久久国模美| 精品不卡国产一区二区三区| 能在线免费观看的黄片| 只有这里有精品99| 免费黄网站久久成人精品| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 91精品一卡2卡3卡4卡| 久久99热这里只频精品6学生| 波野结衣二区三区在线| 国产淫片久久久久久久久| 免费高清在线观看视频在线观看| 亚洲av电影不卡..在线观看| 亚洲自偷自拍三级| 免费av观看视频| 女人十人毛片免费观看3o分钟| h日本视频在线播放| 午夜精品在线福利| 日韩精品青青久久久久久| 国产精品久久久久久久久免| 高清日韩中文字幕在线| 少妇的逼水好多| 国产av国产精品国产| 插逼视频在线观看| 久久精品夜色国产| 少妇丰满av| 久久久a久久爽久久v久久| 久久久久久国产a免费观看| 亚洲在久久综合| 日韩欧美精品免费久久| av网站免费在线观看视频 | 国产爱豆传媒在线观看| 国产v大片淫在线免费观看| 一级毛片 在线播放| 亚洲av中文字字幕乱码综合| 免费观看无遮挡的男女| 精品熟女少妇av免费看| 国产视频内射| 亚洲综合色惰| 人人妻人人看人人澡| 久久99热这里只频精品6学生| 九九在线视频观看精品| 三级毛片av免费| 久久精品国产自在天天线| 又黄又爽又刺激的免费视频.| 国产亚洲最大av| 伦理电影大哥的女人| 高清欧美精品videossex| 日韩av不卡免费在线播放| 久久久久久久大尺度免费视频| 久久久久久久久久黄片| 国产午夜精品论理片| 18禁在线无遮挡免费观看视频| 国产女主播在线喷水免费视频网站 | 欧美三级亚洲精品| 国产精品一二三区在线看| 日韩制服骚丝袜av| 成人欧美大片| 成人午夜高清在线视频| 两个人的视频大全免费| 晚上一个人看的免费电影| 亚洲精品色激情综合| 国产伦理片在线播放av一区| 一区二区三区乱码不卡18| 日韩精品青青久久久久久| 免费黄网站久久成人精品| 九色成人免费人妻av| 欧美日本视频| 一区二区三区免费毛片| 亚洲av国产av综合av卡| 国产在线一区二区三区精| 97在线视频观看| 少妇高潮的动态图| 人妻少妇偷人精品九色| 免费播放大片免费观看视频在线观看| 可以在线观看毛片的网站| 国精品久久久久久国模美| 国产一级毛片七仙女欲春2| 亚洲一区高清亚洲精品| av免费在线看不卡| 色5月婷婷丁香| 日韩欧美精品免费久久| 只有这里有精品99| 国产综合懂色| 激情 狠狠 欧美| 人人妻人人澡人人爽人人夜夜 | 午夜福利在线观看免费完整高清在| 99热6这里只有精品| 亚洲伊人久久精品综合| 亚洲精品一二三| 最后的刺客免费高清国语| 汤姆久久久久久久影院中文字幕 | 亚洲成人中文字幕在线播放| 大又大粗又爽又黄少妇毛片口| 两个人视频免费观看高清| 永久免费av网站大全| 成人亚洲精品av一区二区| 国产黄a三级三级三级人| 国产精品一区二区三区四区免费观看| 乱系列少妇在线播放| 亚洲欧洲国产日韩| 午夜免费男女啪啪视频观看| 午夜老司机福利剧场| 午夜日本视频在线| 极品少妇高潮喷水抽搐| 免费看美女性在线毛片视频| 寂寞人妻少妇视频99o| 欧美 日韩 精品 国产| 一边亲一边摸免费视频| 国产毛片a区久久久久| 精品人妻熟女av久视频| 日韩不卡一区二区三区视频在线| 人人妻人人看人人澡| 精品久久久久久久末码| 国产极品天堂在线| 中文在线观看免费www的网站| 偷拍熟女少妇极品色| 亚洲欧美精品自产自拍| 中文字幕亚洲精品专区| 色哟哟·www| 欧美最新免费一区二区三区| 国产精品一二三区在线看| 爱豆传媒免费全集在线观看| 久久鲁丝午夜福利片| 成人性生交大片免费视频hd| 九九久久精品国产亚洲av麻豆| kizo精华| 99久国产av精品国产电影| 3wmmmm亚洲av在线观看| 午夜视频国产福利| 亚洲高清免费不卡视频| freevideosex欧美| 久久久a久久爽久久v久久| 亚洲精品久久午夜乱码| 1000部很黄的大片| 丝瓜视频免费看黄片| 听说在线观看完整版免费高清| 日本-黄色视频高清免费观看| 又爽又黄无遮挡网站| 亚洲内射少妇av| av福利片在线观看| 三级国产精品片| 国产精品国产三级国产av玫瑰| 国产成人freesex在线| 亚洲成人久久爱视频| 国产精品一及| 日韩,欧美,国产一区二区三区| 亚洲人与动物交配视频| 日韩大片免费观看网站| 人妻制服诱惑在线中文字幕| 少妇人妻精品综合一区二区| 日韩大片免费观看网站| 久久鲁丝午夜福利片| 午夜激情福利司机影院| 久久久久久久久久黄片| 99久久九九国产精品国产免费| 99热全是精品| 国产男人的电影天堂91| 亚洲成色77777| 亚洲一区高清亚洲精品| 日本猛色少妇xxxxx猛交久久| 欧美日韩综合久久久久久| 久久草成人影院| 精品一区在线观看国产| 日本-黄色视频高清免费观看| 精品久久久久久成人av| 一级二级三级毛片免费看| 永久免费av网站大全| 高清在线视频一区二区三区| 久久99热这里只频精品6学生| 男女边摸边吃奶| 一级爰片在线观看|