aic python_python机器学习案例（五）：逻辑回归电信客户流失建模

最新推荐文章于 2024-06-06 01:00:00 发布

望止洋

最新推荐文章于 2024-06-06 01:00:00 发布

阅读量924

点赞数

文章标签： aic python

本文链接：https://blog.csdn.net/weixin_36237278/article/details/112441102

版权

本文介绍了使用Python进行逻辑回归分析的电信客户流失建模案例，涉及数据集获取、字段解释、混淆矩阵、ROC曲线、特征选择与多重共线性处理。通过传统建模和机器学习方法，得出最佳模型并剔除相关变量，如'nrProm', 'avgplan'和'posTrend'。最终模型表现良好，展示了机器学习在客户流失预测中的应用。" 113508635,10548454,Python假设检验详解：U检验、T检验、卡方检验与F检验,"['Python', '统计学', '数据分析', '假设检验']

摘要由CSDN通过智能技术生成

**本逻辑回归电信客户流失建模案例代码及数据集链接及下载密码：关注公众号书豪创投笔记并回复python数据科学即可获取**

关注公众号书豪创投笔记并回复python数据科学即可获取数据集

本案例来源乃是学习朋友常国珍老师的python数据科学书，内容非常实用！希望深入浅出学习机器学习的同学可学习北大博士常国珍建模实战课程系列常国珍录制的Python 数据科学机器学习教学视频，可以由浅入深的学习机器学习！

逻辑回归案例字段解释
#subscriberID="个人客户的ID"
#churn="是否流失：1=流失";
#Age="年龄"
#incomeCode="用户居住区域平均收入的代码"
#duration="在网时长"
#peakMinAv="统计期间内最高单月通话时长"
#peakMinDiff="统计期间结束月份与开始月份相比通话时长增加数量"
#posTrend="该用户通话时长是否呈现出上升态势：是=1"
#negTrend="该用户通话时长是否呈现出下降态势：是=1"
#nrProm="电话公司营销的数量"
#prom="最近一个月是否被营销过：是=1"
#curPlan="统计时间开始时套餐类型：1=最高通过200分钟；2=300分钟；3=350分钟；4=500分钟"
#avPlan="统计期间内平均套餐类型"
#planChange="统计结束时和开始时套餐的变化：正值代表套餐档次提升，负值代表下降，0代表不变"
#posPlanChange="统计期间是否提高套餐：1=是"
#negPlanChange="统计期间是否降低套餐：1=是"
#call_10086="拨打10086的次数"

import os
import numpy as np
from scipy import stats
import pandas as pd
import statsmodels.api as sm
import statsmodels.formula.api as smf
import matplotlib.pyplot as plt

churn = pd.read_csv(r'telecom_churn.csv', skipinitialspace=True)
churn.head()

#1两变量分析：检验该用户通话时长是否呈现出上升态势(posTrend)对流失(churn) 是否有预测价值
# ##  分类变量的相关关系
#
# 交叉表

cross_table = pd.crosstab(churn.posTrend, 
                         churn.churn, margins=True)
cross_table

交叉表

# 列联表

def percConvert(ser):
    return ser/float(ser[-1])

cross_table.apply(percConvert, axis=1)


#卡方检验走一波
print('''chisq = %6.4f 
p-value = %6.4f
dof = %i 
expected_freq = %s'''  %stats.chi2_contingency(cross_table.iloc[:2, :2]))

卡方检验走一波告诉俺们显著的！

#2首先将原始数据拆分为训练和测试数据集，使用训练数据集建立在网时长对流失的逻辑回归，使用测试数据集制作混淆矩阵
#（阈值为0.5），提供准确性、召回率指标，提供ROC曲线和AUC。
# ## 逻辑回归
# duration  在网时长  这是个连续变量 ！

churn.plot(x='duration', y='churn', kind='scatter')