评分卡建模工具scorecardpy全解读

说到评分卡建模工具，做过评分卡的应该都能想到谢博士的scorecardpy和专为工业界模型开发设计的Python工具包--Toad，两者相比，scorecardpy更加轻量级，且依赖较少，可以满足大多数场景下的评分卡构建。

为了使评分卡建模流程更加便捷，该Python包针对建模中各个关键步骤都提供了现成的函数，如下：

• 数据集划分 (split_df)

• 变量筛选(iv, var_filter)

• 变量分箱(woebin, woebin_plot, woebin_adj, woebin_ply)

• 分数转换(scorecard, scorecard_ply)

• 效果评估(perf_eva, perf_psi)

基于github主页上提供的案例，对每一步操作进行详细解读：

 

数据准备

import scorecardpy as sc # 加载德国信用卡相关数据集 dat = sc.germancredit()

得到的数据集情况如下：

 

变量筛选

然后根据变量的缺失率、IV值、同值性等因素，对20个变量进行筛选，并对目标变量进行标记，该数据集中creditability为目标变量：

dt_s = sc.var_filter(dat, y="creditability")

通过 var_filter() 函数筛选后，数据集保留以下13个变量：

默认的参数配置为：iv_limit=0.02, missing_limit=0.95, identical_limit=0.95，即当某个变量的 IV 值小于0.02，或缺失率大于95%，或同值率（除空值外）大于95%，则剔除掉该变量。

此外，该方法还内置了除上述以外的其他参数：

def var_filter(dt, y, x=None, iv_limit=0.02, missing_limit=0.95, identical_limit=0.95, var_rm=None, var_kp=None, return_rm_reason=False, positive='bad|1')

其中各参数含义如下：

• varrm可设置强制保留的变量，默认为空；

• varkp可设置强制剔除的变量，默认为空；

• return_rm_reason可设置是否返回剔除原因，默认为不返回（False）；

• positive可设置坏样本对应的值，默认为“bad|1”。

 

数据集划分

将数据集划分成训练集和测试集：

train, test = sc.split_df(dt_s, 'creditability').values()

split_df()函数包括以下参数：

def split_df(dt, y=None, ratio=0.7, seed=186)

默认的切分比例为7:3，也可自行修改参数ratio设置不同的切分比例。

 

变量分箱

可通过woebin()函数对全部变量进行自动分箱，并基于woe_bin的结果，使用woebin_plot对各变量分箱的count distribution和bad probability进行可视化，可观察是否存在单调性：

bins = sc.woebin(dt_s, y="creditability") sc.woebin_plot(bins)

bins为字典形式，以变量名作为key，可获得单个变量的分箱结果详情：

woebin_plot()得到的可视化结果如下：

woebin()函数包括如下参数：

