基于逻辑回归的银行对公贷款违约判别方法

摘要

本发明提供了基于逻辑回归的银行对公贷款违约判别方法，属于计算机技术领域。本发明采用灰色关联法计算各行业的违约倾向系数；采用单因素分析提取对被解释变量影响最大的几个解释变量，将从单因素分析中提取出来的原始变量，取方差累积贡献率L达到预设阈值的前m个指标为违约预测指标，每指标都是各原始变量的线性组合，各指标之间互不相关，通过二项逻辑回归算法建立二项逻辑回归模型进行违约概率预测，实现了对存在潜在违约风险的金融客户提前预警识别。

权利要求

1.基于逻辑回归的银行对公贷款违约判别方法，其特征在于，包括如下步骤：违约预测指标预提取步骤，采用单因素方差分析法或多重共线性验证方法，根据数据库中的企业基本信息和企业财务指标，提取对违约率影响最大的N个指标；

所述企业基本信息包括：企业名称、数据年度、是否违约、成立日期、在职员工数、总资产和借款人组织机构代码；

所述企业财务指标：主营业务收入、流动比率、速动比率、现金比率、资产负债率、负债权益比率、利息保障倍数、存货周转率、应收账款周转率、营业周期、流动资产周转率、总资产周转率、固定资产周转率、总资产报酬率、成本费用利润率、销售净利率、净资产收益率、销售毛利率、盈余现金保障倍数、现金流动负债比、营业总成本、现金流动负债比、管理费用、销售费用、财务费用、无形资产、长期借款、短期借款、净利润和营业收入；

违约预测指标确定步骤，

根据预提取的N个指标，取方差累积贡献率达到预设阈值的指标作为违约预测的指标；

具体包括如下步骤：

由预提取的N个指标组成指标变量矩阵；

计算指标变量矩阵的协方差，得到指标变量协方差矩阵Σ；

计算指标变量协方差矩阵Σ的特征向量λ＝(λ1,λ2,Λ,λN)；

计算指标变量协方差矩阵Σ的特征向量λ中各特征值λi对应的特征向量γi；

计算特征值λi对应的违约预测指标Yi，采用如下公式计算：Yi＝γ′iX；

其中，

Yi为特征值λi对应的违约预测指标；

X为待判别公司的指标变量矩阵；

按N个预提取指标的违约率影响的由大到小的顺序，依次计算N个预提取指标的方差累加贡献率；

取方差累积贡献率L达到预设阈值的前m个指标为违约预测指标，各指标的累积贡献率按如下公式计算：

其中，

N为预提取的指标数；

m为指标方差累积贡献率L达到预设阈值时的指标数；

λi为指标变量协方差矩阵Σ的特征向量λ中第i个指标的特征值；

二项逻辑回归模型建立步骤，将违约与否作为二项逻辑回归中的被解释变量，建立二项逻辑回归模型；

违约判断步骤，

根据二项逻辑回归模型及企业数据，得到违约概率；

根据违约概率和预设的违约概率阈值，判断该企业是否违约。

2.根据权利要求1所述的基于逻辑回归的银行对公贷款违约判别方法，其特征在于，采用单因素方差分析法提取对违约率影响最大的N个指标具体包括如下步骤：每个指标单独与违约与否进行单因素方差分析，取设定年份内的数据，根据如下公式得到检验统计量F的值：

其中，

k为设定年份数；

n为数据库中该企业同行业内所有公司的数量；

为第i年数据库中该企业同行业内各公司对应该指标变量的均值；

为设定年份内数据库中该企业同行业内各公司对应该指标变量的总均值；

ni为第i年的样本量；

xij第i年数据库中该企业同行业内第j个公司对应该指标变量的取值；

根据得到的检验统计量F的数值，查询F分布表，得到检验p值；

选择N个p值最小的指标，作为对违约率影响最大的N个指标。

3.根据权利要求1所述的基于逻辑回归的银行对公贷款违约判别方法，其特征在于，通过多重共线性验证方法提取对违约率影响最大的N个指标具体包括如下步骤：对各待检验指标进行多重共线性检验，检验公式为：c1X1+c2X2+ΛciXiΛ+cnXn＝c0，若存在c0,c1,c2,K,cn使得该公式成立，则待检验指标存在多重共线性；

其中，

Xi为第i个待检验指标；

ci为第i个待检验指标的系数；

c0为预设常数；

计算各待检验指标的拟合优度：其中，

Rj为第j个待检验指标的拟合优度；

Yj为第j个待检验指标的真实值；

j为解释变量的序号；

为第j个待检验指标以其余指标为解释变量进行回归时得到的回归模型估计量；

为待检验指标真实值的均值；

对各待检验指标的拟合优度进行F检验：其中，

为第j.个解释变量对其他解释变量的回归方程的判定系数；

k为解释变量的个数；

n为样本量；

j为解释变量的序号；

j.个用于F检验的解释变量的序号；

Fj为第j个待检验指标经过F检验后的拟合优度；

F(k‑1,n‑k)是指自由度从k‑1到n‑k的F分布；

根据对各指标的拟合优度进行F检验得到的检验结果，查询F分布表，得到检验p值；

选择N个p值最大的指标，作为对违约率影响最大的N个指标。

4.根据权利要求1所述的基于逻辑回归的银行对公贷款违约判别方法，其特征在于，所述二项逻辑回归模型为：

其中：

pi表示待判别公司在第i年的违约概率，取值在0‑1之间，违约为1，非违约为0；

