航空公司客户价值分析
目标：企业针对不同价值的客户制定个性化的服务，将有限的资源集中于高价值客户。
1、借助航空公司的数据进行客户分类
2、比较不同类客户的价值并制定销策略
识别客户价值的最广泛的模型是通过RFM模型来识别出高价值的客户：
Recency： 最近消费时间间隔
Frequency： 消费频率
Monetary： 消费金额
问题：同样的消费金额的不同旅客对航空公司的价值不同，例如买长航线、低等仓的旅客和买短航线、高等仓的旅客消费金额相同
但是价值确实不同的。显然后者更有价值。因此这个指标可能不合适，故选择客户在一定时间内的飞行里程M和乘坐舱位所对应的折扣洗漱C。
航空公司客户价值分析模型：
L：会员入会时间距观测窗口结束的时间
R：客户最近一次乘坐公司分级距观测窗口结束的时间（月数）
F：客户在观测窗口内乘坐公司飞机的次数
M：客户在观测窗口内累计的飞行里程
C：客户在观测窗口内乘坐舱位所对应的折扣系数的平均值
模型概述：
通过使用k-means聚类算法没使用LRFMC模型进行测试
对模型进行分析，分析客户的价值
import pandas as pddatafile = 'data/air_data.csv'
resultfile = 'data/tmp.xls'
data = pd.read_csv(datafile, encoding='utf-8')
explore = data.describe().T
explore['null'] = len(data) - explore['count'] #计算空值
explore = explore[['null', 'max', 'min']] #得到每一列的空值，最大值，最小值
explore.to_excel(resultfile) # 是结果写入文件
根据上面的数据统计，丢弃所有不符合的数据，
1、票价为空的
2、票价为0，但是折扣不是0，而且飞行里程大于0，
这样的数据是错误数据，直接删除。
data = pd.read_csv(datafile, encoding='utf-8')
cleanedfile = 'data/cleanfile.csv' #存清理后的数据
# 去掉票价是空的元素
data = data[data['SUM_YR_1'].notnull()]
data = data[data['SUM_YR_2'].notnull()]
# 去掉票价为0，但是折扣不是0，而且飞行里程大于0
data = data.drop(data['SUM_YR_1'] ==0 & (data['SEG_KM_SUM'] != 0) | (data['avg_discount'] > 0))
print(len(data))
# print(data.columns)
# 属性规约：去掉不相管的属性，只留下与LRFMC模型相关的属性
# FFP_DATE 入会时间
# LOAD_TIME 观测窗口结束时间
# FLIGHT_COUNT 飞行频率
# avg_discount 平均折扣
# SEG_KM_SUM 观测窗口总飞行公里数
# LAST_TO_END 最后一次乘机时间至观察窗口末端时长
data = data[['FFP_DATE','LOAD_TIME', 'FLIGHT_COUNT', 'avg_discount', 'SEG_KM_SUM','LAST_TO_END']]
data.to_csv(cleanedfile) #写入文件数据变化的LRFMC数据：
L = LOAD_TIME - FFP_DATE （观测窗口时间 - 入会时间）
R = LOAD_TIME - LAST_TO_END （观测窗口时间 - 最后一次乘机时间）
F = FLIGHT_COUNT
M = SEG_KM_SUM
C = avg_discount
from datetime import datetime
import time
def normal_time(date):
'''
格式化数据
'''
return datetime.strptime(date,"%Y/%m/%d")
def interval_time(dd):
'''
计算时间间隔，以月为单位
'''
return dd.days / 30# data_LRFMC数据
data_LRFMC = pd.DataFrame()
# data_LRFMC.columns = ['L', 'R', 'F','M', 'C']
data_LRFMC['L'] = (data['LOAD_TIME'].apply(normal_time) - data['FFP_DATE'].apply(normal_time)).apply(interval_time)
data_LRFMC['R'] = data['LAST_TO_END']
data_LRFMC['F'] = data['FLIGHT_COUNT']
data_LRFMC['M'] = data['SEG_KM_SUM']
data_LRFMC['C'] = data['avg_discount']
# 显示数据的描述，最大值和最小值
data_LRFMC_describe = data_LRFMC.describe().T
data_LRFMC_describe = data_LRFMC_describe[['max','min']].T
data_LRFMC.to_csv('data/LRFMC.csv') #数据写入文件注意到所有的数据值，最大值和最小值间隔较大，需要对数据进行标准化
这里使用标注差标准化,可以使用认识方法进行处理
# 标准化、重命名、写入文件
data_normal = (data_LRFMC - data_LRFMC.mean()) / (data_LRFMC.std())
data_normal.columns = ['Z'+i for i in data_normal.columns]
data_normal.to_csv('data/data_normal')
数据处理完毕，下面进行模型的构建，
1、使用聚类算法，将数据生成5类用户
2、针对聚类结果进行特征分析
from sklearn.cluster import KMeans
k = 5
kmodel = KMeans(k) #得到模型
kmodel.fit(data_normal) #训练模型
# 查看聚类中心和对应的类别
print(kmodel.cluster_centers_)
print(kmodel.labels_)
[[ 1.16443082 -0.3775831 -0.08518766 -0.09323092 -0.15906231]
[-0.70014625 -0.41584827 -0.15935193 -0.15955595 -0.26734394]
[ 0.01594116 0.00784598 -0.24843215 -0.25750534 2.08918975]
[ 0.48503896 -0.80134848 2.48281809 2.43536635 0.31393202]
[-0.31252471 1.68746582 -0.57426021 -0.53886494 -0.18594193]]
[3 3 3 ..., 1 0 4]
import matplotlib.pyplot as plt
clu = kmodel.cluster_centers_
x = [1,2,3,4,5]
colors = ['red','green','yellow','blue','black']
for i in range(5):
plt.plot(x,clu[i],label='clustre '+str(i),linewidth=6-i,color=colors[i],marker='o')
plt.xlabel('L R F M C')
plt.ylabel('values')
plt.show()
客户价值分析：
注意kmeans每次运行的时候得到的类会有差别，簇号也会相应的改变，但是中间点基本不会改变：
L R F M C
客户群1：[ 1.16443082 -0.3775831 -0.08518766 -0.09323092 -0.15906231] red
客户群2：[-0.70014625 -0.41584827 -0.15935193 -0.15955595 -0.26734394] green
客户群3：[ 0.01594116 0.00784598 -0.24843215 -0.25750534 2.08918975] yellow
客户群4：[ 0.48503896 -0.80134848 2.48281809 2.43536635 0.31393202] blue
客户群5：[-0.31252471 1.68746582 -0.57426021 -0.53886494 -0.18594193] black
我们重点关注的是L，F，M，从图中可以看到：
1、客户群4[blue] 的F,M很高，L也不低，可以看做是重要保持的客户；
2、客户群3[yellow] 重要发展客户
3、客户群1[red] 重要挽留客户，原因：入会时间长，但是F,M较低
4、客户群2[green] 一般客户
5、客户群5[black] 低价值客户