三个特征变量怎么聚类python

Posted 2023-04-24

参考技术A 使用聚类代码。三个特征变量是可以使用聚类代码进行聚类的，Python由荷兰数学和计算机科学研究学会的吉多·范罗苏姆于1990年代初设计，作为一门叫做ABC语言的替代品。

Python使用RMF聚类分析客户价值

投资机构或电商企业等积累的客户交易数据繁杂。需要根据用户的以往消费记录分析出不同用户群体的特征与价值，再针对不同群体提供不同的营销策略。

用户分析指标

根据美国数据库营销研究所Arthur Hughes的研究，客户数据库中有三个神奇的要素，这三个要素构成了数据分析最好的指标

    R-最近一次消费(Recency)

    F-消费频率(Frequency)

    M-消费金额(Monetary)

通过该图将用户进行分类：

    R、F、M都很高，重要价值客户(VIP客户)

    F、M很高，R不高，重要保持客户

    R、F、M都很低，流失客户

    M很高，R、F不高，重要挽留客户

 

根据这8个类别的R、F、M指标，对用户进行标注，哪些是重要价值客户，哪些是重要保持客户，哪些是重要发展客户，哪些是流失客户等

 

流程介绍

以R、F、M这三个核心指标为维度进行聚类分析

利用K-means聚类分析将用户分类

根据R、F、M指标，对用户进行标注

 

准备工作：
数据：

     某电商企业客户近期购买的数据。包含客户注册日期，最后购买日期以及购买消费总金额

 

参数：

R-求出最近一次投资时间距提数日天数

F-月均投资次数

M-月均投资金额

目标：分析客户交易数据，用户群体的特征与价值，进行精准营销，降低营销成本，提高销售业绩。

1  分析数据获取RFM

R-求出最近一次投资时间距提数日天数

    确定一个提现日，减去用户的最新投资日期

F-月均投资次数

    总投资次数/总月数

M-月均投资金额

    投资总金额/总月数

 

 

处理数据获取R-F-M

def dataChange(data):
    deadline_time = datetime(2016,7,20)
    print(deadline_time)

    # 时间相减 得到天数查 timedelta64类型
    diff_R = deadline_time - data["最近一次投资时间"]

    # 渠道具体天数
    # days = diff_R[0].days
    R = []
    for i in diff_R:
        days = i.days
        R.append(days)

    print(R)
    \'\'\'
    用户在投时长(月
    Python没有直接获取月数差的函数
    1、获取用户在投天数
    2、月=在投天数/30，向上取整
    \'\'\'
    diff = deadline_time - data["首次投资时间"]
    print(diff)

    # 利用向上取整函数
    months = []
    for i in diff:
        month = ceil(i.days/30)
        months.append(month)

    print(months)

    # 月均投资次数
    month_ave = data["总计投标总次数"]/months
    F = month_ave.values
    print(F)

    # 月均投资金额
    M = (data["总计投资总金额"]/months).values
    print(M)

    return R, M, F

 2 训练KMeans模型

先对数据进行转换，然后通过K—Means模型训练，生产模型

def analy_data(data, R, M, F):
    cdata = DataFrame([R, list(F), list(M)]).T
    # 指定cdata的index和colums
    cdata.index = data.index
    cdata.columns = ["R-最近一次投资时间距提数日的天数", "F-月均投资次数", "月均投资金额"]
    print("cdata_info:\\n", cdata)

    print("cdata:\\n", cdata.describe())

    # K-Means聚类分析

    # 01 数据标准化  均值：cdata.mean()   标准差：cdata.std()
    # 对应位置分别先相减 再相除
    zcdata = (cdata-cdata.mean())/cdata.std()
    print("zcdata:\\n", zcdata)

    # n_clusters：分类种数  n_jobs：计算的分配资源  max_iter：最大迭代次数  random_state：随机数种子，种子相同，参数固定
    kModel = KMeans(n_clusters=4, n_jobs=4, max_iter=100, random_state=0)
    kModel.fit(zcdata)
    print(kModel.labels_)

 

3  通过模型对用户标注

   # 统计每个类别的频率
    value_counts = Series(kModel.labels_).value_counts()
    print(value_counts)

    # 将类别标签赋回原来的数据
    cdata_rst = pd.concat([cdata, Series(kModel.labels_, index=cdata.index)], axis=1)
    print(cdata_rst)

    # 命名最后一列为类别
    cdata_rst.columns = list(cdata.columns) + ["类别"]
    print(cdata_rst)

    # 按照类别分组统计R, F, M的指标均值
    user_ret = cdata_rst.groupby(cdata_rst["类别"]).mean()
    print(user_ret)

 

\'\'\'

        R-最近一次投资时间距提数日的天数   F-月均投资次数         月均投资金额
 类别
  0            27.859375               2.820312          21906.754297
  1            20.684211               4.552632          115842.105263
  2            10.568182               5.579545          26984.313636
  3            12.111111               17.277778         107986.000000
  
  结论：
  类别3：R、F、M都比较高，属于重要价值客户 或 超级用户
  类别0：R、F、M都比较低，属于低价值客户
  类别1：R一般、F一般、M很高，也属于重要价值客户
  
  \'\'\'

通过模型对新用户标注

1、获取新用户数据

2、通过和原数据处理获取RFM

3、通过训练模型得出用户类型

def user_classes(cdata, user_info):
    \'\'\'
    # 模拟一条用户数据
    1、获取当前时间表示为截止时间
    2.计算出: R F M

    \'\'\'
    R, M, F = user_info_change(user_info)
    user_data_info = DataFrame([[R], [F], [M]]).T
    print(user_data_info)

    # user_data_info = DataFrame([[12.5], [18.0], [20000.0]]).T
    user_data_info.index = ["lily"]
    user_data_info.columns = cdata.columns
    print("cdata_info:\\n", user_data_info)

    new_zcdata = (user_data_info-cdata.mean())/cdata.std()
    print("new_zcdata", new_zcdata)

    kModel = load_model("user_classes.pkl")
    ret = kModel.predict(new_zcdata)
    print("new_zcdata_ret:", ret)
    # new_zcdata_ret: [3]