Xi为待判别公司在第i年的指标变量矩阵；

xki表示待判别公司公司第i年在第k个指标上的取值；

β0为回归常数；

βk为回归模型中第k个指标的回归系数。

5.根据权利要求1所述的基于逻辑回归的银行对公贷款违约判别方法，其特征在于，违约概率判断步骤中，还包括在得到违约概率后，结合待判别公司所属行业的违约倾向系数，对违约概率采用如下公式进行修正：ci＝roi×pi

其中：

ci为企业i某一段时期内根据其所属行业的违约倾向系数修正后的违约概率；

roi为行业i某一段时期内的违约倾向系数；

pi为企业i某一段时期内的违约概率。

6.根据权利要求5所述的基于逻辑回归的银行对公贷款违约判别方法，其特征在于，待判别公司所属行业的违约倾向系数采用灰度关联分析法确定。

7.根据权利要求6所述的基于逻辑回归的银行对公贷款违约判别方法，其特征在于，确定行业的违约倾向系数的方法包括如下步骤：选取数据库中所有行业的公司数据的指标变量，组成指标变量矩阵，其中选取的指标变量包括所述企业基本信息和企业财务指标；

分别计算各行业包括的各公司指标变量矩阵与参考序列对应指标取值的灰度关联系数，所述参考序列由该行业中各指标的均值数据组成；

分别计算各行业包括的各公司的各指标变量与参考序列对应元素的灰关联度；

提取各行业包括的各公司的灰关联度，取各行业包括的各公司灰关联度的均值为该行业的违约倾向系数。

8.根据权利要求7所述的基于逻辑回归的银行对公贷款违约判别方法，其特征在于，各公司各个指标变量与参考序列对应元素的灰关联度的计算包括如下步骤：对各公司各个指标变量矩阵进行灰色无量纲化处理，其中，

Xi(k)为第i个公司指标变量矩阵的对应第k个指标变量的取值；

Xi(1)为第i个公司指标变量矩阵的对应第1个指标变量的取值；

计算差异信息序列Δoi(k)，Δoi(k)＝|x0(k)‑xi(k)|其中，

x0(k)为参考序列进行灰色无量纲化处理过后的对应第k个指标变量的取值，所述参考序列由该行业中各指标的均值数据组成；

xi(k)为第i个公司指标变量矩阵灰色无量纲化后对应第k个指标变量的取值；

确定环境参数：

其中，

Δ0i(max)为差异信息序列Δoi(k)中的最大值；

Δ0i(min)为差异信息序列Δoi(k)中的最小值；

分别采用如下公式计算各公司各指标变量与参考序列对应元素的灰关联系数ξi(k)，其中，

ρ为分辨系数，0<ρ<1；

Δ0i(max)为差异信息序列Δoi(k)中的最大值；

Δ0i(min)为差异信息序列Δoi(k)中的最小值；

Δoi(k)为差异信息序列；

分别采用如下公式计算各公司各指标变量与参考序列对应元素的灰关联度r0i，其中，

ξi(k)为第i个公司的第k个指标与参考序列对应元素的灰关联系数；

m为企业数量信息表和所有行业违约数量信息表中包括的指标总数。

说明书

基于逻辑回归的银行对公贷款违约判别方法

技术领域

[0001] 本发明涉及计算机技术领域，尤其涉及一种基于逻辑回归的银行对公贷款违约判别方法。

背景技术

[0002] 随着近年来金融经济的高速发展，贷款成为企业筹资的一大重要支撑手段，不可避免的企业违约现象的发生与之而来，而且越来越频繁的发生。银行业务结构和经营环境
的变化以及金融创新的迅速发展，促使越来越多的银行意识到需要构建更加全面和系统的
信用风险管理系统来应对金融风险日趋复杂化的现实，所以贷款企业违约风险测度具有极
大的现实意义。

[0003] 目前国外学术界和金融界对违约概率的测度与评估的研究概括起来有：(1)巴塞尔委员会新资本协议与国际知名金融机构推出的高级信用风险模型，主要针对企业违约概
率数值如何测算问题的研究，这些模型也是巴塞尔委员会推荐的在资本金及经济资本方面
计算的依据。(2)学术界对违约率的评估研究，集中于影响违约率的关键变量探寻，并利用
企业历史样本建立企业违约率分类判别模型来估计企业违约的可能性。这两个方法相互补
充，相互促进，但也存在过于依靠人工经验等弊端。

[0004] 在中国专利申请文献CN111192140A中，公开了一种客户违约概率预测的方法，根据待识别客户标识从预先保存的数据集中选取对应的指标数据；指标数据包括以下至少之
一：企业基本信息、企业信贷数据、企业财务数据、企业工商数据、企业社保数据；将指标数
据作为预先训练的违约概率预测模型的输入，预测待识别客户的违约概率。该方法还包括：
利用模型训练算法对样本数据进行模型训练得到违约概率预测模型，包括对样本训练集中
的样本数据进行预处理；利用特征相关性分析方法对预处理后的样本数据进行违约变量的
相关性分析得到特征候选集；利用模型训练算法对特征候选集进行模型训练得到违约概率
预测模型。根据违约业务定义对样本数据构造违约特征；利用特征相关性分析方法分析各
个违约特征与违约变量的相关性程度，选择相关性高的违约特征；将选择的违约特征作为
特征候选集。利用Xgboost算法对特征候选集进行模型训练得到违约概率预测模型。在得到
违约概率预测模型之后，根据所述违约概率预测模型对预先准备好的样本测试集中的样本
数据进行违约概率预测；利用ROC曲线、AUC或KS值对违约概率预测模型进行评估。该方法在
预测违约概率时，没有考虑行业平均违约率，也未对行业平均违约率进行分布和结构的分
析，特征分析选择与违约变量相关性高的违约特征进行分析，使得预测准确率、实时性、有
效程度等无法满足实际工作需要。

