SDU项目实训记录3.3——信用评分（7.15）

客官°小女子只卖身不卖艺 2023-01-23 08:55 17阅读 0赞

### SDU项目实训记录3.3——信用评分 ###

*  一、基础知识
 *  二、实现过程

# 一、基础知识 #

评分卡设定的分数刻度值通过将分值表示为比率对数的线性表达式来定义  
score = A –Blog(Odds) , 其中Odds = p/(1-p) ,p为违约概率,常数A B可以通过两个已知假设分值带入到公式中，通常需要以下两个假设：

*  在某个特定比率设定特定的预期分值
 *  指定比率翻倍的分数PDO

**规定：** odds=1:60时 P0=600 PDO=20  
![在这里插入图片描述][watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3dlaXhpbl80Mzk1OTQyMQ_size_16_color_FFFFFF_t_70]  
推导公式可得：

*  B = PDO / log(2)
 *  A = P0 + B \* log(odds)

参考：[信用评分卡模型的理论准备][Link 1]3.1评分卡刻度

# 二、实现过程 #

根据之前训练得到的逻辑回归模型得到对应截距和权重参数：

print("截距项:", model.intercept_[0])
    coes = model.coef_[0]
    print("权重参数:", coes)

![在这里插入图片描述][20210713234704301.png]  
设P0为600，PDO为20（每高20分好坏比翻一倍），好坏比取60的前提下，我们计算得基础分数为558分。

odds = 1/60
    P0 = 600
    PDO = 20
    #计算A、B
    B = PDO / np.log(2)
    A = P0 + B * np.log(odds)
    BaseScore = round(A - B * model.intercept_[0], 0)
    print('BaseScore =', BaseScore)
    print('A =', A)  # <A认为会违约

![在这里插入图片描述][20210713234806281.png]  
定义由权重参数和woe计算区间分值的函数，调用得到各分箱区间对应的分值如下：

def get_score(coe,woe,factor):
        scores=[]
        for w in woe:
            score=round(coe * w * factor,0)
            scores.append(score)
        return scores
    i = 0
    scores = {
        }
    for f in features_saving:
        score = get_score(coes[i], woes[f], B)
        scores[f] = score
        print(f, ':', score)
        i += 1

![在这里插入图片描述][20210713234852253.png]  
对应边界，获得列表：  
![在这里插入图片描述][watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3dlaXhpbl80Mzk1OTQyMQ_size_16_color_FFFFFF_t_70 1]

[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3dlaXhpbl80Mzk1OTQyMQ_size_16_color_FFFFFF_t_70]: /images/20221004/73d66ea130dd4a7583c7639fa23b3ab1.png
[Link 1]: https://www.cnblogs.com/shaocf/p/9147135.html#_label3_0
[20210713234704301.png]: /images/20221004/fc8c4f0b3a4840248094126827c115b6.png
[20210713234806281.png]: /images/20221004/c9c45c998b5248efbd9206cae062e4ab.png
[20210713234852253.png]: /images/20221004/1b04540e1f7e4f7faad16c69efc01b16.png
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3dlaXhpbl80Mzk1OTQyMQ_size_16_color_FFFFFF_t_70 1]: /images/20221004/fe64a6f936fe4eee8b2aa4b5e4c1483e.png