金融风控实战——特征工程下

特征选择

Filter

1. 移除低方差的特征 (Removing features with low variance)
2. 单变量特征选择 (Univariate feature selection)

Wrapper

1. 递归特征消除 (Recursive Feature Elimination)

Embedded

1. 使用SelectFromModel选择特征 (Feature selection using SelectFromModel)
2. 将特征选择过程融入pipeline (Feature selection as part of a pipeline)

当数据预处理完成后，我们需要选择有意义的特征输入机器学习的算法和模型进行训练。

注意：一般在特征衍生前不做特征筛选，除了缺失值特别多又没啥用的变量。

通常来说，从两个方面考虑来选择特征:

特征是否发散

如果一个特征不发散，例如方差接近于0，也就是说样本在这个特征上基本上没有差异，这个特征对于样本的区分并没有什么用。（注意是要在归一化过程后）

特征与目标的相关性

这点比较显见，与目标相关性高的特征，应当优选选择。除移除低方差法外，可从相关性考虑根据特征选择的形式又可以将特征选择方法分为3种：

• Filter：过滤法，按照发散性或者相关性对各个特征进行评分，设定阈值或者待选择阈值的个数，选择特征。
• Wrapper：包装法，根据目标函数（通常是预测效果评分），每次选择若干特征，或者排除若干特征，看模型效果是否提升。
• Embedded：嵌入法，先使用某些机器学习的算法和模型进行训练，得到各个特征的权值系数，根据系数从大到小选择特征。类似于Filter方法，但是是通过训练来确定特征的优劣。（=数模随机森林判断特征重要性）

特征选择主要有两个目的：

• 减少特征数量、降维，使模型泛化能力更强，减少过拟合（主要方式：粗箱、筛选变量）；
• 增强对特征和特征值之间的理解。

拿到数据集，一个特征选择方法，往往很难同时完成这两个目的。通常情况下，选择一种自己最熟悉或者最方便的特征选择方法（往往目的是降维，而忽略了对特征和数据理解的目的）。接下来将结合Scikit-learn提供的例子介绍几种常用的特征选择方法，它们各自的优缺点和问题。

Filter

（1）移除低方差的特征 (Removing features with low variance)

假设某特征的特征值只有0和1，并且在所有输入样本中，95%的实例的该特征取值都是1，那就可以认为这个特征作用不大。如果100%都是1，那这个特征就没意义了。当特征值都是离散型变量的时候这种方法才能用，如果是连续型变量，就需要将连续变量离散化之后才能用。而且实际当中，一般不太会有95%以上都取某个值的特征存在，所以这种方法虽然简单但是不太好用。可以把它作为特征选择的预处理，先去掉那些取值变化小的特征，然后再从接下来提到的的特征选择方法中选择合适的进行进一步的特征选择。

import numpy as np
import pandas as pd
from sklearn.feature_selection import VarianceThreshold
X = np.array([[0, 0, 1], [0, 1, 0], [1, 0, 0], [0, 1, 1], [0, 1, 0], [0, 1, 1]])
sel = VarianceThreshold(threshold=(.8 * (1 - .8)))
sel.fit_transform(X)
#array([[0, 1],
#       [1, 0],
#       [0, 0],
#       [1, 1],
#       [1, 0],
#       [1, 1]])

果然, VarianceThreshold 移除了第一列特征，第一列中特征值为0的概率达到了5/6。

（2）单变量特征选择 (Univariate feature selection)

单变量特征选择的原理是分别单独的计算每个变量的某个统计指标，根据该指标来判断哪些变量重要，剔除那些不重要的变量。

对于分类问题(y离散)，可采用：

• 卡方检验
• f_classif
• mutual_info_classif
• 互信息

对于回归问题(y连续)，可采用：

• 皮尔森相关系数
• f_regression,
• mutual_info_regression
• 最大信息系数

这种方法比较简单，易于运行，易于理解，通常对于理解数据有较好的效果（但对特征优化、提高泛化能力来说不一定有效）。

• SelectKBest 移除得分前k名以外的所有特征(取top k)
• SelectPercentile 移除得分在用户指定百分比以后的特征(取top k%)
• 对每个特征使用通用的单变量统计检验： 假正率(false positive rate) SelectFpr, 伪发现率(false discovery rate) SelectFdr, 或族系误差率 SelectFwe.
• GenericUnivariateSelect 可以设置不同的策略来进行单变量特征选择。同时不同的选择策略也能够使用超参数寻优，从而让我们找到最佳的单变量特征选择策略。

