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

    時(shí)變參數(shù)模型的最優(yōu)滾動(dòng)窗寬選擇標(biāo)準(zhǔn)及應(yīng)用

    2018-09-10 07:15:04張興敏
    中國(guó)管理科學(xué) 2018年8期
    關(guān)鍵詞:估計(jì)量均方斷點(diǎn)

    傅 強(qiáng),張興敏

    (重慶大學(xué)經(jīng)濟(jì)與工商管理學(xué)院,重慶 400030)

    1 引言

    結(jié)構(gòu)不穩(wěn)定性是研究領(lǐng)域中不可忽視的問(wèn)題,靜態(tài)模型難以捕捉樣本關(guān)系的結(jié)構(gòu)性變化[1-2],對(duì)結(jié)構(gòu)不穩(wěn)定性的實(shí)證研究在金融風(fēng)險(xiǎn)測(cè)度、預(yù)警[3-8],格蘭杰因果關(guān)系檢驗(yàn)[2,9]、宏觀經(jīng)濟(jì)預(yù)測(cè)[10-12]等領(lǐng)域已有較大發(fā)展。為捕捉結(jié)構(gòu)關(guān)系中的不穩(wěn)定性,使用時(shí)變參數(shù)模型刻畫(huà)經(jīng)濟(jì)序列間關(guān)系的方法被廣泛使用,而參數(shù)的時(shí)變性常通過(guò)使用最近的觀測(cè)值估計(jì)模型,也即滾動(dòng)樣本估計(jì)法實(shí)現(xiàn)。使用滾動(dòng)樣本技術(shù)的原因有:第一,該方法認(rèn)為各變量間的關(guān)系隨時(shí)間而變化;第二,經(jīng)濟(jì)序列中存在可能的結(jié)構(gòu)突變,而滾動(dòng)樣本法能夠捕捉到這種變化過(guò)程。

    滾動(dòng)樣本技術(shù)是指在整個(gè)樣本期內(nèi)使用固定長(zhǎng)度的滾動(dòng)樣本序列估計(jì)模型,得到一系列隨時(shí)間變化的參數(shù)估計(jì)量,進(jìn)一步地,可根據(jù)不同時(shí)點(diǎn)的估計(jì)模型進(jìn)行樣本外預(yù)測(cè)。選擇多大規(guī)模的窗口寬度進(jìn)行估計(jì),是實(shí)踐過(guò)程中需要解決的關(guān)鍵問(wèn)題。

    盡管關(guān)于滾動(dòng)窗寬的選擇問(wèn)題類似非參估計(jì)中的帶寬(bandwidth)選擇問(wèn)題,但對(duì)于使用滾動(dòng)樣本技術(shù)捕捉時(shí)變參數(shù)或進(jìn)行樣本外預(yù)測(cè)的窗寬選擇問(wèn)題卻受到較少關(guān)注,且在應(yīng)用過(guò)程中由于研究對(duì)象和目標(biāo)有所差異。預(yù)測(cè)模型中,滾動(dòng)樣本技術(shù)的目標(biāo)是使預(yù)測(cè)精度更高,而較少關(guān)注模型估計(jì)量的準(zhǔn)確性,以及模型參數(shù)的時(shí)變性。時(shí)變參數(shù)模型中,滾動(dòng)樣本技術(shù)是為了捕捉參數(shù)隨時(shí)間的變化趨勢(shì),更關(guān)注模型估計(jì)量的精確度及參數(shù)在不同子樣本序列中的變化情況。當(dāng)樣本序列中出現(xiàn)一個(gè)或多個(gè)離散斷點(diǎn)時(shí),Pesaran和Timmermann[13]以提高預(yù)測(cè)精度為目的,提出了五種方法選擇滾動(dòng)窗寬。他們的方法中需要首先運(yùn)用斷點(diǎn)檢驗(yàn)法(如,chow檢驗(yàn),累積求和統(tǒng)計(jì)量檢驗(yàn)法(CUSUM and CUSUM-of-squares test)等)識(shí)別出未知斷點(diǎn)的位置,并對(duì)斷點(diǎn)前后的樣本量進(jìn)行取舍,進(jìn)而估計(jì)出預(yù)測(cè)模型。當(dāng)樣本序列中出現(xiàn)多個(gè)連續(xù)和離散斷點(diǎn)時(shí),Pesaran等[14]推導(dǎo)了最優(yōu)權(quán)重法確定預(yù)測(cè)模型所用樣本量;Giraitis等[15]則基于交叉驗(yàn)證思想選擇平滑參數(shù),以減少早期信息的權(quán)重,進(jìn)而得到更精確的預(yù)測(cè)量。而當(dāng)斷點(diǎn)信息未知時(shí),Inoue等[16]提出以最小化近似預(yù)測(cè)均方誤差選擇最優(yōu)滾動(dòng)窗寬的方法。這些方法僅從預(yù)測(cè)視角出發(fā),以預(yù)測(cè)精度為研究目的,且多以線性回歸模型為研究對(duì)象。對(duì)于以捕捉參數(shù)時(shí)變性為目的的模型,它們選擇的滾動(dòng)窗寬無(wú)法達(dá)到提高模型估計(jì)量的準(zhǔn)確性,捕捉參數(shù)結(jié)構(gòu)突變性的目標(biāo)。因此,預(yù)測(cè)模型下的最優(yōu)窗寬選擇標(biāo)準(zhǔn)對(duì)于時(shí)變參數(shù)模型可能并不適用。

    時(shí)變參數(shù)模型對(duì)滾動(dòng)樣本技術(shù)的應(yīng)用仍處于發(fā)展階段,滾動(dòng)窗寬的選擇問(wèn)題少有文獻(xiàn)探討。現(xiàn)有文獻(xiàn)中,實(shí)驗(yàn)者或隨意選擇,或基于經(jīng)驗(yàn),或以年、季、月等時(shí)間截點(diǎn)選擇滾動(dòng)窗寬[4,6-7,17-19],并沒(méi)有形成統(tǒng)一明確的準(zhǔn)則。實(shí)證表明模型的估計(jì)結(jié)果對(duì)于窗寬是敏感的,Nyakabawo等[2]認(rèn)為滾動(dòng)窗寬選擇需權(quán)衡“模型估計(jì)的準(zhǔn)確性和模型在各子樣本期間的代表性”兩個(gè)目標(biāo)?;诖怂枷?,他們選擇了多個(gè)滾動(dòng)窗寬,比較分析不同窗寬下變量的顯著性水平,發(fā)現(xiàn)結(jié)果有顯著差異,但沒(méi)有建立最優(yōu)窗寬的統(tǒng)計(jì)準(zhǔn)則。Pesaran和Timmermann[20]證明了,當(dāng)子樣本的規(guī)模較小時(shí),在子樣本內(nèi)包含多個(gè)模型的風(fēng)險(xiǎn)會(huì)隨之減小。但較小的子樣本規(guī)模卻意味著參數(shù)估計(jì)的準(zhǔn)確性隨著標(biāo)準(zhǔn)誤增加而下降。Khediri和Charfeddine[17]使用滾動(dòng)樣本法分析了現(xiàn)貨和期貨能源市場(chǎng)有效性的時(shí)變性,他們認(rèn)為通過(guò)結(jié)構(gòu)斷點(diǎn)檢驗(yàn)方法識(shí)別極端事件日期的做法并不合適,且多有文獻(xiàn)進(jìn)行批評(píng)。與此相反,Charfeddine和Benlagha[19]運(yùn)用滾動(dòng)樣本法分析了copula模型中相依參數(shù)的時(shí)變行為,并使用結(jié)構(gòu)斷點(diǎn)檢驗(yàn)法識(shí)別出極端政治、金融事件。他們認(rèn)為,當(dāng)數(shù)據(jù)集中出現(xiàn)一些斷點(diǎn)時(shí),短的滾動(dòng)窗寬相比長(zhǎng)的滾動(dòng)窗寬能更好地捕捉參數(shù)時(shí)變性。Anagnostidis等[21]也認(rèn)為滾動(dòng)窗寬的選擇對(duì)于模型的估計(jì)結(jié)果會(huì)產(chǎn)生較大影響,但他們沒(méi)有建立相應(yīng)的選擇標(biāo)準(zhǔn),而是選擇多個(gè)時(shí)間窗口對(duì)比分析股票市場(chǎng)效率隨時(shí)間的變化情況。可見(jiàn),對(duì)于滾動(dòng)樣本的研究逐漸受到重視,但關(guān)于如何選擇滾動(dòng)窗寬仍缺乏深入的研究。

    與現(xiàn)有研究的不同之處在于:第一,滾動(dòng)樣本技術(shù)是當(dāng)前文獻(xiàn)中常用的動(dòng)態(tài)方法,但罕見(jiàn)文獻(xiàn)探討滾動(dòng)窗寬的選擇準(zhǔn)則,本文以捕捉參數(shù)結(jié)構(gòu)突變性為目的,建立了選擇滾動(dòng)窗寬的統(tǒng)計(jì)準(zhǔn)則,彌補(bǔ)了當(dāng)前該領(lǐng)域的空白。第二,在模型設(shè)定上,并不假定某一具體模型形式。將數(shù)據(jù)生成過(guò)程設(shè)置為單指標(biāo)半?yún)?shù)模型,既可捕捉變量間的線性關(guān)系,也可捕捉非線性關(guān)系,假設(shè)條件較弱,適用性強(qiáng)。同時(shí),當(dāng)指標(biāo)函數(shù)(鏈接函數(shù))已知時(shí),可拓展至通常情況下的線性回歸模型和廣義線性回歸模型,此時(shí)只需使用對(duì)應(yīng)模型的估計(jì)方法即可。因此,本文方法在應(yīng)用范疇上較為廣泛。第三,通過(guò)最小化估計(jì)量的近似二次損失函數(shù)及最大化各子樣本估計(jì)量間的曼哈頓距離選擇窗寬大小,同時(shí)考慮了模型估計(jì)量的準(zhǔn)確性和時(shí)變性兩個(gè)相悖目標(biāo)。與傳統(tǒng)上,僅根據(jù)主觀意愿以時(shí)間截點(diǎn)選擇滾動(dòng)窗寬,以及與以預(yù)測(cè)精度為目標(biāo)選擇的滾動(dòng)窗寬相比,在模型估計(jì)效果上有顯著改善。第四,蒙特卡羅模擬實(shí)驗(yàn)和參數(shù)敏感性分析表明,本文方法能夠捕捉到參數(shù)的結(jié)構(gòu)突變性,且能滿足參數(shù)估計(jì)量的準(zhǔn)確性目標(biāo),證明了本文方法的有效性和穩(wěn)健性。第五,將滾動(dòng)窗寬選擇標(biāo)準(zhǔn)運(yùn)用到我國(guó)金融系統(tǒng)網(wǎng)絡(luò)結(jié)構(gòu)的構(gòu)建中,捕捉金融網(wǎng)絡(luò)的結(jié)構(gòu)突變性,優(yōu)化了動(dòng)態(tài)網(wǎng)絡(luò)模型。

    2 研究設(shè)計(jì)

    2.1 模型設(shè)定

    假設(shè)數(shù)據(jù)生成過(guò)程(DGP)如下:

    (1)

    下文僅以單指標(biāo)半?yún)?shù)模型為例,介紹局部線性迭代估計(jì)過(guò)程[22-24]。注意,以下參數(shù)中沒(méi)有時(shí)間符號(hào),原因在于這里的模型固定在一個(gè)時(shí)間窗口內(nèi),動(dòng)態(tài)性由滾動(dòng)樣本技術(shù)捕捉,單個(gè)時(shí)間窗口的樣本量為T1。首先對(duì)單指標(biāo)半?yún)?shù)模型有以下假設(shè):f(·)存在連續(xù)的二階微分,對(duì)u=θTXt的某鄰域內(nèi)的ν,利用局部線性函數(shù)f(ν)可近似表示為f(ν)≈f(u)+f′(u)(ν-u)≡a+b(ν-u),其中a=f(u),b=f′(u)。并記Kh(·)=K(·/h)h-1,K(·)是核函數(shù),滿足條件:非負(fù)、有界、關(guān)于0對(duì)稱且緊支撐以及Lipschitz連續(xù),h>0為非參估計(jì)中的窗寬。則關(guān)于指標(biāo)函數(shù)和指標(biāo)參數(shù)的估計(jì)過(guò)程如下:

    (2)

    (3)

    2.2 滾動(dòng)窗寬選擇標(biāo)準(zhǔn)

    滾動(dòng)樣本技術(shù)中選擇滾動(dòng)窗寬時(shí),需權(quán)衡兩個(gè)相悖目標(biāo)[2]:增加模型估計(jì)量的準(zhǔn)確性,同時(shí)使各子樣本期間模型的異構(gòu)性最高。一方面,較小的子樣本規(guī)模減少了子樣本期內(nèi)模型的異構(gòu)性,提高了子樣本內(nèi)模型估計(jì)量的代表性,但增加了估計(jì)量的標(biāo)準(zhǔn)誤差,從而降低估計(jì)量的精度。另一方面,較大的子樣本規(guī)模提高了估計(jì)量的準(zhǔn)確性,但降低了子樣本期內(nèi)所估計(jì)模型的代表性,特別是在存在結(jié)構(gòu)突變的情況下。基于此,本文從理論上建立滾動(dòng)窗寬的選擇標(biāo)準(zhǔn)。以下為具體準(zhǔn)則:

    均方誤差MSE表征了模型的估計(jì)優(yōu)劣,可用于度量估計(jì)量的準(zhǔn)確性;曼哈頓距離表征了變量之間的差異,可用于度量各子樣本期間模型的異構(gòu)性。在單指標(biāo)半?yún)?shù)模型框架(式(1))下,指標(biāo)函數(shù)和指標(biāo)參數(shù)均未知,需要對(duì)其進(jìn)行估計(jì)。因此,本文選擇這兩類估計(jì)量的均方誤差之和作為模型估計(jì)優(yōu)劣的判斷標(biāo)準(zhǔn)。在線性回歸框架或廣義線性回歸框架下,則只包含參數(shù)估計(jì)量的均方誤差。設(shè)總樣本量為T,滾動(dòng)窗寬為l,每個(gè)子窗口的樣本量為N=l+1,子窗口數(shù)目為T-l,在第m窗口下的兩類估計(jì)量的均方誤差定義如下:

    指標(biāo)函數(shù)f(·)的估計(jì)優(yōu)劣采用均方根誤差(RASE)表征,定義為:

    RASE(m)=

    (4)

    指標(biāo)參數(shù)θ的估計(jì)優(yōu)劣采用參數(shù)均方根誤差(PMSE)表征,定義為:

    (5)

    則最優(yōu)窗寬準(zhǔn)則定義為:

    (6)

    (7)

    第五步:重復(fù)第二步至第四步K次,得到指標(biāo)函數(shù)和參數(shù)估計(jì)量的Bootstrap值。

    則均方根誤差(RASE)的Bootstrap估計(jì)為:

    (8)

    參數(shù)均方誤差(PMSE)的Bootstrap估計(jì)為:

    (9)

    則模型的Bootstrap均方誤差(MSE)為:

    (10)

    本文中最優(yōu)滾動(dòng)窗寬的算法流程圖見(jiàn)圖1所示。本文算法的整體思路是,給定一個(gè)初始滾動(dòng)窗寬值,采用局部線性迭代估計(jì)法估計(jì)模型,進(jìn)而對(duì)樣本使用Bootstrap法,得到模型的Bootstrap估計(jì),求得近似均方誤差值和曼哈頓距離值。對(duì)所有可行滾動(dòng)窗寬執(zhí)行相同的過(guò)程。最后再使用隨機(jī)權(quán)重法,求得不同滾動(dòng)窗寬下的評(píng)價(jià)函數(shù),選擇最小評(píng)價(jià)函數(shù)對(duì)應(yīng)的滾動(dòng)窗寬作為最終的最優(yōu)滾動(dòng)窗寬。

    圖1 整體算法流程圖

    根據(jù)上述算法可知,第一步的時(shí)間復(fù)雜度為O((T-l))(局部線性迭代估計(jì),一般只關(guān)注收斂速率,不關(guān)注時(shí)間復(fù)雜度),第二步的時(shí)間復(fù)雜度為O((T-l)*K),第三步的時(shí)間復(fù)雜度為O(1),所有可行滾動(dòng)窗寬的循環(huán)步驟的時(shí)間復(fù)雜度為O(widths*(T-l)*K),第四步的時(shí)間復(fù)雜度為O(widths*(T-l)*K*W),則最終求出最優(yōu)滾動(dòng)窗寬的算法的時(shí)間復(fù)雜度為O(widths*(T-l)*K*W)。其中,widths為可行的滾動(dòng)窗寬個(gè)數(shù),(T-l)為選擇滾動(dòng)窗寬為l時(shí)的窗口個(gè)數(shù),K為Bootstrap法中重復(fù)次數(shù),W為隨機(jī)權(quán)重個(gè)數(shù)。可以發(fā)現(xiàn),時(shí)間復(fù)雜度O(nk)型,因此為有效算法。

    3 蒙特卡羅模擬

    本部分進(jìn)行蒙特卡羅模擬實(shí)驗(yàn),以檢驗(yàn)前文中提出的最優(yōu)滾動(dòng)窗寬選擇標(biāo)準(zhǔn)的表現(xiàn)??紤]了兩種變量關(guān)系:線性關(guān)系(線性回歸)和非線性關(guān)系(單指標(biāo)半?yún)?shù)模型),從均方誤差MSE和曼哈頓距離分析本文方法(下文簡(jiǎn)稱最優(yōu)法)在捕捉參數(shù)準(zhǔn)確性和時(shí)變性上的表現(xiàn)。由于現(xiàn)有文獻(xiàn)對(duì)滾動(dòng)窗寬的研究,均從預(yù)測(cè)角度并僅針對(duì)線性模型[13-16]。因此,在線性關(guān)系框架下(鏈接函數(shù)已知),將最優(yōu)法與預(yù)測(cè)均方誤差法的結(jié)果進(jìn)行對(duì)比。在非線性關(guān)系框架下(鏈接函數(shù)未知),僅列出本文方法的估計(jì)結(jié)果。

    預(yù)測(cè)均方誤差法[13-16]通過(guò)最小化近似預(yù)測(cè)均方誤差(MSFE)選擇最優(yōu)滾動(dòng)窗寬。但與以預(yù)測(cè)精度為目的的方法有所區(qū)別的是,本文以估計(jì)時(shí)變參數(shù)為主要目的,故根據(jù)預(yù)測(cè)法[13-16]原理,結(jié)合前文中模型設(shè)定(1)式,具體準(zhǔn)則如下:

    (11)

    3.1 數(shù)據(jù)生成過(guò)程(DGP)

    本文選擇了兩類數(shù)據(jù)生成過(guò)程,線性模型和非線性模型,以說(shuō)明本文的方法對(duì)于線性關(guān)系和非線性關(guān)系均適用。數(shù)據(jù)生成過(guò)程如下:

    yt=at*x1t+bt*x2t+εt

    (12)

    yt=exp(at*x1t+bt*x2t)2+εt

    (13)

    3.2 模擬結(jié)果

    本文選擇在樣本量T=100和T=200時(shí)執(zhí)行5000次蒙特卡羅模擬實(shí)驗(yàn),評(píng)估線性模型和非線性模型的估計(jì)表現(xiàn)。表1和表2給出了模擬結(jié)果,包括各最優(yōu)滾動(dòng)窗寬對(duì)應(yīng)的均方誤差與全樣本均方誤差的相對(duì)比值、各子窗口間系數(shù)的曼哈頓距離。

    均方誤差相對(duì)值及曼哈頓距離值的具體計(jì)算公式如下:

    (14)

    (15)

    表1 線性模型的蒙特卡羅模擬實(shí)驗(yàn)結(jié)果

    表2 非線性模型的蒙特卡羅模擬實(shí)驗(yàn)結(jié)果

    根據(jù)理論分析可知,對(duì)于常參數(shù)模型,全樣本估計(jì)的效果應(yīng)該比滾動(dòng)樣本技術(shù)估計(jì)更優(yōu)。在式(14)中,相對(duì)均方誤差值小于1,說(shuō)明滾動(dòng)樣本技術(shù)的估計(jì)表現(xiàn)優(yōu)于全樣本估計(jì)。根據(jù)表1和表2的結(jié)果,可以得到以下結(jié)論:

    第一,真實(shí)均方誤差(也是不可行均方誤差)在所有數(shù)據(jù)生成過(guò)程中都得到了最小的相對(duì)均方誤差值。在線性模型和非線性模型中,Bootstrap法得到的近似均方誤差與真實(shí)均方誤差比較接近。說(shuō)明Bootstrap法能夠用來(lái)估計(jì)線性模型和非線性模型的均方誤差。

    第二,當(dāng)參數(shù)為常數(shù)時(shí)(DGP1,DGP14),滾動(dòng)樣本技術(shù)相對(duì)全樣本的均方誤差比值較大。在線性模型框架下(DGP1),最優(yōu)法的估計(jì)表現(xiàn)優(yōu)于預(yù)測(cè)法??傮w而言,當(dāng)參數(shù)不存在斷點(diǎn)時(shí),滾動(dòng)樣本技術(shù)的估計(jì)表現(xiàn)較差。因?yàn)榧僭O(shè)樣本中存在斷點(diǎn)是錯(cuò)誤的。但隨著樣本量增加,滾動(dòng)法的表現(xiàn)存在改善。

    第三,當(dāng)參數(shù)出現(xiàn)離散斷點(diǎn)時(shí)(DGP2-7,DGP15-20),滾動(dòng)樣本技術(shù)的表現(xiàn)較常參數(shù)情況下好,且非線性模型較線性模型的表現(xiàn)好,但Bootstrap近似均方誤差與真實(shí)均方誤差仍有較大差異。在線性模型框架下(DGP2-7),預(yù)測(cè)法的估計(jì)表現(xiàn)仍比最優(yōu)法的估計(jì)表現(xiàn)差。

    第四,當(dāng)參數(shù)出現(xiàn)平滑連續(xù)斷點(diǎn)時(shí)(DGP8-11,DGP21-24),滾動(dòng)樣本技術(shù)的表現(xiàn)較好,且在線性模型和非線性模型下差別不大,最優(yōu)法和預(yù)測(cè)法的估計(jì)表現(xiàn)差異較小,Bootstrap均方誤差與真實(shí)均方誤差比較接近,但仍較真實(shí)均方誤差的表現(xiàn)差,因?yàn)锽ootstrap法又增加了額外的誤差。隨著樣本量增加,Bootstrap法的表現(xiàn)傾向于改善。

    第五,當(dāng)參數(shù)出現(xiàn)隨機(jī)游走連續(xù)斷點(diǎn)時(shí)(DGP12-13,DGP25-26),Bootstrap法在線性模型和非線性模型框架下的表現(xiàn)較好,最優(yōu)法和預(yù)測(cè)法的估計(jì)表現(xiàn)較好且差異較小。

    第六,隨著參數(shù)斷點(diǎn)數(shù)目的增加,曼哈頓距離值呈增加趨勢(shì),說(shuō)明該指標(biāo)能夠測(cè)度不同子窗口估計(jì)值之間的差異性。

    以上分析表明,當(dāng)參數(shù)出現(xiàn)連續(xù)斷點(diǎn)時(shí),本文提出的采用最小化Bootstrap近似均方誤差,最大化曼哈頓距離確定滾動(dòng)窗寬,得到的模型估計(jì)量表現(xiàn)較好。即使當(dāng)參數(shù)不是連續(xù)斷點(diǎn)時(shí),Bootstrap近似均方誤差的表現(xiàn)也是合理的,在缺乏更優(yōu)準(zhǔn)則情況時(shí),也可以采用。在線性框架下,最優(yōu)法的估計(jì)表現(xiàn)優(yōu)于預(yù)測(cè)法,也說(shuō)明了不能直接將以預(yù)測(cè)精度為目的的滾動(dòng)窗寬選擇標(biāo)準(zhǔn)用于以捕捉參數(shù)結(jié)構(gòu)突變?yōu)槟康牡哪P椭?。因此,通過(guò)求解最小化近似均方誤差,最大化曼哈頓距離下的雙目標(biāo)優(yōu)化問(wèn)題,從而選擇滾動(dòng)窗寬的做法能夠更好地捕捉參數(shù)的結(jié)構(gòu)突變性,具有理論和實(shí)踐意義。

    3.3 參數(shù)敏感性分析

    同一數(shù)據(jù)生成過(guò)程中,不同的參數(shù)設(shè)置可能會(huì)影響到最優(yōu)滾動(dòng)窗寬的選擇和模型估計(jì)表現(xiàn)的穩(wěn)健性,故對(duì)各數(shù)據(jù)生成過(guò)程中主要參數(shù)進(jìn)行敏感性分析。分別分析線性模型和非線性模型在樣本量T=100,幾類典型數(shù)據(jù)生成過(guò)程:DGP1和DGP14(無(wú)斷點(diǎn)),DGP2-4和DGP15-17(一個(gè)離散斷點(diǎn)),DGP8和DGP21(平滑時(shí)變連續(xù)斷點(diǎn)),DGP10和DGP23(二次平滑時(shí)變連續(xù)斷點(diǎn)),DGP12和DGP25(隨機(jī)游走連續(xù)斷點(diǎn))情形下,bt參數(shù)變化(變小或變大兩種情況)對(duì)模擬結(jié)果(最優(yōu)滾動(dòng)窗寬對(duì)應(yīng)的模型估計(jì)表現(xiàn))的影響,分析結(jié)果見(jiàn)表3。

    表3 主要參數(shù)敏感性分析結(jié)果

    將表3與表1,表2對(duì)比可知,在各數(shù)據(jù)生成過(guò)程中,參數(shù)bt在某區(qū)間內(nèi)變小或變大,對(duì)模擬結(jié)果(最優(yōu)滾動(dòng)窗寬對(duì)應(yīng)的模型估計(jì)表現(xiàn))影響較小,仍能得到與前文中仿真模擬相同的結(jié)論。在所有數(shù)據(jù)生成過(guò)程中,真實(shí)均方誤差比值均是最小的,且Bootstrap近似均方誤差值與真實(shí)值比較接近;隨著斷點(diǎn)數(shù)目的增加和連續(xù)性的增強(qiáng),最優(yōu)法表現(xiàn)更優(yōu)。說(shuō)明,本文的最優(yōu)滾動(dòng)窗寬選擇標(biāo)準(zhǔn)對(duì)同類數(shù)據(jù)生成過(guò)程中參數(shù)變化不敏感,是穩(wěn)健的。

    4 實(shí)證檢驗(yàn)

    宏觀經(jīng)濟(jì)系列通常受到外來(lái)沖擊或社會(huì)經(jīng)濟(jì)制度改變的影響,從而處于長(zhǎng)期不穩(wěn)定狀態(tài),使用靜態(tài)模型剖析金融經(jīng)濟(jì)序列結(jié)構(gòu)顯然不合理。Betz等[4], H?rdle等[7], Paltalidis等[32]的研究已表明金融系統(tǒng)的網(wǎng)絡(luò)結(jié)構(gòu)呈現(xiàn)出了較強(qiáng)的結(jié)構(gòu)突變性。本部分采用滾動(dòng)樣本技術(shù)研究我國(guó)金融系統(tǒng)的網(wǎng)絡(luò)結(jié)構(gòu)隨時(shí)間的變化情況,以檢驗(yàn)第二部分提出的最優(yōu)滾動(dòng)窗寬選擇標(biāo)準(zhǔn)的實(shí)踐價(jià)值。

    4.1 模型簡(jiǎn)介

    運(yùn)用單指標(biāo)分位數(shù)模型[7]構(gòu)建金融網(wǎng)絡(luò)結(jié)構(gòu):

    (16)

    (17)

    (18)

    (19)

    表4 第m窗口下金融系統(tǒng)的網(wǎng)絡(luò)結(jié)構(gòu)

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

    選取2011年前已上市的30家金融機(jī)構(gòu),其中銀行16家,保險(xiǎn)公司3家和證券公司11家,樣本周期為2010年10月16日至2015年9月26日,覆蓋了2012年希臘債務(wù)危機(jī),2013年銀行業(yè)“錢荒”以及2015年股市崩盤事件。周頻股票收盤價(jià)數(shù)據(jù)和財(cái)務(wù)數(shù)據(jù)來(lái)自于國(guó)泰君安數(shù)據(jù)庫(kù),宏觀經(jīng)濟(jì)數(shù)據(jù)來(lái)自于Wind數(shù)據(jù)庫(kù)。低頻控制變量(季、月)采用三次樣條插值技術(shù)轉(zhuǎn)化為周頻數(shù)據(jù),高頻控制變量采用均值法轉(zhuǎn)化為周頻數(shù)據(jù)。各控制變量定義[1,33-35]見(jiàn)表5。下文的實(shí)證過(guò)程在R軟件中實(shí)現(xiàn)。

    4.3 結(jié)果分析

    金融系統(tǒng)的網(wǎng)絡(luò)結(jié)構(gòu)模型的估計(jì)表現(xiàn)及網(wǎng)絡(luò)傳染度由表6和圖2給出。表6中給出了使用不同時(shí)間截點(diǎn)[4,7,36-37]選擇窗寬、本文的最優(yōu)法選擇窗寬的模型估計(jì)表現(xiàn),對(duì)比分析了不同方法的結(jié)果,以說(shuō)明本文提出的滾動(dòng)窗寬選擇標(biāo)準(zhǔn)的優(yōu)良性和可靠性。圖2給出了整個(gè)金融系統(tǒng)的網(wǎng)絡(luò)傳染度隨時(shí)間的變化趨勢(shì),根據(jù)式(19)得到。

    表5 控制變量定義

    表6 金融網(wǎng)絡(luò)結(jié)構(gòu)模型的估計(jì)表現(xiàn)

    比較表6中不同窗寬選擇方法下模型的估計(jì)表現(xiàn),得到的結(jié)論有:第一,當(dāng)采用不同方法選擇窗寬時(shí),我國(guó)金融網(wǎng)絡(luò)結(jié)構(gòu)模型的估計(jì)表現(xiàn)差異較大。與最優(yōu)法相比,以一年、兩年等時(shí)間截點(diǎn)選擇窗寬時(shí),Bootstrap均方誤差較大,說(shuō)明模型的估計(jì)量表現(xiàn)差;曼哈頓距離值較小,說(shuō)明參數(shù)隨時(shí)間的變化情況相對(duì)較小。這可能導(dǎo)致的后果是無(wú)法捕捉到關(guān)鍵結(jié)構(gòu)突變點(diǎn);第二,采用最優(yōu)法選擇滾動(dòng)窗寬時(shí),Bootstrap均方誤差值較小,曼哈頓距離值較高。這說(shuō)明,本文方法能夠兼顧兩相悖目標(biāo):既能捕捉網(wǎng)絡(luò)結(jié)構(gòu)的結(jié)構(gòu)突變性,也能使參數(shù)估計(jì)量的表現(xiàn)較優(yōu),優(yōu)化了金融網(wǎng)絡(luò)的結(jié)構(gòu)突變識(shí)別過(guò)程。

    圖2 整個(gè)金融系統(tǒng)的網(wǎng)絡(luò)傳染度的時(shí)間變化趨勢(shì)圖

    由圖2可以看出,網(wǎng)絡(luò)傳染度隨時(shí)間的變化趨勢(shì)明顯,2011年至2012年呈現(xiàn)出小幅上升后回落的趨勢(shì),這與2012年上半年希臘危機(jī),我國(guó)加快利率、匯率等市場(chǎng)化改革等有關(guān)。在2013年間,利率市場(chǎng)化進(jìn)程加快,各大銀行激烈競(jìng)爭(zhēng)增加了中小銀行的競(jìng)爭(zhēng)成本,隨之出現(xiàn)的“錢荒”問(wèn)題也突顯了銀行系統(tǒng)面臨的短期流動(dòng)性問(wèn)題,使得這期間金融系統(tǒng)的網(wǎng)絡(luò)傳染度水平攀升。隨后由于宏觀政策的調(diào)控,網(wǎng)絡(luò)傳染度有一定回落。而2015年上半年的股市出現(xiàn)較大波動(dòng),產(chǎn)能過(guò)剩問(wèn)題不斷升溫,導(dǎo)致2015年間的金融系統(tǒng)網(wǎng)絡(luò)傳染度不斷升高,經(jīng)政府的不斷注資和出臺(tái)相應(yīng)的調(diào)控政策,金融系統(tǒng)網(wǎng)絡(luò)傳染度逐漸回落但仍然處于一個(gè)相對(duì)高位。可見(jiàn),當(dāng)出現(xiàn)政策沖擊,市場(chǎng)沖擊,風(fēng)險(xiǎn)蔓延等事件時(shí),金融系統(tǒng)的網(wǎng)絡(luò)傳染度會(huì)出現(xiàn)拐點(diǎn),并呈現(xiàn)劇烈波動(dòng)現(xiàn)象,說(shuō)明本文的最優(yōu)窗寬選擇標(biāo)準(zhǔn)能夠捕捉到金融網(wǎng)絡(luò)的結(jié)構(gòu)突變,從實(shí)證上證明了本文方法的有效性。

    5 結(jié)語(yǔ)

    本文針對(duì)時(shí)變參數(shù)模型,提供了一種可行的窗寬選擇方法,具有理論意義和實(shí)踐意義。蒙特卡羅模擬實(shí)驗(yàn)表明:第一,Bootstrap近似均方誤差能夠替代真實(shí)均方誤差;第二,使用本文提出的最優(yōu)窗寬選擇法得到的模型估計(jì)量?jī)?yōu)于全樣本序列得到的模型估計(jì)量,當(dāng)參數(shù)出現(xiàn)連續(xù)斷點(diǎn)時(shí),本文的方法更有效;第三,在線性框架下,最優(yōu)窗寬選擇方法的估計(jì)表現(xiàn)優(yōu)于以預(yù)測(cè)為目的的窗寬選擇方法,說(shuō)明了不能直接將以預(yù)測(cè)精度為目的的窗寬選擇標(biāo)準(zhǔn)用于以捕捉參數(shù)結(jié)構(gòu)突變?yōu)槟康牡哪P椭小5谒?,本文的方法適用于線性關(guān)系和非線性關(guān)系下的時(shí)變參數(shù)模型,且能拓展至常見(jiàn)的線性回歸模型和廣義線性回歸模型中。實(shí)證檢驗(yàn)表明,本文的方法能夠優(yōu)化金融網(wǎng)絡(luò)的結(jié)構(gòu)突變識(shí)別過(guò)程,彌補(bǔ)了當(dāng)前在該領(lǐng)域中缺乏統(tǒng)計(jì)標(biāo)準(zhǔn)選擇滾動(dòng)窗寬的空白。

    猜你喜歡
    估計(jì)量均方斷點(diǎn)
    一類隨機(jī)積分微分方程的均方漸近概周期解
    Beidou, le système de navigation par satellite compatible et interopérable
    一類無(wú)限可能問(wèn)題的解法
    主導(dǎo)電回路發(fā)生斷點(diǎn)故障判斷方法探討
    淺談估計(jì)量的優(yōu)良性標(biāo)準(zhǔn)
    基于抗差最小均方估計(jì)的輸電線路參數(shù)辨識(shí)
    基于配網(wǎng)先驗(yàn)信息的諧波狀態(tài)估計(jì)量測(cè)點(diǎn)最優(yōu)配置
    基于隨機(jī)牽制控制的復(fù)雜網(wǎng)絡(luò)均方簇同步
    負(fù)極值指標(biāo)估計(jì)量的漸近性質(zhì)
    使用變異系數(shù)和Kurtosis系數(shù)的雙輔助變量的比估計(jì)
    偷拍熟女少妇极品色| 日本免费a在线| 最后的刺客免费高清国语| 综合色av麻豆| 亚洲国产精品sss在线观看| 中文字幕人妻熟人妻熟丝袜美| 亚洲天堂国产精品一区在线| 国产成人精品婷婷| 又爽又黄无遮挡网站| 卡戴珊不雅视频在线播放| 国产亚洲5aaaaa淫片| 少妇的逼水好多| 国产中年淑女户外野战色| 国精品久久久久久国模美| 亚洲av成人av| 婷婷色综合www| 亚洲精品一区蜜桃| av.在线天堂| 九九久久精品国产亚洲av麻豆| 舔av片在线| 亚洲欧美一区二区三区国产| 久久亚洲国产成人精品v| 91久久精品电影网| 最近视频中文字幕2019在线8| 成人av在线播放网站| 真实男女啪啪啪动态图| 国产精品久久久久久久电影| 免费电影在线观看免费观看| 久久久色成人| 80岁老熟妇乱子伦牲交| 亚洲精品aⅴ在线观看| 亚洲自偷自拍三级| 久久久色成人| 天堂av国产一区二区熟女人妻| 亚洲欧美精品自产自拍| 久久人人爽人人爽人人片va| 国产人妻一区二区三区在| 精品久久久噜噜| 狂野欧美白嫩少妇大欣赏| 国产成人福利小说| 亚洲欧美成人综合另类久久久| 精品国内亚洲2022精品成人| 一本一本综合久久| 精品久久久久久电影网| 国产伦在线观看视频一区| 91精品国产九色| 黑人高潮一二区| 日韩,欧美,国产一区二区三区| 国产一级毛片在线| 亚洲精品日本国产第一区| 亚洲国产精品成人久久小说| 观看免费一级毛片| 亚洲成人av在线免费| 精品久久久久久久久亚洲| 日韩成人av中文字幕在线观看| 久久99热6这里只有精品| 欧美高清成人免费视频www| 91午夜精品亚洲一区二区三区| 国产黄a三级三级三级人| 在线免费十八禁| 国产淫语在线视频| 2021天堂中文幕一二区在线观| 听说在线观看完整版免费高清| 国产精品一区www在线观看| 亚洲丝袜综合中文字幕| 日本与韩国留学比较| 免费看美女性在线毛片视频| 中文字幕亚洲精品专区| 色综合色国产| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 久久这里只有精品中国| 亚洲va在线va天堂va国产| 99热这里只有精品一区| 亚洲国产精品国产精品| 国产精品蜜桃在线观看| 久久综合国产亚洲精品| 亚洲第一区二区三区不卡| 国产老妇伦熟女老妇高清| 日韩三级伦理在线观看| 我的老师免费观看完整版| 国产久久久一区二区三区| 一区二区三区免费毛片| 国产成人91sexporn| 中文欧美无线码| 亚洲av电影在线观看一区二区三区 | 爱豆传媒免费全集在线观看| 亚洲综合精品二区| 中文在线观看免费www的网站| 人人妻人人澡人人爽人人夜夜 | 精品一区二区三区视频在线| 99热这里只有是精品50| 永久网站在线| 亚洲国产欧美在线一区| 男人和女人高潮做爰伦理| 国产成人a区在线观看| 成年女人看的毛片在线观看| 亚洲欧美成人综合另类久久久| 日本与韩国留学比较| 婷婷六月久久综合丁香| 亚洲av免费在线观看| 成人特级av手机在线观看| 欧美成人一区二区免费高清观看| 精品欧美国产一区二区三| 18禁动态无遮挡网站| 少妇的逼水好多| 成年版毛片免费区| 在线观看av片永久免费下载| 久久精品久久久久久噜噜老黄| 久久久国产一区二区| 国精品久久久久久国模美| 男人爽女人下面视频在线观看| 三级国产精品片| 一个人看视频在线观看www免费| 精品久久久久久久久亚洲| 大片免费播放器 马上看| 在线观看人妻少妇| 男人舔奶头视频| 美女xxoo啪啪120秒动态图| 91精品国产九色| 国产精品一区二区在线观看99 | 欧美日韩亚洲高清精品| 日本免费在线观看一区| 中国国产av一级| 国产精品伦人一区二区| 大话2 男鬼变身卡| 99九九线精品视频在线观看视频| 亚洲天堂国产精品一区在线| 22中文网久久字幕| 91久久精品国产一区二区成人| 99久久人妻综合| 欧美三级亚洲精品| 久久综合国产亚洲精品| 亚洲天堂国产精品一区在线| 亚洲欧美精品专区久久| 最近手机中文字幕大全| 中国美白少妇内射xxxbb| 久久精品国产自在天天线| 一本久久精品| 欧美激情国产日韩精品一区| 精品午夜福利在线看| 久久久久精品性色| 亚洲av在线观看美女高潮| 亚洲精品一区蜜桃| 色播亚洲综合网| 国产av国产精品国产| 国产成人a∨麻豆精品| 美女国产视频在线观看| 国产老妇女一区| 不卡视频在线观看欧美| 免费黄网站久久成人精品| 免费不卡的大黄色大毛片视频在线观看 | eeuss影院久久| 男的添女的下面高潮视频| 国产乱人视频| 亚洲熟妇中文字幕五十中出| 中文字幕制服av| 国产有黄有色有爽视频| 天美传媒精品一区二区| 成人一区二区视频在线观看| 美女大奶头视频| 97热精品久久久久久| 日韩一区二区视频免费看| 亚洲精品日韩av片在线观看| 精品久久久精品久久久| 欧美+日韩+精品| 亚洲国产精品专区欧美| 国产中年淑女户外野战色| 99久久九九国产精品国产免费| 亚洲自偷自拍三级| 99热这里只有是精品在线观看| 免费无遮挡裸体视频| 亚洲av福利一区| 干丝袜人妻中文字幕| 成人亚洲欧美一区二区av| 国产亚洲午夜精品一区二区久久 | 亚洲精品一二三| 国产午夜精品久久久久久一区二区三区| 国产精品久久视频播放| 久久亚洲国产成人精品v| 亚洲国产av新网站| 天堂俺去俺来也www色官网 | 久久韩国三级中文字幕| 男人狂女人下面高潮的视频| 亚洲一区高清亚洲精品| 国产69精品久久久久777片| 狂野欧美激情性xxxx在线观看| 久久久久久伊人网av| 我的老师免费观看完整版| 亚洲无线观看免费| 色尼玛亚洲综合影院| 国模一区二区三区四区视频| 听说在线观看完整版免费高清| 亚洲自偷自拍三级| 日韩av免费高清视频| 久久国内精品自在自线图片| 免费看光身美女| 国产男人的电影天堂91| 狠狠精品人妻久久久久久综合| 国内精品美女久久久久久| 深爱激情五月婷婷| 欧美潮喷喷水| 高清日韩中文字幕在线| 91精品国产九色| 三级国产精品欧美在线观看| 好男人视频免费观看在线| 国产午夜精品一二区理论片| 亚州av有码| 美女高潮的动态| 亚洲第一区二区三区不卡| 亚洲av一区综合| 国产成人a区在线观看| 亚洲欧美日韩无卡精品| 麻豆成人午夜福利视频| 在线a可以看的网站| av一本久久久久| 久久久久久久久久久丰满| 男人舔奶头视频| av国产久精品久网站免费入址| 国产在线男女| 日韩欧美三级三区| 蜜桃久久精品国产亚洲av| 22中文网久久字幕| 大香蕉97超碰在线| 欧美最新免费一区二区三区| 久久99热这里只频精品6学生| 天堂网av新在线| 日韩中字成人| 日韩一区二区三区影片| 欧美xxⅹ黑人| videossex国产| 亚洲精品一二三| 男人舔女人下体高潮全视频| 国产一区亚洲一区在线观看| 国语对白做爰xxxⅹ性视频网站| 听说在线观看完整版免费高清| 最新中文字幕久久久久| 国内精品一区二区在线观看| 欧美日本视频| 国产探花极品一区二区| 亚洲精华国产精华液的使用体验| 精品一区二区三卡| 成人美女网站在线观看视频| 亚洲无线观看免费| 纵有疾风起免费观看全集完整版 | 晚上一个人看的免费电影| 中文天堂在线官网| 校园人妻丝袜中文字幕| 日韩制服骚丝袜av| 日韩人妻高清精品专区| 午夜激情久久久久久久| 国产精品久久久久久av不卡| 高清毛片免费看| 黄色配什么色好看| 亚洲欧美精品专区久久| 听说在线观看完整版免费高清| 日日啪夜夜爽| 精品国产三级普通话版| 久久草成人影院| 久久久久久九九精品二区国产| 午夜福利视频精品| 成人高潮视频无遮挡免费网站| 一级毛片电影观看| 亚洲人成网站高清观看| 七月丁香在线播放| 伊人久久国产一区二区| 亚洲欧洲国产日韩| 晚上一个人看的免费电影| 亚洲精品日韩在线中文字幕| 亚洲国产精品成人久久小说| 精品一区二区三区人妻视频| 国产黄片视频在线免费观看| 欧美3d第一页| 国产淫语在线视频| 1000部很黄的大片| 亚洲精品乱久久久久久| 欧美成人a在线观看| 真实男女啪啪啪动态图| 插逼视频在线观看| 成人国产麻豆网| 午夜精品在线福利| 国产精品.久久久| 国产亚洲最大av| 国产成人a区在线观看| 国产视频首页在线观看| 精品一区二区三卡| 国内揄拍国产精品人妻在线| 18禁动态无遮挡网站| 日产精品乱码卡一卡2卡三| 精品欧美国产一区二区三| 国产高清不卡午夜福利| 好男人视频免费观看在线| 国产探花极品一区二区| 乱人视频在线观看| 小蜜桃在线观看免费完整版高清| 精品久久久久久久久av| 日韩视频在线欧美| 亚洲av成人av| 国产探花在线观看一区二区| 热99在线观看视频| 亚洲国产成人一精品久久久| 日本午夜av视频| 99热全是精品| 神马国产精品三级电影在线观看| 三级毛片av免费| av在线蜜桃| 一本久久精品| 全区人妻精品视频| 亚洲精品国产av成人精品| 九草在线视频观看| 一二三四中文在线观看免费高清| 久久鲁丝午夜福利片| 中国国产av一级| 亚洲av中文字字幕乱码综合| 久久久精品94久久精品| 国产在线男女| 精品久久久久久久久av| 国产探花在线观看一区二区| 一边亲一边摸免费视频| 天堂网av新在线| 最近最新中文字幕免费大全7| 蜜桃久久精品国产亚洲av| 亚洲国产精品sss在线观看| 秋霞在线观看毛片| 午夜日本视频在线| 三级国产精品片| 午夜日本视频在线| 国产乱来视频区| 男女边摸边吃奶| 成人性生交大片免费视频hd| 男的添女的下面高潮视频| 亚洲最大成人手机在线| 精品一区二区三区人妻视频| 国产老妇伦熟女老妇高清| 欧美三级亚洲精品| 精品欧美国产一区二区三| 夜夜看夜夜爽夜夜摸| 大片免费播放器 马上看| 99久久九九国产精品国产免费| 一夜夜www| 色5月婷婷丁香| or卡值多少钱| 国产伦理片在线播放av一区| 一级毛片久久久久久久久女| 伊人久久精品亚洲午夜| 亚洲精品久久午夜乱码| 国产伦一二天堂av在线观看| 日韩成人av中文字幕在线观看| 色5月婷婷丁香| 听说在线观看完整版免费高清| 午夜免费男女啪啪视频观看| 狂野欧美激情性xxxx在线观看| 精品久久久久久久人妻蜜臀av| 精品久久久久久久久av| 少妇的逼水好多| 一级毛片电影观看| 亚洲国产日韩欧美精品在线观看| 一级二级三级毛片免费看| 午夜免费男女啪啪视频观看| 一级毛片电影观看| 色5月婷婷丁香| 99九九线精品视频在线观看视频| 精品少妇黑人巨大在线播放| 69av精品久久久久久| 岛国毛片在线播放| 欧美xxxx黑人xx丫x性爽| 欧美最新免费一区二区三区| 春色校园在线视频观看| 国产一级毛片七仙女欲春2| 国内精品宾馆在线| 国产永久视频网站| 国产美女午夜福利| 国产永久视频网站| 建设人人有责人人尽责人人享有的 | 国产高清国产精品国产三级 | 狠狠精品人妻久久久久久综合| 中文字幕av成人在线电影| 一个人免费在线观看电影| 韩国av在线不卡| 国产午夜福利久久久久久| 亚洲av成人av| 直男gayav资源| 亚洲精品日韩在线中文字幕| 大片免费播放器 马上看| 赤兔流量卡办理| 亚洲成人中文字幕在线播放| 久久韩国三级中文字幕| 日韩电影二区| 成人国产麻豆网| 少妇人妻精品综合一区二区| 亚洲经典国产精华液单| 搡老乐熟女国产| 亚洲人成网站高清观看| 精华霜和精华液先用哪个| 国产美女午夜福利| 国产老妇女一区| 一夜夜www| 国产伦精品一区二区三区四那| 久久久久久久久大av| 国产一级毛片在线| 国产乱人视频| 免费观看的影片在线观看| 亚洲欧美日韩东京热| 欧美zozozo另类| 久久久久久久久久久免费av| 最后的刺客免费高清国语| 成年女人看的毛片在线观看| 亚洲欧美一区二区三区国产| 国产精品无大码| 久久久久久久大尺度免费视频| 热99在线观看视频| 国产成人a区在线观看| 一边亲一边摸免费视频| 尾随美女入室| 午夜福利视频精品| 天堂中文最新版在线下载 | 九九爱精品视频在线观看| 久久久久久久亚洲中文字幕| 波多野结衣巨乳人妻| 欧美精品一区二区大全| 2021天堂中文幕一二区在线观| 久久久久久国产a免费观看| 99久国产av精品国产电影| 亚洲高清免费不卡视频| 成年免费大片在线观看| 中文在线观看免费www的网站| 少妇熟女aⅴ在线视频| 日韩中字成人| 91久久精品国产一区二区成人| 色视频www国产| 国产伦一二天堂av在线观看| 汤姆久久久久久久影院中文字幕 | 精品久久久精品久久久| 国产成人精品一,二区| 欧美bdsm另类| 麻豆av噜噜一区二区三区| 两个人的视频大全免费| 亚洲精品国产av成人精品| 欧美不卡视频在线免费观看| 一级二级三级毛片免费看| 国产伦精品一区二区三区四那| 日韩欧美三级三区| 直男gayav资源| 五月伊人婷婷丁香| 久久99热6这里只有精品| 久久精品人妻少妇| 最近手机中文字幕大全| 白带黄色成豆腐渣| 在线观看美女被高潮喷水网站| 波多野结衣巨乳人妻| 国产成人精品婷婷| 又黄又爽又刺激的免费视频.| 亚洲美女视频黄频| 午夜精品国产一区二区电影 | 日韩强制内射视频| 十八禁网站网址无遮挡 | 国内少妇人妻偷人精品xxx网站| 麻豆精品久久久久久蜜桃| 男人舔女人下体高潮全视频| 国产高清三级在线| or卡值多少钱| 看十八女毛片水多多多| 亚洲国产精品成人久久小说| 2021天堂中文幕一二区在线观| 亚洲人成网站在线播| 国产精品久久久久久精品电影| 狂野欧美白嫩少妇大欣赏| 国产久久久一区二区三区| 日本-黄色视频高清免费观看| 夫妻午夜视频| 中文字幕制服av| 麻豆久久精品国产亚洲av| 综合色丁香网| 欧美97在线视频| 婷婷色综合大香蕉| 国产黄色视频一区二区在线观看| 免费无遮挡裸体视频| 美女xxoo啪啪120秒动态图| 欧美变态另类bdsm刘玥| 天美传媒精品一区二区| 欧美日韩一区二区视频在线观看视频在线 | av一本久久久久| 九草在线视频观看| 国产老妇伦熟女老妇高清| av播播在线观看一区| 亚洲精品一区蜜桃| 成人特级av手机在线观看| 成人高潮视频无遮挡免费网站| 日本一二三区视频观看| 免费看a级黄色片| 亚洲av在线观看美女高潮| 国产黄a三级三级三级人| 久久精品国产亚洲av涩爱| 一区二区三区免费毛片| 日韩av在线免费看完整版不卡| 成人亚洲精品av一区二区| 亚洲自偷自拍三级| 熟妇人妻久久中文字幕3abv| 亚洲精品国产av成人精品| 91午夜精品亚洲一区二区三区| 中文字幕制服av| 欧美精品国产亚洲| 国产探花极品一区二区| 国产精品熟女久久久久浪| 国产av在哪里看| 久久久久久久久久久丰满| 国产亚洲午夜精品一区二区久久 | 亚洲欧美成人综合另类久久久| 大又大粗又爽又黄少妇毛片口| 一个人免费在线观看电影| 欧美日韩一区二区视频在线观看视频在线 | 欧美日韩国产mv在线观看视频 | 亚洲欧美精品自产自拍| 超碰97精品在线观看| 免费不卡的大黄色大毛片视频在线观看 | 国产成人freesex在线| 噜噜噜噜噜久久久久久91| 亚洲图色成人| av在线蜜桃| 成人无遮挡网站| 在线播放无遮挡| 亚洲国产欧美在线一区| 精品99又大又爽又粗少妇毛片| 成人毛片a级毛片在线播放| 久久韩国三级中文字幕| 国产在线男女| 精品午夜福利在线看| 午夜精品国产一区二区电影 | 精品人妻一区二区三区麻豆| 免费高清在线观看视频在线观看| 五月天丁香电影| 日韩欧美 国产精品| 简卡轻食公司| 久久精品夜夜夜夜夜久久蜜豆| 免费人成在线观看视频色| 成人午夜高清在线视频| 国产爱豆传媒在线观看| 欧美极品一区二区三区四区| 国产亚洲91精品色在线| 国产毛片a区久久久久| 人妻少妇偷人精品九色| 亚洲欧美成人综合另类久久久| 亚洲av中文字字幕乱码综合| 久久午夜福利片| 日产精品乱码卡一卡2卡三| 中文在线观看免费www的网站| av在线天堂中文字幕| 大陆偷拍与自拍| 色视频www国产| 国产一区亚洲一区在线观看| 久久鲁丝午夜福利片| 成人av在线播放网站| 又大又黄又爽视频免费| 国产在线男女| 一区二区三区乱码不卡18| 亚洲精品,欧美精品| 日本一本二区三区精品| 久久99热这里只频精品6学生| 永久免费av网站大全| 蜜桃久久精品国产亚洲av| 亚洲精品国产成人久久av| 欧美日韩综合久久久久久| 天堂俺去俺来也www色官网 | 亚洲图色成人| 99九九线精品视频在线观看视频| 日本-黄色视频高清免费观看| 99久国产av精品| 男女边吃奶边做爰视频| 亚洲色图av天堂| 在线免费观看不下载黄p国产| 高清日韩中文字幕在线| 国产探花极品一区二区| 午夜爱爱视频在线播放| 亚洲国产精品专区欧美| 国产精品99久久久久久久久| 亚洲国产色片| 亚洲欧美成人综合另类久久久| kizo精华| 午夜精品在线福利| 国内揄拍国产精品人妻在线| 欧美高清成人免费视频www| 美女cb高潮喷水在线观看| 五月伊人婷婷丁香| 亚洲人成网站在线观看播放| 性色avwww在线观看| 91av网一区二区| 国产激情偷乱视频一区二区| 亚洲国产最新在线播放| 亚洲精品日韩av片在线观看| 亚洲av中文字字幕乱码综合| 插阴视频在线观看视频| 亚洲成人中文字幕在线播放| 赤兔流量卡办理| 免费观看的影片在线观看| 国产不卡一卡二| 亚洲精华国产精华液的使用体验| 国产亚洲精品久久久com| 久久精品久久久久久噜噜老黄| 免费av观看视频| 男人和女人高潮做爰伦理| 国产精品av视频在线免费观看| av在线亚洲专区| 简卡轻食公司| 国产精品一及| 岛国毛片在线播放| 国产精品国产三级国产专区5o| 777米奇影视久久| 夫妻性生交免费视频一级片| 成人二区视频| 男女国产视频网站| 欧美三级亚洲精品| 99热这里只有精品一区| 欧美一区二区亚洲| 日韩欧美 国产精品| 中文在线观看免费www的网站| 国产高清有码在线观看视频|