php网站开发作业,做注册会计师网站,上海小程序定制公司,厦门网站建设报#x1f490;大家好#xff01;我是码银~#xff0c;欢迎关注#x1f490;#xff1a; CSDN#xff1a;码银 公众号#xff1a;码银学编程 实验目的和要求
会用Python创建Kmeans聚类分析模型使用KMeans模型对航空公司客户价值进行聚类分析会对聚类结果进行分析评价
实… 大家好我是码银~欢迎关注 CSDN码银 公众号码银学编程 实验目的和要求
会用Python创建Kmeans聚类分析模型使用KMeans模型对航空公司客户价值进行聚类分析会对聚类结果进行分析评价
实验环境
pycharm2020Win11Python3.7Anaconda2019
KMeans聚类算法简介
KMeans聚类算法是一种基于中心点的聚类方法其目标是将数据点划分为K个簇使得每个簇内的数据点与簇中心的距离之和最小。算法的基本步骤包括
初始化随机选择K个数据点作为初始簇中心。分配将每个数据点分配到最近的簇中心形成K个簇。更新重新计算每个簇的中心点。迭代重复步骤2和3直到簇中心不再变化或达到最大迭代次数。
数据的加载和分析
数据集的获取搜索微信公众号“码银学编程”。回复航空数据集 ZL入会至当前时长反映客户的活跃时间。 ZR最近消费时间间隔反映客户的最近活跃程度。 ZF消费频次反映客户的忠诚度。 ZM消费里程总额反映客户对航空公司服务的依赖程度。 ZC舱位等级对应折扣系数通常舱位等级越高折扣系数越大。
首先使用Pandas库加载CSV格式的环境监测数据文件。
def load_data(filepath):加载CSV数据文件return pd.read_csv(filepath, header0)聚类分析
接着使用Scikit-learn库中的KMeans模型对数据进行聚类分析。通过设置不同的参数如最大迭代次数、簇的数量等可以对模型进行调整以适应不同的数据集。
def perform_kmeans(data, n_clusters):执行KMeans聚类分析model KMeans(max_iter300, n_clustersn_clusters, random_stateNone, tol0.0001)model.fit(data)return model结果可视化
为了直观展示聚类结果使用Matplotlib库绘制聚类图。通过将数据点和簇中心在二维平面上表示可以清晰地观察到数据的分布和簇的划分情况。
def plot_clusters(model, data):绘制聚类结果plt.figure(figsize(10, 6)) # 设置图表大小plt.xlabel(ZL-ZR-ZF-ZM-ZC) # 假设环境指标plt.ylabel(Cluster-center-value)plt.title(聚类分析结果图)colors [r, g, y, b, k]for i in range(model.n_clusters):plt.plot(data.columns, model.cluster_centers_[i], labelfCluster {i}, colorcolors[i], markero)plt.legend()plt.grid(True)plt.show()主函数
def main():# 加载数据data load_data(air_data.csv)print(数据形状, data.shape)print(数据前五行)print(data.head())# 聚类分析kmodel perform_kmeans(data, 5)print(聚类类别数目统计)print(pd.Series(kmodel.labels_).value_counts())# 聚类中心cluster_centers pd.DataFrame(kmodel.cluster_centers_)print(聚类中心)print(cluster_centers)# 聚类中心与类别数目cluster_info pd.concat([cluster_centers, pd.Series(kmodel.labels_).value_counts()], axis1)cluster_info.columns list(data.columns) [Cluster Count]print(聚类中心与类别数目)print(cluster_info)# 绘制聚类结果图plot_clusters(kmodel, data)if __name__ __main__:main()分析与讨论
这个结果展示了使用K-Means聚类算法对航空公司客户数据进行分析后得到的聚类中心和每个聚类的样本数量。每一列ZL、ZR、ZF、ZM、ZC代表数据集中的一个特征这些特征分别表示
ZL入会至当前时长反映客户的活跃时间。ZR最近消费时间间隔反映客户的最近活跃程度。ZF消费频次反映客户的忠诚度。ZM消费里程总额反映客户对航空公司服务的依赖程度。ZC舱位等级对应折扣系数通常舱位等级越高折扣系数越大。
聚类中心Cluster Centers是每个聚类中所有点的均值可以看作是该聚类的“代表”或“典型”客户。在这个例子中我们有5个聚类中心和它们的统计数据 第一个聚类中心Cluster 0的ZL值较低ZR值较高ZF和ZM值较低ZC值也较低。这可能代表一群活跃时间较短、最近消费间隔较长、消费频次和里程较低的客户他们可能对航空公司的忠诚度和依赖程度不高。 第二个聚类中心Cluster 1的ZL值较高ZR值较低ZF值较低ZM值较低ZC值较低。这可能代表一群活跃时间较长但最近不太活跃的客户他们的消费频次和里程也较低。 第三个聚类中心Cluster 2的ZL和ZR值都较低ZF值较低ZM值较低ZC值较高。这可能代表一群活跃时间较短且最近消费间隔较长的客户他们的消费频次和里程较低但可能购买了较高舱位等级的机票。 第四个聚类中心Cluster 3的ZL和ZR值都较高ZF和ZM值较高ZC值也较高。这可能代表一群活跃时间较长、最近消费频繁、消费里程高且购买了较高舱位等级机票的客户他们对航空公司的忠诚度和依赖程度很高。 第五个聚类中心Cluster 4的ZL值较低ZR值较低ZF值较低ZM值较低ZC值较高。这可能代表一群最近活跃且购买了较高舱位等级机票的客户但他们的总体消费频次和里程较低。
完整代码
import pandas as pd
from matplotlib import pyplot as plt
from sklearn.cluster import KMeansdef load_data(filepath):加载CSV数据文件return pd.read_csv(filepath, header0)def perform_kmeans(data, n_clusters):执行KMeans聚类分析model KMeans(max_iter300, n_clustersn_clusters, random_stateNone, tol0.0001)model.fit(data)return modeldef plot_clusters(model, data):绘制聚类结果plt.figure(figsize(10, 6)) # 设置图表大小plt.xlabel(ZL-ZR-ZF-ZM-ZC)plt.ylabel(Cluster-center-value)plt.rcParams[font.sans-serif] [SimHei] # 确保中文标签正常显示plt.title(聚类分析结果图)cluster_centers model.cluster_centers_colors [r, g, y, b, k]for i in range(len(cluster_centers)):plt.plot(data.columns, cluster_centers[i], labelfCluster {i}, colorcolors[i], markero)plt.legend()plt.grid(True) # 添加网格线plt.show()def main():# 加载数据data load_data(air_data.csv)print(数据形状, data.shape)print(数据前五行)print(data.head())# 聚类分析kmodel perform_kmeans(data, 5)print(聚类类别数目统计)print(pd.Series(kmodel.labels_).value_counts())# 聚类中心cluster_centers pd.DataFrame(kmodel.cluster_centers_)print(聚类中心)print(cluster_centers)# 聚类中心与类别数目cluster_info pd.concat([cluster_centers, pd.Series(kmodel.labels_).value_counts()], axis1)cluster_info.columns list(data.columns) [Cluster Count]print(聚类中心与类别数目)print(cluster_info)# 绘制聚类结果图plot_clusters(kmodel, data)if __name__ __main__:main()