Notice:
The methods based on F-test estimate the degree of linear dependency between two random variables. (F检验用于评估两个随机变量的线性相关性) On the other hand, mutual information methods can capture any kind of statistical dependency, but being nonparametric, they require more samples for accurate estimation. (另一方面，互信息的方法可以捕获任何类型的统计依赖关系，但是作为一个非参数方法，估计准确需要更多的样本)

卡方(Chi2)检验（和IV值差不多）

经典的卡方检验是检验定性自变量对定性因变量的相关性（离散与离散变量之间的相关性）。比如，我们可以对样本进行一次chi2测试来选择最佳的两项特征：

from sklearn.datasets import load_iris
from sklearn.feature_selection import SelectKBest
from sklearn.feature_selection import chi2
iris = load_iris()
X, y = iris.data, iris.target
print(X.shape)
X_new = SelectKBest(chi2, k=2).fit_transform(X, y)
print(X_new.shape)
#(150, 4)
#(150, 2)

Pearson相关系数 (Pearson Correlation)

皮尔森相关系数是一种最简单的，能帮助理解特征和响应变量之间关系的方法，该方法衡量的是变量之间的线性相关性，结果的取值区间为[-1，1]，-1表示完全的负相关，+1表示完全的正相关，0表示没有线性相关

import numpy as np
from scipy.stats import pearsonr
 
np.random.seed(0)
size = 300
x = np.random.normal(0, 1, size)
 
""" pearsonr(x, y)的输入为特征矩阵和目标向量，能够同时计算相关系数和p-value. """
print("Lower noise", pearsonr(x, x + np.random.normal(0, 1, size)))
print("Higher noise", pearsonr(x, x + np.random.normal(0, 10, size)))
""" 比较了变量在加入噪音之前和之后的差异。当噪音比较小的时候，相关性很强，p-value很低 """
""" 使用Pearson相关系数主要是为了看特征之间的相关性，而不是和因变量之间的。 """
#Lower noise (0.7182483686213834, 7.324017313000586e-49)
#Higher noise (0.05796429207933808, 0.3170099388532581)

Wrapper

递归特征消除 (Recursive Feature Elimination)

递归消除特征法使用一个基模型来进行多轮训练，每轮训练后，移除若干权值系数的特征，再基于新的特征集进行下一轮训练。

对特征含有权重的预测模型(例如，线性模型对应参数coefficients)，RFE通过递归减少考察的特征集规模来选择特征。首先，预测模型在原始特征上训练，每个特征指定一个权重。之后，那些拥有最小绝对值权重的特征被踢出特征集（做了woe就不能这么判断了，做了woe权重不能说明重要程度）。如此往复递归，直至剩余的特征数量达到所需的特征数量。

RFECV 通过交叉验证的方式执行RFE，以此来选择最佳数量的特征：对于一个数量为d的feature的集合，他的所有的子集的个数是2的d次方减1(包含空集)。指定一个外部的学习算法，比如SVM之类的。通过该算法计算所有子集的validation error。选择error最小的那个子集作为所挑选的特征。

from sklearn.feature_selection import RFE
from sklearn.ensemble import RandomForestClassifier
from sklearn.datasets import load_iris
rf = RandomForestClassifier()
iris=load_iris()
X,y=iris.data,iris.target
rfe = RFE(estimator=rf, n_features_to_select=3)
X_rfe = rfe.fit_transform(X,y)
X_rfe.shape
#(150, 3)

Embedded

使用SelectFromModel选择特征 (Feature selection using SelectFromModel)

基于L1的特征选择 (L1-based feature selection)

使用L1范数作为惩罚项的线性模型(Linear models)会得到稀疏解：大部分特征对应的系数为0。当你希望减少特征的维度以用于其它分类器时，可以通过feature_selection.SelectFromModel来选择不为0的系数。

常用于此目的的稀疏预测模型有linear_model.Lasso（回归），linear_model.Logistic Regression和svm.LinearSVC（分类）

from sklearn.feature_selection import SelectFromModel
from sklearn.svm import LinearSVC
 
lsvc = LinearSVC(C=0.01, penalty="l1", dual=False).fit(X,y)  #dual:布尔值。默认为False。如果等于True，则求解其对偶形式。
model = SelectFromModel(lsvc, prefit=True) #prefit ：布尔，默认为False，是否为训练完的模型，如果是False的话则先fit，再transform。
X_embed = model.transform(X)
X_embed.shape
#(150, 3)

实际业务

首先来回顾一下我们在业务中的模型会遇到什么问题。

