航空公司客户价值分析
data_explore
#-*- coding: utf-8 -*-
# 代码7-1
# 对数据进行基本的探索
# 返回缺失值个数以及最大最小值
import pandas as pd
datafile= 'E:/Third work/数据代码/chapter7/demo/data/air_data.csv' # 航空原始数据,第一行为属性标签
resultfile = 'E:/Third work/数据代码/chapter7/demo/tmp/explore_2.csv' # 数据探索结果表
# 读取原始数据,指定UTF-8编码(需要用文本编辑器将数据装换为UTF-8编码)
data = pd.read_csv(datafile, encoding = 'utf-8')
# 包括对数据的基本描述,percentiles参数是指定计算多少的分位数表(如1/4分位数、中位数等)
explore = data.describe(percentiles = [], include = 'all').T # T是转置,转置后更方便查阅
explore
explore['null'] = len(data)-explore['count'] # describe()函数自动计算非空值数,需要手动计算空值数
explore
explore = explore[['null', 'max', 'min']]
explore.columns = ['空值数', '最大值', '最小值'] # 表头重命名
explore

explore.to_csv(resultfile) # 导出结果
data_distribution
#-*- coding: utf-8 -*-
# 代码7-2
# 对数据的分布分析
import pandas as pd
import matplotlib.pyplot as plt
datafile= 'E:/Third work/数据代码/chapter7/demo/data/air_data.csv' # 航空原始数据,第一行为属性标签
# 读取原始数据,指定UTF-8编码(需要用文本编辑器将数据装换为UTF-8编码)
data = pd.read_csv(datafile, encoding = 'utf-8')
# 客户信息类别
# 提取会员入会年份
from datetime import datetime
ffp = data['FFP_DATE'].apply(lambda x:datetime.strptime(x,'%Y/%m/%d'))
ffp_year = ffp.map(lambda x : x.year)
ffp

ffp_year

# 绘制各年份会员入会人数直方图
fig = plt.figure(figsize = (8 ,5)) # 设置画布大小
plt.rcParams['font.sans-serif'] = 'SimHei' # 设置中文显示
plt.rcParams['axes.unicode_minus'] = False
plt.hist(ffp_year, bins='auto', color='red')
plt.xlabel('年份')
plt.ylabel('入会人数')
plt.title('各年份会员入会人数(20信计1班黄美静2019310143014)')
plt.savefig('E:/Third work/数据代码/chapter7/demo/imag/各年份会员入会人数直方图.jpg',dpi=800)
plt.show()
# plt.close

# 提取会员不同性别人数
male = pd.value_counts(data['GENDER'])['男']
female = pd.value_counts(data['GENDER'])['女']
# 绘制会员性别比例饼图
fig = plt.figure(figsize = (7 ,4)) # 设置画布大小
plt.pie([ male, female], labels=['男','女'], colors=['lightskyblue', 'lightcoral'],
autopct='%1.1f%%')
plt.title('会员性别比例(20信计1班黄美静2019310143014)')
plt.savefig('E:/Third work/数据代码/chapter7/demo/imag/性别饼图.jpg',dpi=800)
plt.show()
# plt.close

# 提取不同级别会员的人数
lv_four = pd.value_counts(data['FFP_TIER'])[4]
lv_five = pd.value_counts(data['FFP_TIER'])[5]
lv_six = pd.value_counts(data['FFP_TIER'])[6]
# 绘制会员各级别人数条形图
fig = plt.figure(figsize = (8 ,5)) # 设置画布大小
plt.bar(x=range(3), height=[lv_four,lv_five,lv_six], width=0.4, alpha=0.8, color='green')
plt.xticks([index for index in range(3)], ['4','5','6'])
plt.xlabel('会员等级')
plt.ylabel('会员人数')
plt.title('会员各级别人数(20信计1班黄美静2019310143014)')
plt.savefig('E:/Third work/数据代码/chapter7/demo/imag/会员级别柱状图.jpg',dpi=800)
plt.show()
# plt.close()

# 提取会员年龄
age = data['AGE'].dropna()
age = age.astype('int64')
# 绘制会员年龄分布箱型图
fig = plt.figure(figsize = (5 ,10))
plt.boxplot(age,
patch_artist=True,
labels = ['会员年龄'], # 设置x轴标题
boxprops = {'facecolor':'yellow'}) # 设置填充颜色
plt.title('会员年龄分布箱线图(20信计1班黄美静2019310143014)')
# 显示y坐标轴的底线
plt.grid(axis='y')
plt.savefig('E:/Third work/数据代码/chapter7/demo/imag/会员年龄箱型图.jpg',dpi=800)
plt.show()
# plt.close

