GBDT与xgboost :流失预测 shap解释 调参 保存调参好的模型

发布时间:2023-06-14 08:30

集成学习

集成学习的方式分为两类:
个体的学习器之间存在强依赖关系,必须串行生成序列化方法,代表
Boosting;
个体学习器之间不存在强依赖关系,可同时生成并行化方法,代表是Bagging和随机森林。

GBDT与xgboost :流失预测 shap解释 调参 保存调参好的模型_第1张图片
GBDT与xgboost :流失预测 shap解释 调参 保存调参好的模型_第2张图片

bagging

GBDT与xgboost :流失预测 shap解释 调参 保存调参好的模型_第3张图片

boosting

GBDT与xgboost :流失预测 shap解释 调参 保存调参好的模型_第4张图片
GBDT与xgboost :流失预测 shap解释 调参 保存调参好的模型_第5张图片

stacking

GBDT与xgboost :流失预测 shap解释 调参 保存调参好的模型_第6张图片
2)点击率预估
使用GBDT+LR进行点击率(CTR)预估。
https://blog.csdn.net/Snoopy_Yuan/article/details/80703175?depth_1-utm_source=distribute.pc_r
elevant.none-task&utm_source=distribute.pc_relevant.none-task

GBDT与xgboost :流失预测 shap解释 调参 保存调参好的模型_第7张图片
GBDT与xgboost :流失预测 shap解释 调参 保存调参好的模型_第8张图片

XGBOOST 的推导过程

GBDT与xgboost :流失预测 shap解释 调参 保存调参好的模型_第9张图片
实例:

在这里插入图片描述
GBDT与xgboost :流失预测 shap解释 调参 保存调参好的模型_第10张图片
XGBOOST 是一个加法模型
GBDT与xgboost :流失预测 shap解释 调参 保存调参好的模型_第11张图片
GBDT与xgboost :流失预测 shap解释 调参 保存调参好的模型_第12张图片
T为叶子节点的个数,Wj为叶子节点的权重。

GBDT与xgboost :流失预测 shap解释 调参 保存调参好的模型_第13张图片
GBDT与xgboost :流失预测 shap解释 调参 保存调参好的模型_第14张图片
GBDT与xgboost :流失预测 shap解释 调参 保存调参好的模型_第15张图片
GBDT与xgboost :流失预测 shap解释 调参 保存调参好的模型_第16张图片
GBDT与xgboost :流失预测 shap解释 调参 保存调参好的模型_第17张图片

Xgboost-重要参数

GBDT与xgboost :流失预测 shap解释 调参 保存调参好的模型_第18张图片
GBDT与xgboost :流失预测 shap解释 调参 保存调参好的模型_第19张图片
对于General paramaters,表示使用何种基本模型,通常是决策树或线性模型,一般不需要调整,采
用默认值即可,参数调优不包括这部分

主要参数是以下几类:
GBDT与xgboost :流失预测 shap解释 调参 保存调参好的模型_第20张图片
GBDT与xgboost :流失预测 shap解释 调参 保存调参好的模型_第21张图片
对于Booster paramaters,表示若学习器相关参数,需要仔细调整,会影响模型性能。.
eta 学习率
gamma
GBDT与xgboost :流失预测 shap解释 调参 保存调参好的模型_第22张图片
GBDT与xgboost :流失预测 shap解释 调参 保存调参好的模型_第23张图片
GBDT与xgboost :流失预测 shap解释 调参 保存调参好的模型_第24张图片
GBDT与xgboost :流失预测 shap解释 调参 保存调参好的模型_第25张图片
对于Learning task parameters,与任务有关,定下来后通常不需要调整

GBDT与xgboost :流失预测 shap解释 调参 保存调参好的模型_第26张图片

XGBoost调参技巧

GBDT与xgboost :流失预测 shap解释 调参 保存调参好的模型_第27张图片
XGBOOST 比GBDT的特点
GBDT与xgboost :流失预测 shap解释 调参 保存调参好的模型_第28张图片

选用哪一种xgboost


*

代码实操* :

## 获取数据
import xgboost as xgb # 调用陈天奇的库
from xgboost.sklearn import XGBClassifier # 使用sklearn api
import numpy as np 
import pandas as pd 
import matplotlib.pyplot as plt

from sklearn import metrics
data =  pd.read_table('C:/Users/lb/Desktop/test/gouwu.txt',sep='\t',engine="python",encoding = 'utf-8')
data.columns.values
data.head()