def woebin(dt, y, x=None, var_skip=None, breaks_list=None, special_values=None, stop_limit=0.1, count_distr_limit=0.05, bin_num_limit=8, # min_perc_fine_bin=0.02, min_perc_coarse_bin=0.05, max_num_bin=8, positive="bad|1", no_cores=None, print_step=0, method="tree", ignore_const_cols=True, ignore_datetime_cols=True, check_cate_num=True, replace_blank=True, save_breaks_list=None, **kwargs)

woebin()可针对数值型和类别型变量生成最优分箱结果，方法可选择决策树分箱、卡方分箱或自定义分箱。其他各参数的含义如下：

• var_skip: 设置需要跳过分箱操作的变量；

• breaks_list: 切分点列表，默认为空。如果非空，则按设置的切分点进行分箱处理。

• special_values: 设置需要单独分箱的值，默认为空。

• count_distr_limit: 设置分箱占比的最小值，一般可接受范围为0.01-0.2，默认值为0.05；

• stop_limit: 当IV值的增长率小于所设置的stop_limit，或卡方值小于qchisq(1-stoplimit, 1)时，停止分箱。一般可接受范围为0-0.5，默认值为0.1；

• bin_num_limit: 该参数为整数，代表最大分箱数。

• positive: 指定样本中正样本对应的标签，默认为"bad|1"；

• no_cores: 设置用于并行计算的 CPU 数目；

• print_step: 该参数为非负数，默认值为1。若print_step>0，每次迭代会输出变量名。若iteration=0或no_cores>1，不会输出任何信息；

• method: 设置分箱方法，可设置"tree"（决策树）或"chimerge"（卡方），默认值为"tree";

• ignore_const_cols: 是否忽略常数列，默认值为True，即忽略常数列；

• ignore_datetime_cols: 是否忽略日期列，默认值为True，即忽略日期列；

• check_cate_num: 检查类别变量中枚举值数目是否大于50，默认值为True，即自动进行检查。若枚举值过多，会影响分箱过程的速度；

• replace_blank: 设置是否将空值填为None，默认为True。

若对自动分箱结果不满意，还可手动自定义分箱：

breaks_adj = { 'age.in.years': [26, 35, 40], 'other.debtors.or.guarantors': ["none", "co-applicant%,%guarantor"] } bins_adj = sc.woebin(dt_s, y="creditability", breaks_list=breaks_adj)

于是年龄变量的分箱结果变成：

分箱完成后，使用woebin_ply()函数对变量进行woe变换，之后需要把所得到的woe值作为模型的输入：

train_woe = sc.woebin_ply(train, bins_adj) test_woe = sc.woebin_ply(test, bins_adj)

 

模型训练

y_train = train_woe.loc[:,'creditability'] X_train = train_woe.loc[:,train_woe.columns != 'creditability'] y_test = test_woe.loc[:,'creditability'] X_test = test_woe.loc[:,train_woe.columns != 'creditability'] # 逻辑回归 ------ from sklearn.linear_model import LogisticRegression lr = LogisticRegression(penalty='l1', C=0.9, solver='saga', n_jobs=-1) lr.fit(X_train, y_train) lr.coef_ lr.intercept_ # 预测 train_pred = lr.predict_proba(X_train)[:,1] test_pred = lr.predict_proba(X_test)[:,1]

 

效果评估

可使用perf_eva()函数对模型效果进行计算及可视化，基于预测的概率值和label值，提供KS(kolmogorov-smirnow), ROC, lift以及precision-recall四种评估指标：

train_perf = sc.perf_eva(y_train, train_pred, title = "train") test_perf = sc.perf_eva(y_test, test_pred, title = "test")

该函数还包括以下参数：

def perf_eva(label, pred, title=None, groupnum=None, plot_type=["ks", "roc"], show_plot=True, positive="bad|1", seed=186)

参数plot_type可设置为："ks", "lift", "roc", "pr"，默认为["ks", "roc"]。

 

分数转换

然后基于’woebin'的结果和sklearn.linear_model的LogisticRegression，创建scorecard()函数，用于构建评分卡，只需一行代码：

card = sc.scorecard(bins_adj, lr, X_train.columns)

生成的结果为各变量名及其分箱、对应得分组成的字典，例如：

scorecard()包括以下参数：

def scorecard(bins, model, xcolumns, points0=600, odds0=1/19, pdo=50, basepoints_eq0=False, digits=0)

各参数含义如下：

• bins：由`woebin`得到的分箱信息；

• model：LogisticRegression模型对象；

• points0：基准分数，默认值为600；

• odds0: 基准 Odds（好坏比），与真实违约概率对应，可换算得到违约概率，Odds = p/(1-p)。默认值为 1/19；

• pdo: Points toDouble theOdds，即Odds变成2倍时，所增加的信用分。默认值为50。

然后基于`scorecard`的结果，用scorecard_ply()函数计算train和test数据集的信用分数：

train_score = sc.scorecard_ply(train, card, print_step=0) test_score = sc.scorecard_ply(test, card, print_step=0)

直接生成各样本的评分：

最后用perf_psi()得到该评分卡在测试数据集上的表现：

sc.perf_psi( score = {'train':train_score, 'test':test_score}, label = {'train':y_train, 'test':y_test} )

perf_psi()可计算群体稳定性指标（PSI），并提供信用分数的分布图：

 

参考资料

scorecardpy github主页：https://github.com/ShichenXie/scorecardpy

---

⬇️ 扫描下方二维码关注公众号【小数据研究院】 ⬇️

回复【算法】，获取最全面的机器学习算法网络图：