• <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)控難題
    成年av动漫网址| 国产成人欧美| 日韩欧美精品免费久久| 99香蕉大伊视频| 婷婷色麻豆天堂久久| 欧美日韩一区二区视频在线观看视频在线| 在线精品无人区一区二区三| 天天影视国产精品| 蜜桃在线观看..| 精品福利永久在线观看| 亚洲成国产人片在线观看| 啦啦啦中文免费视频观看日本| 成年人午夜在线观看视频| 乱码一卡2卡4卡精品| 国产精品.久久久| 国产精品一国产av| 男女高潮啪啪啪动态图| 欧美少妇被猛烈插入视频| www.熟女人妻精品国产 | 天堂中文最新版在线下载| 久久这里只有精品19| 天堂8中文在线网| 少妇 在线观看| 成年女人在线观看亚洲视频| 国产熟女欧美一区二区| xxx大片免费视频| 性色avwww在线观看| 熟妇人妻不卡中文字幕| 国产男女超爽视频在线观看| 一本—道久久a久久精品蜜桃钙片| 热re99久久国产66热| 国产一级毛片在线| 欧美激情极品国产一区二区三区 | 伊人久久国产一区二区| 熟女电影av网| 我的女老师完整版在线观看| 女性被躁到高潮视频| 久久久国产欧美日韩av| 国产熟女欧美一区二区| 国产男女内射视频| 久久人人爽人人片av| 日韩一区二区视频免费看| 国产免费一级a男人的天堂| 日韩大片免费观看网站| 亚洲精品国产av成人精品| 老熟女久久久| 久热这里只有精品99| 国产精品麻豆人妻色哟哟久久| 久久久久网色| 精品国产国语对白av| 亚洲av电影在线观看一区二区三区| 大片电影免费在线观看免费| 美女国产视频在线观看| 日韩人妻精品一区2区三区| 少妇被粗大的猛进出69影院 | 亚洲av日韩在线播放| 99国产综合亚洲精品| 欧美最新免费一区二区三区| 欧美人与善性xxx| 亚洲精华国产精华液的使用体验| 伦精品一区二区三区| 国产免费福利视频在线观看| 国产精品三级大全| 激情视频va一区二区三区| 久久亚洲国产成人精品v| 久久99精品国语久久久| 卡戴珊不雅视频在线播放| 欧美成人精品欧美一级黄| 久久精品国产亚洲av涩爱| 又粗又硬又长又爽又黄的视频| 国产精品一区二区在线不卡| 国产亚洲一区二区精品| 看十八女毛片水多多多| 日韩三级伦理在线观看| 曰老女人黄片| 色5月婷婷丁香| 亚洲av中文av极速乱| 搡老乐熟女国产| 韩国av在线不卡| 日韩成人伦理影院| 少妇被粗大猛烈的视频| 嫩草影院入口| 在线观看人妻少妇| 男女下面插进去视频免费观看 | 日韩在线高清观看一区二区三区| 日产精品乱码卡一卡2卡三| 久久影院123| 婷婷色av中文字幕| 看免费成人av毛片| 少妇人妻久久综合中文| www日本在线高清视频| 国产色爽女视频免费观看| 制服丝袜香蕉在线| 中文精品一卡2卡3卡4更新| 久久久久网色| 99re6热这里在线精品视频| 亚洲色图 男人天堂 中文字幕 | 欧美变态另类bdsm刘玥| 美女主播在线视频| 国产亚洲精品久久久com| 午夜影院在线不卡| 汤姆久久久久久久影院中文字幕| 午夜免费鲁丝| 熟妇人妻不卡中文字幕| 边亲边吃奶的免费视频| 精品一区在线观看国产| 色5月婷婷丁香| 边亲边吃奶的免费视频| av国产精品久久久久影院| 国产无遮挡羞羞视频在线观看| 免费大片黄手机在线观看| 亚洲欧美成人综合另类久久久| av国产久精品久网站免费入址| 丝袜喷水一区| 蜜臀久久99精品久久宅男| 日本av手机在线免费观看| 国产精品久久久久久久久免| 在线天堂最新版资源| 国产男人的电影天堂91| 九色亚洲精品在线播放| 黄色配什么色好看| av又黄又爽大尺度在线免费看| av女优亚洲男人天堂| 午夜免费男女啪啪视频观看| 久久久久精品久久久久真实原创| 亚洲精品国产av成人精品| 久久精品熟女亚洲av麻豆精品| 91午夜精品亚洲一区二区三区| 黄片播放在线免费| tube8黄色片| 亚洲一码二码三码区别大吗| 国产精品女同一区二区软件| 新久久久久国产一级毛片| 国产精品久久久久久久电影| 精品久久久久久电影网| 最近中文字幕高清免费大全6| 国产精品久久久久久久久免| 在线天堂最新版资源| 搡女人真爽免费视频火全软件| 欧美精品av麻豆av| 中文字幕制服av| 成年动漫av网址| 永久免费av网站大全| 狂野欧美激情性xxxx在线观看| 亚洲精品中文字幕在线视频| 欧美日韩一区二区视频在线观看视频在线| 久久热在线av| 少妇精品久久久久久久| 大香蕉久久成人网| 人妻人人澡人人爽人人| 亚洲国产精品999| 久久精品国产自在天天线| 99九九在线精品视频| 国产色婷婷99| 中文欧美无线码| 一二三四中文在线观看免费高清| 免费人成在线观看视频色| 蜜臀久久99精品久久宅男| 久久韩国三级中文字幕| 日韩电影二区| 看非洲黑人一级黄片| 亚洲欧美清纯卡通| 人妻一区二区av| 亚洲国产av新网站| 亚洲五月色婷婷综合| 久久鲁丝午夜福利片| 亚洲成国产人片在线观看| 免费黄色在线免费观看| 国产一区二区三区综合在线观看 | 国产av一区二区精品久久| 亚洲熟女精品中文字幕| 在现免费观看毛片| 日韩熟女老妇一区二区性免费视频| 精品久久蜜臀av无| 香蕉精品网在线| 插逼视频在线观看| 黑人猛操日本美女一级片| 女性被躁到高潮视频| 亚洲国产精品国产精品| 老女人水多毛片| 秋霞伦理黄片| 精品国产一区二区三区久久久樱花| 日本与韩国留学比较| 国产亚洲欧美精品永久| 亚洲精品久久午夜乱码| 老司机影院成人| 女的被弄到高潮叫床怎么办| 日韩精品免费视频一区二区三区 | 一区二区三区精品91| 中文精品一卡2卡3卡4更新| 91精品国产国语对白视频| 黄网站色视频无遮挡免费观看| 国产精品麻豆人妻色哟哟久久| 免费女性裸体啪啪无遮挡网站| 国国产精品蜜臀av免费| 伊人久久国产一区二区| 亚洲色图 男人天堂 中文字幕 | 黄色配什么色好看| 日本猛色少妇xxxxx猛交久久| 国产精品久久久久久久久免| 久久99一区二区三区| 欧美精品av麻豆av| 天天影视国产精品| 亚洲国产av新网站| 99热这里只有是精品在线观看| 人体艺术视频欧美日本| 国产综合精华液| 少妇精品久久久久久久| 亚洲欧洲国产日韩| 韩国av在线不卡| 国产男女超爽视频在线观看| 亚洲精品久久午夜乱码| 少妇的逼水好多| 免费在线观看黄色视频的| 亚洲国产色片| 欧美激情极品国产一区二区三区 | 人人妻人人爽人人添夜夜欢视频| 日韩三级伦理在线观看| 亚洲图色成人| 国产黄频视频在线观看| 中文乱码字字幕精品一区二区三区| 欧美日本中文国产一区发布| 国产激情久久老熟女| 女人被躁到高潮嗷嗷叫费观| 在线看a的网站| 日本av手机在线免费观看| 少妇精品久久久久久久| 精品国产一区二区久久| 啦啦啦啦在线视频资源| 日本免费在线观看一区| 综合色丁香网| 精品国产一区二区三区四区第35| 亚洲欧美精品自产自拍| 久久人人爽人人爽人人片va| 女性被躁到高潮视频| 在线观看人妻少妇| 久久精品人人爽人人爽视色| 久久免费观看电影| 男女高潮啪啪啪动态图| 9热在线视频观看99| 男女午夜视频在线观看 | 久久精品久久精品一区二区三区| 美女国产高潮福利片在线看| 亚洲第一av免费看| 九九在线视频观看精品| 久久这里只有精品19| 夜夜爽夜夜爽视频| 欧美+日韩+精品| 伊人亚洲综合成人网| 有码 亚洲区| 日韩制服丝袜自拍偷拍| 波多野结衣一区麻豆| 国产精品久久久久成人av| 中国美白少妇内射xxxbb| 日韩大片免费观看网站| 欧美亚洲日本最大视频资源| 亚洲av国产av综合av卡| 美女国产高潮福利片在线看| 亚洲第一区二区三区不卡| 国产精品女同一区二区软件| 在线观看免费视频网站a站| 69精品国产乱码久久久| 热re99久久精品国产66热6| 午夜福利在线观看免费完整高清在| 久久久久久久精品精品| 欧美激情 高清一区二区三区| 精品久久久久久电影网| 国产爽快片一区二区三区| 日产精品乱码卡一卡2卡三| 国产精品熟女久久久久浪| 性高湖久久久久久久久免费观看| 亚洲综合色网址| 亚洲精品日本国产第一区| 日韩一区二区视频免费看| 久久99精品国语久久久| 狂野欧美激情性xxxx在线观看| 在线免费观看不下载黄p国产| 免费黄网站久久成人精品| 一级爰片在线观看| 日韩人妻精品一区2区三区| 国产免费一区二区三区四区乱码| www.av在线官网国产| 亚洲国产精品专区欧美| 午夜激情av网站| 亚洲国产日韩一区二区| 99久久中文字幕三级久久日本| 夜夜骑夜夜射夜夜干| 欧美激情极品国产一区二区三区 | 十八禁网站网址无遮挡| 国产高清国产精品国产三级| 久久久a久久爽久久v久久| 少妇被粗大猛烈的视频| 亚洲精品中文字幕在线视频| 国产精品成人在线| 精品一区二区三区视频在线| 在线天堂最新版资源| 久久国产亚洲av麻豆专区| 日韩欧美一区视频在线观看| 在线观看三级黄色| 精品亚洲乱码少妇综合久久| 一区在线观看完整版| 亚洲精品成人av观看孕妇| 国产成人精品久久久久久| 久久韩国三级中文字幕| 伦理电影大哥的女人| 青青草视频在线视频观看| 两个人免费观看高清视频| 少妇的丰满在线观看| 99久久人妻综合| 欧美日韩视频精品一区| 97精品久久久久久久久久精品| 大话2 男鬼变身卡| 女性被躁到高潮视频| 美女脱内裤让男人舔精品视频| 日日摸夜夜添夜夜爱| 久久亚洲国产成人精品v| 国产精品三级大全| 久久这里有精品视频免费| 亚洲精品色激情综合| 自线自在国产av| 国产精品国产三级专区第一集| 久久久久精品性色| 母亲3免费完整高清在线观看 | 最近2019中文字幕mv第一页| 日韩熟女老妇一区二区性免费视频| 亚洲av.av天堂| 2021少妇久久久久久久久久久| 亚洲欧美色中文字幕在线| 最近中文字幕2019免费版| 黑人高潮一二区| 97精品久久久久久久久久精品| 不卡视频在线观看欧美| av免费在线看不卡| 亚洲五月色婷婷综合| videosex国产| 一级片'在线观看视频| 久久久精品区二区三区| 女人久久www免费人成看片| 国产精品人妻久久久影院| 天堂俺去俺来也www色官网| 人人妻人人澡人人爽人人夜夜| av又黄又爽大尺度在线免费看| 国产精品国产三级国产专区5o| 国产成人精品福利久久| 国产成人精品婷婷| 国产精品99久久99久久久不卡 | 亚洲精品一二三| 黑人猛操日本美女一级片| 男女无遮挡免费网站观看| 黑人高潮一二区| 91成人精品电影| 亚洲国产av影院在线观看| 午夜精品国产一区二区电影| 视频区图区小说| 在线亚洲精品国产二区图片欧美| 视频区图区小说| 如何舔出高潮| 欧美日韩精品成人综合77777| 热re99久久精品国产66热6| 欧美亚洲 丝袜 人妻 在线| 水蜜桃什么品种好| 久久久亚洲精品成人影院| 国产精品久久久久成人av| 国产男女超爽视频在线观看| 一本久久精品| 亚洲欧美一区二区三区黑人 | 精品亚洲成a人片在线观看| 91aial.com中文字幕在线观看| 久久99热6这里只有精品| 2022亚洲国产成人精品| 97精品久久久久久久久久精品| 国产精品 国内视频| 中文乱码字字幕精品一区二区三区| 寂寞人妻少妇视频99o| av不卡在线播放| 精品国产一区二区三区久久久樱花| 精品卡一卡二卡四卡免费| 久久久国产精品麻豆| 男男h啪啪无遮挡| 性高湖久久久久久久久免费观看| 国产av国产精品国产| 少妇被粗大的猛进出69影院 | 亚洲国产毛片av蜜桃av| 日韩中字成人| 亚洲国产成人一精品久久久| 久久久久精品久久久久真实原创| 纯流量卡能插随身wifi吗| 99热国产这里只有精品6| 天天躁夜夜躁狠狠躁躁| 久久久久久伊人网av| 男女午夜视频在线观看 | 国产不卡av网站在线观看| 久久久久人妻精品一区果冻| 精品少妇黑人巨大在线播放| 亚洲精品久久午夜乱码| 色94色欧美一区二区| 精品卡一卡二卡四卡免费| 欧美日韩综合久久久久久| 久久午夜福利片| 美女脱内裤让男人舔精品视频| 亚洲综合精品二区| 男女免费视频国产| 欧美xxxx性猛交bbbb| 久久久亚洲精品成人影院| 一级黄片播放器| 色94色欧美一区二区| 久久精品国产综合久久久 | 水蜜桃什么品种好| 久久99蜜桃精品久久| 黄色 视频免费看| 欧美bdsm另类| 伊人久久国产一区二区| 十八禁网站网址无遮挡| 人人妻人人添人人爽欧美一区卜| 欧美激情极品国产一区二区三区 | 亚洲av欧美aⅴ国产| 国产黄色视频一区二区在线观看| 18在线观看网站| 在线免费观看不下载黄p国产| 丝袜人妻中文字幕| 亚洲成人一二三区av| 日韩中字成人| 国产亚洲精品第一综合不卡 | 亚洲国产精品999| 日本av手机在线免费观看| 五月开心婷婷网| 王馨瑶露胸无遮挡在线观看| 在线看a的网站| 亚洲人与动物交配视频| 在线观看免费日韩欧美大片| 欧美成人精品欧美一级黄| 国产精品秋霞免费鲁丝片| 亚洲国产精品国产精品| 大片电影免费在线观看免费| 日韩成人伦理影院| 国产1区2区3区精品| 人人妻人人澡人人爽人人夜夜| 插逼视频在线观看| av卡一久久| 国产av精品麻豆| 免费日韩欧美在线观看| 男女边摸边吃奶| 欧美丝袜亚洲另类| 国产一区二区在线观看日韩| 亚洲伊人色综图| 最近最新中文字幕免费大全7| 亚洲欧美色中文字幕在线| 男女边吃奶边做爰视频| 91精品三级在线观看| 国产成人a∨麻豆精品| 熟女电影av网| 精品国产国语对白av| 国产在线免费精品| 久久精品久久精品一区二区三区| 蜜桃在线观看..| 国产精品国产三级国产av玫瑰| 晚上一个人看的免费电影| 大香蕉久久成人网| 国产成人aa在线观看| av.在线天堂| 午夜91福利影院| 少妇高潮的动态图| 天堂8中文在线网| 亚洲精品一二三| 免费av中文字幕在线| 久久久久久久亚洲中文字幕| 日韩不卡一区二区三区视频在线| 国产69精品久久久久777片| 亚洲av欧美aⅴ国产| av国产久精品久网站免费入址| av有码第一页| 看免费av毛片| 熟女av电影| 男女高潮啪啪啪动态图| 国产毛片在线视频| 人人澡人人妻人| 成人免费观看视频高清| 亚洲欧美清纯卡通| 五月开心婷婷网| 久久人人爽人人爽人人片va| 久久99热这里只频精品6学生| 一级爰片在线观看| 一级,二级,三级黄色视频| 伦理电影免费视频| 老女人水多毛片| 国产一区二区在线观看av| 曰老女人黄片| 美女脱内裤让男人舔精品视频| 久久综合国产亚洲精品| 日本欧美视频一区| 草草在线视频免费看| 亚洲成人一二三区av| 国产黄色免费在线视频| 亚洲四区av| 国产成人精品福利久久| 亚洲国产精品一区三区| 久久精品aⅴ一区二区三区四区 | 久久久久久人人人人人| 汤姆久久久久久久影院中文字幕| 51国产日韩欧美| 中文天堂在线官网| 国产片特级美女逼逼视频| 日本黄大片高清| 丰满少妇做爰视频| 久久久国产精品麻豆| 少妇猛男粗大的猛烈进出视频| 久久久久国产网址| 黄片播放在线免费| 在线观看一区二区三区激情| 香蕉丝袜av| 久久久久久久精品精品| 亚洲欧美成人精品一区二区| 99久久中文字幕三级久久日本| 欧美日韩亚洲高清精品| 街头女战士在线观看网站| 在线免费观看不下载黄p国产| 国产av码专区亚洲av| 青青草视频在线视频观看| 亚洲国产精品999| 午夜久久久在线观看| 欧美97在线视频| 又粗又硬又长又爽又黄的视频| 久久狼人影院| 日本猛色少妇xxxxx猛交久久| 亚洲丝袜综合中文字幕| 啦啦啦在线观看免费高清www| 精品一区二区三卡| 欧美精品国产亚洲| 精品一品国产午夜福利视频| 在线观看三级黄色| 成人综合一区亚洲| 日本猛色少妇xxxxx猛交久久| 日韩欧美一区视频在线观看| 欧美精品国产亚洲| 日韩精品免费视频一区二区三区 | 成人无遮挡网站| 99视频精品全部免费 在线| 男男h啪啪无遮挡| 亚洲在久久综合| 国精品久久久久久国模美| 精品少妇黑人巨大在线播放| 多毛熟女@视频| 在线观看免费视频网站a站| 久久人人爽人人爽人人片va| 在线观看人妻少妇| av免费在线看不卡| 久久人人爽人人片av| 午夜日本视频在线| 国产成人aa在线观看| 亚洲av欧美aⅴ国产| 我要看黄色一级片免费的| 国产一区二区激情短视频 | 伊人久久国产一区二区| 22中文网久久字幕| 中文字幕精品免费在线观看视频 | 大片免费播放器 马上看| videosex国产| 在线观看免费高清a一片| 国产成人午夜福利电影在线观看| 免费不卡的大黄色大毛片视频在线观看| 一级爰片在线观看| 亚洲综合色惰| 在线观看www视频免费| 日本免费在线观看一区| 久久韩国三级中文字幕| 91精品三级在线观看| 亚洲情色 制服丝袜| 亚洲精品一二三| 久久午夜综合久久蜜桃| 久久亚洲国产成人精品v| 久久热在线av| 欧美老熟妇乱子伦牲交| 中国国产av一级| 午夜免费男女啪啪视频观看| 搡老乐熟女国产| 人人妻人人爽人人添夜夜欢视频| xxxhd国产人妻xxx| 免费av中文字幕在线| 日韩成人伦理影院| 久久精品久久精品一区二区三区| 国产一区二区激情短视频 | 欧美日韩视频精品一区| 免费黄网站久久成人精品| av免费在线看不卡| 精品99又大又爽又粗少妇毛片| 日日摸夜夜添夜夜爱| 亚洲丝袜综合中文字幕| 亚洲国产精品专区欧美| 久久久久久久大尺度免费视频| 波多野结衣一区麻豆| 97在线视频观看| 国产成人精品一,二区| av又黄又爽大尺度在线免费看| 久久国产精品男人的天堂亚洲 | 国产亚洲精品久久久com| 高清黄色对白视频在线免费看| 日韩 亚洲 欧美在线| av播播在线观看一区| 久久久久久久国产电影| 草草在线视频免费看| 国产精品久久久av美女十八| 多毛熟女@视频| 欧美+日韩+精品| 毛片一级片免费看久久久久| 99视频精品全部免费 在线| 色婷婷久久久亚洲欧美| 交换朋友夫妻互换小说| 欧美xxxx性猛交bbbb| 亚洲精品久久久久久婷婷小说| 在线免费观看不下载黄p国产| 国产精品偷伦视频观看了| 欧美 日韩 精品 国产| 日韩精品有码人妻一区| 亚洲国产精品国产精品| 777米奇影视久久|