[0005] 现有技术至少存在以下不足：

[0006] 1.对公贷款违约存在样本量小、考虑周期短和较少考虑地区因素等缺陷，使得违约概率预测不准，算法的不精确导致业务人员工作效率降低。

[0007] 2.所使用的算法单一，无法满足业务部门实时、动态进行风险排查的需求。

[0008] 3.模型与算法阈值的设定不合理，频繁预警导致风险排查频率过高，加重了金融行业客户在还款期间带来应对排查的负担。

发明内容

[0009] 为解决现有技术中存在的问题，本发明提供了一种基于逻辑回归的银行对公贷款违约判别方法，采用灰色关联法计算各行业的违约倾向系数；采用单因素分析提取对被解
释变量影响最大的几个解释变量，将从单因素分析中提取出来的原始变量，取方差累积贡
献率L达到预设阈值的前m个指标为违约预测指标，每指标都是各原始变量的线性组合，各
指标之间互不相关，并通过财务指标、工商登记信息指标、行业贷款计划指标等信息根据二
项逻辑回归算法建立二项逻辑回归模型进行违约概率预测，实现了对金融客户还款期间所
蕴含的经营风险、还款能力进行实时、动态分析和预警。

[0010] 本发明提供了一种基于逻辑回归的银行对公贷款违约判别方法，包括如下步骤：

[0011] 违约预测指标预提取步骤，

[0012] 采用单因素方差分析法或多重共线性验证方法，根据数据库中的企业基本信息和企业财务指标，提取对违约率影响最大的N个指标；

[0013] 所述企业基本信息包括：企业名称、数据年度、是否违约、成立日期、在职员工数、总资产和借款人组织机构代码；是否违约、成立日期和在职员工数参与到指标变量矩阵中；

[0014] 所述企业财务指标：主营业务收入、流动比率、速动比率、现金比率、资产负债率、负债权益比率、利息保障倍数、存货周转率、应收账款周转率、营业周期、流动资产周转率、
总资产周转率、固定资产周转率、总资产报酬率、成本费用利润率、销售净利率、净资产收益
率、销售毛利率、盈余现金保障倍数、现金流动负债比、营业总成本、现金流动负债比、管理
费用、销售费用、财务费用、无形资产、长期借款、短期借款、净利润和营业收入；

[0015] 违约预测指标确定步骤，

[0016] 根据预提取的N个指标，取方差累积贡献率达到预设阈值的指标作为违约预测的指标；

[0017] 二项逻辑回归模型建立步骤，

[0018] 将违约与否作为二项逻辑回归中的被解释变量，建立二项逻辑回归模型；

[0019] 违约判断步骤，

[0020] 根据二项逻辑回归模型及企业数据，得到违约概率；

[0021] 根据违约概率和预设的违约概率阈值，判断该企业是否违约。

[0022] 优选地，采用单因素方差分析法提取对违约率影响最大的N个指标具体包括如下步骤：

[0023] 每个指标单独与违约与否进行单因素方差分析，取设定年份内的数据，根据如下公式得到检验统计量F的值：

[0024]

[0025] 其中，

[0026] k为设定年份数；

[0027] n为数据库中该企业同行业内所有公司的数量；

[0028] 为第i年数据库中该企业同行业内各公司对应该指标变量的均值；

[0029] 为设定年份内数据库中该企业同行业内各公司对应该指标变量的总均值；

[0030] ni为第i年的样本量；

[0031] xij第i年数据库中该企业同行业内第j个公司对应该指标变量的取值；

[0032] 根据得到的检验统计量F的数值，查询F分布表，得到检验p值；

[0033] 选择N个p值最小的指标，作为对违约率影响最大的N个指标。

[0034] 优选地，通过多重共线性验证方法提取对违约率影响最大的N个指标具体包括如下步骤：

[0035] 对各待检验指标进行多重共线性检验，检验公式为：c1X1+c2X2+…ciXi…+cnXn＝c0，若存在c0,c1,c2,…,cn使得该公式成立，则待检验指标存在多重共线性；

[0036] 其中，

[0037] Xi为第i个待检验指标；

[0038] ci为第i个待检验指标的系数；

[0039] c0为预设常数；

[0040] 计算各待检验指标的拟合优度：

[0041]

[0042] 其中，

[0043] Rj为第j个待检验指标的拟合优度；

[0044] Yj为第j个待检验指标的真实值；

[0045] 为第j个待检验指标以其余指标为解释变量进行回归时得到的回归模型估计量；

[0046] 为待检验指标真实值的均值；

[0047] 对各待检验指标的拟合优度进行F检验：

[0048]

[0049] 其中，

[0050] 为第j.个解释变量对其他解释变量的回归方程的判定系数；

[0051] k为解释变量的个数；

[0052] n为样本量；

[0053] j为解释变量的序号；

[0054] j.个用于F检验的解释变量的序号；

[0055] Fj为第j个待检验指标经过F检验后的拟合优度；

