人工智能-聚类算法、误差评估、算法优化、特征降维

1、聚类算法简介

典型的无监督算法,主要用于将相似的样本自动归到一个类别中。
根据样本之间的相似性,将样本划分到不同的类别中,对于不同的相似度计算方法,会得到不同的聚类结果。常用的相似度计算方法是欧式距离法

聚类算法与分类算法的区别:
聚类算法是无监督的学习算法,而分类算法属于监督学习。

2、聚类算法API

API:sklearn.cluster.KMeans(n_clusters=8)
参数:n_clusters要聚类的中心数,默认值8
调用方法:estimator.fit(x)、estimator.predict(x)、estimator.fit_predect(x)-指先fit再预测

2.1案例

import matplotlib.pyplot as plt
#from sklearn.datasets.samples_generator import make_blobs
from sklearn.cluster import KMeans
#from sklearn.metrics import calinski_harabaz_score

#创建数据
x,y = make_bloks(n_samples=1000,n_features=2,centers=[[-1,-1],,[0,0],[1,1],[2,2]],cluster_std=[0.4,0.2,0.2,0.2],random_state=9)
#生成1000个数,两个特征值,4个中心点,标准差越大越离散,
plt.scatter(x[:,0],x[:,1],marker="o")
plt.show()

在这里插入图片描述

#kmeans训练且可视化
y_pre = KMeans(n_clusters=2,random_stats=9).fit_predict(x)

#可视化展示
plt.scatter(x[:,0],x[:,1],c=y_pre)
plt.show()
#用cn_sccie查看最后效果
print(calinski_harabaz_score[x,y_pre])  #评估函数越大越好,输出:3116

在这里插入图片描述

3、聚类算法实现流程

3.1思路

K-初始中心点个数
means - 中心点到其他数据点距离的平均值

流程:

  1. 随机设置k个特征空间内的点作为初始的聚类中心
  2. 对于其他每个点,计算到这k个中心的距离,最近的中心作为标记类别
  3. 对已经标记好类别后,重新计算出每个聚类的心中心点(平均值)
  4. 如果计算得到的新中心点与原中心点一样则结束,否则重新回到第二步
    在这里插入图片描述
    在大规模数据集上,K-Means的收敛速度比较慢
    K-Means一点会停止迭代的

4、误差评估

4.1误差平方和

K-Means中的SSE:
在这里插入图片描述
在这里插入图片描述
SSE越小越好
缺点:

4.2 肘方法–确定K值

横轴:k值,纵轴:方差
取拐点的k值
在这里插入图片描述

4.3 轮廓系数法

在这里插入图片描述
S越接近1,效果越好。越接近-1,效果越不好
使得内部距离最小化,外部距离最大化

4.4 CH系数

迹,对角线元素之和=秩的和
在这里插入图片描述
S(K)越大聚类越好

要达到的目的:用尽量少的类别聚类尽量多的样本,同时获得好的聚类效果

5、算法优化

5.1KMeans算法小结

优点

  1. 原理简单,实现容易
  2. 聚类效果中上(依赖K的选择)
  3. 空间复杂度o(N),时间复杂度o(lkN)(N为样本点个数,K中心点个数,I迭代次数)
    缺点
  4. 对离群点,噪声敏感
  5. 结果不一定全局最优,只能保证局部最优(与K值和初始点的选择有关)
  6. 很难发现大小差别很大的簇及进行增量计算

5.2 改进1-Canopy算法

初始中心点的选择:避免离的太近
确定初始点的选择步骤:
先随机一个,然后在外围再选一个
在这里插入图片描述

5.3 改进2-K-means++

目的也是让初始中心点尽可能远
在这里插入图片描述

5.4 改进3- 二分k-means

步骤:
1,将所有点作为一个簇
2,将该簇一分为二
3,选择能最大限度降低聚类代价函数(即误差平方和)的簇划分为两个簇
4,依次进行下去,直到簇的数目等于用户给定的数目K为止

5.5 改进4-K-medoids

与k-means的区别在于:确定新点时,不是取的平均值,而是取的到其他点的距离之和最小的点
因此,计算量剧增

在这里插入图片描述

6、特征降维

6.1 降维定义

降维:在某些特定限定条件下,降低随机变量个数,得到一组"不相关"的主变量的过程

降维的两种方式:特征选择、主成分分析

6.2 特征选择

数据中包含冗余或无关变量(或称特征、属性、指标等),从原特征中找出主要特征
方法:

  • Filter过滤式:主要探究特征本身特点、特征与特征和目标值之间关联(方差选择)
  • Embedded(嵌入式):算法自动选择特征

API :
sklearn.feature_selection.VarianceThreshold(threshold = 0.0)
删除所有低方差特征

#删除方差小的列
import pandas as pd
from sklearn.feature_selection import VarianceThreshold
def var_thr():
	data = pd.read_csv("./data/factor_returns.csv")
	print(data.shape)
	transfer = VarianceThreshold(threshold=1)
	transfer_data = transfer.fit_transform(data.iloc[:,1:10])
	print(transfer_data)
	print(data.iloc[:,1:10].shape)
	print(transfer_data.shape)

if __name__ == '__main__':
	var_thr()

相关系数

皮尔逊相关系数

from scipy.stats import pearsonr
def pea_demo():
	x1 = [12.5,13,11.5,11.9,33.5]
	x2 = [12,13,11,11,32]
	ret = pearsonr(x1,x2)
	print("皮尔逊相关系数结果是:\n",ret)

if __name__ == '__main__':
	pea_demo()

斯皮尔曼相关系数

from scipy.stats import pearsonr
def spea_demo():
	x1 = [12.5,13,11.5,11.9,33.5]
	x2 = [12,13,11,11,32]
	ret = spearmanr(x1,x2)
	print("斯皮尔曼相关系数结果是:\n",ret)

if __name__ == '__main__':
	pea_demo()

6.3 主成分分析

高维数据转化为低维数据

def pca_demo():
	data = [[2,8,4,5],[6,3,0,8],[5,4,9,1]]
	transfer = PCA(n_components=0.9)
	trans_data = transfer.fit_transform(data)
	print('保留0.9的数据最后维度为:\n',trans_data)

	transfer2 = PCA(n_components=3)
	trans_data = transfer2.fit_transform(data)
	print('降到3维度的结果为:\n',trans_data)

if __name__ == '__main__':
	pca_demo()

7、探究用户对物品类别的喜好

import pandas as pd
from sklearn.decomposition import PCA
from sklearn.cluster import KMeans
from sklearn.metrics import silhouette_score
#获取数据
order_product = pd.read_csv('./order_products_prior.csv')
products = pd.read_csv('./products.csv')
orders = pd.read_csv('./order_products_prior.csv')
aisles = pd.read_csv('./aisles.csv')
order_product.head()
# 合并表格
table1 = pd.merge(order_product,products,on=['product_id','product_id'])
table2 = pd.merge(table1,orders,on=['order_id','order_id'])
table = pd.merge(table2,aisles,on=['aisles_id','aisles_id'])
table.shape

#交叉表合并
pd.crosstab(table['user_id'],table['aisle'])
data.shape
#数据截取
new_data= data[:1000]
#特征工程-pca
transfer = PCA(n_components=0.9)
trans_data = transfer.fit_transform(new_data)
trans_data.shape
#k-means
estimator = KMeans(n_clusters=5)
y_pre = estimator.fit_predict(trans_data)

#模型评估
silhouette_score(trans_data,y_pre)

  • 2
    点赞
  • 7
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值