GBDT与xgboost :流失预测 shap解释 调参 保存调参好的模型_第29张图片

data.shape

在这里插入图片描述
分割数据集

from sklearn.model_selection import train_test_split
data_x = data.drop(['留存标签'],axis=1)
data_y = data['留存标签']
train_x,test_x,train_y,test_y = train_test_split(data_x,data_y,test_size=0.3,random_state=5)

行标签的恢复以后注意添加这个步骤

for i in [train_x,test_x]:
    i.index = range(i.shape[0])
train_x.head()

GBDT与xgboost :流失预测 shap解释 调参 保存调参好的模型_第30张图片
模型的训练

from sklearn.metrics import accuracy_score

model = XGBClassifier( ) 
# 训练模型
model.fit(train_x,train_y) 
# 预测模型
pred_y = model.predict(test_x)
metrics.accuracy_score(test_y, pred_y)

分类的准确率

在这里插入图片描述

xgboost 进阶用法

GBDT与xgboost :流失预测 shap解释 调参 保存调参好的模型_第31张图片
监控每一步的表现

# 监控模型每一步的表现
# eval_set=[(x_test,y_test)]  评估数据集,list类型
# eval_metric 评估标准(多分类问题,使用mlogloss作为损失函数),二分类 auc
# early_stopping_rounds= 10   如果模型的loss十次内没有减小,则提前结束模型训练
# verbose = True  True显示,False不显示
eval_set = [(test_x, test_y)] 
model.fit(train_x, train_y, early_stopping_rounds=10, eval_metric="auc",
          eval_set=eval_set, verbose=True)

GBDT与xgboost :流失预测 shap解释 调参 保存调参好的模型_第32张图片

显示特征的重要性

#显示各个特征的重要性  ----  基于训练的特征树输出
from xgboost import plot_importance
from matplotlib import pyplot
plt.rcParams['font.sans-serif']=['SimHei'] #显示中文标签
plt.rcParams['axes.unicode_minus']=False #用来正常显示负号
plot_importance(model)
pyplot.show()

GBDT与xgboost :流失预测 shap解释 调参 保存调参好的模型_第33张图片
为了得到影响因素正向还是逆向的影响,用shap

shap 来解释xgboost模型

import shap
shap.initjs()
explainer = shap.TreeExplainer(model)
shap_values = explainer.shap_values(train_x)
print(shap_values)

GBDT与xgboost :流失预测 shap解释 调参 保存调参好的模型_第34张图片

shap_values.shape 
train_x.shape 

GBDT与xgboost :流失预测 shap解释 调参 保存调参好的模型_第35张图片

plt.rcParams['font.sans-serif']=['SimHei'] #显示中文标签
plt.rcParams['axes.unicode_minus']=False #用来正常显示负号
shap.summary_plot(shap_values,train_x,max_display=30)
#红色高 蓝色低

GBDT与xgboost :流失预测 shap解释 调参 保存调参好的模型_第36张图片

shap.summary_plot(shap_values,train_x, plot_type="bar")

GBDT与xgboost :流失预测 shap解释 调参 保存调参好的模型_第37张图片
做一个简单的验证

pd.DataFrame(abs(shap_values),columns=train_x.columns).mean().sort_values(ascending=False)

GBDT与xgboost :流失预测 shap解释 调参 保存调参好的模型_第38张图片

# 多个变量交互
shap_interaction_values = shap.TreeExplainer(model).shap_interaction_values(train_x)
shap.summary_plot(shap_interaction_values, train_x, max_display=4)

GBDT与xgboost :流失预测 shap解释 调参 保存调参好的模型_第39张图片

依赖图

依赖图
同样可以通过单变量的依赖图看出单变量对目标的影响,横坐标是特征值纵坐标是shap值:

shap.dependence_plot('购物距离注册时长', 
                     shap_values,train_x, 
                     interaction_index=None,
                     show=False)

GBDT与xgboost :流失预测 shap解释 调参 保存调参好的模型_第40张图片
可以看出购物距离注册市场对于用户的留存存在正相关的作用
类似案例的参考解释:
GBDT与xgboost :流失预测 shap解释 调参 保存调参好的模型_第41张图片
GBDT与xgboost :流失预测 shap解释 调参 保存调参好的模型_第42张图片
那么双变量呢? 参考案例如下:
其他别人案例代码:

shap.dependence_plot('Pclass', shap_values[1], X_test, interaction_index='Sex_encoding', show=False)

GBDT与xgboost :流失预测 shap解释 调参 保存调参好的模型_第43张图片
得出一个经验:不要看太多的颜色而是方向
[添加链接描述] 可再去参考这个教程(https://blog.csdn.net/weixin_43615654/article/details/103436632?utm_medium=distribute.pc_relevant.none-task-blog-BlogCommendFromMachineLearnPai2-1.channel_param&depth_1-utm_source=distribute.pc_relevant.none-task-blog-BlogCommendFromMachineLearnPai2-1.channel_param)

观察某一行各个特征对结果的影响

shap.initjs()
shap.force_plot(explainer.expected_value, shap_values[5], train_x.iloc[5])

GBDT与xgboost :流失预测 shap解释 调参 保存调参好的模型_第44张图片

模型调参

from sklearn.model_selection import GridSearchCV
#这里交叉验证也可以  K折
from sklearn.model_selection import StratifiedKFold 
#设定网格搜索的xgboost参数搜索范围,值搜索XGBoost的主要6个参数.
#注意看这里是个字典
param_dist = {
        'n_estimators':range(80,150,5),
        'max_depth':range(3,10,1),
        'learning_rate':np.linspace(0.01,0.1,20),
        'subsample':np.linspace(0.7,0.9,20),
        'colsample_bytree':np.linspace(0.5,0.98,10),
        'min_child_weight':range(1,9,1)
        }

网格搜索

grid = GridSearchCV(model,param_dist,cv = 3,scoring = 'neg_log_loss',
                    n_jobs = -1)
#在训练集上训练
grid.fit(train_x, train_y)
#返回最优的训练器
grid.best_params_

model_best = XGBClassifier(grid.best_params_)

速度慢下3.

#上面这个速度太慢
model1 = XGBClassifier()
learning_rate = [0.0001, 0.001, 0.01, 0.1,0.12,0.15,0.2,0.3]
#转成字典
param_grid = dict(learning_rate=learning_rate)```

```python
param_grid 

在这里插入图片描述
score 在不同场景情况下的取值
Python-sklearn包中StratifiedKFold和KFold生成交叉验证数据集的区别
利用KFOLD 实现交叉验证
GBDT与xgboost :流失预测 shap解释 调参 保存调参好的模型_第45张图片

GBDT与xgboost :流失预测 shap解释 调参 保存调参好的模型_第46张图片

model1 = XGBClassifier()
learning_rate = [0.0001, 0.001, 0.01, 0.1,0.12,0.15,0.2,0.3]
param_grid = dict(learning_rate=learning_rate)
kfold = StratifiedKFold(n_splits=10, shuffle=True, random_state=7)
grid_search = GridSearchCV(model1, param_grid, scoring="neg_log_loss",
                           n_jobs=-1, cv=kfold)
grid_result = grid_search.fit(train_x, train_y)
print("Best: %f using %s" % (grid_result.best_score_, grid_result.best_params_))

在这里插入图片描述

model2 = XGBClassifier()
n_estimators = range(80,200,5)
param_grid = dict(n_estimators = n_estimators)
kfold = StratifiedKFold(n_splits=10, shuffle=True, random_state=7)
grid_search = GridSearchCV(model2, param_grid, scoring="neg_log_loss",
                           n_jobs=-1, cv=kfold)
grid_result = grid_search.fit(train_x, train_y)
print("Best: %f using %s" % (grid_result.best_score_, grid_result.best_params_))

在这里插入图片描述

# 调参
model3 = XGBClassifier(n_estimators=110,learning_rate= 0.12)
model3.fit(train_x,train_y) 
# 预测模型
pred_y = model3.predict(test_x)
metrics.accuracy_score(test_y, pred_y)

0.7625979843225084

保存模型

# 保存模型
import pickle
pickle.dump(model,open("C:/Users/lb/Desktop/test/tl.dat","wb")) 
# 可以进行读
# 读取模型,重新打开jupyter
import pickle
load_file = open("C:/Users/lb/Desktop/test/tl.dat","rb")
load_game_data = pickle.load(load_file)
load_file.close()
load_game_data

在这里插入图片描述
接着预测

load_game_data.predict(train_x)

ItVuer - 免责声明 - 关于我们 - 联系我们

本网站信息来源于互联网,如有侵权请联系:561261067@qq.com

桂ICP备16001015号