# 代码7-3
# 乘机信息类别
lte = data['LAST_TO_END']
fc = data['FLIGHT_COUNT']
sks = data['SEG_KM_SUM']
# 绘制最后乘机至结束时长箱线图
fig = plt.figure(figsize = (5 ,8))
plt.boxplot(lte,
patch_artist=True,
labels = ['时长'], # 设置x轴标题
boxprops = {'facecolor':'pink'}) # 设置填充颜色
plt.title('会员最后乘机至结束时长分布箱线图(20信计1班黄美静2019310143014)')
# 显示y坐标轴的底线
plt.grid(axis='y')
plt.savefig('E:/Third work/数据代码/chapter7/demo/imag/最后乘机至结束时长箱线图.jpg',dpi=800)
plt.show()
# plt.close

# 绘制客户飞行次数箱线图
fig = plt.figure(figsize = (5 ,8))
plt.boxplot(fc,
patch_artist=True,
labels = ['飞行次数'], # 设置x轴标题
boxprops = {'facecolor':'purple'}) # 设置填充颜色
plt.title('会员飞行次数分布箱线图(20信计1班黄美静2019310143014)')
# 显示y坐标轴的底线
plt.grid(axis='y')
plt.savefig('E:/Third work/数据代码/chapter7/demo/imag/客户飞行次数箱线图.jpg',dpi=800)
plt.show()
# plt.close

# 绘制客户总飞行公里数箱线图
fig = plt.figure(figsize = (5 ,10))
plt.boxplot(sks,
patch_artist=True,
labels = ['总飞行公里数'], # 设置x轴标题
boxprops = {'facecolor':'orange'}) # 设置填充颜色
plt.title('客户总飞行公里数箱线图(20信计1班黄美静2019310143014)')
# 显示y坐标轴的底线
plt.grid(axis='y')
plt.savefig('E:/Third work/数据代码/chapter7/demo/imag/客户总飞行公里数箱线图.jpg',dpi=800)
plt.show()
# plt.close

# 代码7-4
# 积分信息类别
# 提取会员积分兑换次数
ec = data['EXCHANGE_COUNT']
# 绘制会员兑换积分次数直方图
fig = plt.figure(figsize = (8 ,5)) # 设置画布大小
plt.hist(ec, bins=5, color='#0504ab')
plt.xlabel('兑换次数')
plt.ylabel('会员人数')
plt.title('会员兑换积分次数分布直方图(20信计1班黄美静2019310143014)')
plt.savefig('E:/Third work/数据代码/chapter7/demo/imag/会员兑换积分次数直方图.jpg',dpi=800)
plt.show()
# plt.close

# 提取会员总累计积分
ps = data['Points_Sum']
# 绘制会员总累计积分箱线图
fig = plt.figure(figsize = (5 ,8))
plt.boxplot(ps,
patch_artist=True,
labels = ['总累计积分'], # 设置x轴标题
boxprops = {'facecolor':'lightgreen'}) # 设置填充颜色
plt.title('客户总累计积分箱线图(20信计1班黄美静2019310143014)')
# 显示y坐标轴的底线
plt.grid(axis='y')
plt.savefig('E:/Third work/数据代码/chapter7/demo/imag/会员总累计积分箱线图.jpg',dpi=800)
plt.show()
# plt.close

# 代码7-5
# 提取属性并合并为新数据集
data_corr = data[['FFP_TIER','FLIGHT_COUNT','LAST_TO_END',
'SEG_KM_SUM','EXCHANGE_COUNT','Points_Sum']]
age1 = data['AGE'].fillna(0)
data_corr['AGE'] = age1.astype('int64')
data_corr['ffp_year'] = ffp_year
# 计算相关性矩阵
dt_corr = data_corr.corr(method = 'pearson')
print('相关性矩阵为:\n',dt_corr)
# 绘制热力图
import seaborn as sns
plt.subplots(figsize=(10, 10)) # 设置画面大小
sns.heatmap(dt_corr, annot=True, vmax=1, square=True, cmap='pink')
plt.title("热力图(20信计1班黄美静2019310143014)")
plt.savefig('E:/Third work/数据代码/chapter7/demo/imag/热力图.jpg',dpi=600)
plt.show()
# plt.close


浙公网安备 33010602011771号