1. 模型效果不好：大概率数据有问题
2. 训练集效果好，跨时间测试(一般测试样本是训练数据的1/10)效果不好：测试数据分布与训练数据不太一样导致的，说明选入特征变量有问题波动比较大，查看分析比较波动的特征变量
   跨时间测试效果也好，上线之后效果不好：线下和线上变量的逻辑出了问题，线下特征信息可能包含未来变量
3. 上线之后效果还好，几周之后分数分布开始下滑：说明模型效果不行，说明一两个变量在跨时间上效果比较差
4. 一两个月内都比较稳定，突然分数分布骤降：可能是外部因素，如运营部门一些操作或国家政策导致
5. 没有明显问题，但模型每个月逐步失效

然后我们来考虑一下业务所需要的变量是什么。

1. 变量必须对模型有贡献，也就是说必须能对客群加以区分
2. 逻辑回归要求变量之间线性无关
3. 逻辑回归评分卡也希望变量呈现单调趋势（有一部分也是业务原因，但从模型角度来看，单调变量未必一定比有转折的变量好）
4. 客群在每个变量上的分布稳定，分布迁移无可避免，但不能波动太大

为此我们从上述方法中找到最贴合当前使用场景的几种方法。

import pandas as pd
import numpy as np
df_train = pd.read_csv('/Users/zhucan/Desktop/金融风控实战/第三课资料/train.csv')
df_train.head()

​​​​​​

1）变量重要性

• IV值
• 卡方检验
• 模型筛选

这里我们使用IV值或者模型筛选多一点

IV其实就是在WOE前面加上一项。

• $p_{y_i}=\frac{y_i}{y_T}$
• $p_{n_i}=\frac{n_i}{n_T}$
• $wo{e}_i=ln(\frac{p_{y_i}}{p_{n_i}})$
• $i{v}_i=(p_{y_i}-p_{n_i})\times wo{e}_i$

最后只需要将每一个区间的iv加起来就得到总的iv值：
$$IV=\sum i{v}_i$$

a = 0.4  
b = 0.6
iv = (a - b) * math.log(a / b)

分箱、WOE、IV

import numpy as np
import pandas as pd
from scipy import stats

def mono_bin(Y,X,n=20):
    r=0
    good = Y.sum()
    bad = Y.count()-good
    while np.abs(r)< 1:
        d1=pd.DataFrame({"X":X,"Y":Y,"Bucket":pd.qcut(X,n)})
        d2=d1.groupby('Bucket',as_index=True)
        r,p=stats.spearmanr(d2.mean().X,d2.mean().Y)
        n=n-1
    d3=pd.DataFrame(d2.X.min(),columns=['min'])
    d3['min']=d2.min().X
    d3['max']=d2.max().X
    d3['sum']=d2.sum().Y
    d3['total']=d2.count().Y
    d3['rate']=d2.mean().Y
    d3['woe']=np.log((d3['rate']/(1-d3['rate']))/(good/bad))
    d3['iv']=(d3['rate']/(1-d3['rate']) - (good/bad)) * np.log((d3['rate']/(1-d3['rate']))/(good/bad))
    d4=(d3.sort_values(by='min')).reset_index(drop=True)
    print("="*60)
    print(d4)
    return d4
mono_bin(df_train["label"],df_train["Age"],n = 20)

​​​​​​
或者集成模型输出特征重要性：

#lightGBM中的特征重要性
feature = pd.DataFrame(
            {'name' : model.booster_.feature_name(),
            'importance' : model.feature_importances_
          }).sort_values(by =  ['importance'],ascending = False)

2）共线性

• 相关系数 COR
• 方差膨胀系数 VIF

在做很多基于空间划分思想的模型的时候，我们必须关注变量之间的相关性。单独看两个变量的时候我们会使用皮尔逊相关系数。

df_train.corr()

import seaborn as sns
sns.set(color_codes=True)
np.random.seed(sum(map(ord, "distributions")))
# 在数据集中绘制成对关系
sns.pairplot(df_train) # 对角线上是单维度分布

在多元回归中，我们可以通过计算方差膨胀系数VIF来检验回归模型是否存在严重的多重共线性问题。定义：
$$VIF=\frac{1}{1-R_i^2}$$

其中，$R_i$为自变量 对其余自变量作回归分析的负相关系数。方差膨胀系数是容忍度$1-R^2$的倒数。

方差膨胀系数VIF越大，说明自变量之间存在共线性的可能性越大。一般来讲，如果方差膨胀因子超过10，则回归模型存在严重的多重共线性。又根据Hair(1995)的共线性诊断标准，当自变量的容忍度大于0.1，方差膨胀系数小于10的范围是可以接受的，表明自变量之间没有共线性问题存在。