[0056] F(k‑1,n‑k)是指自由度从k‑1到n‑k的F分布；

[0057] 根据对各指标的拟合优度进行F检验得到的检验结果，查询F分布表，得到检验p值；

[0058] 选择N个p值最大的指标，作为对违约率影响最大的N个指标。

[0059] 优选地，违约预测指标确定步骤具体包括如下步骤：

[0060] 由预提取的N个指标组成指标变量矩阵；

[0061] 计算指标变量矩阵的协方差，得到指标变量协方差矩阵Σ；

[0062] 计算指标变量协方差矩阵Σ的特征向量λ＝(λ1,λ2,…,λN)；

[0063] 计算指标变量协方差矩阵Σ的特征向量λ中各特征值λi对应的特征向量γi；

[0064] 计算特征值λi对应的违约预测指标Yi，采用如下公式计算：

[0065] Yi＝γ′iX；

[0066] 其中，

[0067] Yi为特征值λi对应的违约预测指标；

[0068] X为待判别公司的指标变量矩阵；

[0069] 按N个预提取指标的违约率影响的由大到小的顺序，依次计算N个预提取指标的方差累加贡献率；

[0070] 取方差累积贡献率L达到预设阈值的前m个指标为违约预测指标，各指标的累积贡献率按如下公式计算：

[0071]

[0072] 其中，

[0073] N为预提取的指标数；

[0074] m为指标方差累积贡献率L达到预设阈值时的指标数；

[0075] λi为指标变量协方差矩阵Σ的特征向量λ中第i个指标的特征值。

[0076] 优选地，所述二项逻辑回归模型为：

[0077]

[0078] 其中：

[0079] pi表示待判别公司在第i年的违约概率，取值在0‑1之间，违约为1，非违约为0；

[0080] Xi为待判别公司在第i年的指标变量矩阵；

[0081] xki表示待判别公司公司第i年在第k个指标上的取值；

[0082] β0为回归常数；

[0083] βk为回归模型中第k个指标的回归系数。

[0084] 优选地，违约概率判断步骤中，还包括在得到违约概率后，结合待判别公司所属行业的违约倾向系数，对违约概率采用如下公式进行修正：

[0085] ci＝roi×pi

[0086] 其中：

[0087] ci为企业i某一段时期内根据其所属行业的违约倾向系数修正后的违约概率；

[0088] roi为行业i某一段时期内的违约倾向系数；

[0089] pi为企业i某一段时期内的违约概率。

[0090] 优选地，待判别公司所属行业的违约倾向系数采用灰度关联分析法确定。

[0091] 优选地，确定行业的违约倾向系数的方法包括如下步骤：

[0092] 选取数据库中所有行业的公司数据的指标变量，组成指标变量矩阵，其中选取的指标变量包括所述企业基本信息和企业财务指标；

[0093] 分别计算各行业包括的各公司指标变量矩阵与参考序列对应指标取值的灰度关联系数，所述参考序列由该行业中各指标的均值数据组成；

[0094] 分别计算各行业包括的各公司的各指标变量与参考序列对应元素的灰关联度；

[0095] 提取各行业包括的各公司的灰关联度，取各行业包括的各公司灰关联度的均值为该行业的违约倾向系数。

[0096] 优选地，各公司各个指标变量与参考序列对应元素的灰关联度的计算包括如下步骤：

[0097] 对各公司各个指标变量矩阵进行灰色无量纲化处理，

[0098]

[0099] 其中，

[0100] Xi(k)为第i个公司指标变量矩阵的对应第k个指标变量的取值；

[0101] Xi(1)为第i个公司指标变量矩阵的对应第1个指标变量的取值；

[0102] 计算差异信息序列Δoi(k)，

[0103] Δoi(k)＝|x0(k)‑xi(k)|

[0104] 其中，

[0105] x0(k)为参考序列进行灰色无量纲化处理过后的对应第k个指标变量的取值，所述参考序列由该行业中各指标的均值数据组成；

[0106] xi(k)为第i个公司指标变量矩阵灰色无量纲化后对应第k个指标变量的取值；

[0107] 确定环境参数：

[0108]

[0109]

[0110] 其中，

[0111] Δ0i(max)为差异信息序列Δoi(k)中的最大值；

[0112] Δ0i(min)为差异信息序列Δoi(k)中的最小值；

[0113] 分别采用如下公式计算各公司各指标变量与参考序列对应元素的灰关联系数ξi(k)，

[0114]

[0115] 其中，

[0116] ρ为分辨系数，0＜ρ＜1；

[0117] Δ0i(max)为差异信息序列Δoi(k)中的最大值；

[0118] Δ0i(min)为差异信息序列Δoi(k)中的最小值；

[0119] Δoi(k)为差异信息序列；

[0120] 分别采用如下公式计算各公司各指标变量与参考序列对应元素的灰关联度r0i，

[0121]

[0122] 其中，

[0123] ξi(k)为第i个公司的第k个指标与参考序列对应元素的灰关联系数；

[0124] m为企业数量信息表和所有行业违约数量信息表中包括的指标总数。

[0125] 与现有技术相对比，本发明的有益效果如下：

[0126] (1)本发明采用单因素分析或多重共线性验证法提取对被解释变量影响最大的几个解释变量，从而达到在多个基础变量中筛选出有效解释变量，达到初步降维的作用，同时
单因素分析选取出来的解释变量为对违约与否影响最大的变量。

