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

    基于投資混頻數(shù)據(jù)的中國經(jīng)濟(jì)波動預(yù)測

    2024-05-26 01:21:54岳金秀
    統(tǒng)計與決策 2024年9期
    關(guān)鍵詞:模型

    岳金秀,陳 堯

    (湖北民族大學(xué)數(shù)學(xué)與統(tǒng)計學(xué)院,湖北 恩施 445000)

    0 引言

    一直以來,經(jīng)濟(jì)波動預(yù)測對國家制定經(jīng)濟(jì)政策和企業(yè)、部門采取相應(yīng)措施至關(guān)重要。而經(jīng)濟(jì)波動與各宏觀經(jīng)濟(jì)指標(biāo)密不可分,其中,固定資產(chǎn)投資與經(jīng)濟(jì)增長之間具有顯著的雙向影響關(guān)系[1]。在現(xiàn)代經(jīng)濟(jì)周期理論中,投資被認(rèn)為是導(dǎo)致經(jīng)濟(jì)波動的主要原因。消費(fèi)需求基本保持比較平滑的變化,隨著國際競爭日漸激烈,以及在擴(kuò)大內(nèi)需的大背景下,凈出口對經(jīng)濟(jì)的驅(qū)動作用也在逐漸降低。因此,投資需求是對經(jīng)濟(jì)波動影響較大的成分。投資的變動將直接或間接影響經(jīng)濟(jì)波動:一方面,投資可以作為資本,直接使GDP增加;另一方面,投資還具有乘數(shù)效應(yīng),投資的獲利會刺激更多的投資,加速消費(fèi)的刺激作用,促進(jìn)居民消費(fèi)。

    目前僅根據(jù)投資預(yù)測經(jīng)濟(jì)波動的研究較少,只有少量文獻(xiàn)涉及相關(guān)內(nèi)容。如戴卓爾等(2023)[2]引入外商直接投資(FDI),基于“三駕馬車”對中國GDP 增長率進(jìn)行預(yù)測,發(fā)現(xiàn)引入FDI后不僅預(yù)測精度更高,而且實(shí)時預(yù)測也比未引入FDI的預(yù)測更精準(zhǔn)。潘雅茹和羅良文(2020)[3]實(shí)證檢驗(yàn)了固定資產(chǎn)投資對經(jīng)濟(jì)高質(zhì)量發(fā)展的推動作用。現(xiàn)有文獻(xiàn)已經(jīng)證實(shí)了投資與經(jīng)濟(jì)波動之間有顯著的相關(guān)性?;诖?,本文從投資角度構(gòu)建中國的宏觀經(jīng)濟(jì)波動預(yù)測指標(biāo)并進(jìn)行實(shí)證分析。隨著網(wǎng)絡(luò)技術(shù)的迅速發(fā)展,使用網(wǎng)絡(luò)數(shù)據(jù)進(jìn)行預(yù)測模型修正成為研究熱點(diǎn)。陳衛(wèi)華和徐國祥(2018)[4]利用股票論壇發(fā)帖數(shù)增長率預(yù)測股票指數(shù)波動,發(fā)現(xiàn)股票論壇數(shù)據(jù)對預(yù)測精度的提高有所貢獻(xiàn)。徐映梅和高一銘(2017)[5]根據(jù)百度取詞構(gòu)建了高頻和低頻輿情指數(shù),提高了對CPI 的預(yù)測精度。Hochreiter 和Schmidhuber(1997)[6]從經(jīng)濟(jì)學(xué)角度給出了個體關(guān)注的定義,他認(rèn)為投資者對資源的認(rèn)知能力是有限的,個體在單位時間內(nèi)只能進(jìn)行有限的信息關(guān)注和處理。通常投資者更容易關(guān)注表現(xiàn)突出和自己感興趣的信息。投資者在單位時間內(nèi)的有限關(guān)注行為,體現(xiàn)了對當(dāng)前項(xiàng)目的關(guān)注。當(dāng)投資者通過網(wǎng)絡(luò)搜索行為產(chǎn)生了正向、積極的心理預(yù)期時,就會轉(zhuǎn)化為投資決策。在整個過程中,網(wǎng)絡(luò)搜索數(shù)據(jù)起到了建議和修正的作用。隨著深度學(xué)習(xí)技術(shù)的發(fā)展,長短期記憶(Long-Short Time Memory,LSTM)神經(jīng)網(wǎng)絡(luò)模型因可以更好地學(xué)習(xí)時間序列之間的長期依賴關(guān)系,開始被廣泛應(yīng)用于預(yù)測模型中。徐映梅和陳堯(2021)[7]將傳統(tǒng)的ARIMA模型與LSTM 模型在我國經(jīng)濟(jì)波動預(yù)測上比較,發(fā)現(xiàn)LSTM 模型在預(yù)測復(fù)雜的非線性時間序列數(shù)據(jù)時,能夠取得更好的效果。梁龍躍和陳玉霞(2023)[8]結(jié)合LSTM模型和小波分析對季度GDP進(jìn)行預(yù)測,結(jié)果顯示,加入小波分析的LSTM 預(yù)測模型具有較好的泛化能力,預(yù)測精度更高。然而傳統(tǒng)的LSTM神經(jīng)網(wǎng)絡(luò)模型在處理混頻數(shù)據(jù)時,模型的學(xué)習(xí)性能也會受到影響。對此,Neil 等(2016)[9]提出了一種改進(jìn)的長短期記憶(Phased-LSTM)網(wǎng)絡(luò)模型,適用于處理異步混頻數(shù)據(jù)的建模。

    綜上所述,宏觀經(jīng)濟(jì)波動預(yù)測多采用傳統(tǒng)計量模型,使用混頻數(shù)據(jù)并利用神經(jīng)網(wǎng)絡(luò)模型進(jìn)行經(jīng)濟(jì)波動預(yù)測的研究還比較少。本文旨在更好地利用高頻數(shù)據(jù)的原始信息,同時實(shí)現(xiàn)變量之間非線性關(guān)系的探討,構(gòu)建多源混頻LSTM模型(Multi Source Mixing LSTM,MM-LSTM),采用政府統(tǒng)計月度數(shù)據(jù)、網(wǎng)絡(luò)搜索日度數(shù)據(jù)以及GDP 季度數(shù)據(jù)多指標(biāo)聯(lián)合,從投資視角研究和預(yù)測宏觀經(jīng)濟(jì)波動。為了驗(yàn)證本文提出的混頻數(shù)據(jù)預(yù)測宏觀經(jīng)濟(jì)波動的有效性,以及構(gòu)建的MM-LSTM 模型的預(yù)測效果,將其與傳統(tǒng)的AR模型、MIDAS模型進(jìn)行對比研究,以評估不同數(shù)據(jù)以及不同方法的研究結(jié)果。

    1 研究設(shè)計

    1.1 統(tǒng)計數(shù)據(jù)收集與處理

    本文基于《中國統(tǒng)計年鑒》中公布的相關(guān)數(shù)據(jù),同時參考與投資相關(guān)的文獻(xiàn),并且兼顧數(shù)據(jù)的可得性,初步挑選出10個候選指標(biāo)。包括工業(yè)增加值累計增長率[10]、發(fā)電量累計增長率[11]、房地產(chǎn)投資累計增長率[12]、制造業(yè)采購經(jīng)理指數(shù)(PMI 指數(shù))[13]、固定資產(chǎn)投資額累計增長率[14]等。收集并整理得到上述10 個指標(biāo)2011 年1 月至2022 年12 月的月度統(tǒng)計數(shù)據(jù),并收集整理2011—2022 年國內(nèi)生產(chǎn)總值指數(shù)累計值(上年同期=100)季度數(shù)據(jù)。對于政府統(tǒng)計數(shù)據(jù),進(jìn)行如下處理:

    (1)GDP增長率的構(gòu)建。本文選用反映一定時期內(nèi)國內(nèi)生產(chǎn)總值變動趨勢的國內(nèi)生產(chǎn)總值指數(shù)累計值(上年同期=100)作為預(yù)測變量。為了消除各指標(biāo)間的數(shù)量差異,把國內(nèi)生產(chǎn)總值指數(shù)累計值(上年同期=100)數(shù)據(jù)轉(zhuǎn)換為增長率,即GDP 增長率=國內(nèi)生產(chǎn)總值指數(shù)累計值(上年同期=100)-100,記為GDP_IG。

    (2)缺失數(shù)據(jù)處理。根據(jù)官方公布的數(shù)據(jù),工業(yè)增加值累計增長率、固定資產(chǎn)投資額累計增長率、水泥產(chǎn)量累計增長率等指標(biāo)缺少各年份1月的數(shù)據(jù),為了數(shù)據(jù)的連續(xù)性,取2月和3月的均值補(bǔ)全1月的數(shù)據(jù)。

    (3)投資統(tǒng)計指標(biāo)篩選。計算各指標(biāo)數(shù)據(jù)與GDP_IG的皮爾遜相關(guān)系數(shù),篩選出與經(jīng)濟(jì)波動相關(guān)性較高的8個指標(biāo),結(jié)果如表1所示。

    (4)投資統(tǒng)計指數(shù)的構(gòu)建。前文統(tǒng)計指標(biāo)篩選已經(jīng)確定了各指標(biāo)與因變量的強(qiáng)相關(guān)關(guān)系。但皮爾遜相關(guān)系數(shù)只能衡量兩個變量間的線性關(guān)系,而隨機(jī)森林算法可以有效處理非線性數(shù)據(jù),因此采用隨機(jī)森林算法計算各指標(biāo)的權(quán)重。利用Python中Sklearn機(jī)器學(xué)習(xí)庫建立隨機(jī)森林回歸模型,得到各特征的重要性得分。模型輸出特征的重要性得分如表2 所示。把每個特征的重要性得分作為每個特征的權(quán)重系數(shù),合成一個月度時間序列,記為投資統(tǒng)計指數(shù)(Investment Statistics Index,ISI)。

    表2 各特征重要性得分

    1.2 網(wǎng)絡(luò)搜索數(shù)據(jù)指標(biāo)收集及處理

    本文選用百度指數(shù)[15]作為網(wǎng)絡(luò)搜索指數(shù)的代表。百度指數(shù)的計算是以網(wǎng)民的搜索量為基礎(chǔ),以關(guān)鍵詞為統(tǒng)計對象。關(guān)鍵詞的選取是合成投資網(wǎng)絡(luò)搜索指數(shù)的重要內(nèi)容。本文選取“投資”作為關(guān)鍵詞,利用文本檢索和篩選,在中國知網(wǎng)檢索CSSCI期刊,得到5389條記錄。下載這些文獻(xiàn),然后利用中文分詞技術(shù)對標(biāo)題進(jìn)行分詞,篩選出“投資”“創(chuàng)新”“經(jīng)濟(jì)”等100個關(guān)鍵詞作為候選關(guān)鍵詞。同時利用百度指數(shù)需求圖譜相關(guān)詞推薦功能,搜集并整理得到2011 年1 月1 日至2022 年12 月31 日與“投資”“房產(chǎn)”“固定資產(chǎn)”等相關(guān)的28 個關(guān)鍵詞,將其加入候選關(guān)鍵詞中。另外,根據(jù)經(jīng)驗(yàn)取詞,選取“發(fā)電量”“金融”“GDP”等8 個關(guān)鍵詞加入候選關(guān)鍵詞中。選取的關(guān)鍵詞還需被百度指數(shù)收錄,因此剔除沒有被收錄的詞,最終整理得到108 個關(guān)鍵詞的實(shí)時搜索量數(shù)據(jù)集。對于投資網(wǎng)絡(luò)搜索數(shù)據(jù),數(shù)據(jù)處理過程如下:

    (1)數(shù)據(jù)預(yù)處理。百度指數(shù)數(shù)據(jù)是以日為頻率的時間序列,但每月天數(shù)并不相等,將日度數(shù)據(jù)換算成月度、季度等時間序列會產(chǎn)生誤差和波動。因此,本文參照徐映梅和高一銘(2017)[5]的方法,將各月份轉(zhuǎn)化為相同的天數(shù)。即刪除2月29日及5月、7月、8月、10月、12月中每月31日的數(shù)據(jù),把1月最后一天的數(shù)據(jù)作為2月第一天的數(shù)據(jù),把3月1日的數(shù)據(jù)作為2月最后一天的數(shù)據(jù)。處理后,每年和每月的天數(shù)相等。即每月30天,每年均為360天。

    (2)消除網(wǎng)絡(luò)發(fā)展趨勢。網(wǎng)絡(luò)搜索數(shù)據(jù)不但會受到網(wǎng)民關(guān)注的影響,而且會受到互聯(lián)網(wǎng)發(fā)展的影響,例如2013年智能手機(jī)的逐漸普及,致使手機(jī)客戶端搜索數(shù)量增加。因此,在數(shù)據(jù)分析之前還需消除網(wǎng)絡(luò)發(fā)展對搜索量波動的影響。借鑒徐映梅和高一銘(2017)[5]的方法,采用與投資關(guān)鍵詞相關(guān)性較低的中性詞來消除網(wǎng)絡(luò)社會發(fā)展對搜索量的影響。

    其中,anum為關(guān)鍵詞搜索量數(shù)據(jù),bnum為中性詞“百度”的搜索量數(shù)據(jù),Anew為處理后的搜索量數(shù)據(jù),記為搜索熱度數(shù)據(jù)。為了防止取對數(shù)時出現(xiàn)異常,對所有搜索量值進(jìn)行加1的平移處理。

    (3)篩選核心關(guān)鍵詞。獲得各關(guān)鍵詞的搜索熱度數(shù)據(jù)后,利用時差相關(guān)分析,求出關(guān)鍵詞搜索熱度序列和房地產(chǎn)投資累計增長率的時差和相關(guān)系數(shù),進(jìn)一步篩選重要的關(guān)鍵詞①根據(jù)各關(guān)鍵詞搜索熱度序列與選取的8個政府統(tǒng)計指標(biāo)選擇最優(yōu)時差相關(guān)分析結(jié)果確定。。利用Python軟件編程計算,最終選取23個相關(guān)性大于0.5的關(guān)鍵詞,其中,相同時差的關(guān)鍵詞有18個,具有領(lǐng)先和滯后階數(shù)的關(guān)鍵詞有5個,將他們調(diào)整為相同時差。這23個關(guān)鍵詞將用于投資網(wǎng)絡(luò)搜索指數(shù)的合成。選取關(guān)鍵詞的時差階數(shù)和最大相關(guān)系數(shù)如表3所示。

    表3 投資網(wǎng)絡(luò)搜索指數(shù)關(guān)鍵詞

    (4)投資網(wǎng)絡(luò)搜索指數(shù)的構(gòu)建。由于選取的關(guān)鍵詞數(shù)量較多且部分關(guān)鍵詞序列存在較強(qiáng)的相關(guān)性,因此本文利用因子分析法,使用少數(shù)的幾個變量來解釋眾多變量,達(dá)到降維的目的,同時又保持較高的解釋力。在對關(guān)鍵詞序列進(jìn)行因子分析之前,需要進(jìn)行KMO 值和Bartlett 球形度檢驗(yàn),檢驗(yàn)結(jié)果中KMO值為0.96>0.6,Bartlett球形度檢驗(yàn)的P值小于0.05,兩個檢驗(yàn)均通過。利用SPSS分析軟件進(jìn)行因子分析,前兩個因子的累計方差貢獻(xiàn)率為85.15%。其中,因子1在商務(wù)部、長期投資、發(fā)改委、服務(wù)業(yè)、管理者等12 個成分上有較大載荷,這些關(guān)鍵詞主要體現(xiàn)了投資者對實(shí)體投資的關(guān)注。因子2 在理財、投資公司、風(fēng)險投資、融資等11個成分上有較大載荷,這些關(guān)鍵詞大部分反映了投資者對金融投資的關(guān)注。

    (5)合成投資網(wǎng)絡(luò)搜索指數(shù)。首先將各因子的成分得分系數(shù)作為關(guān)鍵詞的權(quán)重系數(shù),合成兩個因子的月度數(shù)據(jù)。然后以兩個因子的月度數(shù)據(jù)為特征,以房地產(chǎn)投資累計增長率月度數(shù)據(jù)為目標(biāo)變量,利用Python建立隨機(jī)森林回歸模型,得到兩個因子的特征重要性。模型輸出特征的重要性如表4 所示。最后把每個特征因子的重要性得分作為每個特征的權(quán)重系數(shù),合成了月度投資網(wǎng)絡(luò)搜索指數(shù)(Investment Network Search Index,INSI)。

    表4 各因子重要性得分

    1.3 模型建立

    LSTM 模型最早由Hochreiter 和Schmidhuber 在1997年提出。由于能更好地記憶時間序列的長期依賴關(guān)系,因此被廣泛用于處理序列信息。Phased-LSTM 模型是在基礎(chǔ)模型上添加了新的時間門kt來擴(kuò)展LSTM,該門的開啟和關(guān)閉由具有一定頻率范圍的參數(shù)振蕩控制。模型結(jié)構(gòu)如圖1所示,原文中的定義公式如下:

    圖1 MM-LSTM模型結(jié)構(gòu)

    其中,it、ft、ot分別為輸入門、遺忘門、輸出門的計算公式;°表示矩陣按元素相乘運(yùn)算;xt和ht分別為輸入特征向量和隱藏輸出向量;σi、σf、σo、σc為激活函數(shù);σh與權(quán)重矩陣W連接不同的輸入門和記憶單元的輸出;b為偏置項(xiàng),是模型訓(xùn)練時要學(xué)習(xí)的參數(shù)。

    與傳統(tǒng)的LSTM 模型相比,Phased-LSTM 模型的更新可以在不規(guī)則的時間點(diǎn)tj上執(zhí)行。本文用簡寫的cj=ctj表示儲存記憶單元在tj時刻的狀態(tài),用cj-1表示儲存記憶單元先前tj-1時刻的狀態(tài),則式(2)中ct和ht更新的方程為:

    其中,cj只有在時間門打開的時候才允許更新,從而實(shí)現(xiàn)一個周期的一部分區(qū)間產(chǎn)生儲存記憶單元的更新,這使得模型具備了處理不同頻率數(shù)據(jù)的能力。本文合成的INSI為日度數(shù)據(jù),ISI為月度數(shù)據(jù),GDP_IG為季度數(shù)據(jù),不同的模型需要的數(shù)據(jù)和數(shù)據(jù)頻率存在差異。AR 模型僅利用GDP_IG自身的滯后性建模。MIDAS 模型雖然可以用多元高頻數(shù)據(jù)來預(yù)測低頻數(shù)據(jù),但多元高頻數(shù)據(jù)由于倍差較大,需要調(diào)整為同頻的,本文把INSI降頻為月度數(shù)據(jù),記為INSI_M。MM-LSTM 模型能直接利用月度數(shù)據(jù)、日度數(shù)據(jù)對季度數(shù)據(jù)進(jìn)行預(yù)測,有效地利用了混頻數(shù)據(jù)的信息。

    2 實(shí)證研究

    2.1 基于AR模型的實(shí)證

    AR 模型通過利用時間序列自身的滯后性,對未來指標(biāo)的變動進(jìn)行預(yù)測。利用該模型進(jìn)行時間序列分析時,需要檢驗(yàn)被研究的序列是否為平穩(wěn)序列,以防止偽回歸。本文采用EViews 10.0 對GDP_IG數(shù)據(jù)進(jìn)行建模,根據(jù)AIC與SC 準(zhǔn)則,其中添加截距項(xiàng)的未進(jìn)行差分的序列通過平穩(wěn)性檢驗(yàn),拒絕了存在單位根的原假設(shè),序列平穩(wěn)。同時,根據(jù)序列的自相關(guān)圖和偏自相關(guān)圖(圖略)可知,自相關(guān)圖呈現(xiàn)拖尾特征,偏自相關(guān)圖呈現(xiàn)一階截尾特征,因此應(yīng)建立AR 模型。通過各種模型參數(shù)建模比較,根據(jù)AIC、SC、HQ值和回歸方程各系數(shù)的顯著性,最終確定應(yīng)建立添加截距項(xiàng)的AR(1)模型。為了進(jìn)行擴(kuò)展預(yù)測,將樣本內(nèi)估計區(qū)間設(shè)置為2011 年第一季度至2020 年第四季度,樣本外預(yù)測區(qū)間擴(kuò)展到2022年第四季度①本文采用靜態(tài)預(yù)測,每次向前預(yù)測一個數(shù)據(jù),并且假設(shè)上一期的預(yù)測值等于真實(shí)值。。其中,樣本內(nèi)均方根誤差為2.139,樣本外均方根誤差為5.708。

    2.2 基于AR-MIDAS模型的實(shí)證

    MIDAS 模型最初是由Ghysels 等在2004 年提出的,MIDAS 模型的顯著特點(diǎn)是可以直接使用不同頻率的數(shù)據(jù)進(jìn)行建模?;A(chǔ)的MIDAS模型主要包括單變量和多變量MIDAS模型,以及加入自回歸項(xiàng)的AR-MIDAS模型和h步向前預(yù)測的MIDAS 模型。因宏觀經(jīng)濟(jì)運(yùn)行的慣性作用,經(jīng)濟(jì)時間序列普遍存在自相關(guān)性,因變量的若干期會對當(dāng)期產(chǎn)生影響,因此此處采用含有自回歸的AR-MIDAS 模型。為了研究不同類型數(shù)據(jù)和不同MIDAS 模型預(yù)測差異,本文分別使用ISI和INSI_M數(shù)據(jù)建立單因素AR-MIDAS 模型,并同時使用以上兩類數(shù)據(jù)建立雙因素AR-MIDAS模型進(jìn)行對比研究。

    前文基于AR 模型的分析已經(jīng)討論了GDP_IG的滯后階數(shù),這里直接令A(yù)R-MIDAS模型的自回歸項(xiàng)滯后階數(shù)p=1。利用EViews 軟件建模時,還需要確定高頻回歸器抽取的樣本數(shù)和高頻數(shù)據(jù)的滯后階數(shù)。本文選取阿爾蒙多項(xiàng)式,軟件可以根據(jù)AIC準(zhǔn)則自動尋找高頻回歸器的最佳滯后階數(shù),這里最大滯后階數(shù)設(shè)置為12,即高頻回歸器最多使用前12 個月的數(shù)據(jù)來擬合當(dāng)前的GDP_IG數(shù)據(jù)。關(guān)于滯后階數(shù)的選擇,通過各種模型建模比較,根據(jù)回歸方程各系數(shù)的顯著性和樣本區(qū)間內(nèi)外的預(yù)測誤差等參數(shù)最終確定。根據(jù)ISI建模得到的模型為AR(1)-MIDAS(4,2),模型的擬合優(yōu)度為0.976,樣本內(nèi)均方根誤差為0.910,樣本外均方根誤差為1.243。根據(jù)INSI_M建模得到的模型為AR(1)-MIDAS(3,0),樣本內(nèi)均方根誤差為2.332,樣本外均方根誤差為5.350。結(jié)合ISI、INSI_M構(gòu)建AR-MIDAS模型,為了保證各模型擬合效果和預(yù)測結(jié)果具有可比性,此處直接引用上述單因素AR-MIDAS建模的各種參數(shù)設(shè)置進(jìn)行建模,記為AR(1)-MIDAS_tn。模型的樣本內(nèi)均方根誤差為0.763,樣本外均方根誤差為1.212,其樣本內(nèi)外預(yù)測誤差與單因素AR-MIDAS 模型相比均有所降低。

    2.3 基于MM-LSTM模型的實(shí)證

    通過上述AR模型、MIDAS模型可知,合成的ISI、INSI數(shù)據(jù)可以用于中國經(jīng)濟(jì)波動預(yù)測,這使得構(gòu)建的MM-LSTM 模型具有理論基礎(chǔ)。由于GDP_IG與ISI和INSI數(shù)據(jù)存在較大差異,對模型進(jìn)行訓(xùn)練和測試之前,先對數(shù)據(jù)進(jìn)行歸一化處理。歸一化公式如下:

    其中,xi表示序列X中的第i個值,max(X)、min(X)分別表示序列X的最大值、最小值。與前面模型的樣本內(nèi)區(qū)間和樣本外區(qū)間設(shè)置相同,2011—2020 年的數(shù)據(jù)作為訓(xùn)練集,2021—2022年的數(shù)據(jù)作為測試集。

    (1)MM-LSTM模型參數(shù)設(shè)置

    通過設(shè)置合適的參數(shù),可以獲得良好的深度學(xué)習(xí)性能,本文的目標(biāo)是預(yù)測經(jīng)濟(jì)波動,故選用均方根誤差(RMSE)作為損失函數(shù)。優(yōu)化器方面,本文采用Adam 優(yōu)化器進(jìn)行優(yōu)化訓(xùn)練。Adam 優(yōu)化器是目前最常用的算法,與其他優(yōu)化器相比,Adam算法收斂速度最快,學(xué)習(xí)效率更高。各網(wǎng)絡(luò)層的激活函數(shù)均采用ReLu 函數(shù),模型的最后一層采用線性linear 函數(shù)。為防止出現(xiàn)過擬合現(xiàn)象,對INSI數(shù)據(jù)設(shè)置Dropout層,網(wǎng)絡(luò)節(jié)點(diǎn)的舍棄率設(shè)為0.2。

    (2)MM-LSTM模型實(shí)證

    為了與上文的MIDAS 模型進(jìn)行對比研究,分別使用ISI、使用INSI以及同時使用ISI和INSI數(shù)據(jù)分開建模。三種模型分別記為MM-LSTM_t、MM-LSTM_n、MM-LSTM_tn。反復(fù)修改參數(shù)后,根據(jù)模型的訓(xùn)練集誤差和測試集誤差得到最優(yōu)參數(shù)設(shè)置值,見表5。

    表5 最優(yōu)參數(shù)設(shè)置

    表5 中,n_test表示測試集使用的月份數(shù),epochs和batch_size分別代表模型訓(xùn)練時的迭代次數(shù)和一次訓(xùn)練所取的樣本數(shù),nodes代表隱藏神經(jīng)元個數(shù),month_t和month_p分別代表訓(xùn)練時使用當(dāng)季和上一季的月份數(shù)量,gdp_p代表GDP_IG數(shù)據(jù)的滯后期數(shù)。以MM-LSTM_t模型為例,從表5 中的結(jié)果來看,其在預(yù)測當(dāng)季的GDP_IG時,使用當(dāng)季三個月和上季最后一個月的數(shù)據(jù)將取得相對較好的結(jié)果。

    根據(jù)表5 還可以發(fā)現(xiàn),使用ISI數(shù)據(jù)建立的模型訓(xùn)練集、測試集誤差與使用INSI數(shù)據(jù)建立的模型誤差都相對較小,這與MIDAS 模型的結(jié)果一致。原因可能是,ISI數(shù)據(jù)是由各宏觀經(jīng)濟(jì)指標(biāo)合成的,與經(jīng)濟(jì)波動直接相關(guān);而INSI數(shù)據(jù)是根據(jù)與投資相關(guān)指標(biāo)的搜索數(shù)據(jù)建立的聯(lián)系,具有一定的不確定性。

    由上述結(jié)果可知,使用混頻數(shù)據(jù)的MM-LSTM模型訓(xùn)練集和預(yù)測集誤差都最小。為了驗(yàn)證該模型是否具有提前預(yù)測GDP_IG數(shù)據(jù)的功能,將模型參數(shù)month_t分別設(shè)置為1、2、3。通過調(diào)整表5中其他參數(shù)得到模型的訓(xùn)練和測試誤差,如表6所示。

    表6 MM-LSTM_tn模型訓(xùn)練和預(yù)測誤差比較

    表6中,30+1表示使用每季度前30天的INSI數(shù)據(jù)和1 個月的ISI數(shù)據(jù)預(yù)測該季度的GDP_IG數(shù)據(jù)。對比發(fā)現(xiàn),僅使用當(dāng)季前30天和1個月數(shù)據(jù)就可以比國家統(tǒng)計局提前約60天預(yù)測該GDP_IG數(shù)據(jù)。而使用當(dāng)季前60天和2 個月數(shù)據(jù)預(yù)測時其預(yù)測結(jié)果與AR-MIDAS 模型精度相當(dāng),但其可提前30天預(yù)測該GDP_IG數(shù)據(jù)。因此,可以利用這一預(yù)測的先行性,為政府部門制定經(jīng)濟(jì)政策提供決策依據(jù)。

    3 模型預(yù)測效果評估

    3.1 各模型預(yù)測性能比較

    利用實(shí)證分析部分各模型的估計結(jié)果,計算各模型樣本內(nèi)擬合誤差和樣本外預(yù)測誤差,如表7所示。

    表7 各模型樣本內(nèi)外預(yù)測性能比較

    第一,根據(jù)投資合成的INSI和ISI數(shù)據(jù)與GDP_IG序列之間存在正向關(guān)系,說明自2011年以來,較高的投資率是引起經(jīng)濟(jì)波動的重要因素,若想經(jīng)濟(jì)有較大的波動,則需要用適度的投資刺激我國經(jīng)濟(jì),但同時應(yīng)警惕過高的投資率所帶來的GDP大幅度波動。

    第二,從所有模型的估計結(jié)果來看,MIDAS 模型與MM-LSTM模型比僅利用GDP_IG數(shù)據(jù)的AR模型預(yù)測效果有所提高,說明本文合成的投資統(tǒng)計指數(shù)和投資網(wǎng)絡(luò)搜索指數(shù)對中國經(jīng)濟(jì)波動的預(yù)測具有優(yōu)化作用。其中,投資統(tǒng)計指數(shù)數(shù)據(jù)對經(jīng)濟(jì)波動的預(yù)測優(yōu)化效果更好。但綜合來看,預(yù)測經(jīng)濟(jì)波動取得了最好的效果。

    第三,從使用混頻數(shù)據(jù)估計結(jié)果來看,MM-LSTM 模型的樣本內(nèi)擬合誤差和樣本外預(yù)測誤差比AR(1)-MIDAS模型的預(yù)測精度分別提高了3.5和38.4個百分點(diǎn)。究其原因,AR(1)-MIDAS 模型直接使用了月度投資統(tǒng)計指數(shù)數(shù)據(jù),但把日度投資網(wǎng)絡(luò)搜索指數(shù)數(shù)據(jù)降頻為月度數(shù)據(jù);而MM-LSTM 模型直接使用了月度數(shù)據(jù)和日度數(shù)據(jù)。這說明MM-LSTM模型更為有效地利用了樣本信息,得到了更加精確的預(yù)測結(jié)果,同時也說明了MM-LSTM模型可以更好地擬合我國宏觀經(jīng)濟(jì)波動。

    第四,從使用投資網(wǎng)絡(luò)搜索指數(shù)數(shù)據(jù)來看,MIDAS 模型和MM-LSTM模型在添加了投資網(wǎng)絡(luò)搜索指數(shù)數(shù)據(jù)后,樣本外預(yù)測誤差分別比僅使用投資統(tǒng)計數(shù)據(jù)的樣本外預(yù)測誤差降低了3.1和52.2個百分點(diǎn)。這說明本文添加網(wǎng)絡(luò)搜索數(shù)據(jù)提高了模型的估計精度?;谶@一特點(diǎn),本文認(rèn)為網(wǎng)絡(luò)搜索數(shù)據(jù)作為反映用戶關(guān)注變化的即時指標(biāo),為傳統(tǒng)的監(jiān)測方法提供了有效補(bǔ)充,可以用來預(yù)測經(jīng)濟(jì)波動。

    綜上,從MM-LSTM模型和對比模型對中國宏觀經(jīng)濟(jì)波動的預(yù)測結(jié)果可以看出,MM-LSTM 模型在中國宏觀經(jīng)濟(jì)波動預(yù)測中的有效性得到了證實(shí)。

    3.2 模型預(yù)測結(jié)果比較

    從圖2 來看,AR 模型對GDP_IG的估計結(jié)果類似于將原始值往后平移1 期,當(dāng)GDP_IG出現(xiàn)較大波動時,不能及時預(yù)測。例如,2020年第一季度、2021年第一季度和2022年第二季度,預(yù)測結(jié)果出現(xiàn)了很大偏差。從下頁圖3來看,與AR模型相比,加入投資統(tǒng)計、網(wǎng)絡(luò)搜索指數(shù)數(shù)據(jù)后,AR(1)-MIDAS-tn模型準(zhǔn)確預(yù)測了2020年第一季度、2021 年第一季度和2022 年第二季度的波動,同時樣本外短期預(yù)測效果較好,但2022 年第二季度波動預(yù)測誤差較大。從下頁圖4 來看,MM-LSTM_tn 模型短期預(yù)測比MIDAS 模型更加準(zhǔn)確,中期預(yù)測誤差也相對較小,這說明本文提出的MM-LSTM 模型提高了樣本外短期和中期的預(yù)測精度。

    圖2 AR模型樣本內(nèi)擬合和樣本外預(yù)測結(jié)果對比

    圖3 AR(1)-MIDAS_tn模型樣本內(nèi)擬合和樣本外預(yù)測結(jié)果對比

    圖4 MM-LSTM_tn模型樣本內(nèi)擬合和樣本外預(yù)測結(jié)果對比

    4 結(jié)論

    實(shí)時、準(zhǔn)確地預(yù)測我國經(jīng)濟(jì)波動是宏觀經(jīng)濟(jì)研究中的重要課題。本文從宏觀和微觀兩個層面分析了投資對中國經(jīng)濟(jì)增長的影響??紤]到宏觀經(jīng)濟(jì)指標(biāo)間并不都是同頻數(shù)據(jù),以及大量的網(wǎng)絡(luò)搜索數(shù)據(jù)的有效性已得到學(xué)者們證實(shí),本文構(gòu)建了能夠利用季度數(shù)據(jù)、月度數(shù)據(jù)和日度數(shù)據(jù)的多源MM-LSTM 模型。通過選擇國內(nèi)生產(chǎn)總值指數(shù)累計值(上年同期=100)數(shù)據(jù)和月度統(tǒng)計指標(biāo)與網(wǎng)絡(luò)搜索日度指標(biāo),本文從投資角度對我國2011年1月1日至2022年12月31日的經(jīng)濟(jì)波動進(jìn)行了實(shí)證分析。本文的主要結(jié)論如下:

    (1)投資相關(guān)指標(biāo)與我國GDP_IG之間存在正向關(guān)系。加入投資統(tǒng)計指數(shù)或投資網(wǎng)絡(luò)搜索指數(shù)的兩種模型比基準(zhǔn)模型的預(yù)測誤差更小。這說明了用投資指標(biāo)預(yù)測我國經(jīng)濟(jì)波動的可行性。

    (2)添加投資網(wǎng)絡(luò)搜索數(shù)據(jù)提高了模型的估計精度。MIDAS 模型和MM-LSTM 模型在加入投資網(wǎng)絡(luò)搜索指數(shù)后,樣本外預(yù)測誤差分別降低了3.1和52.2個百分點(diǎn),說明網(wǎng)絡(luò)搜索數(shù)據(jù)對經(jīng)濟(jì)預(yù)測具有積極影響。

    (3)MM-LSTM 模型具有更小的預(yù)測誤差,提高了短期和中期的預(yù)測精度。本文構(gòu)建的MM-LSTM 模型充分利用了混頻數(shù)據(jù)的信息,與AR 模型、MIDAS 模型的最優(yōu)結(jié)果相比,能夠使樣本內(nèi)擬合誤差降低3.5個百分點(diǎn),樣本外預(yù)測誤差降低38.4 個百分點(diǎn),并且準(zhǔn)確預(yù)測了2020 年第一季度、2021 年第一季度、2022 年第二季度的波動,從而驗(yàn)證了模型的有效性。

    (4)本文利用混頻數(shù)據(jù)構(gòu)建MM-LSTM模型對經(jīng)濟(jì)波動進(jìn)行預(yù)測時,僅使用當(dāng)季前30天和1個月數(shù)據(jù)就可以比國家統(tǒng)計局提前約60 天預(yù)測該季度GDP_IG數(shù)據(jù),其預(yù)測精度高于AR模型。而使用當(dāng)季前60天和2個月數(shù)據(jù)預(yù)測時,其預(yù)測結(jié)果與AR-MIDAS 模型精度相當(dāng),但其預(yù)測結(jié)果可提前30天預(yù)測該季度GDP_IG數(shù)據(jù)。可以利用這一預(yù)測的先行性,為政府部門制定經(jīng)濟(jì)政策提供決策依據(jù)。

    猜你喜歡
    模型
    一半模型
    一種去中心化的域名服務(wù)本地化模型
    適用于BDS-3 PPP的隨機(jī)模型
    提煉模型 突破難點(diǎn)
    函數(shù)模型及應(yīng)用
    p150Glued在帕金森病模型中的表達(dá)及分布
    函數(shù)模型及應(yīng)用
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計的漸近分布
    3D打印中的模型分割與打包
    草草在线视频免费看| 丰满的人妻完整版| 嫩草影院精品99| 动漫黄色视频在线观看| 久久久久九九精品影院| 欧美成人a在线观看| 天天躁日日操中文字幕| 久久天躁狠狠躁夜夜2o2o| 欧美成人a在线观看| 久久午夜亚洲精品久久| 亚洲最大成人手机在线| 99精品久久久久人妻精品| 亚洲人与动物交配视频| 国产色爽女视频免费观看| 欧美不卡视频在线免费观看| 国产一区二区三区av在线 | 国产一区二区三区av在线 | 免费电影在线观看免费观看| 亚洲黑人精品在线| 内地一区二区视频在线| 亚洲av二区三区四区| 国产 一区精品| 丰满乱子伦码专区| 久久九九热精品免费| 亚洲不卡免费看| 国产免费男女视频| 日韩欧美精品v在线| 久久久久久大精品| 亚洲成人久久爱视频| 亚洲在线自拍视频| 国产熟女欧美一区二区| 欧美精品国产亚洲| 久9热在线精品视频| 国产综合懂色| 亚洲第一区二区三区不卡| 男女下面进入的视频免费午夜| 偷拍熟女少妇极品色| 日韩 亚洲 欧美在线| 精品久久久噜噜| 在线观看免费视频日本深夜| 国产主播在线观看一区二区| 欧美激情在线99| 欧美3d第一页| 亚洲中文日韩欧美视频| 2021天堂中文幕一二区在线观| 成年女人永久免费观看视频| 简卡轻食公司| 18+在线观看网站| 一区二区三区高清视频在线| 国产精品人妻久久久影院| 色精品久久人妻99蜜桃| 色av中文字幕| 三级男女做爰猛烈吃奶摸视频| 午夜福利欧美成人| 日韩av在线大香蕉| 99热这里只有精品一区| 一本久久中文字幕| 国产色爽女视频免费观看| 欧美最新免费一区二区三区| 搡女人真爽免费视频火全软件 | www.www免费av| 别揉我奶头 嗯啊视频| 黄色一级大片看看| 午夜a级毛片| 欧美日韩乱码在线| 亚洲avbb在线观看| 亚洲人成网站在线播放欧美日韩| 国产蜜桃级精品一区二区三区| 少妇猛男粗大的猛烈进出视频 | 在现免费观看毛片| 免费一级毛片在线播放高清视频| 成人毛片a级毛片在线播放| 色哟哟·www| 久久国内精品自在自线图片| 亚洲人成网站在线播| 少妇熟女aⅴ在线视频| 亚洲欧美日韩无卡精品| 欧美zozozo另类| 国产精品伦人一区二区| 欧美激情国产日韩精品一区| 亚洲国产精品sss在线观看| 超碰av人人做人人爽久久| 麻豆精品久久久久久蜜桃| 少妇的逼水好多| 国产乱人视频| 丰满乱子伦码专区| 亚洲av不卡在线观看| 成人国产综合亚洲| 一进一出好大好爽视频| 免费观看的影片在线观看| 久久精品国产鲁丝片午夜精品 | 1024手机看黄色片| 在线国产一区二区在线| 久久人人精品亚洲av| 搡老妇女老女人老熟妇| 久久精品国产鲁丝片午夜精品 | 国产精品美女特级片免费视频播放器| 成人一区二区视频在线观看| 在线观看av片永久免费下载| 精品久久久久久久久av| 乱系列少妇在线播放| 少妇的逼水好多| 日本五十路高清| 欧美性猛交╳xxx乱大交人| 亚洲成a人片在线一区二区| 亚洲国产日韩欧美精品在线观看| av天堂在线播放| 我要看日韩黄色一级片| 一卡2卡三卡四卡精品乱码亚洲| 国内精品一区二区在线观看| 草草在线视频免费看| 国产成人影院久久av| 国产又黄又爽又无遮挡在线| 身体一侧抽搐| 俄罗斯特黄特色一大片| 欧美国产日韩亚洲一区| 亚洲人成网站在线播放欧美日韩| 久久精品影院6| 午夜爱爱视频在线播放| 一个人看视频在线观看www免费| 婷婷色综合大香蕉| 午夜免费激情av| 国产成人a区在线观看| 麻豆久久精品国产亚洲av| 变态另类丝袜制服| 91麻豆av在线| 欧美最新免费一区二区三区| 色吧在线观看| 99在线视频只有这里精品首页| 日日夜夜操网爽| 99国产极品粉嫩在线观看| 俄罗斯特黄特色一大片| 中文在线观看免费www的网站| 变态另类成人亚洲欧美熟女| 亚洲av免费在线观看| 麻豆成人午夜福利视频| 国产成人一区二区在线| 99久久精品热视频| 春色校园在线视频观看| 免费av观看视频| 国产精品福利在线免费观看| 丝袜美腿在线中文| 黄色欧美视频在线观看| 亚洲乱码一区二区免费版| 国产伦人伦偷精品视频| 一本精品99久久精品77| 成人国产一区最新在线观看| 国产黄色小视频在线观看| 极品教师在线视频| 亚洲第一区二区三区不卡| 亚洲av日韩精品久久久久久密| 毛片女人毛片| 国产真实乱freesex| www.色视频.com| 免费看美女性在线毛片视频| 最近最新免费中文字幕在线| 国产美女午夜福利| 精品国内亚洲2022精品成人| 午夜福利在线观看免费完整高清在 | 婷婷精品国产亚洲av| 亚洲精品亚洲一区二区| 国产国拍精品亚洲av在线观看| 亚洲精品久久国产高清桃花| 欧美中文日本在线观看视频| 亚洲三级黄色毛片| 欧美日韩综合久久久久久 | 亚洲七黄色美女视频| 久久精品国产亚洲网站| 99久久精品一区二区三区| 小说图片视频综合网站| 亚洲中文字幕一区二区三区有码在线看| 亚洲中文字幕日韩| 亚洲精品在线观看二区| 亚洲精品在线观看二区| 亚洲成人中文字幕在线播放| 亚洲av日韩精品久久久久久密| 一区二区三区免费毛片| xxxwww97欧美| 国内久久婷婷六月综合欲色啪| 伊人久久精品亚洲午夜| 美女免费视频网站| 欧美极品一区二区三区四区| 波多野结衣巨乳人妻| 日本-黄色视频高清免费观看| 国产免费av片在线观看野外av| 亚洲黑人精品在线| 中文字幕精品亚洲无线码一区| 91在线观看av| 人人妻,人人澡人人爽秒播| 亚洲av成人av| 色精品久久人妻99蜜桃| 欧美日韩亚洲国产一区二区在线观看| 成人国产综合亚洲| 亚洲美女视频黄频| 午夜亚洲福利在线播放| 欧美中文日本在线观看视频| 黄色配什么色好看| 欧美bdsm另类| 一个人观看的视频www高清免费观看| 亚洲久久久久久中文字幕| 亚洲欧美日韩东京热| 在线观看一区二区三区| 老熟妇仑乱视频hdxx| 日本免费一区二区三区高清不卡| 在线观看av片永久免费下载| 欧美日韩黄片免| 永久网站在线| 伦理电影大哥的女人| 国产精品人妻久久久久久| 人妻制服诱惑在线中文字幕| 日韩在线高清观看一区二区三区 | 极品教师在线视频| 可以在线观看的亚洲视频| 国产真实乱freesex| 乱人视频在线观看| 久久6这里有精品| 成年女人永久免费观看视频| 露出奶头的视频| 91久久精品国产一区二区三区| 欧美黑人欧美精品刺激| 久久精品国产亚洲av香蕉五月| 在线观看一区二区三区| eeuss影院久久| 国产主播在线观看一区二区| 日本免费a在线| videossex国产| 欧美另类亚洲清纯唯美| 国产精品乱码一区二三区的特点| 国产女主播在线喷水免费视频网站 | 久久欧美精品欧美久久欧美| 日本五十路高清| 天堂动漫精品| 禁无遮挡网站| 亚洲av不卡在线观看| 欧美日韩亚洲国产一区二区在线观看| 国产大屁股一区二区在线视频| 国产成人一区二区在线| 国产精品福利在线免费观看| 国国产精品蜜臀av免费| 午夜福利在线观看吧| 一级黄片播放器| 99久久精品一区二区三区| 免费在线观看影片大全网站| 亚洲一级一片aⅴ在线观看| 欧美性猛交╳xxx乱大交人| 少妇猛男粗大的猛烈进出视频 | 色噜噜av男人的天堂激情| 欧美性猛交黑人性爽| 男人舔奶头视频| 精品午夜福利在线看| 男人和女人高潮做爰伦理| 国产aⅴ精品一区二区三区波| 精品一区二区三区视频在线观看免费| 免费不卡的大黄色大毛片视频在线观看 | 国产精品国产高清国产av| 美女大奶头视频| 欧美zozozo另类| 天堂影院成人在线观看| 国产精品98久久久久久宅男小说| 婷婷六月久久综合丁香| 最后的刺客免费高清国语| 精品日产1卡2卡| 国产乱人伦免费视频| 久久久久久久久久黄片| 精品久久国产蜜桃| www.www免费av| 免费一级毛片在线播放高清视频| 最好的美女福利视频网| 久久久久久久午夜电影| 亚洲人成网站在线播| 全区人妻精品视频| 亚洲精品色激情综合| 国产淫片久久久久久久久| 91在线观看av| 久久精品国产自在天天线| 91久久精品国产一区二区三区| x7x7x7水蜜桃| 中国美白少妇内射xxxbb| 久久这里只有精品中国| 国产精品一及| 欧美在线一区亚洲| 真实男女啪啪啪动态图| 夜夜看夜夜爽夜夜摸| 国产私拍福利视频在线观看| videossex国产| 亚洲欧美激情综合另类| 精品无人区乱码1区二区| 国内精品宾馆在线| 亚洲一区高清亚洲精品| av黄色大香蕉| 免费看光身美女| 国产午夜福利久久久久久| 99久久久亚洲精品蜜臀av| 久久人人精品亚洲av| 亚洲精品在线观看二区| 精品人妻1区二区| 日本爱情动作片www.在线观看 | 赤兔流量卡办理| 成人高潮视频无遮挡免费网站| 亚洲美女黄片视频| 99九九线精品视频在线观看视频| 午夜爱爱视频在线播放| 真人做人爱边吃奶动态| 午夜福利在线在线| 欧美成人性av电影在线观看| 国产色爽女视频免费观看| 成人国产一区最新在线观看| 成人午夜高清在线视频| 中文字幕免费在线视频6| 亚洲成人免费电影在线观看| 精品久久久久久久久亚洲 | 日韩大尺度精品在线看网址| 国产一区二区三区视频了| 18禁在线播放成人免费| 成人二区视频| .国产精品久久| 精品人妻视频免费看| 国产精品精品国产色婷婷| 乱系列少妇在线播放| 久久久久免费精品人妻一区二区| 亚洲电影在线观看av| 黄色女人牲交| 真人一进一出gif抽搐免费| 精品福利观看| 亚洲精品影视一区二区三区av| 亚洲最大成人中文| 岛国在线免费视频观看| 国产精品国产高清国产av| 日本免费a在线| 国产精品不卡视频一区二区| 色哟哟·www| 人妻丰满熟妇av一区二区三区| 丰满的人妻完整版| 国产av一区在线观看免费| www日本黄色视频网| 97超级碰碰碰精品色视频在线观看| 中文字幕av成人在线电影| 日韩亚洲欧美综合| 国产女主播在线喷水免费视频网站 | 午夜福利视频1000在线观看| 一级黄色大片毛片| 18禁裸乳无遮挡免费网站照片| 一个人看的www免费观看视频| 校园人妻丝袜中文字幕| 人妻制服诱惑在线中文字幕| 日本免费a在线| 国产午夜精品论理片| 色综合站精品国产| 少妇的逼水好多| 有码 亚洲区| 久久久久久久亚洲中文字幕| 99热只有精品国产| 人妻丰满熟妇av一区二区三区| 亚洲精品色激情综合| 国产又黄又爽又无遮挡在线| 91麻豆精品激情在线观看国产| 99热这里只有精品一区| 九九久久精品国产亚洲av麻豆| 亚洲av不卡在线观看| 亚洲精品成人久久久久久| 成年免费大片在线观看| 亚洲国产精品sss在线观看| 免费看av在线观看网站| 亚洲欧美日韩无卡精品| 精品日产1卡2卡| 国产三级在线视频| 男插女下体视频免费在线播放| 亚洲中文日韩欧美视频| 亚洲av中文av极速乱 | 色av中文字幕| 欧美最新免费一区二区三区| 国产精品三级大全| 亚洲av免费高清在线观看| 99热这里只有是精品在线观看| 欧美最黄视频在线播放免费| 国产精品1区2区在线观看.| 久久国产乱子免费精品| 狂野欧美激情性xxxx在线观看| 男女视频在线观看网站免费| 国产亚洲欧美98| 九九爱精品视频在线观看| 免费av不卡在线播放| 99热网站在线观看| 午夜福利视频1000在线观看| 精品国内亚洲2022精品成人| 在线观看av片永久免费下载| 听说在线观看完整版免费高清| 亚洲精品456在线播放app | 免费搜索国产男女视频| av黄色大香蕉| 国产三级在线视频| 午夜免费激情av| 特级一级黄色大片| 99在线人妻在线中文字幕| 亚洲avbb在线观看| 国产男靠女视频免费网站| 男女视频在线观看网站免费| a级一级毛片免费在线观看| 亚洲精品一区av在线观看| 国模一区二区三区四区视频| 日日夜夜操网爽| av国产免费在线观看| 夜夜爽天天搞| 男女那种视频在线观看| 男人舔奶头视频| 欧美zozozo另类| 亚洲中文日韩欧美视频| 久久欧美精品欧美久久欧美| 丝袜美腿在线中文| 在线播放国产精品三级| 欧美最新免费一区二区三区| 免费高清视频大片| aaaaa片日本免费| av福利片在线观看| 91久久精品国产一区二区三区| 精品日产1卡2卡| 日韩欧美 国产精品| 国产精品国产高清国产av| 欧美精品国产亚洲| x7x7x7水蜜桃| 精品福利观看| 欧美性猛交黑人性爽| 欧美一级a爱片免费观看看| 丰满的人妻完整版| 国产中年淑女户外野战色| 亚洲av五月六月丁香网| 极品教师在线免费播放| 国产成人a区在线观看| 亚洲成人中文字幕在线播放| 国产成人aa在线观看| 亚洲电影在线观看av| 亚洲av第一区精品v没综合| 天堂av国产一区二区熟女人妻| 美女黄网站色视频| 国产成人av教育| 搞女人的毛片| 免费av毛片视频| 午夜激情欧美在线| 少妇高潮的动态图| 国产精品国产三级国产av玫瑰| 免费观看的影片在线观看| 小蜜桃在线观看免费完整版高清| 国内揄拍国产精品人妻在线| 日本一二三区视频观看| 国产午夜精品论理片| 日韩精品中文字幕看吧| 国产精品亚洲一级av第二区| 天美传媒精品一区二区| 国产精品一区二区免费欧美| 国产精品久久久久久亚洲av鲁大| 国产高清有码在线观看视频| 国内精品久久久久精免费| 国产一区二区在线观看日韩| 成年人黄色毛片网站| 少妇高潮的动态图| 嫩草影院精品99| 久久久久久久亚洲中文字幕| 精品久久久久久久久亚洲 | 一区二区三区激情视频| 偷拍熟女少妇极品色| 欧美潮喷喷水| 国内揄拍国产精品人妻在线| 精品人妻熟女av久视频| 成年版毛片免费区| 亚洲无线在线观看| 久久久久精品国产欧美久久久| 亚洲乱码一区二区免费版| 国产精品嫩草影院av在线观看 | 久久国产乱子免费精品| 好男人在线观看高清免费视频| 国内精品一区二区在线观看| 欧美性猛交╳xxx乱大交人| 老师上课跳d突然被开到最大视频| 国产精品福利在线免费观看| 日本欧美国产在线视频| 成人国产综合亚洲| 有码 亚洲区| 精品久久久久久久久久久久久| 亚洲精品日韩av片在线观看| 国产亚洲av嫩草精品影院| 国语自产精品视频在线第100页| 久久久久久久午夜电影| 中文在线观看免费www的网站| 欧美bdsm另类| 最新在线观看一区二区三区| 村上凉子中文字幕在线| 三级毛片av免费| 国产爱豆传媒在线观看| 久久久国产成人免费| 一级黄色大片毛片| 最近中文字幕高清免费大全6 | 日本免费一区二区三区高清不卡| 毛片一级片免费看久久久久 | av国产免费在线观看| x7x7x7水蜜桃| 国产爱豆传媒在线观看| 成人无遮挡网站| 成人永久免费在线观看视频| 色综合亚洲欧美另类图片| 国产美女午夜福利| 99视频精品全部免费 在线| 欧美激情在线99| 韩国av一区二区三区四区| 精品乱码久久久久久99久播| 亚洲一区高清亚洲精品| 男人舔奶头视频| 成人av在线播放网站| 欧美日韩黄片免| 最后的刺客免费高清国语| 少妇熟女aⅴ在线视频| 国产精品日韩av在线免费观看| 黄色丝袜av网址大全| 午夜福利高清视频| 男女之事视频高清在线观看| 日韩欧美在线二视频| 日韩欧美国产一区二区入口| 久久午夜福利片| 一个人看视频在线观看www免费| 精品久久国产蜜桃| 国产精品一区二区免费欧美| 国内揄拍国产精品人妻在线| 男人狂女人下面高潮的视频| 国产精品久久久久久亚洲av鲁大| 男人和女人高潮做爰伦理| 精品福利观看| 美女xxoo啪啪120秒动态图| 日日摸夜夜添夜夜添小说| 国产精品av视频在线免费观看| 人妻久久中文字幕网| 欧美bdsm另类| 亚洲久久久久久中文字幕| 国产精品久久久久久精品电影| 亚洲不卡免费看| 亚洲综合色惰| 日本一本二区三区精品| 波多野结衣高清作品| 亚洲av中文字字幕乱码综合| 亚洲最大成人av| 尾随美女入室| 最近中文字幕高清免费大全6 | 欧美日韩精品成人综合77777| 免费看日本二区| 午夜a级毛片| av黄色大香蕉| 97超视频在线观看视频| 国产在线精品亚洲第一网站| 99在线人妻在线中文字幕| a级一级毛片免费在线观看| 久久久精品欧美日韩精品| 日本熟妇午夜| 中亚洲国语对白在线视频| 国产精品久久久久久久电影| 亚洲四区av| 国产女主播在线喷水免费视频网站 | 中文亚洲av片在线观看爽| 国产成人影院久久av| 色哟哟哟哟哟哟| 少妇人妻精品综合一区二区 | 我的女老师完整版在线观看| 国产aⅴ精品一区二区三区波| 99热只有精品国产| 精品久久久久久久久久免费视频| 哪里可以看免费的av片| 窝窝影院91人妻| 极品教师在线免费播放| 欧美激情国产日韩精品一区| 在线观看美女被高潮喷水网站| 老司机深夜福利视频在线观看| 国产一区二区三区在线臀色熟女| 看片在线看免费视频| 免费av毛片视频| 99热这里只有精品一区| 丰满人妻一区二区三区视频av| 久久国产精品人妻蜜桃| 国产精品三级大全| 老司机深夜福利视频在线观看| 国产精品爽爽va在线观看网站| 亚洲第一电影网av| 国产主播在线观看一区二区| 亚洲人成网站在线播| 欧美日本亚洲视频在线播放| 久久久国产成人免费| 嫩草影视91久久| 麻豆av噜噜一区二区三区| 男人舔女人下体高潮全视频| 亚洲人成网站高清观看| 久久久久久国产a免费观看| 国产伦一二天堂av在线观看| 嫩草影院入口| 全区人妻精品视频| bbb黄色大片| 亚洲精品日韩av片在线观看| 久久久久精品国产欧美久久久| 国产精品一区二区三区四区久久| 国产精品人妻久久久久久| 久久6这里有精品| 精品一区二区三区视频在线观看免费| 无人区码免费观看不卡| 18禁黄网站禁片免费观看直播| 一本精品99久久精品77| 亚洲精品色激情综合| 免费观看精品视频网站| 乱人视频在线观看| 免费无遮挡裸体视频| 国模一区二区三区四区视频| 91久久精品电影网| 欧美丝袜亚洲另类 | 一本久久中文字幕| 九九在线视频观看精品| 18禁裸乳无遮挡免费网站照片| 亚洲av二区三区四区| 好男人在线观看高清免费视频| 成人国产一区最新在线观看| 中出人妻视频一区二区| www.色视频.com|