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

    我國上市公司債券信用評級指標體系研究

    2024-12-31 00:00:00姜昱汐高倫雨宋旭
    海南金融 2024年11期

    摘" "要:債券信用評級對債券信用風險防范具有重要作用,而涵蓋信息充分、能有效區(qū)分債券信息差異的債券信用評級指標體系是識別債券信用影響因素、提高債券違約鑒別精度、提升債券信用評級可靠性的重要基準。本文選取2012—2022年我國上市公司債券為樣本,引入風險披露指標增大信息充分性、構(gòu)建宏觀指標與債券其他特征的克羅內(nèi)克積指標增大信息差異性,并基于AdaBoost算法的預測結(jié)果反推最優(yōu)指標體系。通過計算加入風險披露等指標和宏觀指標與債券其他特征克羅內(nèi)克積指標前后的多個模型性能評估指標,驗證指標體系的有效性;將AdaBoost模型與其他常用機器學習模型預測結(jié)果對比,驗證根據(jù)AdaBoost方法確定的最終指標體系可靠性;通過計算指標在模型預測中的重要性程度,提高指標體系的可解釋性。研究發(fā)現(xiàn),在多個預測性能評估指標上,增加指標體系的信息充分性和信息差異性,可顯著優(yōu)于未加入此類指標的預測結(jié)果;基于宏觀指標與債券其他特征克羅內(nèi)克積的多個指標對債券違約預測有重大影響;上市公司風險信息披露情況等指標是債券違約預測的重要指標。

    關(guān)鍵詞:債券信用評級指標體系;信息充分性;信息差異性;債券違約鑒別;AdaBoost

    DOI:10.3969/j.issn.1003-9031.2024.11.005

    中圖分類號:F224.9;F812.5" " " " 文獻標識碼:A" " "文章編號:1003-9031(2024)11-0049-15

    一、引言及文獻綜述

    1981年國務(wù)院決定恢復國債發(fā)行,拉開了我國債券市場從無到有、長足發(fā)展的序幕。截至2023年12月,我國債券市場托管余額157.9萬億元,同比增長9.1%,實現(xiàn)了債券市場規(guī)模的穩(wěn)定增長,我國債券市場規(guī)模穩(wěn)居全球第二①。債券市場蓬勃發(fā)展的背后是債券違約事件的頻發(fā)和各種極端風險的挑戰(zhàn)(王培輝和張猛,2023),2023年,全國共有41只債券主體發(fā)生違約,違約債券余額高達308.19億元②。為了盡可能地規(guī)避風險,減少損失,實現(xiàn)我國債券市場的穩(wěn)健發(fā)展,對債券進行合理的違約風險評估、建立有效的信用評級指標體系顯得愈發(fā)重要。目前,在債券信用評級領(lǐng)域中,國內(nèi)外學者的研究多聚焦于以下兩個方面。

    一是債券信用評級影響因素研究。鄒昆侖和陸萍(2018)基于債券發(fā)行人的盈利能力、資本結(jié)構(gòu)、償債能力、成長能力、每股指標五個方面,考慮了一系列財務(wù)指標對債券信用評級的影響。Ma C et al.(2021)認為債券信息、債券發(fā)行主體信息、債券信用評級和宏觀經(jīng)濟指標等會影響債券信用水平。姚瀟等(2022)將債券違約事件視為一種外生事件,同時考慮發(fā)債主體基本信息和財務(wù)指標對于債券信用評級的影響。生柳榮等(2019)在篩選信用評級影響因素時,綜合考慮了涵蓋企業(yè)屬性、企業(yè)償債能力和杠桿水平等六個方面的定性因素和定量因素。蔣敏等(2021)基于上市公司債券違約樣本,從公司經(jīng)營狀況等內(nèi)部經(jīng)營問題和外部經(jīng)濟環(huán)境兩方面出發(fā),探究內(nèi)部及外部因素對債券違約的“聯(lián)合效應”。肖艷麗和向有濤(2021)以上市公司發(fā)行的企業(yè)債券作為研究樣本,實證發(fā)現(xiàn)發(fā)債主體的財務(wù)指標和債券自身屬性等非財務(wù)指標都會對債券信用評級產(chǎn)生一定的影響。陳學彬等(2021)基于債券發(fā)行公司的財務(wù)指標、債券市場指標、債券評級指標、宏觀經(jīng)濟環(huán)境指標、行業(yè)景氣指標、地區(qū)經(jīng)濟景氣指標等六類指標進行信用債違約風險預測。Zhou Lu et al.(2021)采用半年期形式的數(shù)據(jù),衡量上市狀態(tài)、所有權(quán)結(jié)構(gòu)、行業(yè)分類和信用評級降級信息四個方面對債券信用水平的影響。近年來,部分學者發(fā)現(xiàn)債券發(fā)行人的風險信息披露情況也會在一定程度上影響債券的信用水平,如吳育輝和唐浩博(2021)發(fā)現(xiàn)債券募集說明書中披露的風險信息對債券信用水平有顯著影響;吳建華等(2017)推導了信息偏誤下資產(chǎn)價值的條件分布、違約概率和信用價差的解析表達式,基于該模型發(fā)現(xiàn)信息披露下的財務(wù)信息扭曲在一定程度上會影響債券信用水平。

    二是債券信用評級指標篩選方法研究。由于確定的初始指標體系往往含有較多的備選指標,為了篩選出信息利用度高的指標,學者們采用了不同的方法對初始指標體系進行降維處理?,F(xiàn)有文獻多在構(gòu)建信用評級模型前便對指標體系進行篩選,如高秀存等(2022)采用因子分析法對初始指標體系進行篩選,最終確定了包含流動比率、短期債務(wù)比總債務(wù)等14個指標的最優(yōu)指標體系。王玉龍等(2022)分別依據(jù)相關(guān)性和冗余性對初始指標體系進行第一次和第二次指標篩選,最終選定營業(yè)凈利率、凈利潤、現(xiàn)金比率、財務(wù)費用和資產(chǎn)負債率等5個財務(wù)指標共同組成最優(yōu)指標體系。龐春潮和羅苑瑋(2023)綜合考慮了各指標間的顯著性、相關(guān)性和多重共線性,最終構(gòu)建了包含20個定量指標、24個定性指標的最優(yōu)指標體系。謝玲玲和范龍振(2020)首先將債券數(shù)據(jù)分為違約組和非違約組,并針對初始指標體系中的指標依次進行均值差異比較和相關(guān)性分析,最終建立了包括14個指標的最優(yōu)指標體系。徐舒玥和曹艷華(2023)利用Pearson相關(guān)分析法對全部特征指標進行篩選,最終保留20個變量組成最優(yōu)指標體系。

    通過回顧現(xiàn)有文獻,發(fā)現(xiàn)學者們已經(jīng)對債券信用評級指標體系構(gòu)建做出了諸多探索,也取得了眾多成果,但仍存在進一步研究的空間。

    一是債券信用評級海選指標信息涵蓋不全面。債券信用評級情況受多方面因素影響,在指標海選時,現(xiàn)有研究雖然是從債券信用評級影響因素出發(fā)去選擇指標,但多是考慮了發(fā)債主體的財務(wù)指標、非財務(wù)指標、宏觀指標和債券基本信息等結(jié)構(gòu)化指標,對債券發(fā)行人的風險披露信息并未考慮,導致構(gòu)建的債券信用評級指標體系涵蓋信息不全面,影響信用評級的可靠性。

    二是指標篩選方法不科學帶來信息損失。很多研究在構(gòu)建信用評級模型前便通過相關(guān)性等統(tǒng)計量對指標進行篩選,但指標篩選時無法保證保留的指標體系組合對債券違約的預測精度,沒有體現(xiàn)出指標體系應保證債券違約可能性大、其信用評分低的合理思路,這種篩選方法雖然可以降低信用評級指標體系規(guī)模,但可能帶來信息損失,進而影響最終信用評級可靠性。

    三是不同債券部分指標取值相同,指標差異性不強。現(xiàn)有文獻中所考慮的宏觀指標多為債券發(fā)行人所在省級層面的相關(guān)數(shù)據(jù),這些指標對相同省份公司發(fā)行的債券取值相同,當相同地區(qū)的債券數(shù)量較多時,指標體系中這些債券的宏觀指標取值便會相同,這會導致該類指標不能反映出宏觀因素對債券信用影響的差別,導致錯失重要影響因素。

    本文的研究目標是確定一個違約鑒別能力強、復雜度低且具備經(jīng)濟可解釋性的最優(yōu)債券信用評級指標體系。基于此目標,本文從信息充分性和信息差異性視角出發(fā),在指標海選、指標篩選和宏觀指標差異化方面進行改進。

    根據(jù)上述分析,本文的主要創(chuàng)新點如下:第一,在債券信用評級指標海選時,將發(fā)債主體財務(wù)指標、非財務(wù)指標、宏觀指標、債券基本信息指標和根據(jù)上市公司發(fā)布的文本信息處理得到的風險披露方面的指標(如管理層凈語調(diào)等)等均納入海選指標體系,保證指標體系涵蓋信息更充分;第二,本文基于AdaBoost的嵌入式特征選擇方法對指標體系進行降維處理,在模型訓練的同時完成最優(yōu)指標體系篩選構(gòu)建,保證用篩選后的指標體系進行信用評級可以取得較好的債券違約鑒別精度,減少因只考慮指標相關(guān)性等進行指標體系降維而導致的信息損失;第三,考慮到同省份不同債券宏觀指標取值相同導致該部分指標信息差異小的弊端,一方面采用城市層面宏觀指標,增大同省份債券宏觀指標信息差異,另一方面引入克羅內(nèi)克積構(gòu)建宏觀指標與債券其他特征指標的相互作用項,并將其加入指標體系,保證指標體系中該類指標能體現(xiàn)不同債券的信息差異性,有助于識別宏觀因素是否影響債券信用評級。

    二、我國上市公司債券信用評級指標體系構(gòu)建原理

    (一)樣本數(shù)據(jù)處理

    1.缺失值處理

    樣本數(shù)據(jù)的大量缺失可能會導致模型精度的下降,為了在不損失大量信息的基礎(chǔ)上盡可能地減小分析結(jié)果與真實情況之間的誤差,本文在數(shù)據(jù)預處理階段對樣本存在的缺失值進行處理(白鈺銘和姜昱汐,2024)。

    缺失值處理的整體思路是利用樣本數(shù)據(jù)中的最差值來填補缺失值。對于流動比率、速動比率等正向定量指標,考慮用式(1)來填補缺失值:

    其中,ν'代表用于填充缺失值的指標值,i代表債券(i=1,2,…,N),j代表指標,ν代表現(xiàn)有的原始指標值,N代表債券總數(shù)。若式(1)得到的ν'偏離了該指標的合理范圍,則用現(xiàn)有指標值中的最小值來填補。

    對于有形凈值債務(wù)率、長期資本負債率等負向定量指標,考慮用式(2)來填補缺失值:

    式(2)中各參數(shù)含義同式(1),若式(2)得到的ν'偏離了該指標的合理范圍,則用現(xiàn)有指標值中的最大值來填補。

    對于董事會規(guī)模、管理層平均年齡等區(qū)間指標,考慮用式(3)來填補缺失值:

    其中,q代表區(qū)間指標的最小值,q代表區(qū)間指標的最大值,(q+ q)/2代表區(qū)間指標的中間值。

    式(1)—(3)為定量指標的缺失值填補方法,對于股權(quán)性質(zhì)等定性指標,考慮按照類別進行等級劃分,并在[0,1]范圍內(nèi)為其賦分,然后用最低分來填補缺失值。

    2.指標數(shù)據(jù)異常值處理和標準化

    在缺失值填補完成后,對獲得的完整數(shù)據(jù)進行前后1%縮尾的異常值處理(白鈺銘和姜昱汐,2024)。同時,考慮到樣本數(shù)據(jù)間的單位和量綱不同,本文在數(shù)據(jù)預處理階段對樣本數(shù)據(jù)進行標準化處理,即將所有定量指標數(shù)據(jù)歸于[0,1]之間,便于后續(xù)模型的構(gòu)建與評估,具體標準化公式如下(遲國泰等,2016):

    式(4)—(6)為定量指標的標準化公式,其中,式(4)為正向指標的歸一化過程,式(5)為負向指標的歸一化過程,式(6)為區(qū)間指標的歸一化過程。公式中各參數(shù)含義為:xij代表第i只債券第j個指標標準化后的數(shù)值,其它參數(shù)含義同前文。

    3.基于SMOTE方法的非均衡樣本處理

    對于我國上市公司債券信用評級問題而言,違約債券的數(shù)量遠少于非違約債券,因此我國上市公司債券數(shù)據(jù)為非均衡樣本。為了避免非均衡樣本數(shù)據(jù)可能導致的分類預測精度降低等問題,本文采用SMOTE方法對樣本進行均衡化處理。

    不同于隨機過采樣技術(shù)直接復制合成新的少數(shù)類樣本即違約樣本來解決樣本非均衡問題,欠采樣技術(shù)直接刪除多數(shù)類樣本即非違約樣本使樣本數(shù)據(jù)達到均衡,本文采用SMOTE方法通過對違約樣本進行插值來創(chuàng)建新的違約樣本,在實現(xiàn)樣本均衡的同時,幫助提高模型的泛化能力和準確性。

    用SMOTE方法實現(xiàn)上市公司債券樣本均衡的過程如下:首先,從違約債券中隨機選擇一個樣本計算它與其k個最近鄰樣本之間的距離,將計算得到的差值乘以0到1之間的隨機數(shù),然后將結(jié)果添加到違約樣本中作為新的合成樣本。重復上述過程,直到樣本數(shù)據(jù)達到均衡。

    (二)基于AdaBoost的上市公司債券信用評級指標體系構(gòu)建與模型精度評估

    信用評級是一個全面評估債券違約風險的過程,合理有效的信用評級要求能夠充分反映債券的違約可能性,即能夠滿足債券違約可能性大、其信用評級低的合理條件。因此,本文以債券違約預測精度高和指標體系規(guī)模合理為指標體系建立的原則,確立債券信用評級指標體系。

    1.指標體系的構(gòu)建

    (1)初始指標體系構(gòu)建

    參照惠譽、標普、穆迪、中誠信等國內(nèi)外權(quán)威評級機構(gòu)的信用評級指標體系,同時結(jié)合已有的經(jīng)典文獻,從發(fā)債主體內(nèi)部財務(wù)指標、發(fā)債主體內(nèi)部非財務(wù)指標、債券基本信息指標和宏觀指標等多個方面,選取債券信用評級可能的影響因素,建立本文的初始指標體系,具體包括各種結(jié)構(gòu)化指標和根據(jù)上市公司發(fā)布的文本信息處理得到的風險披露方面的指標(如管理層凈語調(diào)等)。

    (2)基于數(shù)據(jù)可獲性的指標第一次篩選

    由于部分指標數(shù)據(jù)無法獲取或存在大量缺失,若不在模型訓練前對指標進行處理可能會導致模型精度的下降,因此,本文基于指標數(shù)據(jù)的可獲得性對初始指標體系進行第一次篩選。

    (3)宏觀指標差異化處理

    為保證每個指標都能反映不同債券的差異,更好地識別債券信用風險因素,本文構(gòu)造了宏觀指標與債券其他特征指標的相互作用項,即通過計算宏觀指標與債券其他特征指標的克羅內(nèi)克積,得到具有債券差異性的新的宏觀指標,并將得到的新的特征指標納入指標池,參與指標的第二次篩選。

    宏觀指標差異化處理的目的:通過計算宏觀指標與債券其他特征指標的克羅內(nèi)克積來構(gòu)造相互作用項,可以保證每一只債券的宏觀指標都能在最大程度上反映債券的真實信息,且在不同債券間具有差異性,進而能更準確識別信用風險關(guān)鍵因素;考慮到初始指標池中可能含有大量無法獲取數(shù)據(jù)的指標,若直接計算宏觀指標與債券其他特征指標的克羅內(nèi)克積,可能會產(chǎn)生大量數(shù)據(jù)缺失嚴重的指標,造成信息噪聲增加,同時增加后續(xù)計算工作量,因此本文對第一次篩選后的宏觀指標做債券差異化處理。宏觀指標差異化處理的具體公式如下(Leippold M et al.,2022):

    (4)基于AdaBoost的指標第二次篩選

    為了進一步提高模型的預測精度,同時降低模型的復雜度,本文利用基于AdaBoost的嵌入式特征選擇方法對指標體系進行降維處理,在模型訓練的同時完成指標第二次篩選。

    對模型訓練得到的結(jié)果進行精度和指標體系精簡性的對比,依據(jù)結(jié)果反推同時滿足預測精度高且模型更精簡的最優(yōu)指標體系。

    AdaBoost算法作為一種重要的集成學習算法,主要通過不斷將多個弱學習器進行訓練,最終加權(quán)組合為一個強學習器,實現(xiàn)模型精度的提高,并對樣本準確分類(Freund Y and Schapire R E,1997)。具體計算步驟如下:

    假設(shè)債券樣本的訓練集為S=(x,y),(x,y),…,(x,y),其中,xR為q維債券特征空間的樣本,y0,1(i=1,2,…,N)為債券違約狀態(tài)(其中yi=1代表債券違約,否則為不違約)。

    初始化債券樣本數(shù)據(jù)的權(quán)值分布,對每一個訓練樣本賦予相同的權(quán)值w:

    進行m(m=1,2,…,M)輪迭代,依次訓練每個弱分類器Cm(x):

    選取弱分類器決策樹Cm(x)訓練債券樣本數(shù)據(jù)的權(quán)值,計算每一次迭代時弱分類器的錯誤率em:

    計算弱分類器Cm(x)在最終分類器中的權(quán)重α:

    更新債券樣本數(shù)據(jù)的權(quán)重w:

    其中,Z為:

    對債券樣本進行反復迭代,直至迭代次數(shù)m等于設(shè)置的迭代次數(shù)M,將多次迭代得到的各個弱分類器進行加權(quán)組合,得到f(x):

    通過符號函數(shù)sign(·)的作用,得到最終的強分類器G(x):

    2.模型評估

    為了更好地評估模型的分類性能,本文引入了5個經(jīng)典的判別精度指標對模型進行精度檢驗,分別為準確率(Accuracy)、AUC(ROC曲線下面積)、第一類錯誤(Type I Error)、第二類錯誤(Type Ⅱ Error)以及幾何平均值(G-mean),各精度指標的含義及計算參照經(jīng)典文獻(遲國泰等,2020),本文不再贅述。選取的上述5個判別精度指標既涵蓋了對模型整體性能的評估,又考慮到了模型的分類性能,可以較好地衡量模型的綜合精度。

    三、基于AdaBoost的上市公司債券信用評級指標體系構(gòu)建

    (一)樣本選取及數(shù)據(jù)來源

    1.樣本選取

    考慮到信用評級應該對債券違約有預警作用,因此本文通過構(gòu)建t-1年時間窗對我國上市公司發(fā)行的38233只在市債券進行違約預測,并依據(jù)違約預測的精度結(jié)果反推確定債券信用評級指標體系,即所使用的指標數(shù)據(jù)為2012—2021年,對應的違約狀態(tài)為2013—2022年。

    2.樣本數(shù)據(jù)來源

    本文的債券樣本數(shù)據(jù)來自多個渠道,債券基本信息、發(fā)債主體內(nèi)部財務(wù)指標數(shù)據(jù)和非財務(wù)指標數(shù)據(jù)從國泰安數(shù)據(jù)庫獲取;債券違約情況從銳思數(shù)據(jù)庫獲取;宏觀指標數(shù)據(jù)則從《中國城市統(tǒng)計年鑒》《中國城市建設(shè)統(tǒng)計年鑒》及各地級市統(tǒng)計公報獲取。

    需要特別指出的是,本文在發(fā)債主體內(nèi)部非財務(wù)指標中補充考慮了債券發(fā)行人的風險信息披露情況,并采用管理層凈語調(diào)對其進行衡量。管理層凈語調(diào)無法從國泰安數(shù)據(jù)庫中直接獲取,需要依據(jù)正面詞匯數(shù)量和負面詞匯數(shù)量進行進一步計算,具體計算如式(16)所示(謝德仁和林樂,2015)。

    其中,pos代表管理層信息披露中的正面詞匯數(shù)量,neg代表管理層信息披露中的負面詞匯數(shù)量,tone則代表管理層凈語調(diào),tone越大,證明管理層信息披露的語調(diào)越積極。

    (二)樣本數(shù)據(jù)處理

    1.指標數(shù)據(jù)預處理

    依據(jù)前文所述,對于獲取的原始指標數(shù)據(jù),首先利用已有數(shù)據(jù)中的最差值對缺失值進行填補(其原因是本文通過預測債券違約來對債券進行信用評級,信息缺失時用最差值填補缺失值,即對信息披露不完全的債券,將缺失的信息視為該方面表現(xiàn)差,這樣處理的目的是為了更好地防控債券違約風險),進一步地,對無缺失值的完整數(shù)據(jù)進行前后1%縮尾的異常值處理,并將所有指標數(shù)據(jù)歸于[0,1]之間。

    2.樣本劃分和基于SMOTE方法的非均衡樣本處理

    對于全部債券樣本,本文首先按照8:2的比例進行分層抽樣,將樣本隨機劃分為訓練集和測試集。劃分后的訓練集樣本中,違約樣本數(shù)為200,非違約樣本數(shù)為30386,顯然本文選用的樣本數(shù)據(jù)為類別非均衡數(shù)據(jù)。根據(jù)已有研究可知,對類別非均衡樣本,模型對小類別的特征學習相對不足,因此通常需要對非均衡樣本進行處理。

    本文對訓練集樣本數(shù)據(jù)按照前文所述SMOTE方法步驟進行均衡化處理,來保證訓練集中兩類樣本數(shù)量相同,降低樣本類別不均衡對預測結(jié)果產(chǎn)生的影響。通過利用SMOTE方法處理上述非均衡樣本數(shù)據(jù),得到違約樣本數(shù)與非違約樣本數(shù)達到1:1的均衡樣本,SMOTE非均衡樣本處理前后對比如表1所示。

    為了提高模型的泛化性能,對于均衡化的訓練集樣本,本文利用十折交叉驗證方法將數(shù)據(jù)按照9:1的比例劃分為新的訓練集和驗證集,參與后續(xù)模型的訓練,占樣本總數(shù)20%的測試集樣本則用于檢驗訓練完成的模型精度。

    (三)指標體系構(gòu)建

    1.初始指標體系構(gòu)建

    本研究所涉及的初始指標共377個,其中包含“資產(chǎn)負債率”“長期借款與總資產(chǎn)比”等157個發(fā)債主體內(nèi)部財務(wù)指標,“股權(quán)制衡度”“管理層凈語調(diào)”等125個發(fā)債主體內(nèi)部非財務(wù)指標,“目前狀態(tài)”“品種類別”等50個債券基本信息指標以及“城市產(chǎn)業(yè)結(jié)構(gòu)高級化”“城市人均GDP”等45個宏觀指標。

    2.第一次指標篩選

    對于指標海選得到的初始指標體系,刪掉無法獲取數(shù)據(jù)或數(shù)據(jù)存在大量缺失的指標。刪除數(shù)據(jù)獲取困難的指標、保留數(shù)據(jù)可獲取的指標可以在最大程度上保證構(gòu)建的指標體系信息的有效性。

    經(jīng)過篩選后保留了“資產(chǎn)負債率”“資產(chǎn)報酬率”“城市數(shù)字普惠金融指數(shù)”等200個指標,第一次篩選后的指標體系如表2所示。

    3.宏觀指標差異化處理

    計算第一次指標篩選所得指標體系的宏觀指標與債券其他特征指標的克羅內(nèi)克積,并將得到的相互作用項作為新的宏觀指標納入后續(xù)模型訓練及最優(yōu)債券信用評級指標體系的構(gòu)建中,最終得到6951個特征指標。

    (四)上市公司債券信用評級指標體系構(gòu)建與精度檢驗

    1.基于AdaBoost的第二次指標篩選

    對于第一次篩選完成的指標體系,在宏觀指標完成債券差異化處理的基礎(chǔ)上,通過基于AdaBoost的嵌入式特征選擇方法對指標體系進行第二次指標篩選,并依據(jù)模型的精度及降維效果,遴選出對債券違約狀態(tài)有顯著影響的特征指標,最終確定本文的最優(yōu)指標體系如表3所示。

    2.債券信用評級違約預測精度測算

    模型精度的測算離不開混淆矩陣中的相關(guān)指標,本文表4、表5分別給出了模型訓練集和測試集的混淆矩陣。依據(jù)表4,可以得到訓練集的鑒別精度結(jié)果,其中,acc=98.76%、AUC=99.88%、Type I Error=1.74%、TypeⅡError=0.74%、G-mean=98.76%。

    依據(jù)表5,可以得到測試集的鑒別精度結(jié)果,其中,acc=98.35%、AUC=99.07%、Type I Error=1.62%、Type Ⅱ Error=6.00%、G-mean=96.17%。可以看出,無論是訓練集還是測試集,基于AdaBoost算法的模型鑒別精度均在94%以上,表明用篩選后的指標體系建立模型具有一定的泛化性能,同時對上市公司債券有較強的違約判別能力。

    3.不同模型精度對比

    為了進一步論證基于AdaBoost算法對指標體系進行篩選的合理性及有效性,本文選取了4個常用模型與AdaBoost模型得到的結(jié)果進行精度及指標體系規(guī)模對比,4個模型分別為梯度提升樹(Gradient Boosting Decision Tree,GBDT)、邏輯回歸(Logistic Regression,LR)、隨機森林(Random Forest,RF)以及隨機梯度下降(Stochastic Gradient Descent,SGD)。

    表6給出了不同違約預測模型的指標體系規(guī)模,可以看出,AdaBoost的最優(yōu)指標體系包含45個特征,指標個數(shù)明顯少于其他模型的最優(yōu)指標體系,因此,對于本文而言,AdaBoost模型降維效果最好。

    較訓練集而言,測試集精度可以更好地衡量模型的泛化性能,因此在模型精度對比時,本文僅列出不同模型的測試集精度結(jié)果,不同模型測試集精度結(jié)果對比如表7所示。

    由表7可知,AdaBoost模型的Type Ⅱ Error指標明顯小于LR模型等其他模型。對于債券違約預警問題而言,Type Ⅱ Error是投資者更為關(guān)注的精度指標,Type Ⅱ Error越大意味著將違約樣本誤判為非違約樣本的可能性越大,帶來的損失也會越大。同時,AdaBoost模型的其他精度指標雖然不全為最好,但同其他模型中該指標最好值相差很小。因此,在Type Ⅱ Error最小的基礎(chǔ)上,綜合考慮各模型的精度指標,確定AdaBoost模型擁有最好的違約鑒別精度以及泛化能力。

    綜合考慮各模型的指標體系規(guī)模和鑒別精度,不難看出AdaBoost模型的指標體系規(guī)模最小、鑒別精度綜合表現(xiàn)最好,表明本文基于AdaBoost算法篩選得到的信用評級指標體系可信度高。

    4.引入風險信息披露情況指標前后的精度對比

    為了驗證信息充分性在債券信用評級中的重要作用,表8列出了引入風險信息披露情況指標前后測試集預測精度的結(jié)果對比。

    由表8可知,在引入風險信息披露情況指標后,模型的綜合精度和類別精度指標均得到了一定的提升,明顯優(yōu)于未考慮風險信息披露情況指標的樣本模型,其中在債券違約鑒別中重點關(guān)注的指標Type Ⅱ Error精度提升了2%,且對非違約的債券鑒別精度也提升了3.26%。顯然,加入風險信息披露情況指標來增大指標體系的信息充分性,的確可以得到更高的債券違約鑒別精度,同時降維效果也更優(yōu),即指標信息充分性提高,會提高債券違約的預測精度,因此得到的信用評級指標體系對衡量債券信用評級更可信。

    5.宏觀指標差異化處理前后的精度對比

    為了進一步驗證使用具有債券差異性的宏觀指標進行模型訓練能夠取得更好的綜合鑒別精度和降維效果,即信息差異對識別債券信用評級關(guān)鍵因素具有積極作用,本文計算了宏觀指標差異化處理前后模型的預測精度,表9列出了宏觀指標差異化處理前后測試集的預測精度對比。

    由表9可知,當宏觀指標具有債券差異性時,模型的各個精度指標得到了較顯著的提升,其中綜合精度指標acc、G-mean和非違約債券的鑒別指標Type I Error、違約債券的鑒別精度指標Type Ⅱ Error提升較為顯著,最大的提升幅度達到2.12%,且綜合精度指標AUC也有小幅提升。同時,對于使用具有債券差異性的宏觀指標樣本的模型而言,其最優(yōu)指標體系中指標個數(shù)占第一次篩選后的指標體系的0.65%,明顯優(yōu)于未使用具有債券差異性的宏觀指標樣本模型的15.50%,表明使用具有債券差異性的宏觀指標可以更好地提升指標體系預測精度和降維效率。

    綜合考慮模型的精度及降維效率,發(fā)現(xiàn)采用克羅內(nèi)克積構(gòu)造具有債券差異性的宏觀指標方法有效,得到的債券違約鑒別精度更高,即提升債券信用評級指標體系信息差異的思路可行,指標信息差異性增大,會提高債券違約的預測精度,提升債券信用評級指標體系可靠性。

    6.信用評級指標重要性排序

    前文已經(jīng)驗證本文構(gòu)建的信用評級模型具有不同信用的區(qū)分能力,且判別精度較高。在此基礎(chǔ)上,計算該模型對應的各特征指標權(quán)重,并對其重要性進行排序,最優(yōu)信用評級指標體系中各特征指標權(quán)重及重要性排序如表10所示。

    由表10可知,在包含45個指標的信用評級指標體系中,城市數(shù)字普惠金融指數(shù)*管理層平均年齡、城市產(chǎn)業(yè)結(jié)構(gòu)*監(jiān)事會規(guī)模、綜合收益增長率等指標的重要性得分較高,表明這些指標對于債券違約鑒別、信用評級而言經(jīng)濟意義更顯著。最終指標體系中包含債券發(fā)行人的風險信息披露情況指標(如管理層凈語調(diào))和多個經(jīng)過差異化處理后的宏觀指標,進一步驗證了具有信息充分性和差異性的指標對模型預測有影響。

    四、結(jié)論與政策建議

    (一)結(jié)論

    本文通過對我國上市公司發(fā)行的債券進行違約鑒別,并依據(jù)違約預測精度結(jié)果反推最優(yōu)指標體系,得到主要結(jié)論如下:

    一是本文構(gòu)建的最終指標體系包含管理層風險披露信息的管理層凈語調(diào)等指標,表明上市公司管理層信息披露對債券信用評級有影響,引入該類指標,可提高債券信用評級指標體系涵蓋信息的充分性。通過對比引入風險信息披露情況指標前后的預測結(jié)果,發(fā)現(xiàn)指標體系信息充分性的提高能提升預測精度,進而增強債券信用評級指標體系的可靠性。

    二是本文利用基于AdaBoost的嵌入式特征選擇方法對指標體系進行降維處理,在指標篩選時能更好地保證指標體系信息充分性。該方法在模型訓練的同時完成指標篩選,相比在預測前進行指標篩選的處理方法,本文的方法在預測時輸入更多特征,因此涵蓋的信息更充分。且根據(jù)預測精度結(jié)果可知,本文模型的預測精度AUC=99.07%gt;70.00%,且降維后只包含第一次篩選后指標數(shù)的0.65%(即從原有的6951個指標降低到45個指標),表明了本文篩選方法的可靠性和有效性。

    三是根據(jù)最終確定的最優(yōu)指標體系可知,基于克羅內(nèi)克積構(gòu)建的多個差異性指標對債券信用評級有影響,表明信息差異對識別債券信用評級關(guān)鍵因素有積極作用。

    四是根據(jù)宏觀指標的差異性處理前后預測結(jié)果對比,引入進行差異性處理的宏觀指標構(gòu)建的指標體系預測精度更高,表明提高指標體系信息差異性能改善預測精度,提高債券信用評級指標體系可靠性。

    (二)政策建議

    對于發(fā)債主體而言,應充分了解可能影響信用水平的各因素,加強企業(yè)自身財務(wù)建設(shè)、非財務(wù)建設(shè),進一步提高企業(yè)及其發(fā)行債券的透明度和風險信息披露。同時在企業(yè)內(nèi)部建立完善的債券管理和監(jiān)控機制,依據(jù)信用評級影響因素做好自身違約風險評估與防范,提升其發(fā)行債券的信用水平。

    對于相關(guān)機構(gòu)而言,應及時更新信用評級指標體系,根據(jù)企業(yè)的發(fā)展情況不斷進行自我革新,同時要清晰地披露信用評級標準和方法。此外,相關(guān)機構(gòu)要加強對企業(yè)綜合收益增長率等指標的關(guān)注,幫助企業(yè)有的放矢地提升自身及發(fā)行債券的信用水平,為企業(yè)健康良性發(fā)展營造一個好的信用環(huán)境。

    對于投資者而言,要對債券信用評級的影響因素有一定的了解,不能盲目投資。不僅要關(guān)注企業(yè)公布的重要財務(wù)指標,還應該關(guān)注企業(yè)的風險披露信息,同時,要考慮發(fā)債主體所在城市的宏觀指標和微觀指標的綜合狀況,對債券信用水平關(guān)鍵影響因素表現(xiàn)不佳的債券要謹慎投資,要及時關(guān)注相關(guān)機構(gòu)公布的債券信用影響因素變化情況,盡可能地提前規(guī)避風險,減少自身損失。

    參考文獻:

    [1]王培輝,張猛.中國債券市場信息溢出效應研究[J].海南金融,2023(9):3-20.

    [2]鄒昆侖,陸萍.基于Z模型對中國企業(yè)債信用評級的實證研究[J].宏觀經(jīng)濟研究,2018(5):72-83+175.

    [3]Chi M,Hongyan S,Shaofan W,et al.Bond default prediction based on deep learning and knowledge graph technology[J].IEEE Access,2021(9):12750-12761.

    [4]姚瀟,李可,余樂安.非平衡樣本下基于生成對抗網(wǎng)絡(luò)過抽樣技術(shù)的公司債券違約風險預測研究[J].系統(tǒng)工程理論與實踐,2022(10):2617-2634.

    [5]生柳榮,陳海華,胡施聰,等.企業(yè)債券信用風險預警模型及其運用[J].投資研究,2019(6):25-35.

    [6]蔣敏,周煒,史濟川,等.基于fsQCA的上市企業(yè)債券違約影響因素研究[J].管理學報,2021(7):1076-1085.

    [7]肖艷麗,向有濤.企業(yè)債券違約風險預警——基于GWO-XGBoost方法[J].上海金融,2021(10):44-54.

    [8]陳學彬,武靖,徐明東.我國信用債個體違約風險測度與防范——基于LSTM深度學習模型[J].復旦學報(社會科學版),2021(3):159-173.

    [9]Lu Z,Zhuo Z.Modelling of Chinese corporate bond default–A machine learning approach[J].Accounting amp; Finance,2021(5):6147-6191.

    [10]吳育輝,唐浩博.風險信息披露能預測債券違約風險嗎?——來自債券募集說明書的經(jīng)驗證據(jù)[J].財務(wù)研究,2021(5):15-28.

    [11]吳建華,張穎,王新軍.信息披露扭曲下企業(yè)債券違約風險量化研究[J].數(shù)理統(tǒng)計與管理,2017(1):175-190.

    [12]高秀存,趙明航,楊倩芳.基于因子分析和SVM的信用債違約風險預警實證研究——以河北省為例[J].征信,2022(8):87-92.

    [13]王玉龍,周榴,張滌霏.企業(yè)債務(wù)違約風險預測——基于機器學習的視角[J].財政科學,2022(6):62-74.

    [14]龐春潮,羅苑瑋.基于投資者視角的信用類債券違約風險預警[J].統(tǒng)計與決策,2023(2):152-157.

    [15]謝玲玲,范龍振.上市公司信用債違約風險評估[J].管理現(xiàn)代化,2020(2):104-108.

    [16]徐舒玥,曹艷華.基于可解釋機器學習的信用債違約研究[J].科學決策,2023(5):190-200.

    [17]白鈺銘,姜昱汐.基于特征選擇和違約鑒別的我國上市公司債券違約預警模型研究[J].工程數(shù)學學報,2024(1):67-87.

    [18]遲國泰,張亞京,石寶峰.基于Probit回歸的小企業(yè)債信評級模型及實證[J].管理科學學報,2016(6):136-156.

    [19]Leippold M,Wang Q,Zhou W.Machine learning in the Chinese stock market[J].Journal of Financial Economics,2022(2):64-82.

    [20]Freund Y,Schapire R E.A Decision-Theoretic Generalization of On-Line Learning and an Application to Boosting[J].Journal of Computer and System Sciences,1997(1):119-139.

    [21]遲國泰,章彤,張志鵬.基于非平衡數(shù)據(jù)處理的上市公司ST預警混合模型[J].管理評論,2020(3):3-20.

    [22]謝德仁,林樂.管理層語調(diào)能預示公司未來業(yè)績嗎?——基于我國上市公司年度業(yè)績說明會的文本分析[J].會計研究,2015(2):20-27+93.

    精品国产乱码久久久久久小说| 啦啦啦在线免费观看视频4| 激情视频va一区二区三区| 亚洲欧美清纯卡通| 777久久人妻少妇嫩草av网站| 如何舔出高潮| av线在线观看网站| 熟妇人妻不卡中文字幕| 99久久人妻综合| 国产色婷婷99| 蜜桃在线观看..| 午夜福利,免费看| 中文字幕人妻丝袜一区二区 | 涩涩av久久男人的天堂| 国产淫语在线视频| 人人妻人人爽人人添夜夜欢视频| 老司机影院毛片| 亚洲一区二区三区欧美精品| 成人国产麻豆网| 亚洲内射少妇av| 高清欧美精品videossex| 97在线人人人人妻| 国产熟女欧美一区二区| 又大又黄又爽视频免费| 久久久久久久精品精品| 国产无遮挡羞羞视频在线观看| 久久毛片免费看一区二区三区| 国产成人精品久久二区二区91 | 国产精品 欧美亚洲| 老司机亚洲免费影院| 乱人伦中国视频| 久久久久久久久久久免费av| 精品一区二区免费观看| 亚洲精品国产av成人精品| 日日啪夜夜爽| 制服丝袜香蕉在线| 97在线视频观看| 人妻少妇偷人精品九色| 1024香蕉在线观看| 国产熟女午夜一区二区三区| 热re99久久国产66热| 999精品在线视频| 人妻一区二区av| 只有这里有精品99| 免费观看无遮挡的男女| 高清在线视频一区二区三区| 亚洲成人av在线免费| 欧美日韩一级在线毛片| 亚洲婷婷狠狠爱综合网| 精品亚洲成a人片在线观看| videosex国产| 中国三级夫妇交换| 成年人午夜在线观看视频| 国产xxxxx性猛交| 大陆偷拍与自拍| 国产极品粉嫩免费观看在线| 少妇猛男粗大的猛烈进出视频| 久久精品久久久久久久性| 久久av网站| 国产黄频视频在线观看| 国产成人a∨麻豆精品| 久久人人97超碰香蕉20202| 亚洲成av片中文字幕在线观看 | av国产精品久久久久影院| 国产无遮挡羞羞视频在线观看| 欧美+日韩+精品| 欧美精品一区二区免费开放| 国产成人a∨麻豆精品| 国产深夜福利视频在线观看| 永久网站在线| 国产激情久久老熟女| 亚洲国产色片| 考比视频在线观看| freevideosex欧美| 日日摸夜夜添夜夜爱| 天天躁夜夜躁狠狠躁躁| 亚洲av日韩在线播放| 777米奇影视久久| 成年人午夜在线观看视频| 精品一品国产午夜福利视频| 街头女战士在线观看网站| 丰满迷人的少妇在线观看| 久久精品久久久久久久性| 十八禁网站网址无遮挡| 一区在线观看完整版| av在线播放精品| 999久久久国产精品视频| 国产精品.久久久| 一区二区av电影网| 亚洲国产最新在线播放| 久热这里只有精品99| 日韩一本色道免费dvd| 欧美少妇被猛烈插入视频| 人妻系列 视频| 18禁国产床啪视频网站| 成人手机av| 人人妻人人添人人爽欧美一区卜| 精品人妻一区二区三区麻豆| 日韩人妻精品一区2区三区| 黑人猛操日本美女一级片| 国产乱来视频区| 国产成人一区二区在线| av福利片在线| 在线 av 中文字幕| 久久亚洲国产成人精品v| 欧美日韩视频精品一区| 天天躁日日躁夜夜躁夜夜| 少妇精品久久久久久久| freevideosex欧美| av国产久精品久网站免费入址| 一个人免费看片子| 亚洲av欧美aⅴ国产| 亚洲精品国产色婷婷电影| 激情视频va一区二区三区| 成年人免费黄色播放视频| 视频区图区小说| 日韩欧美精品免费久久| 国产精品国产av在线观看| 一区二区三区精品91| 最近中文字幕高清免费大全6| 18禁观看日本| 满18在线观看网站| xxxhd国产人妻xxx| 午夜影院在线不卡| 日本av免费视频播放| 国产av国产精品国产| 亚洲伊人色综图| 肉色欧美久久久久久久蜜桃| 男女高潮啪啪啪动态图| 久久久久久久国产电影| 最近最新中文字幕大全免费视频 | 91成人精品电影| 精品国产乱码久久久久久男人| 亚洲国产欧美在线一区| 晚上一个人看的免费电影| 国产男女超爽视频在线观看| 午夜久久久在线观看| 麻豆精品久久久久久蜜桃| 不卡视频在线观看欧美| 亚洲av.av天堂| 色播在线永久视频| 乱人伦中国视频| 欧美xxⅹ黑人| 啦啦啦啦在线视频资源| 午夜免费观看性视频| 亚洲久久久国产精品| av在线观看视频网站免费| av一本久久久久| 97精品久久久久久久久久精品| 天天影视国产精品| 精品国产露脸久久av麻豆| 一区在线观看完整版| 国产人伦9x9x在线观看 | 人妻一区二区av| 国产成人精品久久二区二区91 | 国产黄色免费在线视频| 各种免费的搞黄视频| 国精品久久久久久国模美| 天堂俺去俺来也www色官网| 涩涩av久久男人的天堂| 国产一区二区三区综合在线观看| 美女大奶头黄色视频| 少妇的丰满在线观看| 午夜福利视频精品| 在线观看www视频免费| 久久婷婷青草| 亚洲欧美一区二区三区黑人 | 97在线人人人人妻| 成年人免费黄色播放视频| 麻豆av在线久日| 自拍欧美九色日韩亚洲蝌蚪91| 国产免费视频播放在线视频| 一级a爱视频在线免费观看| 波多野结衣av一区二区av| 伦理电影大哥的女人| 18在线观看网站| 国产亚洲午夜精品一区二区久久| 日韩免费高清中文字幕av| 狠狠精品人妻久久久久久综合| 亚洲精华国产精华液的使用体验| 亚洲第一区二区三区不卡| 国产成人免费观看mmmm| 毛片一级片免费看久久久久| 国产精品久久久久久久久免| 亚洲,欧美精品.| 久久久久久久久久久久大奶| 有码 亚洲区| 97精品久久久久久久久久精品| 18+在线观看网站| 亚洲欧美一区二区三区黑人 | 国产精品香港三级国产av潘金莲 | 久久久久久久国产电影| 国产精品二区激情视频| 91国产中文字幕| 男女国产视频网站| 亚洲精品日本国产第一区| 亚洲av在线观看美女高潮| 免费大片黄手机在线观看| 这个男人来自地球电影免费观看 | 国产人伦9x9x在线观看 | 欧美日韩成人在线一区二区| 中文字幕最新亚洲高清| 一级黄片播放器| 18禁动态无遮挡网站| 一区福利在线观看| 肉色欧美久久久久久久蜜桃| 美女国产视频在线观看| 热re99久久国产66热| 18禁观看日本| 深夜精品福利| 亚洲av欧美aⅴ国产| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲四区av| 国产成人精品一,二区| 国产精品av久久久久免费| 国产欧美日韩综合在线一区二区| 两个人看的免费小视频| 亚洲五月色婷婷综合| 成人午夜精彩视频在线观看| 国产男女超爽视频在线观看| 日本色播在线视频| 两个人看的免费小视频| 97精品久久久久久久久久精品| 精品一区二区三卡| 叶爱在线成人免费视频播放| 国产一区二区三区av在线| 青青草视频在线视频观看| 亚洲欧美日韩另类电影网站| 美女中出高潮动态图| 亚洲国产精品一区二区三区在线| 亚洲av成人精品一二三区| 啦啦啦在线免费观看视频4| 女人被躁到高潮嗷嗷叫费观| 午夜福利,免费看| 亚洲一区二区三区欧美精品| 国产亚洲午夜精品一区二区久久| 国产老妇伦熟女老妇高清| 精品一区二区三卡| 国产成人精品无人区| 免费久久久久久久精品成人欧美视频| 亚洲欧美清纯卡通| 欧美精品高潮呻吟av久久| 国产高清不卡午夜福利| 久久久精品国产亚洲av高清涩受| av国产精品久久久久影院| 欧美日韩亚洲高清精品| 亚洲精品美女久久av网站| 哪个播放器可以免费观看大片| 91在线精品国自产拍蜜月| 如何舔出高潮| 久久久国产精品麻豆| 婷婷色av中文字幕| 两个人免费观看高清视频| 一级,二级,三级黄色视频| 国产精品欧美亚洲77777| 亚洲国产色片| 日本欧美视频一区| 十八禁高潮呻吟视频| 一个人免费看片子| 永久网站在线| 亚洲国产欧美在线一区| 欧美 亚洲 国产 日韩一| 天天影视国产精品| 亚洲av欧美aⅴ国产| 国产男女内射视频| 美女国产高潮福利片在线看| 成年女人在线观看亚洲视频| 午夜福利网站1000一区二区三区| 免费不卡的大黄色大毛片视频在线观看| 午夜福利视频精品| 高清在线视频一区二区三区| 久久国产精品大桥未久av| 精品国产一区二区久久| 99久久综合免费| 国产成人精品一,二区| 欧美精品人与动牲交sv欧美| 国产亚洲午夜精品一区二区久久| 久久精品夜色国产| 国产又色又爽无遮挡免| 亚洲人成电影观看| 观看美女的网站| 久久人人97超碰香蕉20202| 精品国产露脸久久av麻豆| 可以免费在线观看a视频的电影网站 | 亚洲男人天堂网一区| 免费观看a级毛片全部| 热re99久久国产66热| 一级爰片在线观看| 麻豆av在线久日| 综合色丁香网| 久久人人97超碰香蕉20202| 国产av精品麻豆| 亚洲国产欧美在线一区| 一二三四在线观看免费中文在| 熟女av电影| 亚洲欧洲国产日韩| 一区二区三区乱码不卡18| 免费黄频网站在线观看国产| 制服丝袜香蕉在线| 日韩av免费高清视频| 国产乱来视频区| 成年人免费黄色播放视频| 国产97色在线日韩免费| 热99国产精品久久久久久7| 80岁老熟妇乱子伦牲交| 中文欧美无线码| 日日爽夜夜爽网站| 精品亚洲成国产av| 欧美xxⅹ黑人| www.精华液| 国产精品国产三级国产专区5o| 久久久欧美国产精品| 男人操女人黄网站| 亚洲国产欧美网| 欧美av亚洲av综合av国产av | 日韩人妻精品一区2区三区| 黄片播放在线免费| 精品国产超薄肉色丝袜足j| 国产成人精品在线电影| 国产男人的电影天堂91| 国产精品久久久久久精品电影小说| kizo精华| 中文字幕人妻熟女乱码| 亚洲综合精品二区| 人妻 亚洲 视频| 在线观看免费日韩欧美大片| 免费观看无遮挡的男女| av在线老鸭窝| 久久99蜜桃精品久久| 男男h啪啪无遮挡| 亚洲国产成人一精品久久久| 熟女电影av网| 国产亚洲欧美精品永久| 亚洲精品在线美女| 久久女婷五月综合色啪小说| 亚洲国产欧美在线一区| 亚洲国产欧美日韩在线播放| 国产免费福利视频在线观看| 黑人欧美特级aaaaaa片| 老女人水多毛片| 亚洲人成电影观看| av电影中文网址| 午夜福利一区二区在线看| 日韩制服骚丝袜av| 国产日韩欧美视频二区| 国产高清不卡午夜福利| 国产人伦9x9x在线观看 | 日本-黄色视频高清免费观看| 国产成人av激情在线播放| 欧美97在线视频| 香蕉国产在线看| www日本在线高清视频| 亚洲男人天堂网一区| 日韩大片免费观看网站| 天堂俺去俺来也www色官网| 免费观看在线日韩| 男女高潮啪啪啪动态图| 自线自在国产av| 久久久亚洲精品成人影院| 久久久久久久久免费视频了| 97精品久久久久久久久久精品| 18禁国产床啪视频网站| 亚洲精品,欧美精品| 亚洲欧洲日产国产| 麻豆精品久久久久久蜜桃| 国产精品一区二区在线观看99| 欧美日韩亚洲国产一区二区在线观看 | 久久97久久精品| 秋霞伦理黄片| 亚洲欧洲国产日韩| 久久久久国产精品人妻一区二区| 看免费成人av毛片| 国产精品 欧美亚洲| 精品人妻一区二区三区麻豆| 老鸭窝网址在线观看| 在线观看美女被高潮喷水网站| 日本91视频免费播放| 国产精品 国内视频| 欧美日韩国产mv在线观看视频| 亚洲人成网站在线观看播放| 9色porny在线观看| 久久久国产欧美日韩av| 久久99精品国语久久久| 春色校园在线视频观看| 国产欧美亚洲国产| 一本色道久久久久久精品综合| 黄色 视频免费看| 精品人妻偷拍中文字幕| 日日爽夜夜爽网站| 久久这里有精品视频免费| 一区在线观看完整版| 免费日韩欧美在线观看| 国产毛片在线视频| 一本色道久久久久久精品综合| 少妇熟女欧美另类| 精品酒店卫生间| 男女高潮啪啪啪动态图| 久久久久久久国产电影| 纵有疾风起免费观看全集完整版| 制服人妻中文乱码| 精品99又大又爽又粗少妇毛片| 成人毛片60女人毛片免费| 黄色怎么调成土黄色| av在线老鸭窝| 欧美人与性动交α欧美精品济南到 | 一区二区三区四区激情视频| 久久久久久久亚洲中文字幕| 男人操女人黄网站| 久久久久网色| 波野结衣二区三区在线| av一本久久久久| www.自偷自拍.com| 交换朋友夫妻互换小说| 一二三四中文在线观看免费高清| 午夜影院在线不卡| av福利片在线| 亚洲精品第二区| 99国产综合亚洲精品| 丰满迷人的少妇在线观看| 少妇精品久久久久久久| 成人国产麻豆网| 日韩一区二区三区影片| 大话2 男鬼变身卡| 国产av国产精品国产| 精品亚洲成国产av| av免费观看日本| a级毛片黄视频| 老司机影院毛片| 亚洲情色 制服丝袜| 精品国产一区二区三区四区第35| 少妇被粗大猛烈的视频| 人妻少妇偷人精品九色| www.av在线官网国产| 国产一区亚洲一区在线观看| 麻豆精品久久久久久蜜桃| 国产成人精品在线电影| 天美传媒精品一区二区| 丝袜喷水一区| 少妇人妻 视频| 亚洲精品第二区| 精品国产露脸久久av麻豆| 青青草视频在线视频观看| 久久国产精品男人的天堂亚洲| 视频在线观看一区二区三区| 日日摸夜夜添夜夜爱| 国产午夜精品一二区理论片| 亚洲国产av影院在线观看| 97精品久久久久久久久久精品| 国产高清不卡午夜福利| 我的亚洲天堂| 精品一品国产午夜福利视频| 亚洲色图 男人天堂 中文字幕| 老司机亚洲免费影院| 久久久久国产网址| √禁漫天堂资源中文www| 国产深夜福利视频在线观看| 免费高清在线观看视频在线观看| 一区二区三区乱码不卡18| 性高湖久久久久久久久免费观看| 成人影院久久| 伦精品一区二区三区| 一区二区日韩欧美中文字幕| 久久精品国产自在天天线| 狂野欧美激情性bbbbbb| 大陆偷拍与自拍| 久久综合国产亚洲精品| 国产精品香港三级国产av潘金莲 | 久久精品国产鲁丝片午夜精品| 视频在线观看一区二区三区| 久久午夜福利片| 不卡视频在线观看欧美| 精品久久蜜臀av无| 欧美成人精品欧美一级黄| 最新的欧美精品一区二区| 欧美 亚洲 国产 日韩一| 国产在线免费精品| 亚洲欧洲精品一区二区精品久久久 | 国产成人精品婷婷| 国产精品久久久久久久久免| 午夜日本视频在线| 亚洲av免费高清在线观看| 久久国产精品大桥未久av| 中文字幕av电影在线播放| 十八禁网站网址无遮挡| www.精华液| 一级爰片在线观看| 搡女人真爽免费视频火全软件| 亚洲国产成人一精品久久久| 两性夫妻黄色片| 免费黄色在线免费观看| 人成视频在线观看免费观看| 午夜福利影视在线免费观看| 久久久久久久久久久免费av| 日韩成人av中文字幕在线观看| 亚洲人成网站在线观看播放| 777米奇影视久久| 国产野战对白在线观看| 日韩不卡一区二区三区视频在线| 亚洲国产欧美日韩在线播放| 国产极品天堂在线| 国产 精品1| 寂寞人妻少妇视频99o| 国产免费视频播放在线视频| 欧美人与善性xxx| 久久综合国产亚洲精品| 亚洲视频免费观看视频| 一边亲一边摸免费视频| 久久精品熟女亚洲av麻豆精品| 女人精品久久久久毛片| 久久久久久人妻| 欧美日韩一级在线毛片| 婷婷色av中文字幕| 欧美亚洲 丝袜 人妻 在线| 在线免费观看不下载黄p国产| 男的添女的下面高潮视频| 亚洲四区av| 夜夜骑夜夜射夜夜干| 欧美激情高清一区二区三区 | 97在线视频观看| 成人午夜精彩视频在线观看| 中文字幕色久视频| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲av成人精品一二三区| 亚洲精品国产av蜜桃| 国产国语露脸激情在线看| 晚上一个人看的免费电影| 欧美日韩精品成人综合77777| 中文字幕制服av| 一区二区av电影网| 久久久久久人妻| 亚洲欧美成人综合另类久久久| 国产精品蜜桃在线观看| 亚洲天堂av无毛| 亚洲精品国产av蜜桃| 99香蕉大伊视频| 免费看av在线观看网站| 国产精品偷伦视频观看了| 丰满饥渴人妻一区二区三| 免费高清在线观看视频在线观看| 亚洲三区欧美一区| 欧美日韩视频精品一区| 久热久热在线精品观看| 青草久久国产| 国产一级毛片在线| 国产熟女午夜一区二区三区| 女的被弄到高潮叫床怎么办| 又大又黄又爽视频免费| 三上悠亚av全集在线观看| 亚洲美女视频黄频| 女的被弄到高潮叫床怎么办| 日本午夜av视频| 成人二区视频| 欧美日韩一级在线毛片| 在线亚洲精品国产二区图片欧美| 交换朋友夫妻互换小说| 国产免费一区二区三区四区乱码| 色94色欧美一区二区| 中文字幕最新亚洲高清| 久久久a久久爽久久v久久| 国产麻豆69| 涩涩av久久男人的天堂| 色吧在线观看| 国产成人精品一,二区| 国产成人精品婷婷| 国产av一区二区精品久久| 亚洲欧洲精品一区二区精品久久久 | 9色porny在线观看| 91精品三级在线观看| 啦啦啦在线免费观看视频4| 亚洲欧美成人精品一区二区| av天堂久久9| 久久鲁丝午夜福利片| 亚洲,一卡二卡三卡| 欧美日本中文国产一区发布| 久久99热这里只频精品6学生| 久久鲁丝午夜福利片| 丝袜美足系列| 国产成人精品福利久久| 免费观看无遮挡的男女| 九草在线视频观看| 国产精品久久久久久av不卡| 国产av精品麻豆| 久久午夜福利片| 国产av精品麻豆| 国产在视频线精品| 亚洲精品aⅴ在线观看| 久久人人爽人人片av| 国产黄色免费在线视频| 久久久久久免费高清国产稀缺| 欧美日韩视频精品一区| 另类精品久久| 欧美bdsm另类| 自线自在国产av| 色婷婷av一区二区三区视频| 亚洲国产毛片av蜜桃av| 色吧在线观看| 久久国产精品男人的天堂亚洲| 亚洲内射少妇av| 两个人免费观看高清视频| 国产成人av激情在线播放| 黄色配什么色好看| 大片免费播放器 马上看| 人体艺术视频欧美日本| 一区二区三区乱码不卡18| 国产免费福利视频在线观看| 男人爽女人下面视频在线观看| 日韩一本色道免费dvd| 麻豆乱淫一区二区| 人人妻人人澡人人看| av在线app专区| 精品久久蜜臀av无| 秋霞伦理黄片| 久久久久久久久免费视频了| 精品酒店卫生间| 9191精品国产免费久久|