[0127] (2)本发明取方差累积贡献率L达到预设阈值的前m个指标为违约预测指标，每一个指标都是各原始变量的线性组合，各指标之间互不相关，从而能有效利用大量统计数据
进行定量分析。

[0128] (3)本发明通过二项逻辑回归算法建立违约概率识别模型进行违约概率预测，实现了对未来一个财务周期内发生违约行为概率的识别。

附图说明

[0129] 图1是本发明的一个实施例的银行对公贷款违约判别方法流程图；

[0130] 图2是本发明的又一个实施例的银行对公贷款违约判别方法流程图；

[0131] 图3是本发明的一个实施例的单因素分析法预提取N个指标的流程图；

[0132] 图4是本发明的一个实施例的多重共线性验证法预提取N个指标的流程图；

[0133] 图5是本发明的一个实施例的各行业违约倾向系数计算步骤流程图。

具体实施方式

[0134] 下面结合附图1‑5，对本发明的具体实施方式作详细的说明。

[0135] 本发明提供了一种基于二项逻辑回归的银行对公贷款违约判别方法，包括如下步骤：

[0136] 违约预测指标预提取步骤，

[0137] 采用单因素方差分析法或多重共线性验证方法，根据数据库中的企业基本信息和企业财务指标，提取对违约率影响最大的N个指标；

[0138] 所述企业基本信息包括：企业名称、数据年度、是否违约、成立日期、在职员工数、总资产和借款人组织机构代码；其中是否违约、成立日期和在职员工数会在后续构成指标
变量矩阵时进行考虑；

[0139] 所述企业财务指标：主营业务收入、流动比率、速动比率、现金比率、资产负债率、负债权益比率、利息保障倍数、存货周转率、应收账款周转率、营业周期、流动资产周转率、
总资产周转率、固定资产周转率、总资产报酬率、成本费用利润率、销售净利率、净资产收益
率、销售毛利率、盈余现金保障倍数、现金流动负债比、营业总成本、现金流动负债比、管理
费用、销售费用、财务费用、无形资产、长期借款、短期借款、净利润和营业收入；

[0140] 违约预测指标确定步骤，

[0141] 根据预提取的N个指标，取方差累积贡献率达到预设阈值的指标作为违约预测的指标；

[0142] 二项逻辑回归模型建立步骤，

[0143] 将违约与否作为二项逻辑回归中的被解释变量，建立二项逻辑回归模型；

[0144] 违约判断步骤，

[0145] 根据二项逻辑回归模型及企业数据，得到违约概率；

[0146] 根据违约概率和预设的违约概率阈值，判断该企业是否违约。

[0147] 作为优选实施方式，采用单因素方差分析法提取对违约率影响最大的N个指标具体包括如下步骤：

[0148] 每个指标单独与违约与否进行单因素方差分析，取设定年份内的数据，根据如下公式得到检验统计量F的值：

[0149]

[0150] 其中，

[0151] k为设定年份数；

[0152] n为数据库中该企业同行业内所有公司的数量；

[0153] 为第i年数据库中该企业同行业内各公司对应该指标变量的均值；

[0154] 为设定年份内数据库中该企业同行业内各公司对应该指标变量的总均值；

[0155] ni为第i年的样本量；

[0156] xij第i年数据库中该企业同行业内第j个公司对应该指标变量的取值；

[0157] 根据得到的检验统计量F的数值，查询F分布表，得到检验p值；

[0158] 选择N个p值最小的指标，作为对违约率影响最大的N个指标，p值越小，表示对被解释变量违约率的影响越显著。

[0159] 作为优选实施方式，通过多重共线性验证方法提取对违约率影响最大的N个指标具体包括如下步骤：

[0160] 对各待检验指标进行多重共线性检验，检验公式为：c1X1+c2X2+…ciXi…+cnXn＝c0，若存在c0,c1,c2,…,cn使得该公式成立，则待检验指标存在多重共线性；

[0161] 其中，

[0162] Xi为第i个待检验指标；

[0163] ci为第i个待检验指标的系数；

[0164] c0为预设常数；

[0165] 计算各待检验指标的拟合优度：

[0166]

[0167] 其中，

[0168] Rj为第j个待检验指标的拟合优度；

[0169] Yj为第j个待检验指标的真实值；

[0170] 为第j个待检验指标以其余指标为解释变量进行回归时得到的回归模型估计量；

[0171] 为待检验指标真实值的均值；

[0172] 对各待检验指标的拟合优度进行F检验：

[0173]

[0174] 其中，

[0175] 为第j.个解释变量对其他解释变量的回归方程的判定系数；

[0176] k为解释变量的个数；

[0177] n为样本量；

[0178] j为解释变量的序号；

[0179] j.个用于F检验的解释变量的序号；

[0180] Fj为第j个待检验指标经过F检验后的拟合优度；

[0181] F(k‑1,n‑k)是指自由度从k‑1到n‑k的F分布；

[0182] 根据对各指标的拟合优度进行F检验得到的检验结果，查询F分布表，得到检验p值；

[0183] 选择N个p值最大的指标，作为对违约率影响最大的N个指标。

[0184] 为了得到准确的回归结果需要进行多重共线性检验(因为统计学认为有的指标存在相同的含义，因此需要去掉相同的指标)，为了确保多重共线性检验的准确性需要进行拟
合优度检验。检验p值小于给定显著性水平α，表示被解释变量与其他解释变量间存在显著
的线性关系，且p值越小，这种线性关系越是显著，挑选出p值最大的N个指标变量。