3）单调性

• bivar图

# 等频切分
df_train.loc[:,'fare_qcut'] = pd.qcut(df_train['Fare'], 10)
df_train.head()
df_train = df_train.sort_values('Fare')
alist = list(set(df_train['fare_qcut']))
badrate = {}
for x in alist:
    a = df_train[df_train.fare_qcut == x]
    bad = a[a.label == 1]['label'].count()
    good = a[a.label == 0]['label'].count()
    badrate[x] = bad/(bad+good)
f = zip(badrate.keys(),badrate.values())
f = sorted(f,key = lambda x : x[1],reverse = True )
badrate = pd.DataFrame(f)
badrate.columns = pd.Series(['cut','badrate'])
badrate = badrate.sort_values('cut')
print(badrate)
badrate.plot('cut','badrate')
#                 cut   badrate
#9     (-0.001, 7.55]  0.141304
#6      (7.55, 7.854]  0.298851
#8      (7.854, 8.05]  0.179245
#7       (8.05, 10.5]  0.230769
#3     (10.5, 14.454]  0.428571
#4   (14.454, 21.679]  0.420455
#2     (21.679, 27.0]  0.516854
#5     (27.0, 39.688]  0.373626
#1   (39.688, 77.958]  0.528090
#0  (77.958, 512.329]  0.758621

def binn(x):
    if x <10.5:
        return 0
    elif x <39.688:
        return 1
    else:
        return 2
df_train["fare_cut_new"] = df_train.Fare.map(lambda x:binn(x))
df_train.sample(100)

df_train = df_train.sort_values('Fare')
alist = list(set(df_train['fare_cut_new']))
badrate = {}
for x in alist:
    a = df_train[df_train.fare_cut_new == x]
    bad = a[a.label == 1]['label'].count()
    good = a[a.label == 0]['label'].count()
    badrate[x] = bad/(bad+good)
f = zip(badrate.keys(),badrate.values())
f = sorted(f,key = lambda x : x[1],reverse = True )
badrate = pd.DataFrame(f)
badrate.columns = pd.Series(['cut','badrate'])
badrate = badrate.sort_values('cut')
print(badrate)
badrate.plot('cut','badrate')

4）稳定性

• PSI
• 跨时间交叉检验

跨时间交叉检验

就是将样本按照月份切割，一次作为训练集和测试集来训练模型，取进入模型的变量之间的交集，但是要小心共线特征！

（1）第一个月作为测试集，后面十一个月作为训练集，训练模型，输出变量的重要性
（2）第二个月作为测试集，剩余十一个月作为训练集，训练模型，输出变量的重要性
。。。
（12）最后一个月作为测试集，剩余十一个月作为训练集，训练模型，输出变量的重要性
（13）取交集

解决方法

• 不需要每次都进入模型，大部分都在即可
• 先去除共线性（这也是为什么集成模型我们也会去除共线性）

群体稳定性指标(population stability index)

公式：

$$PSI=\sum (实际占比-预期占比)\ast \ln (\frac{实际占比}{预期占比})$$

来自知乎的例子：
比如训练一个logistic回归模型，预测时候会有个概率输出p。
你测试集上的输出设定为p1吧，将它从小到大排序后10等分，如0-0.1,0.1-0.2,…。
现在你用这个模型去对新的样本进行预测，预测结果叫p2,按p1的区间也划分为10等分。
实际占比就是p2上在各区间的用户占比，预期占比就是p1上各区间的用户占比。
意义就是如果模型跟稳定，那么p1和p2上各区间的用户应该是相近的，占比不会变动很大，也就是预测出来的概率不会差距很大。
一般认为psi小于0.1时候模型稳定性很高，0.1-0.25一般，大于0.25模型稳定性差，建议重做。

def var_PSI(dev_data, val_data):
    dev_cnt, val_cnt = sum(dev_data), sum(val_data)
    if dev_cnt * val_cnt == 0:
        return None
    PSI = 0
    for i in range(len(dev_data)):
        dev_ratio = dev_data[i] / dev_cnt
        val_ratio = val_data[i] / val_cnt + 1e-10
        psi = (dev_ratio - val_ratio) * math.log(dev_ratio/val_ratio)
        PSI += psi
    return PSI

注意分箱的数量将会影响着变量的PSI值。
PSI并不只可以对模型来求，对变量来求也一样。只需要对跨时间分箱的数据分别求PSI即可。

excel计算PSI

PSI评判标准