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

    小微企業(yè)信用評級模型及比較研究

    2017-01-20 11:16:20肖斌卿楊旸余哲沈才勝
    系統(tǒng)工程學(xué)報 2016年6期
    關(guān)鍵詞:識別率評級小微

    肖斌卿楊 旸余 哲沈才勝

    (1.南京大學(xué)工程管理學(xué)院,江蘇南京210093;2.南京大學(xué)商學(xué)院,江蘇南京210093; 3.鄭州商品交易所,河南鄭州450008;4.紫金農(nóng)商銀行,江蘇南京210019)

    小微企業(yè)信用評級模型及比較研究

    肖斌卿1,楊 旸2*,余 哲3,沈才勝1,4

    (1.南京大學(xué)工程管理學(xué)院,江蘇南京210093;2.南京大學(xué)商學(xué)院,江蘇南京210093; 3.鄭州商品交易所,河南鄭州450008;4.紫金農(nóng)商銀行,江蘇南京210019)

    在調(diào)查和文獻(xiàn)基礎(chǔ)上,進(jìn)行信用風(fēng)險來源識別、評級指標(biāo)分類和評級方法選擇,構(gòu)建商業(yè)銀行內(nèi)部信用評級模型,以期在授信審批環(huán)節(jié)提高信用風(fēng)險管理水平.基于某商業(yè)銀行2008–2013年小微企業(yè)實際信貸數(shù)據(jù),運(yùn)用線性判別分析、二項邏輯回歸和十種基于不同學(xué)習(xí)算法的BP神經(jīng)網(wǎng)絡(luò)模型構(gòu)建內(nèi)部信用評級模型,并在評級指標(biāo)體系中加入宏觀經(jīng)濟(jì)變量,使度量風(fēng)險的穩(wěn)健性進(jìn)一步得到提升.最后通過四種方法對不同模型的結(jié)果和評級有效性進(jìn)行了對比分析,認(rèn)為基于Levenbery-Marquardt學(xué)習(xí)算法的NN10模型具有最優(yōu)的評級有效性.

    信用評級;小微企業(yè);模型比較;學(xué)習(xí)算法

    1 引 言

    我國當(dāng)前正處在經(jīng)濟(jì)發(fā)展“提質(zhì)換擋”的“新常態(tài)”發(fā)展階段,經(jīng)濟(jì)結(jié)構(gòu)和產(chǎn)業(yè)結(jié)構(gòu)正在進(jìn)行深度調(diào)整和優(yōu)化.在“新常態(tài)”的經(jīng)濟(jì)氛圍下,新興科技、新興產(chǎn)業(yè)和新興市場的發(fā)展和壯大,為小微企業(yè)的發(fā)展和壯大提供了難得的歷史機(jī)遇.也正是在這一階段,經(jīng)濟(jì)內(nèi)生增長要求更多的人力資本和技術(shù)創(chuàng)新參與社會經(jīng)濟(jì)分工,并最終轉(zhuǎn)化為經(jīng)濟(jì)增長產(chǎn)出.覆蓋廣、種類多的小微企業(yè)在國家經(jīng)濟(jì)結(jié)構(gòu)調(diào)整和升級中起到的促進(jìn)作用,在今后將會更加凸顯.

    小微企業(yè)的經(jīng)營發(fā)展和技術(shù)創(chuàng)新離不開資金的支持,特別是在當(dāng)前信貸資本緊縮投放的時期,小微企業(yè)融資難的問題日益凸顯.Schreiner[1]在研究中認(rèn)為,企業(yè)處在初創(chuàng)和成長階段更加依賴從傳統(tǒng)金融機(jī)構(gòu)那里獲取資金支持.然而,小微企業(yè)在通過銀行進(jìn)行信貸融資的過程中,會遇到諸多的問題.一方面,小微企業(yè)受其經(jīng)營時間、經(jīng)營業(yè)績、公司治理等因素的制約,普遍不能提供類似大型企業(yè)財務(wù)報表的詳細(xì)業(yè)務(wù)信息、財務(wù)信息以及其他融資擔(dān)保證明;另一方面,商業(yè)銀行經(jīng)營強(qiáng)調(diào)貸款的風(fēng)險控制,在開展小微企業(yè)貸款業(yè)務(wù)時,銀行對小微企業(yè)的資信審查和授信工作正逐步趨于嚴(yán)格和規(guī)范,與此相悖的是上市銀行內(nèi)部控制缺陷認(rèn)定存在查找范圍不統(tǒng)一、認(rèn)定標(biāo)準(zhǔn)不完善、認(rèn)定數(shù)量和嚴(yán)重程度明顯不足等問題[2].上述這種銀企之間存在的信息不對稱問題[3],是制約小微企業(yè)從商業(yè)銀行獲取信貸融資的主要因素.

    針對小微企業(yè)當(dāng)前面臨的融資困境,商業(yè)銀行所采用的信用評級卻暴露出諸多問題,例如專家判斷法比重大于模型法比重、評級方法與實際需求存在出入以及小微企業(yè)信用評級體系缺乏針對性.隨著貸款業(yè)務(wù)的不斷精細(xì)化,各家商業(yè)銀行均推出各自的小微企業(yè)貸款方案,更加注重開展專業(yè)化的小微企業(yè)貸前審查工作,同時也對信用評級方法以及評級方法的有效性提出了更高要求.在理論界,學(xué)者們依實務(wù)操作的需要,針對小微企業(yè)貸款業(yè)務(wù)的信貸模式、評級指標(biāo)、評級方法以及違約概率估計開展了廣泛的研究和討論,匯總出豐富的研究成果.其中,對評級指標(biāo)和評級方法的討論尤為突出.

    在研究指標(biāo)方面,郭小波等[4]在小微企業(yè)信用風(fēng)險的識別因子研究中,引入企業(yè)財務(wù)指標(biāo)、企業(yè)定性指標(biāo)以及與企業(yè)主有關(guān)的指標(biāo)作為研究變量.謝平等[5]認(rèn)為信用評級需考慮的因素包括借款企業(yè)經(jīng)營環(huán)境、所有制與經(jīng)營權(quán)、管理水平、營運(yùn)價值、盈利能力、風(fēng)險程度等因素.徐超等[6]基于多智能體仿真方法,從中小企業(yè)組成的聯(lián)保貸款入手,重點(diǎn)考察了聯(lián)保的組織規(guī)模和組織成員初始信用水平對信用風(fēng)險演化的影響.陳其安等[7]認(rèn)為在一定條件下,投融資平臺和銀行之間的均衡貸款數(shù)量將隨著政府信用擔(dān)保履約率增加而增加.同時,李毅等[8]也認(rèn)為小微企業(yè)在融資過程中若能夠擁有信用擔(dān)保,將有益于其獲取貸款. Hajek等[9]對企業(yè)信用評級預(yù)測中的財務(wù)狀況、經(jīng)營狀況等建模指標(biāo)進(jìn)行了研究.Doumpos等[10]在其研究中構(gòu)建了一個包含財務(wù)數(shù)據(jù)和結(jié)構(gòu)化模型的信用評級預(yù)測系統(tǒng).同時,現(xiàn)有文獻(xiàn)亦強(qiáng)調(diào)加入定性指標(biāo)對研究小微企業(yè)信貸業(yè)務(wù)的重要性.Vassiliou[11]在對印度小微企業(yè)貸款案例分析的基礎(chǔ)上,認(rèn)為信用風(fēng)險影響因素包括貸款人經(jīng)營理念、貸款人經(jīng)營水平、貸款人有無違法記錄、貸款人經(jīng)營思路、貸款利率、用途等; Malhotra等[12]通過對首發(fā)小微信貸的孟加拉國發(fā)放小微企業(yè)信貸的風(fēng)險情況進(jìn)行研究分析,了解到影響小微企業(yè)信用風(fēng)險的因素有貸款數(shù)量、貸款人抵押品價值、其所擁有的耐用商品的價值、貸款期限的長短等. Meyer[13]指出銀行可以通過“軟信息”對小企業(yè)信用狀況進(jìn)行評估,這些軟信息包括銀行對借款企業(yè)資信的初步判斷等.Lussier[14]在其模型中加入了企業(yè)管理水平、所處經(jīng)濟(jì)周期、產(chǎn)品生命周期等因素.此外,張良貴等[15]基于DSGE模型研究表明宏觀流動性及其狀態(tài)變化、基準(zhǔn)利率狀態(tài)變化對企業(yè)杠桿有重要影響.由此可見,小微企業(yè)信用評級指標(biāo)需要同時包含企業(yè)內(nèi)部定量指標(biāo)和定性指標(biāo)與企業(yè)外部信息.

    在研究方法方面,張大斌等[16]用差分進(jìn)化算法(DE)優(yōu)化違約點(diǎn)系數(shù),提出一種中國上市公司信用風(fēng)險測度的不確定性DE-KMV模型,來提高上市公司信用風(fēng)險測度的準(zhǔn)確性.于立勇等[17]在其研究中開展基于Logistic回歸分析的違約概率研究,利用正向逐步選擇法建立了信用風(fēng)險評估指標(biāo)體系.黃苒等[18]重點(diǎn)關(guān)注了含跳躍風(fēng)險的公司貸款違約率測度問題,并探索了基于首達(dá)時模型的理論擴(kuò)展,給出了違約概率參數(shù)估計的方法.龐素琳等[19]在銀行個人信用評級中建立C5.0分類算法,在該算法內(nèi)嵌入Boosting算法技術(shù),并構(gòu)造了成本矩陣和參數(shù)調(diào)整后的決策樹.肖進(jìn)等[20]提出了面向缺失數(shù)據(jù)的動態(tài)分類器集成選擇模型DCESM,并認(rèn)為該模型能夠取得更好的客戶信用評估性能.Che等[21]運(yùn)用層次分析法(AHP)和包絡(luò)分析法(DEA)對臺灣地區(qū)商業(yè)銀行小微企業(yè)信貸決策進(jìn)行研究,并得到有效的建模結(jié)果.Angilella等[22]在研究中建立了小微企業(yè)SMMAA-TRI多準(zhǔn)則信用評級模型,并對模型有效性進(jìn)行了驗證.同時,非參數(shù)統(tǒng)計的研究方法也逐步應(yīng)用于信用評級建模.然而,隨著銀行業(yè)可訓(xùn)練數(shù)據(jù)量的大幅提升和數(shù)據(jù)維度的不斷提高,機(jī)器學(xué)習(xí)方法開始嶄露頭角.統(tǒng)計學(xué)和機(jī)器學(xué)習(xí)最本質(zhì)的區(qū)別在于根本目標(biāo)不同,統(tǒng)計學(xué)家更關(guān)心模型的可解釋性,而機(jī)器學(xué)習(xí)專家更關(guān)心模型的預(yù)測能力.由于銀行業(yè)對預(yù)測準(zhǔn)確度要求的越來越高,統(tǒng)計學(xué)方法的不適應(yīng)性開始顯現(xiàn),很多問題不能構(gòu)建出嚴(yán)謹(jǐn)?shù)慕y(tǒng)計模型,而機(jī)器學(xué)習(xí)理論已被證明是此類研究有效的建模方法.王春峰等[23]在其研究中對人工神經(jīng)網(wǎng)絡(luò)在銀行信用風(fēng)險評價中的應(yīng)用展開對比討論,結(jié)果表明,與傳統(tǒng)統(tǒng)計方法相比,神經(jīng)網(wǎng)絡(luò)技術(shù)具有較高的穩(wěn)定性和判別準(zhǔn)確度.Lee等[24]認(rèn)為非參數(shù)統(tǒng)計的方法在信用評級的研究中優(yōu)于經(jīng)典統(tǒng)計模型.Yu等[25]運(yùn)用基于多層神經(jīng)網(wǎng)絡(luò)的機(jī)器學(xué)習(xí)方法研究信用評分模型,并對效果良好的建模結(jié)果進(jìn)行敘述.Zhong等[26]在企業(yè)信用評級中,對BP、ELM、I-ELM和SVM的學(xué)習(xí)訓(xùn)練效果進(jìn)行對比分析,其中ELM和BP神經(jīng)網(wǎng)絡(luò)的建模效果較優(yōu).

    2 信用評級模型及模型有效性衡量方法

    商業(yè)銀行在信貸審批操作中運(yùn)用信用評級模型對申請借款企業(yè)提供的信息進(jìn)行判斷,對企業(yè)在未來是否會出現(xiàn)違約進(jìn)行估計,已經(jīng)是普遍存在的貸前審批流程,因此評級模型的有效性是各家商業(yè)銀行共同關(guān)注的問題.結(jié)合現(xiàn)有文獻(xiàn)中的信用評級方法,本文構(gòu)建包含線性判別分析(linear discriminant analysis, LDA)、二項邏輯回歸(binary logistic regression,BLR)、基于多種學(xué)習(xí)算法的BP神經(jīng)網(wǎng)絡(luò)(BP neural network, BPNN)的評分模型,運(yùn)用建模樣本訓(xùn)練模型,并用測試樣本檢驗?zāi)P陀行院头€(wěn)健性,挑選出較優(yōu)的信用評級模型.LDA是在分類確定的條件下,根據(jù)某一研究對象的各個特征值判別其類型歸屬問題的一種多變量統(tǒng)計分析方法.BLR模型主要是用來對多因素影響的事件進(jìn)行概率預(yù)測,是普通多元線性回歸模型向非線性模型的擴(kuò)展.基于多種學(xué)習(xí)算法的BPNN模型是前饋型神經(jīng)網(wǎng)絡(luò)模型,其將輸入的多個數(shù)據(jù)集映射到單一的輸出的數(shù)據(jù)集上,通過訓(xùn)練與測試將對象進(jìn)行分類.本節(jié)就LDA、BLR、基于多種學(xué)習(xí)算法的BPNN模型以及模型有效性的衡量方法進(jìn)行簡要介紹.

    2.1 線性判別分析(LDA)

    對于k個組別的分類問題,假設(shè)k個組別構(gòu)成的總體分別為G1,G2,...,Gk,于是若要判斷樣本x各來自于哪一個總體,首先必須計算樣本x到每個總體Gi(i=1,2,...,k)的距離d(x,Gi),然后再比較這些距離,其中樣本x到總體Gi的距離d(x,Gi)采用Mahalanobis距離,即其中μi和Σi分別為Gi的均值和協(xié)方差矩陣,為協(xié)方差矩陣Σi的逆矩陣.如果x距某個Gj(j= 1,2,...,k)最近,則認(rèn)為x∈Gj.其判別規(guī)則為

    對于本文組別分類的問題,可通過建立判別分析模型進(jìn)行判別.

    2.2 二項邏輯回歸(BLR)

    對于一個二項分類和n個定量預(yù)測變量x1,x2,...,xn(包含虛擬賦值變量),BLR模型假設(shè)目標(biāo)響應(yīng)的概率為

    其中β0為常數(shù)項,稱β1,β2,...,βn為模型回歸系數(shù),解釋變量可以是連續(xù)變量,也可以是分類變量或啞變量(dummy variable).

    正式的決策框架中通常包含多種評判模型,在多種模型有效性對比過程中,BLR模型可以作為一個判斷依據(jù).BLR與LDA一樣,在多元正態(tài)分布和相同協(xié)方差矩陣的假設(shè)條件下具有最優(yōu)判別能力.BLR要求較大的輸入樣本量以取得較為穩(wěn)定的計算結(jié)果,同時應(yīng)對自變量與因變量的復(fù)雜非線性關(guān)系進(jìn)行適當(dāng)處理.

    2.3 基于多種學(xué)習(xí)算法的BPNN

    BPNN的構(gòu)成包括一個輸入層,一個或多個隱含層和一個輸出層.各層由若干個神經(jīng)元構(gòu)成,每一個節(jié)點(diǎn)的輸出值由輸入值、作用函數(shù)和閾值決定.神經(jīng)網(wǎng)絡(luò)的學(xué)習(xí)過程包含兩個方面:信息正向傳播和誤差反向傳播.在正向傳播過程中,輸入信息從輸入層經(jīng)隱含層傳到輸出層,經(jīng)作用函數(shù)運(yùn)算后得到輸出值與期望值進(jìn)行比較,若有誤差,則誤差反向傳播,沿原先的連接通路返回,通過逐層修改各層神經(jīng)元的權(quán)值,減少誤差,如此循環(huán)直到輸出的結(jié)果符合精度要求為止.具體步驟如下:

    步驟1BPNN的初始化,確定隱含層節(jié)點(diǎn)的個數(shù).將各個權(quán)值和閾值的初始值設(shè)為比較小的隨機(jī)數(shù);

    步驟2輸入樣本和相應(yīng)的輸出進(jìn)行網(wǎng)絡(luò)訓(xùn)練,即對每一個樣本數(shù)據(jù)進(jìn)行步驟3到步驟5的過程;

    步驟3依據(jù)輸入樣本計算實際輸出及其隱含層神經(jīng)元的輸出;

    步驟4計算期望輸出與實際輸出之間的差值,求輸出層和隱含層的誤差;

    步驟5根據(jù)步驟4得出的誤差更新輸入層到隱含層節(jié)點(diǎn)之間以及隱含層到輸出層節(jié)點(diǎn)和之間的連接權(quán)值;

    步驟6求出誤差函數(shù)e,判斷e是否收斂到給定的學(xué)習(xí)精度以內(nèi),即e≤ε,其中ε為擬定誤差,如果滿足則訓(xùn)練結(jié)束,否則轉(zhuǎn)向步驟2.

    運(yùn)用不同的學(xué)習(xí)算法對BPNN進(jìn)行訓(xùn)練,將會得到不同的建模效果[27].MATLAB R2012a軟件中的神經(jīng)網(wǎng)絡(luò)工具箱(neural network toolbox)中包括多種適用于BPNN建模的學(xué)習(xí)算法.本文考慮運(yùn)用10種學(xué)習(xí)算法對BPNN進(jìn)行訓(xùn)練,具體如表1所示.

    表1 基于多種學(xué)習(xí)算法的BPNNTable 1 The BPNN based on multiple learning algorithms

    2.4 模型有效性的衡量方法

    ROC曲線面積(AUC值)是常用的用于評價分類模型有效性的方法.ROC曲線是根據(jù)一系列不同的二分類方式,以真陽性率(敏感度,不良貸款樣本被判正確)為縱坐標(biāo),假陽性率(1–特異性,良性貸款被判錯誤)為橫坐標(biāo)繪制的曲線.傳統(tǒng)的診斷試驗評價方法有一個共同的特點(diǎn),必須將試驗結(jié)果分為兩類,再進(jìn)行統(tǒng)計分析.本文運(yùn)用SPSS軟件對12種模型的檢驗結(jié)果繪制ROC曲線,并比較AUC值,AUC值越大,表明模型判別水平越高.

    均方誤差(mean squared error,MSE)是指參數(shù)估計值與參數(shù)真值之差平方的期望值.在此,MSE是預(yù)測值與期望值之差平方的期望值,即

    其中

    y

    t

    是期望值,

    t

    是預(yù)測值.MSE的值越小,說明預(yù)測模型描述實驗數(shù)據(jù)具有更好的精確度.

    在對模型有效性進(jìn)行衡量時,需要考慮不同錯誤類型的成本,本文依據(jù)Blanco等[27]的兩種錯誤分類構(gòu)造誤判成本.I類錯誤指良性貸款錯評為不良貸款的比率;II類錯誤指不良貸款錯評為良性貸款的比率,相比于I類錯誤,II類錯誤的發(fā)生將會對銀行帶來更大的損失.因此本文將II類錯誤的成本設(shè)為I類錯誤成本的5倍.并在此基礎(chǔ)上計算誤判成本

    其中C21和C12分別表示發(fā)生I類錯誤和II類錯誤的成本,π21和π12分別表示發(fā)生I類錯誤和II類錯誤的概率,p1和p2分別表示樣本在到期時是良性貸款和不良貸款的先驗概率.

    識別率包括不良識別率和總識別率.不良識別率表示檢測樣本中的不良貸款被模型識別出的比例,該指標(biāo)可以反映模型對I類錯誤的規(guī)避能力;總識別率表示全部檢測樣本的二級分類屬性被識別正確的比例,該指標(biāo)可以反映模型的整體判別能力.

    3 數(shù)據(jù)與變量

    本文采集江蘇某商業(yè)銀行2008–2013針對小微企業(yè)信貸的過往數(shù)據(jù),與公開數(shù)據(jù)庫中的宏觀數(shù)據(jù)合并后,總數(shù)據(jù)內(nèi)容主要包括:1)財務(wù)信息(資產(chǎn)總額,經(jīng)營收入等);2)本期貸款基本信息(發(fā)放額度,發(fā)放利率,貸款方式,支付方式等);3)當(dāng)前客戶非財務(wù)信息(客戶信用等級,客戶資信等級,行業(yè)分類,擔(dān)保情況等); 4)與宏觀經(jīng)濟(jì)有關(guān)的變量(用電量,GDP,CPI等);5)信貸的二級分類情況.經(jīng)過剔除缺失和不正常的數(shù)據(jù),得到2 115組有效數(shù)據(jù),其中,良性與不良貸款比例約為9∶1.為了對評分模型進(jìn)行科學(xué)對比,我們把數(shù)據(jù)集隨機(jī)分成兩個互補(bǔ)的子數(shù)據(jù)集:80%數(shù)據(jù)作為建模集,20%數(shù)據(jù)作為測試集.各模型的變量結(jié)構(gòu)都是通過10次交叉驗證法進(jìn)行篩選,交叉驗證法的一個優(yōu)勢就是這樣的信用評分模型是基于較大比例的有效數(shù)據(jù)(80%)開發(fā).

    有關(guān)信用評級指標(biāo)的文獻(xiàn)中,由“硬”信息到“軟”信息的選擇是信用評級指標(biāo)選擇的一個顯著性變化,從最早標(biāo)準(zhǔn)化的財務(wù)指標(biāo)到如今各類非財務(wù)指標(biāo)的不斷加入更新,信用評級系統(tǒng)也因此能夠從規(guī)范化的大公司推廣到小微企業(yè)及個人類客戶.然而,這些指標(biāo)多為微觀指標(biāo),宏觀指標(biāo)鮮有涉及.近年來,有學(xué)者對宏觀經(jīng)濟(jì)變量的研究提出了自己的看法,認(rèn)為借款人的違約與一般經(jīng)濟(jì)狀況密切相關(guān),與宏觀經(jīng)濟(jì)環(huán)境變量也應(yīng)當(dāng)被考慮作為輸入變量[27,28].顏新秀[29]認(rèn)為,不同宏觀經(jīng)濟(jì)環(huán)境下各指標(biāo)變量與個人住房抵押貸款違約率之間存在一定的影響關(guān)系.Kim等[30]指出宏觀經(jīng)濟(jì)環(huán)境是一個關(guān)鍵因素,直接關(guān)系到借款者的償付行為.因此,本文認(rèn)為貸款期間內(nèi)的宏觀經(jīng)濟(jì)變動對貸款主體的還款意愿具有一定影響,應(yīng)當(dāng)作為重要變量參與信用評級研究.因此,本文選取的指標(biāo)類別包括財務(wù)類指標(biāo),非財務(wù)指標(biāo)和宏觀環(huán)境指標(biāo).

    建模自變量和因變量(輸入變量和輸出變量)如表2和表3所示.本文在考慮現(xiàn)有文獻(xiàn)貸款對象財務(wù)信息、非財務(wù)信息、本期貸款基本信息的基礎(chǔ)上,還加入與宏觀經(jīng)濟(jì)有關(guān)的變量作為自變量或輸入變量,主要基于以下考慮:傳統(tǒng)模型的因變量和自變量存在不同期的問題,傳統(tǒng)信用評級模型將貸款期初的相關(guān)變量指標(biāo)作為自變量,將貸款期末的五級分類情況作為因變量來進(jìn)行建模,即自變量和因變量不是同一時間點(diǎn)上的面板數(shù)據(jù),因此不能直接看作面板數(shù)據(jù)處理,需加入時間因素才能解決矛盾,同時,宏觀經(jīng)濟(jì)變量是時間因素的合理體現(xiàn),因為宏觀經(jīng)濟(jì)變量并非某一時間點(diǎn)上的宏觀經(jīng)濟(jì)變量,而是貸款期間相關(guān)宏觀經(jīng)濟(jì)變量變化趨勢的反映,因此不是一個時間點(diǎn)的概念而是某一時間段的變化趨勢的概念,所以能夠反映出時間因素.在實證分析部分,本文將分別運(yùn)用包含宏觀變量和不包含宏觀變量的指標(biāo)體系進(jìn)行建模,并比較模型的有效性.同時,本文借鑒Blanco等[27]的做法,在指標(biāo)體系中加入“客戶經(jīng)理主觀評價”和“客戶資信狀況調(diào)查結(jié)果”兩個指標(biāo),以囊括銀行自有的貸款初期審核信息.

    本文引入貸款期限內(nèi)經(jīng)濟(jì)周期信息的變量,Kim等[30]指出宏觀經(jīng)濟(jì)環(huán)境是一個關(guān)鍵因素,直接關(guān)系到借款者的償付行為.該做法的優(yōu)點(diǎn)還在于考慮到還貸期限內(nèi)省內(nèi)宏觀經(jīng)濟(jì)變化對小微企業(yè)的影響,由于各家商業(yè)銀行在計劃財務(wù)部進(jìn)行財務(wù)預(yù)算的編制過程中均擁有未來至少三年的宏觀經(jīng)濟(jì)變化預(yù)測數(shù)值,因此在實務(wù)操作中,該宏觀變量數(shù)據(jù)采用商業(yè)銀行預(yù)測數(shù)值.

    表2 建模自變量(輸入變量)名稱及解釋Table 2 The names and explanations for independent(input)variables

    表3 建模因變量(輸出變量)名稱及解釋Table 3 The names and explanations for dependent(output)variables

    本文借鑒Blanco等[27]在設(shè)計信用評分模型所使用的宏觀經(jīng)濟(jì)變量的表達(dá)式

    其中ΔVMi,j是宏觀經(jīng)濟(jì)變量的變化率,VM是宏觀經(jīng)濟(jì)變量,i表示放款的時間,j表示貸款的時限.

    本文在對宏觀濟(jì)變量進(jìn)行篩選之后選取用電量、GDP和CPI的季度數(shù)據(jù)作為建模的三個自變量.在實踐中發(fā)現(xiàn),除CPI外用電量和GDP均有很強(qiáng)的周期性,經(jīng)過分析和比較之后,本文決定采用GDP、用電量的同比增量的變化率和CPI增長率來作為自變量,其中CPI增長率為(計算期CPI數(shù)值-基期CPI數(shù)值)×100%/基期CPI數(shù)值,GDP、用電量的同比增量的變化率為

    其中Δ2VMi,j是宏觀經(jīng)濟(jì)變量增量的變化率,i表示放款的時間,j表示貸款的時限.

    4 實證分析

    綜合數(shù)據(jù)可獲取性、數(shù)據(jù)可靠性因素的影響,樣本選取2008–2013年貸款數(shù)據(jù),其中宏觀變量數(shù)據(jù)均進(jìn)行季節(jié)性調(diào)整,并對離散變量進(jìn)行賦值.神經(jīng)網(wǎng)絡(luò)建模時為加快訓(xùn)練網(wǎng)絡(luò)的收斂性的需要,對樣本數(shù)據(jù)中各個指標(biāo)下的數(shù)據(jù)分別進(jìn)行歸一化處理.如第3節(jié)所述,將樣本分為兩部分,其中1 687組數(shù)據(jù)作為訓(xùn)練樣本,428組數(shù)據(jù)作為檢驗樣本,訓(xùn)練樣本與檢驗樣本比例約為4∶1.

    運(yùn)用SPSS 21軟件,以表2中15個變量為判別變量,以表3中B1變量為分組變量,選擇Fisher函數(shù)系數(shù),建立LDA模型,得到分類函數(shù)系數(shù).根據(jù)分類函數(shù)系數(shù),對檢驗樣本中各指標(biāo)進(jìn)行加權(quán)求和,選取得數(shù)較大的類別作為判別結(jié)果.

    運(yùn)用PASW Statistics 18統(tǒng)計分析軟件建立二項邏輯回歸(BLR)通過描述性統(tǒng)計,根據(jù)相關(guān)性、正態(tài)性檢驗和t檢驗,本文最終在Logistic回歸中保留9項指標(biāo),分別為年利率、貸款方式、用途、支付方式、客戶信用等級、客戶資信等級、營業(yè)收入/人、GDP增量增長率和CPI增量增長率,具體分析過程介紹從略.

    運(yùn)用MATLAB R2012a建立BP神經(jīng)網(wǎng)絡(luò)(BPNN),表2中15個變量作為輸入神經(jīng)元,表3中B2變量各分別作為輸出神經(jīng)元.分別運(yùn)用表2中10種學(xué)習(xí)算法進(jìn)行模型訓(xùn)練,各個訓(xùn)練函數(shù)即表1中MATLAB調(diào)用函數(shù).神經(jīng)網(wǎng)絡(luò)參數(shù)設(shè)定:隱含層節(jié)點(diǎn)范圍[7,13];最大迭代次數(shù)=1 000;訓(xùn)練誤差=0.001.

    運(yùn)用訓(xùn)練后的模型對檢驗樣本進(jìn)行測試,即分別將檢驗樣本中的自變量與LDA分類函數(shù)系數(shù)加權(quán)求和取大值,代入BLR模型預(yù)測因變量值,代入訓(xùn)練好的BPNN模型得出預(yù)測輸出變量值.當(dāng)輸出值小于并等于0.5時,判定該組樣本為正常貸款;當(dāng)輸出值大于0.5時,判定該組樣本為不正常貸款.首先,運(yùn)用不包含宏觀變量的指標(biāo)體系建立模型,根據(jù)前文模型有效性的衡量方法,計算衡量各個模型有效性的指標(biāo)數(shù)值,如表4所示,優(yōu)化的ROC曲線如圖1所示.在ROC曲線圖中,縱坐標(biāo)δ代表真陽性率,橫坐標(biāo)1-η代表假陽性率,其中η為特異性,下同.

    表4 模型有效性衡量方法的計算數(shù)值(模型不含宏觀變量)Table 4 The effectiveness measurement value of different models(exclude macroeconomic variables)

    由表4可知,12種模型中,75%的模型的AUC值大于0.8;有11個模型的總識別率超過80%,其中六個模型的總識別率超過90%.但在關(guān)鍵指標(biāo)中,有11個模型的不良識別率低于80%,有11個模型的誤判成本大于0.2.從整體識別效果來看,所構(gòu)建的模型對不良貸款缺乏識別能力是一個共性,因此本文認(rèn)為有必要嘗試通過加入宏觀變量對指標(biāo)體系進(jìn)行改進(jìn)并開展建模.在指標(biāo)體系中加入宏觀變量后,判斷各個模型有效性的指標(biāo)數(shù)值如表5所示,優(yōu)化的ROC曲線如圖2所示.

    1)AUC方面,共有6個模型的AUC值超過0.9,其中,NN6與NN10模型的AUC值超過0.95,分別為0.952和0.959,預(yù)測能力較強(qiáng);2)均方誤差方面,共有8個模型的均方誤差小于0.1,其中BLR的均方誤差為0.071,NN3的均方誤差為0.072,NN6的均方誤差為0.075;3)誤判成本方面,共有5個模型的誤判成本低于0.2,其中NN6誤判成本為0.117,NN7誤判成本為0.110,NN10誤判成本為0.098,誤判成本較低.I類錯誤率方面,NN6、NN7、NN10的II類錯誤率小于0.1,判別效果較優(yōu);4)不良識別率方面,共有4個模型的不良識別率超過0.8,其中NN6的不良識別率為0.915,NN7的不良識別率為1.000,NN10的不良識別率為0.979,說明該三個模型對潛在不良貸款的判別能力較強(qiáng).總識別率方面,共有5個模型的總識別率超過0.9,即表明這5個模型的總體判別正確率超過九成,具有較好的整體判別能力.

    圖1 ROC曲線(模型不含宏觀變量)Fig.1 The ROC curve(the model does not contain macro variables)

    表5 模型有效性衡量方法的計算數(shù)值(模型含宏觀變量)Table 5 The effectiveness measurement value of different models(include macroeconomic variables)

    圖2 ROC曲線(模型含宏觀變量)Fig.2 The ROC curve(the model contains macro variables)

    根據(jù)以上分析,可見LDA和BLR模型在四類模型有效性衡量方法中表現(xiàn)一般,其中BLR模型在AUC和總識別率方面優(yōu)于LDA模型,LDA模型在II類錯誤率方面優(yōu)于BLR模型,但該兩個模型劣于或等于NN6和NN10在AUC、II類錯誤率、誤判成本、不良識別率和總識別率方面的表現(xiàn).當(dāng)將NN6與NN10在模型有效性方面進(jìn)行比較時,NN10在AUC、均方誤差、II類錯誤率、誤判成本和不良識別率的表現(xiàn)均優(yōu)于NN6,而NN6僅在總識別率方面優(yōu)于NN10,因此認(rèn)為NN10具有最優(yōu)的模型有效性,其次為NN6.

    由上述分析結(jié)果可知,1)在整體模型有效性方面,NN10(基于Levenbery-Marquardt學(xué)習(xí)算法)模型表現(xiàn)最優(yōu),NN6(基于Polak-Ribiere共軛梯度學(xué)習(xí)算法)模型表現(xiàn)次優(yōu).2)在各模型有效性衡量方法方面, NN10具有最大的AUC值,BLR具有最小的均方誤差,NN7具有最小的II類錯誤率,NN10具有最小的誤判成本,NN7具有最高的不良識別率,NN3具有最高的總識別率,這些結(jié)果表明在進(jìn)行具體信用評級水平衡量時,需要重視上述模型在信用評級中的應(yīng)用.3)在觀察各個模型檢測結(jié)果中出現(xiàn)兩類錯誤的樣本數(shù)據(jù)時,本文發(fā)現(xiàn)樣本數(shù)據(jù)中的客戶信用等級、客戶資信等級、營業(yè)收入/人指標(biāo)易出現(xiàn)區(qū)別于普通的數(shù)值,因此此類樣本需要在人工貸前審查時給予關(guān)注.

    5 結(jié)束語

    本研究選擇合適的指標(biāo)體系和研究方法,構(gòu)建線性判別分析、二項邏輯回歸和基于多種學(xué)習(xí)算法的BP神經(jīng)網(wǎng)絡(luò)模型對小微企業(yè)進(jìn)行信用評級,并通過四種衡量方法分析模型的有效性,提高研究結(jié)論的可信度.研究結(jié)果表明機(jī)器學(xué)習(xí)方法能夠為商業(yè)銀行授信審批和貸前審查工作提供參考依據(jù),并且在指標(biāo)體系中加入宏觀變量能夠有效地提高模型識別不良貸款的能力,提升模型的穩(wěn)健性.因BP神經(jīng)網(wǎng)絡(luò)具有容易陷入局部極值、“過擬合”等問題,其預(yù)測精度有待優(yōu)化,今后研究中嘗試使用遺傳算法優(yōu)化BP神經(jīng)網(wǎng)絡(luò)強(qiáng)化全局搜索能力.取得更大的數(shù)據(jù)樣本,并嘗試根據(jù)小微企業(yè)資產(chǎn)總量大小對樣本進(jìn)行分層,提高模型對于貸款主體的評級能力.

    [1]Schreiner M.Scoring arrears at a microlender in Bolivia.Journal of Microfnance,2004,6(2):65–88.

    [2]董卉娜,朱志雄.上市銀行內(nèi)部控制缺陷認(rèn)定研究.證券市場導(dǎo)報,2014,24(6):57–62. Dong H N,Zhu Z X.Identifcation of internal control defciencies in listed banks.Securities Market Herald,2014,24(6):57–62.(in Chinese)

    [3]屈文洲,謝雅璐,葉玉妹.信息不對稱,融資約束與投資–現(xiàn)金流敏感性:基于市場微觀結(jié)構(gòu)理論的實證研究.經(jīng)濟(jì)研究,2011, 46(6):105–117. Qu W Z,Xie Y L,Ye Y M.Information asymmetry and investment-cash fow sensitivity:An empirical research based on market microstructure theory.Economic Research Journal,2011,46(6):105–117.(in Chinese)

    [4]郭小波,王婉婷,周 欣.我國中小企業(yè)信貸風(fēng)險識別因子的有效性分析:基于北京地區(qū)中小企業(yè)的信貸數(shù)據(jù).國際金融研究, 2011,28(4):62–67. Guo X B,Wang W T,Zhou X.Effectiveness of identifying factors of SME credit risk:Evidence from SMEs in Beijing.Studies of International Finance,2011,28(4):62–67.(in Chinese)

    [5]謝 平,許國平,李 德.運(yùn)用信用評級原理加強(qiáng)金融監(jiān)管.管理世界,2001,17(1):125–131. Xie P,Xu G P,Li D.Using credit rating theory to strengthen fnancial regulation.Management World,2001,17(1):125–131.(in Chinese)

    [6]徐 超,周宗放.中小企業(yè)聯(lián)保貸款信用行為演化博弈仿真研究.系統(tǒng)工程學(xué)報,2014,29(4):477–486. Xu C,Zhou Z F.The evolutionary game simulation on credit behavior of SMEs’guaranteed loans.Journal of Systems Engineering, 2014,29(4):477–486.(in Chinese)

    [7]陳其安,張紅真,高國婷.基于地方政府擔(dān)保的投融資平臺融資行為模型.系統(tǒng)工程學(xué)報,2015,30(3):319–330. Chen Q A,Zhang H Z,Gao G T.Financing behavior model of the investing and fnancing platform based on local government guarantee.Journal of Systems Engineering,2015,30(3):319–330.(in Chinese)

    [8]李 毅,向 黨.中小企業(yè)信貸融資信用擔(dān)保缺失研究.金融研究,2008,30(12):179–192. Li Y,Xiang D.Research on the lack of credit guarantee in fnancing of SMEs.Journal of Financial Research,2008,30(12):179–192. (in Chinese)

    [9]Hajek P,Michalak K.Feature selection in corporate credit rating prediction.Knowledge-Based Systems,2013,51(3):72–84.

    [10]Doumpos M,Niklis D,Zopounidis C,et al.Combining accounting data and a structural model for predicting credit ratings:Empirical evidence from European listed frms.Journal of Banking&Finance,2015,50(5):599–607.

    [11]Vassiliou P C G.Fuzzy semi-Markov migration process in credit risk.Fuzzy Sets and Systems,2013,223(4):39-58.

    [12]Malhotra R,Malhotra D K.Differentiating between good credits and bad credits using neuro-fuzzy systems.European Journal of Operational Research,2002,136(1):190-211.

    [13]Meyer L H.The present and future roles of banks in small business fnance.Journal of Banking&Finance,1998,22(6):1109-1116.

    [14]Lussier R N,A cross-national prediction model for business success.Journal of Small Business Management.2001,39(3):228–239.

    [15]張良貴,孫久文.信用利差變化對企業(yè)杠桿的動態(tài)影響.證券市場導(dǎo)報,2014,24(6):53–56. Zhang L G,Sun J W.Dynamic effect of changes in credit spreads on corporate leverage.Securities Market Herald,2014,24(6): 53–56.(in Chinese)

    [16]張大斌,周志剛,劉 雯,等.上市公司信用風(fēng)險測度的不確定性DE-KMV模型.系統(tǒng)工程學(xué)報,2015,30(2):165–173. Zhang D B,Zhou Z G,Liu W,et al.Uncertainty DE-KMV model to measure credit risk of listed companies.Journal of Systems Engineering,2015,30(2):165–173.(in Chinese)

    [17]于立勇,詹捷輝.基于Logistic回歸分析的違約概率預(yù)測研究.財經(jīng)研究,2004,30(9):15–23. Yu L Y,Zhan J H.Prediction of default rate based on logistic regression analysis.Study of Finance and Economics,2004,30(9): 15–23.(in Chinese)

    [18]黃 苒,唐齊鳴.含跳躍風(fēng)險的公司貸款違約率測度:基于首達(dá)時模型的理論擴(kuò)展.管理科學(xué)學(xué)報,2015,18(7):93–102. Huang R,Tang Q M.The measurement of loan default rate of frms with jumping risk:Based on the theory of frst reaching time model.Journal of Management Science in China,2015,18(7):93–102.(in Chinese)

    [19]龐素琳,鞏吉璋.C5.0分類算法及在銀行個人信用評級中的應(yīng)用.系統(tǒng)工程理論與實踐,2009,29(12):94–104. Pang S L,Gong J Z.C5.0 classifcation algorithm and its application on individual credit score for banks.Systems Engineering: Theory&Practice,2009,29(12):94–104.(in Chinese)

    [20]肖 進(jìn),劉敦虎,顧 新,等.銀行客戶信用評估動態(tài)分類器集成選擇模型.管理科學(xué)學(xué)報,2015,18(3):114–126. Xiao J,Liu D H,Gu X,Wang S Y.Dynamic classifer ensemble selection model for bank customer’s credit scoring.Journal of Management Science in China,2015,18(3):114–126.(in Chinese)

    [21]Che Z H,Wang H S,Chuang C L.A fuzzy AHP and DEA approach for making bank loan decisions for small and medium enterprises in Taiwan.Expert Systems with Applications,2010,37(10):7189–7199.

    [22]Angilella S,Mazz`u S.The fnancing of innovative SMEs:A multicriteria credit rating model.European Journal of Operational Research,2013,244(2):540–554.

    [23]王春峰,萬海暉,張 維.基于神經(jīng)網(wǎng)絡(luò)技術(shù)的商業(yè)銀行信用風(fēng)險評估.系統(tǒng)工程理論與實踐,1999 19(9):24–32. Wang C F,Wan H H,Zhang W.Credit risk assessment in commercial banks using neural networks.Systems Engineering:Theory& Practice,1999,19(9):24–32.(in Chinese)

    [24]Lee T S,Chen I F.A two-stage hybrid credit scoring model using artifcial neural networks and multivariate adaptive regression splines.Expert Systems with Applications,2005,28(4):743–752.

    [25]Yu L,Wang S,Lai K K.Credit risk assessment with a multistage neural network ensemble learning approach.Expert Systems with Applications,2008,34(2):1434–1444.

    [26]Zhong H,Miao C,Shen Z,et al.Comparing the learning effectiveness of BP,ELM,I-ELM,and SVM for corporate credit ratings. Neurocomputing,2013,128(5):285–295.

    [27]Blanco A,Pino-Mejˊ?as R,Lara J,et al.Credit scoring models for the microfnance industry using neural networks:Evidence from Peru.Expert Systems with Applications,2013,40(1):356–364.

    [28]Tinoco M H,Wilson N.Financial distress and bankruptcy prediction among listed companies using accounting,market and macroeconomic variables.International Review of Financial Analysis,2013,30(4):394–419.

    [29]顏新秀.個貸違約率與宏觀經(jīng)濟(jì)指標(biāo)相關(guān)性研究.國際金融研究,2009,26(10):59–67. Yan X X.Correlation of default rate of personal loan and macroeconomic indicators.Studies of International Finance,2009,26(10): 59–67.(in Chinese)[30]Kim H S,Sohn S Y.Support vector machines for default prediction of SMEs based on technology credit.European Journal of Operational Research,2010,201(3):838–846.

    Comparative study of credit rating models for small and micro enterprises

    Xiao Binqing1,Yang Yang2,Yu Zhe3,Shen Caisheng1,4
    (1.School of Engineering and Management,Nanjing University,Nanjing 210093,China; 2.School of Business,Nanjing University,Nanjing 210093,China; 3.Zhengzhou Commodity Exchange,Zhengzhou 450008,China; 4.Zijin Rural Commercial Bank,Nanjing 210019,China)

    On the basis of investigation and literature research,conducting risk source identifcation,rating indicators classifcation and rating methods assessment,the paper constructs commercial bank’s internal credit rating models to improve the credit risk management in the credit approval procedures.Based on the credit data of small and micro enterprises in a commercial bank from 2008 to 2013,using the linear discriminant analysis, logistic regression and 10 types of BP neural network relying on different learning algorithms,internal credit rating models are constructed with macroeconomic variables,which may further improve the robustness of risk measurement.Finally,results and rating effectiveness of different models are analyzed and compared, and show that the NN10 model based on Levenbery-Marquardt learning algorithm performs optimal rating effectiveness.

    credit rating;small and micro enterprises;model comparison;learning algorithm

    F832.59

    A

    1000-5781(2016)06-0798-10

    10.13383/j.cnki.jse.2016.06.008

    肖斌卿(1979―),男,福建南靖人,博士,副教授,研究方向:金融工程與金融管理,Email:bengking@nju.edu.cn;

    楊 旸(1990―),男,江蘇南京人,博士生,研究方向:金融工程,Email:yangyang68nj@163.com;

    余 哲(1990―),男,河南新野人,碩士,結(jié)算部主管,研究方向:金融工程與金融管理,Email:yuzhe19900915@163.com;

    沈才勝(1982―),男,安徽無為人,博士生,風(fēng)險管理部,研究方向:金融工程,Email:starting_scs@126.com.

    2016-01-07;

    2016-05-19.

    國家自然科學(xué)基金重點(diǎn)資助項目(70932003);國家自然科學(xué)基金資助項目(71271109;71201074;70901037;71271-110;71501131);教育部科技創(chuàng)新工程重大項目培育資金資助項目(708044);教育部人文社會科學(xué)研究青年資助項目(13YJC790174).

    *通信作者

    猜你喜歡
    識別率評級小微
    小微課大應(yīng)用
    基于類圖像處理與向量化的大數(shù)據(jù)腳本攻擊智能檢測
    基于真耳分析的助聽器配戴者言語可懂度指數(shù)與言語識別率的關(guān)系
    提升高速公路MTC二次抓拍車牌識別率方案研究
    小微企業(yè)借款人
    分析師最新給予買入評級的公司
    高速公路機(jī)電日常維護(hù)中車牌識別率分析系統(tǒng)的應(yīng)用
    小微企業(yè)經(jīng)營者思想動態(tài)調(diào)查
    百度遭投行下調(diào)評級
    IT時代周刊(2015年8期)2015-11-11 05:50:22
    解決小微金融機(jī)構(gòu)的風(fēng)控難題
    91国产中文字幕| 老司机靠b影院| 91老司机精品| 亚洲色图 男人天堂 中文字幕| 老司机靠b影院| 成人精品一区二区免费| 亚洲一区高清亚洲精品| 欧美黑人巨大hd| 国产精品九九99| avwww免费| 悠悠久久av| 又紧又爽又黄一区二区| 50天的宝宝边吃奶边哭怎么回事| 精品久久蜜臀av无| 99在线视频只有这里精品首页| 久久久久久久精品吃奶| 国产精品日韩av在线免费观看| 午夜a级毛片| 欧美+亚洲+日韩+国产| 日本黄色视频三级网站网址| 成人永久免费在线观看视频| 怎么达到女性高潮| 人妻久久中文字幕网| 亚洲人成电影免费在线| av视频在线观看入口| 91成年电影在线观看| 91麻豆精品激情在线观看国产| 日韩成人在线观看一区二区三区| 精品一区二区三区视频在线观看免费| 免费在线观看视频国产中文字幕亚洲| 99热这里只有是精品50| 亚洲精品国产一区二区精华液| 日本一区二区免费在线视频| 午夜福利视频1000在线观看| 免费观看人在逋| 免费av毛片视频| 成人三级黄色视频| 麻豆久久精品国产亚洲av| 成人国产综合亚洲| 好看av亚洲va欧美ⅴa在| 国产在线精品亚洲第一网站| 午夜福利成人在线免费观看| 国产成人aa在线观看| 国产成人影院久久av| 一个人免费在线观看的高清视频| 免费高清视频大片| avwww免费| 美女大奶头视频| 日韩精品青青久久久久久| 日本成人三级电影网站| 欧美高清成人免费视频www| xxx96com| 丰满人妻熟妇乱又伦精品不卡| 欧美黄色淫秽网站| 99久久国产精品久久久| 黄片大片在线免费观看| 中亚洲国语对白在线视频| 亚洲性夜色夜夜综合| 日韩三级视频一区二区三区| 欧美黄色片欧美黄色片| 一级a爱片免费观看的视频| 此物有八面人人有两片| 人妻夜夜爽99麻豆av| 黄色成人免费大全| 国产伦在线观看视频一区| 久久久久久久久中文| 日日夜夜操网爽| 精品久久久久久久久久免费视频| 国产激情久久老熟女| 日本撒尿小便嘘嘘汇集6| 99精品久久久久人妻精品| 亚洲五月婷婷丁香| 中国美女看黄片| 狠狠狠狠99中文字幕| 国产精品精品国产色婷婷| 亚洲精品粉嫩美女一区| 久久99热这里只有精品18| 一级a爱片免费观看的视频| 国产av在哪里看| 久久精品人妻少妇| 一进一出抽搐gif免费好疼| 熟女电影av网| 国产精品免费视频内射| 天天一区二区日本电影三级| 一本综合久久免费| 欧美高清成人免费视频www| 国产精品野战在线观看| 黑人欧美特级aaaaaa片| 成人三级做爰电影| 久久久国产成人精品二区| 国产97色在线日韩免费| 中文字幕精品亚洲无线码一区| 欧美人与性动交α欧美精品济南到| 一边摸一边抽搐一进一小说| 久久精品91无色码中文字幕| 国产精品 国内视频| 欧美黑人巨大hd| 老鸭窝网址在线观看| 国产精品美女特级片免费视频播放器 | 亚洲第一欧美日韩一区二区三区| tocl精华| 99在线视频只有这里精品首页| 99在线视频只有这里精品首页| 亚洲色图 男人天堂 中文字幕| 青草久久国产| 99riav亚洲国产免费| 中国美女看黄片| 免费看十八禁软件| 男插女下体视频免费在线播放| 又大又爽又粗| 欧美成人一区二区免费高清观看 | 亚洲av日韩精品久久久久久密| 十八禁网站免费在线| 亚洲男人天堂网一区| 亚洲成人国产一区在线观看| 亚洲美女视频黄频| 国产精品av视频在线免费观看| 日韩精品免费视频一区二区三区| 巨乳人妻的诱惑在线观看| 国产精品久久久av美女十八| 巨乳人妻的诱惑在线观看| 在线观看一区二区三区| 国产主播在线观看一区二区| 99久久精品热视频| 日韩精品免费视频一区二区三区| 久久精品夜夜夜夜夜久久蜜豆 | 深夜精品福利| 男女床上黄色一级片免费看| 亚洲精品粉嫩美女一区| 午夜福利免费观看在线| 午夜福利免费观看在线| 欧美性猛交╳xxx乱大交人| 极品教师在线免费播放| 91av网站免费观看| 精品免费久久久久久久清纯| 51午夜福利影视在线观看| 一个人免费在线观看电影 | 国产亚洲欧美在线一区二区| 淫秽高清视频在线观看| 久久中文字幕人妻熟女| 一夜夜www| 亚洲国产欧美一区二区综合| 日韩大码丰满熟妇| 久久性视频一级片| 久久久久九九精品影院| 亚洲欧美日韩高清专用| 精品熟女少妇八av免费久了| 操出白浆在线播放| 非洲黑人性xxxx精品又粗又长| xxxwww97欧美| 成年女人毛片免费观看观看9| 色av中文字幕| www.精华液| 国产精品一区二区免费欧美| 最新美女视频免费是黄的| 国产99久久九九免费精品| 怎么达到女性高潮| 午夜日韩欧美国产| 又黄又粗又硬又大视频| 日本熟妇午夜| 国内精品一区二区在线观看| 一级黄色大片毛片| 精品国产亚洲在线| 国产精品影院久久| 久久人人精品亚洲av| 丁香欧美五月| 国产主播在线观看一区二区| 亚洲,欧美精品.| 久久久久久久久久黄片| 黄色视频不卡| 国产亚洲精品久久久久5区| 一级片免费观看大全| 两性夫妻黄色片| 18禁黄网站禁片午夜丰满| 桃色一区二区三区在线观看| 精品午夜福利视频在线观看一区| 亚洲欧洲精品一区二区精品久久久| 黄色a级毛片大全视频| 国产亚洲精品av在线| 久久久久国产精品人妻aⅴ院| 亚洲专区国产一区二区| 日韩中文字幕欧美一区二区| 日本免费一区二区三区高清不卡| 国产男靠女视频免费网站| av在线播放免费不卡| 美女黄网站色视频| 久久精品综合一区二区三区| 久久中文字幕人妻熟女| 亚洲精华国产精华精| 欧美色视频一区免费| 成年女人毛片免费观看观看9| 成人18禁在线播放| 99热这里只有是精品50| 91在线观看av| 精品少妇一区二区三区视频日本电影| 国产一级毛片七仙女欲春2| 亚洲国产高清在线一区二区三| 少妇的丰满在线观看| 神马国产精品三级电影在线观看 | 看免费av毛片| 亚洲一码二码三码区别大吗| 国产不卡一卡二| 久久久久久久久免费视频了| 99国产综合亚洲精品| 久久精品国产亚洲av香蕉五月| 制服人妻中文乱码| 午夜福利成人在线免费观看| 在线看三级毛片| av有码第一页| 欧美性长视频在线观看| 久久国产精品人妻蜜桃| 成年免费大片在线观看| 欧美在线一区亚洲| 美女黄网站色视频| 动漫黄色视频在线观看| 日韩欧美在线二视频| 欧美人与性动交α欧美精品济南到| 午夜免费激情av| 色精品久久人妻99蜜桃| 亚洲九九香蕉| 国产在线观看jvid| 国产蜜桃级精品一区二区三区| 岛国视频午夜一区免费看| 亚洲国产精品999在线| 国产激情欧美一区二区| 悠悠久久av| 男人舔奶头视频| 深夜精品福利| 久久午夜亚洲精品久久| 国产aⅴ精品一区二区三区波| 99久久99久久久精品蜜桃| 人妻夜夜爽99麻豆av| 精品乱码久久久久久99久播| 欧美一区二区国产精品久久精品 | 成熟少妇高潮喷水视频| 757午夜福利合集在线观看| 两个人的视频大全免费| 久久精品aⅴ一区二区三区四区| 国产精品久久久人人做人人爽| 国产成年人精品一区二区| 亚洲五月天丁香| 后天国语完整版免费观看| 动漫黄色视频在线观看| 精品少妇一区二区三区视频日本电影| 欧美三级亚洲精品| 日韩免费av在线播放| 亚洲 欧美 日韩 在线 免费| 日本免费一区二区三区高清不卡| 久久久久久人人人人人| 伊人久久大香线蕉亚洲五| 老熟妇仑乱视频hdxx| 日本成人三级电影网站| 全区人妻精品视频| 国产精品99久久99久久久不卡| 日本三级黄在线观看| 老司机午夜十八禁免费视频| 一边摸一边抽搐一进一小说| 久久亚洲真实| 1024视频免费在线观看| 亚洲av日韩精品久久久久久密| 久久99热这里只有精品18| 精品电影一区二区在线| 黄色 视频免费看| 夜夜看夜夜爽夜夜摸| 午夜激情av网站| 日日夜夜操网爽| 中出人妻视频一区二区| 国产欧美日韩一区二区精品| 亚洲av第一区精品v没综合| 欧美黑人欧美精品刺激| 欧美av亚洲av综合av国产av| 亚洲精品在线美女| 9191精品国产免费久久| 日本一区二区免费在线视频| 麻豆国产av国片精品| 久久久国产精品麻豆| 久久精品国产亚洲av香蕉五月| ponron亚洲| 老司机深夜福利视频在线观看| 成人18禁在线播放| av福利片在线| 国产精品亚洲av一区麻豆| 中文字幕高清在线视频| 麻豆一二三区av精品| 国产精品一区二区精品视频观看| 天天躁夜夜躁狠狠躁躁| 两个人视频免费观看高清| 欧美乱妇无乱码| 黑人欧美特级aaaaaa片| 一边摸一边做爽爽视频免费| 精品久久久久久,| 国产黄色小视频在线观看| 亚洲第一欧美日韩一区二区三区| 99国产精品一区二区三区| 亚洲精品一区av在线观看| 美女黄网站色视频| 日韩高清综合在线| 欧美中文日本在线观看视频| 免费在线观看成人毛片| 国产精品 国内视频| 精品日产1卡2卡| 亚洲色图 男人天堂 中文字幕| 久9热在线精品视频| 日本一区二区免费在线视频| 99精品在免费线老司机午夜| av福利片在线观看| 2021天堂中文幕一二区在线观| 视频区欧美日本亚洲| 久久天躁狠狠躁夜夜2o2o| 欧美日韩一级在线毛片| √禁漫天堂资源中文www| 亚洲欧美日韩东京热| 欧美绝顶高潮抽搐喷水| 777久久人妻少妇嫩草av网站| 亚洲国产高清在线一区二区三| 国产精品免费一区二区三区在线| 午夜老司机福利片| 一区二区三区国产精品乱码| 亚洲自拍偷在线| 啪啪无遮挡十八禁网站| 黑人欧美特级aaaaaa片| 人人妻人人看人人澡| 欧美日韩黄片免| 亚洲,欧美精品.| 日本一二三区视频观看| 免费在线观看成人毛片| 久久国产精品人妻蜜桃| 99久久99久久久精品蜜桃| 久久香蕉精品热| 午夜免费观看网址| 国产成人精品久久二区二区91| 狂野欧美白嫩少妇大欣赏| 午夜成年电影在线免费观看| 精品久久蜜臀av无| 两个人看的免费小视频| 韩国av一区二区三区四区| 天天躁狠狠躁夜夜躁狠狠躁| av有码第一页| 国产麻豆成人av免费视频| 午夜精品一区二区三区免费看| 91九色精品人成在线观看| 午夜两性在线视频| 黄色女人牲交| 听说在线观看完整版免费高清| 欧美人与性动交α欧美精品济南到| √禁漫天堂资源中文www| 黑人操中国人逼视频| 一区二区三区高清视频在线| 亚洲中文av在线| 两人在一起打扑克的视频| 女人爽到高潮嗷嗷叫在线视频| av中文乱码字幕在线| 久久久国产精品麻豆| 99精品久久久久人妻精品| 国产片内射在线| 国产精品98久久久久久宅男小说| av国产免费在线观看| 亚洲aⅴ乱码一区二区在线播放 | 久久国产精品人妻蜜桃| 妹子高潮喷水视频| 丰满人妻一区二区三区视频av | 麻豆成人午夜福利视频| 成人三级做爰电影| 91大片在线观看| 国产精品一区二区精品视频观看| avwww免费| 天堂√8在线中文| 欧美成人性av电影在线观看| 十八禁人妻一区二区| 欧美另类亚洲清纯唯美| 亚洲一区中文字幕在线| 嫩草影视91久久| 亚洲熟妇熟女久久| 久久久久久久精品吃奶| 亚洲一码二码三码区别大吗| 又爽又黄无遮挡网站| 中文字幕久久专区| 日本黄大片高清| 一区福利在线观看| 男女之事视频高清在线观看| 首页视频小说图片口味搜索| 久久国产精品影院| 好男人电影高清在线观看| 日韩成人在线观看一区二区三区| 岛国视频午夜一区免费看| 这个男人来自地球电影免费观看| 国产精品久久电影中文字幕| 亚洲一区二区三区不卡视频| 欧美+亚洲+日韩+国产| 国产一区二区激情短视频| 国产熟女午夜一区二区三区| 黄色a级毛片大全视频| 国产99白浆流出| 国产精品一区二区三区四区久久| 久久天躁狠狠躁夜夜2o2o| 老司机福利观看| 国产熟女xx| 嫩草影视91久久| 精品久久久久久久毛片微露脸| 免费看十八禁软件| √禁漫天堂资源中文www| 757午夜福利合集在线观看| 精华霜和精华液先用哪个| 欧美一区二区国产精品久久精品 | 99国产精品一区二区蜜桃av| 校园春色视频在线观看| 国产伦在线观看视频一区| 人人妻,人人澡人人爽秒播| 俄罗斯特黄特色一大片| 久久久久久亚洲精品国产蜜桃av| 国产亚洲精品一区二区www| 午夜影院日韩av| 国产精品精品国产色婷婷| 亚洲第一欧美日韩一区二区三区| 777久久人妻少妇嫩草av网站| 天堂动漫精品| 亚洲欧美激情综合另类| 777久久人妻少妇嫩草av网站| 日韩三级视频一区二区三区| 国产97色在线日韩免费| 欧美性长视频在线观看| 国产蜜桃级精品一区二区三区| 国产av又大| 欧美日韩国产亚洲二区| 日韩欧美一区二区三区在线观看| 亚洲人成电影免费在线| 男男h啪啪无遮挡| 午夜福利在线在线| 18禁美女被吸乳视频| 成人三级做爰电影| 女同久久另类99精品国产91| 91麻豆精品激情在线观看国产| 欧美日韩亚洲国产一区二区在线观看| 91成年电影在线观看| 日韩欧美在线二视频| 国产不卡一卡二| 精品午夜福利视频在线观看一区| 最好的美女福利视频网| 亚洲人成电影免费在线| 国产99白浆流出| 国产亚洲av嫩草精品影院| 熟女电影av网| 久久中文字幕人妻熟女| 欧美一级a爱片免费观看看 | 久99久视频精品免费| 国产精品免费视频内射| 久久久精品欧美日韩精品| 在线播放国产精品三级| 丁香欧美五月| 亚洲一卡2卡3卡4卡5卡精品中文| 全区人妻精品视频| 国产精品久久久久久人妻精品电影| 日韩欧美免费精品| www.熟女人妻精品国产| 2021天堂中文幕一二区在线观| 1024视频免费在线观看| 久久久久久久久中文| 波多野结衣高清作品| 国产成人欧美在线观看| 动漫黄色视频在线观看| 亚洲aⅴ乱码一区二区在线播放 | 特大巨黑吊av在线直播| 在线观看免费视频日本深夜| 午夜福利免费观看在线| 高清在线国产一区| 国产av在哪里看| 亚洲成av人片免费观看| 老司机深夜福利视频在线观看| 精品国产美女av久久久久小说| 亚洲人成电影免费在线| 手机成人av网站| 日韩欧美三级三区| 精品欧美一区二区三区在线| 日韩欧美 国产精品| 国产亚洲精品久久久久5区| 午夜精品一区二区三区免费看| avwww免费| 97超级碰碰碰精品色视频在线观看| 悠悠久久av| 国产成人影院久久av| 国内精品一区二区在线观看| 国产精品日韩av在线免费观看| 久久久久国产一级毛片高清牌| 曰老女人黄片| 一个人免费在线观看的高清视频| 国产又黄又爽又无遮挡在线| 日韩精品免费视频一区二区三区| 亚洲欧美激情综合另类| 国产视频一区二区在线看| 久久久久久国产a免费观看| 色尼玛亚洲综合影院| 亚洲国产看品久久| x7x7x7水蜜桃| 三级国产精品欧美在线观看 | 我要搜黄色片| 高潮久久久久久久久久久不卡| 中文字幕人成人乱码亚洲影| 亚洲成人久久爱视频| 日本成人三级电影网站| 精品免费久久久久久久清纯| 亚洲av中文字字幕乱码综合| 免费av毛片视频| 夜夜夜夜夜久久久久| ponron亚洲| 欧美成人一区二区免费高清观看 | 国产又黄又爽又无遮挡在线| 日韩 欧美 亚洲 中文字幕| 欧美日韩一级在线毛片| 亚洲成人免费电影在线观看| 岛国在线免费视频观看| 淫妇啪啪啪对白视频| 村上凉子中文字幕在线| 午夜福利在线观看吧| 精品高清国产在线一区| 国产精品久久久人人做人人爽| 亚洲免费av在线视频| 国产一级毛片七仙女欲春2| 手机成人av网站| 一二三四社区在线视频社区8| 国产激情久久老熟女| 大型黄色视频在线免费观看| a在线观看视频网站| 又黄又粗又硬又大视频| 五月伊人婷婷丁香| 成人手机av| 欧美乱色亚洲激情| 亚洲精品粉嫩美女一区| 特大巨黑吊av在线直播| 精品福利观看| 亚洲成人国产一区在线观看| 欧美色欧美亚洲另类二区| 女人爽到高潮嗷嗷叫在线视频| 亚洲精品一区av在线观看| 欧美成狂野欧美在线观看| 成人高潮视频无遮挡免费网站| 国产av不卡久久| 黄色丝袜av网址大全| 亚洲国产精品sss在线观看| 国产欧美日韩一区二区精品| 国产精品免费视频内射| 制服丝袜大香蕉在线| 日日摸夜夜添夜夜添小说| 婷婷六月久久综合丁香| 很黄的视频免费| 麻豆久久精品国产亚洲av| 男人的好看免费观看在线视频 | 久久亚洲真实| 别揉我奶头~嗯~啊~动态视频| 亚洲午夜精品一区,二区,三区| 亚洲一码二码三码区别大吗| 成年版毛片免费区| 欧美黑人巨大hd| 女同久久另类99精品国产91| 久久久国产成人免费| 99国产极品粉嫩在线观看| 亚洲欧美精品综合久久99| 国产在线观看jvid| 又爽又黄无遮挡网站| 亚洲精品在线观看二区| 国产精品美女特级片免费视频播放器 | 免费在线观看日本一区| 91麻豆精品激情在线观看国产| 精品久久久久久成人av| 色在线成人网| 日本撒尿小便嘘嘘汇集6| 亚洲精品在线美女| 国产亚洲av嫩草精品影院| 麻豆成人av在线观看| 一级片免费观看大全| 国产精品一区二区三区四区免费观看 | 看黄色毛片网站| 亚洲国产欧洲综合997久久,| 国产精品永久免费网站| 日韩精品免费视频一区二区三区| 成熟少妇高潮喷水视频| 丁香欧美五月| 亚洲一区高清亚洲精品| 久久久久久九九精品二区国产 | 欧美日韩国产亚洲二区| 亚洲人成77777在线视频| 麻豆一二三区av精品| 亚洲中文av在线| 成人精品一区二区免费| 亚洲av电影在线进入| 一本综合久久免费| 精品久久久久久久毛片微露脸| 看黄色毛片网站| 最近视频中文字幕2019在线8| 亚洲国产日韩欧美精品在线观看 | 国内精品久久久久精免费| 久久久久久久久久黄片| 亚洲第一电影网av| 97人妻精品一区二区三区麻豆| 在线观看免费午夜福利视频| 欧美三级亚洲精品| 18禁裸乳无遮挡免费网站照片| 午夜两性在线视频| 一本一本综合久久| 日本成人三级电影网站| 亚洲精品久久国产高清桃花| 全区人妻精品视频| 久久精品国产亚洲av高清一级| 少妇熟女aⅴ在线视频| 成人午夜高清在线视频| 激情在线观看视频在线高清| 夜夜爽天天搞| 国产精品亚洲一级av第二区| 麻豆av在线久日| 全区人妻精品视频| 日本一二三区视频观看| 色哟哟哟哟哟哟| 一个人观看的视频www高清免费观看 | 国产三级黄色录像| 人人妻人人看人人澡| 日本免费a在线|