[0185] 多重共线性检验的任务是：(1)检验多重共线性是否存在；(2)判明存在多重共线性的范围。利用多重共线性，找出对被解释变量影响最大的前N个解释变量，并且经过多重
共线性检验的解释变量之间已经不存在共线性，可以达到跟单因素分析算法同样的效果。

[0186] 作为优选实施方式，违约预测指标确定步骤具体包括如下步骤：

[0187] 由预提取的N个指标组成指标变量矩阵；

[0188] 计算指标变量矩阵的协方差，得到指标变量协方差矩阵Σ；

[0189] 计算指标变量协方差矩阵Σ的特征向量λ＝(λ1,λ2,…,λN)；

[0190] 计算指标变量协方差矩阵Σ的特征向量λ中各特征值λi对应的特征向量γi；

[0191] 计算特征值λi对应的违约预测指标Yi，采用如下公式计算：

[0192] Yi＝γ′iX；

[0193] 其中，

[0194] Yi为特征值λi对应的违约预测指标；

[0195] X为待判别公司的指标变量矩阵；

[0196] 按N个预提取指标的违约率影响的由大到小的顺序，依次计算N个预提取指标的方差累加贡献率；

[0197] 取方差累积贡献率L达到预设阈值的前m个指标为违约预测指标，各指标的累积贡献率按如下公式计算：

[0198]

[0199] 其中，

[0200] N为预提取的指标数；

[0201] m为指标方差累积贡献率L达到预设阈值时的指标数；

[0202] λi为指标变量协方差矩阵Σ的特征向量λ中第i个指标的特征值。

[0203] 作为优选实施方式，所述二项逻辑回归模型为：

[0204]

[0205] 其中：

[0206] pi表示待判别公司在第i年的违约概率，取值在0‑1之间，违约为1，非违约为0；

[0207] Xi为待判别公司在第i年的指标变量矩阵；

[0208] xki表示待判别公司公司第i年在第k个指标上的取值；

[0209] β0为回归常数；

[0210] βk为回归模型中第k个指标的回归系数。

[0211] 上面二项逻辑回归模型的公式还可以变形为

[0212] 作为优选实施方式，违约概率判断步骤中，还包括在得到违约概率后，结合待判别公司所属行业的违约倾向系数，对违约概率采用如下公式进行修正：

[0213] ci＝roi×pi

[0214] 其中：

[0215] ci为企业i某一段时期内根据其所属行业的违约倾向系数修正后的违约概率；

[0216] roi为行业i某一段时期内的违约倾向系数；

[0217] pi为企业i某一段时期内的违约概率。

[0218] 作为优选实施方式，待判别公司所属行业的违约倾向系数采用灰度关联分析法确定。

[0219] 作为优选实施方式，确定行业的违约倾向系数的方法包括如下步骤：

[0220] 选取数据库中所有行业的公司数据的指标变量，组成指标变量矩阵，其中选取的指标变量包括所述企业基本信息和企业财务指标；

[0221] 分别计算各行业包括的各公司指标变量矩阵与参考序列对应指标取值的灰度关联系数，所述参考序列由该行业中各指标的均值数据组成；

[0222] 分别计算各行业包括的各公司的各指标变量与参考序列对应元素的灰关联度；

[0223] 提取各行业包括的各公司的灰关联度，取各行业包括的各公司灰关联度的均值为该行业的违约倾向系数。

[0224] 作为优选实施方式，各公司各个指标变量与参考序列对应元素的灰关联度的计算包括如下步骤：

[0225] 对各公司各个指标变量矩阵进行灰色无量纲化处理，

[0226]

[0227] 其中，

[0228] Xi(k)为第i个公司指标变量矩阵的对应第k个指标变量的取值；

[0229] Xi(1)为第i个公司指标变量矩阵的对应第1个指标变量的取值；

[0230] 计算差异信息序列Δoi(k)，

[0231] Δoi(k)＝|x0(k)‑xi(k)|

[0232] 其中，

[0233] x0(k)为参考序列进行灰色无量纲化处理过后的对应第k个指标变量的取值，所述参考序列由该行业中各指标的均值数据组成；

[0234] xi(k)为第i个公司指标变量矩阵灰色无量纲化后对应第k个指标变量的取值；

[0235] 确定环境参数：

[0236]

[0237]

[0238] 其中，

[0239] Δ0i(max)为差异信息序列Δoi(k)中的最大值；

[0240] Δ0i(min)为差异信息序列Δoi(k)中的最小值；

[0241] 分别采用如下公式计算各公司各指标变量与参考序列对应元素的灰关联系数ξi(k)，

[0242]

[0243] 其中，

[0244] ρ为分辨系数，0＜ρ＜1，ρ越小，关联系数间差异越大，区分能力越强，通常取ρ＝0.5；

[0245] Δ0i(max)为差异信息序列Δoi(k)中的最大值；

[0246] Δ0i(min)为差异信息序列Δoi(k)中的最小值；

[0247] Δoi(k)为差异信息序列；

[0248] 分别采用如下公式计算各公司各指标变量与参考序列对应元素的灰关联度r0i，灰关联度用于反映各公司与参考序列的关联关系，

[0249]

[0250] 其中，

[0251] ξi(k)为第i个公司的第k个指标与参考序列对应元素的灰关联系数；

[0252] m为企业数量信息表和所有行业违约数量信息表中包括的指标总数。

[0253] 实施例1

[0254] 根据本发明的一个具体实施方案，下面对本发明的基于二项逻辑回归的银行对公贷款违约判别方法进行详细说明。

[0255] 本发明提供了一种基于逻辑回归的银行对公贷款违约判别方法，包括如下步骤：

[0256] 违约预测指标预提取步骤，

[0257] 采用单因素方差分析法或多重共线性验证方法，根据数据库中的企业基本信息和企业财务指标，提取对违约率影响最大的N个指标；

[0258] 所述企业基本信息包括：企业名称、数据年度、是否违约、成立日期、在职员工数、总资产和借款人组织机构代码；

[0259] 所述企业财务指标：主营业务收入、流动比率、速动比率、现金比率、资产负债率、负债权益比率、利息保障倍数、存货周转率、应收账款周转率、营业周期、流动资产周转率、
总资产周转率、固定资产周转率、总资产报酬率、成本费用利润率、销售净利率、净资产收益
率、销售毛利率、盈余现金保障倍数、现金流动负债比、营业总成本、现金流动负债比、管理
费用、销售费用、财务费用、无形资产、长期借款、短期借款、净利润和营业收入；

[0260] 违约预测指标确定步骤，

[0261] 根据预提取的N个指标，取方差累积贡献率达到预设阈值的指标作为违约预测的指标；

[0262] 二项逻辑回归模型建立步骤，

[0263] 将违约与否作为二项逻辑回归中的被解释变量，建立二项逻辑回归模型；

[0264] 违约判断步骤，

[0265] 根据二项逻辑回归模型及企业数据，得到违约概率；

[0266] 根据违约概率和预设的违约概率阈值，判断该企业是否违约。

[0267] 实施例2

[0268] 根据本发明的一个具体实施方案，下面对本发明的基于二项逻辑回归的银行对公贷款违约判别方法进行详细说明。

[0269] 本发明提供了一种基于逻辑回归的银行对公贷款违约判别方法，包括如下步骤：

[0270] 违约预测指标预提取步骤，

[0271] 采用单因素方差分析法或多重共线性验证方法，根据数据库中的企业基本信息和企业财务指标，提取对违约率影响最大的N个指标；

[0272] 所述企业基本信息包括：企业名称、数据年度、是否违约、成立日期、在职员工数、总资产和借款人组织机构代码；

[0273] 所述企业财务指标：主营业务收入、流动比率、速动比率、现金比率、资产负债率、负债权益比率、利息保障倍数、存货周转率、应收账款周转率、营业周期、流动资产周转率、
总资产周转率、固定资产周转率、总资产报酬率、成本费用利润率、销售净利率、净资产收益
率、销售毛利率、盈余现金保障倍数、现金流动负债比、营业总成本、现金流动负债比、管理
费用、销售费用、财务费用、无形资产、长期借款、短期借款、净利润和营业收入；

[0274] 违约预测指标确定步骤，

[0275] 根据预提取的N个指标，取方差累积贡献率达到预设阈值的指标作为违约预测的指标；

[0276] 二项逻辑回归模型建立步骤，

[0277] 将违约与否作为二项逻辑回归中的被解释变量，建立二项逻辑回归模型；

[0278] 违约判断步骤，

[0279] 根据二项逻辑回归模型及企业数据，得到违约概率；

[0280] 结合待判别公司所属行业的违约倾向系数，对违约概率采用如下公式进行修正；

[0281] ci＝roi×pi

[0282] 其中：

[0283] ci为企业i某一段时期内根据其所属行业的违约倾向系数修正后的违约概率；

[0284] roi为行业i某一段时期内的违约倾向系数；

[0285] pi为企业i某一段时期内的违约概率；

[0286] 根据违约概率和预设的违约概率阈值，判断该企业是否违约。

[0287] 实施例3

[0288] 根据本发明的一个具体实施方案，下面对本发明的单因素方差分析法提取对违约率影响最大的N个指标的过程进行详细说明。

[0289] 采用单因素方差分析法提取对违约率影响最大的N个指标具体包括如下步骤：

[0290] 每个指标单独与违约与否进行单因素方差分析，取设定年份内的数据，根据如下公式得到检验统计量F的值：

[0291]

[0292] 其中，

[0293] k为设定年份数；

[0294] n为数据库中该企业同行业内所有公司的数量；

[0295] 为第i年数据库中该企业同行业内各公司对应该指标变量的均值；

[0296] 为设定年份内数据库中该企业同行业内各公司对应该指标变量的总均值；

[0297] ni为第i年的样本量；

[0298] xij第i年数据库中该企业同行业内第j个公司对应该指标变量的取值；

[0299] 根据得到的检验统计量F的数值，查询F分布表，得到检验p值；

[0300] 选择N个p值最小的指标，作为对违约率影响最大的N个指标。

[0301] 实施例4

[0302] 根据本发明的一个具体实施方案，下面对本发明的通过多重共线性验证方法提取对违约率影响最大的N个指标的过程进行详细说明。

[0303] 通过多重共线性验证方法提取对违约率影响最大的N个指标具体包括如下步骤：

[0304] 对各待检验指标进行多重共线性检验，检验公式为：c1X1+c2X2+…ciXi…+cnXn＝c0，若存在c0,c1,c2,…,cn使得该公式成立，则待检验指标存在多重共线性；

[0305] 其中，

[0306] Xi为第i个待检验指标；

[0307] ci为第i个待检验指标的系数；

[0308] c0为预设常数；

[0309] 计算各待检验指标的拟合优度：

[0310]

[0311] 其中，

[0312] Rj为第j个待检验指标的拟合优度；

[0313] Yj为第j个待检验指标的真实值；

[0314] 为第j个待检验指标以其余指标为解释变量进行回归时得到的回归模型估计量；

[0315] 为待检验指标真实值的均值；

[0316] 对各待检验指标的拟合优度进行F检验：

[0317]

[0318] 其中，

[0319] 为第j.个解释变量对其他解释变量的回归方程的判定系数；

[0320] k为解释变量的个数；

[0321] n为样本量；

[0322] j为解释变量的序号；

[0323] j.个用于F检验的解释变量的序号；

[0324] Fj为第j个待检验指标经过F检验后的拟合优度；

[0325] F(k‑1,n‑k)是指自由度从k‑1到n‑k的F分布；

[0326] 根据对各指标的拟合优度进行F检验得到的检验结果，查询F分布表，得到检验p值；

[0327] 选择N个p值最大的指标，作为对违约率影响最大的N个指标。

[0328] 实施例5

[0329] 根据本发明的一个具体实施方案，下面对本发明根据预提取的预测指标确定违约预测指标的过程进行详细说明。

[0330] 违约预测指标确定步骤具体包括如下步骤：

[0331] 由预提取的N个指标组成指标变量矩阵；

[0332] 计算指标变量矩阵的协方差，得到指标变量协方差矩阵Σ；

[0333] 计算指标变量协方差矩阵Σ的特征向量λ＝(λ1,λ2,…,λN)；

[0334] 计算指标变量协方差矩阵Σ的特征向量λ中各特征值λi对应的特征向量γi；

[0335] 计算特征值λi对应的违约预测指标Yi，采用如下公式计算：

[0336] Yi＝γ′iX；

[0337] 其中，

[0338] Yi为特征值λi对应的违约预测指标；

[0339] X为待判别公司的指标变量矩阵；

[0340] 按N个预提取指标的违约率影响的由大到小的顺序，依次计算N个预提取指标的方差累加贡献率；

[0341] 取方差累积贡献率L达到预设阈值的前m个指标为违约预测指标，各指标的累积贡献率按如下公式计算：

[0342]

[0343] 其中，

[0344] N为预提取的指标数；

[0345] m为指标方差累积贡献率L达到预设阈值时的指标数；

[0346] λi为指标变量协方差矩阵Σ的特征向量λ中第i个指标的特征值。

[0347] 实施例6

[0348] 根据本发明的一个具体实施方案，下面对本发明确定行业的违约倾向系数的方法进行详细说明。

[0349] 确定行业的违约倾向系数的方法包括如下步骤：

[0350] 选取数据库中所有行业的公司数据的指标变量，组成指标变量矩阵，其中选取的指标变量包括所述企业基本信息和企业财务指标；

[0351] 分别计算各行业包括的各公司指标变量矩阵与参考序列对应指标取值的灰度关联系数，所述参考序列由该行业中各指标的均值数据组成；

[0352] 分别计算各行业包括的各公司的各指标变量与参考序列对应元素的灰关联度；

[0353] 提取各行业包括的各公司的灰关联度，取各行业包括的各公司灰关联度的均值为该行业的违约倾向系数。

[0354] 实施例7

[0355] 根据本发明的一个具体实施方案，下面对本发明计算各公司各个指标变量与参考序列对应元素的灰关联度的过程进行详细说明。

[0356] 各公司各个指标变量与参考序列对应元素的灰关联度的计算包括如下步骤：

[0357] 对各公司各个指标变量矩阵进行灰色无量纲化处理，

[0358]

[0359] 其中，

[0360] Xi(k)为第i个公司指标变量矩阵的对应第k个指标变量的取值；

[0361] Xi(1)为第i个公司指标变量矩阵的对应第1个指标变量的取值；

[0362] 计算差异信息序列Δoi(k)，

[0363] Δoi(k)＝|x0(k)‑xi(k)|

[0364] 其中，

[0365] x0(k)为参考序列进行灰色无量纲化处理过后的对应第k个指标变量的取值，所述参考序列由该行业中各指标的均值数据组成；

[0366] xi(k)为第i个公司指标变量矩阵灰色无量纲化后对应第k个指标变量的取值；

[0367] 确定环境参数：

[0368]

[0369]

[0370] 其中，

[0371] Δ0i(max)为差异信息序列Δoi(k)中的最大值；

[0372] Δ0i(min)为差异信息序列Δoi(k)中的最小值；

[0373] 分别采用如下公式计算各公司各指标变量与参考序列对应元素的灰关联系数ξi(k)，

[0374]

[0375] 其中，

[0376] ρ为分辨系数，0＜ρ＜1；

[0377] Δ0i(max)为差异信息序列Δoi(k)中的最大值；

[0378] Δ0i(min)为差异信息序列Δoi(k)中的最小值；

[0379] Δoi(k)为差异信息序列；

[0380] 分别采用如下公式计算各公司各指标变量与参考序列对应元素的灰关联度r0i，

[0381]

[0382] 其中，

[0383] ξi(k)为第i个公司的第k个指标与参考序列对应元素的灰关联系数；

[0384] m为企业数量信息表和所有行业违约数量信息表中包括的指标总数。

[0385] 以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修
改、等同替换、改进等，均包含在本发明的保护范围之内。