首页 国际 注册加微信送体验金可提现_哈工大硕士生用 Python 实现了 11 种经典数据降维算法,源代码库已开放
注册加微信送体验金可提现_哈工大硕士生用 Python 实现了 11 种经典数据降维算法,源代码库已开放
发布时间: 2020-01-02 08:30:24 | 点击量: 1100
摘要:雷锋网 ai 开发者按:网上关于各种降维算法的资料参差不齐,同时大部分不提供源代码。这里有个 github 项目整理了使用 python 实现了 11 种经典的数据抽取算法,包括:pca、lda、mds、lle、tsne 等,并附有相关资料、展示效果;非常适合机器学习初学者和刚刚入坑数据挖掘的小伙伴。它由 karl pearson 在 1901 年提出,属于线性降维方法。

注册加微信送体验金可提现_哈工大硕士生用 Python 实现了 11 种经典数据降维算法,源代码库已开放

注册加微信送体验金可提现,雷锋网 ai 开发者按:网上关于各种降维算法的资料参差不齐,同时大部分不提供源代码。这里有个 github 项目整理了使用 python 实现了 11 种经典的数据抽取(数据降维)算法,包括:pca、lda、mds、lle、tsne 等,并附有相关资料、展示效果;非常适合机器学习初学者和刚刚入坑数据挖掘的小伙伴。

所谓降维,即用一组个数为 d 的向量 zi 来代表个数为 d 的向量 xi 所包含的有用信息,其中 d

通常,我们会发现大部分数据集的维度都会高达成百乃至上千,而经典的 mnist,其维度都是 64。

mnist 手写数字数据集

但在实际应用中,我们所用到的有用信息却并不需要那么高的维度,而且每增加一维所需的样本个数呈指数级增长,这可能会直接带来极大的「维数灾难」;而数据降维就可以实现:

使得数据集更易使用

确保变量之间彼此独立

降低算法计算运算成本

去除噪音

一旦我们能够正确处理这些信息,正确有效地进行降维,这将大大有助于减少计算量,进而提高机器运作效率。而数据降维,也常应用于文本处理、人脸识别、图片识别、自然语言处理等领域。

往往高维空间的数据会出现分布稀疏的情况,所以在降维处理的过程中,我们通常会做一些数据删减,这些数据包括了冗余的数据、无效信息、重复表达内容等。

例如:现有一张 1024*1024 的图,除去中心 50*50 的区域其它位置均为零值,这些为零的信息就可以归为无用信息;而对于对称图形而言,对称部分的信息则可以归为重复信息。

因此,大部分经典降维技术也是基于这一内容而展开,其中降维方法又分为线性和非线性降维,非线性降维又分为基于核函数和基于特征值的方法。

线性降维方法:

pca 、ica lda、lfa、lpp(le 的线性表示)

非线性降维方法:

基于核函数的非线性降维方法——kpca 、kica、kda

基于特征值的非线性降维方法(流型学习)——isomap、lle、le、lpp、ltsa、mvu

哈尔滨工业大学计算机技术专业的在读硕士生 heucoder 则整理了 pca、kpca、lda、mds、isomap、lle、tsne、autoencoder、fastica、svd、le、lpp 共 12 种经典的降维算法,并提供了相关资料、代码以及展示,下面将主要以 pca 算法为例介绍降维算法具体操作。

pca 是一种基于从高维空间映射到低维空间的映射方法,也是最基础的无监督降维算法,其目标是向数据变化最大的方向投影,或者说向重构误差最小化的方向投影。它由 karl pearson 在 1901 年提出,属于线性降维方法。与 pca 相关的原理通常被称为最大方差理论或最小误差理论。这两者目标一致,但过程侧重点则不同。

最大方差理论降维原理

将一组 n 维向量降为 k 维(k 大于 0,小于 n),其目标是选择 k 个单位正交基,各字段两两间 cov(x,y) 为 0,而字段的方差则尽可能大。因此,最大方差即使得投影数据的方差被最大化,在这过程中,我们需要找到数据集 xmxn 的最佳的投影空间 wnxk、协方差矩阵等,其算法流程为:

算法输入:数据集 xmxn;

按列计算数据集 x 的均值 xmean,然后令 xnew=x−xmean;

求解矩阵 xnew 的协方差矩阵,并将其记为 cov;

计算协方差矩阵 cov 的特征值和相应的特征向量;

将特征值按照从大到小的排序,选择其中最大的 k 个,然后将其对应的 k 个特征向量分别作为列向量组成特征向量矩阵 wnxk;

计算 xneww,即将数据集 xnew 投影到选取的特征向量上,这样就得到了我们需要的已经降维的数据集 xneww。

最小误差理论降维原理

而最小误差则是使得平均投影代价最小的线性投影,这一过程中,我们则需要找到的是平方错误评价函数 j0(x0) 等参数。

详细步骤可参考《从零开始实现主成分分析 (pca) 算法》:

https://blog.csdn.net/u013719780/article/details/78352262

关于 pca 算法的代码如下:

from __future__ import print_function

from sklearn import datasets

import matplotlib.pyplot as plt

import matplotlib.cm as cmx

import matplotlib.colors as colors

import numpy as np

%matplotlib inline

def shuffle_data(x, y, seed=none):

if seed:

np.random.seed(seed)

idx = np.arange(x.shape[0])

np.random.shuffle(idx)

return x[idx], y[idx]

# 正规化数据集 x

def normalize(x, axis=-1, p=2):

lp_norm = np.atleast_1d(np.linalg.norm(x, p, axis))

lp_norm[lp_norm == 0] = 1

return x / np.expand_dims(lp_norm, axis)

# 标准化数据集 x

def standardize(x):

x_std = np.zeros(x.shape)

mean = x.mean(axis=0)

std = x.std(axis=0)

# 做除法运算时请永远记住分母不能等于 0 的情形

# x_std = (x - x.mean(axis=0)) / x.std(axis=0)

for col in range(np.shape(x)[1]):

if std[col]:

x_std[:, col] = (x_std[:, col] - mean[col]) / std[col]

return x_std

# 划分数据集为训练集和测试集

def train_test_split(x, y, test_size=0.2, shuffle=true, seed=none):

if shuffle:

x, y = shuffle_data(x, y, seed)

n_train_samples = int(x.shape[0] * (1-test_size))

x_train, x_test = x[:n_train_samples], x[n_train_samples:]

y_train, y_test = y[:n_train_samples], y[n_train_samples:]

return x_train, x_test, y_train, y_test

# 计算矩阵 x 的协方差矩阵

def calculate_covariance_matrix(x, y=np.empty((0,0))):

if not y.any:

y = x

n_samples = np.shape(x)[0]

covariance_matrix = (1 / (n_samples-1)) * (x - x.mean(axis=0)).t.dot(y - y.mean(axis=0))

return np.array(covariance_matrix, dtype=float)

# 计算数据集 x 每列的方差

def calculate_variance(x):

n_samples = np.shape(x)[0]

variance = (1 / n_samples) * np.diag((x - x.mean(axis=0)).t.dot(x - x.mean(axis=0)))

return variance

# 计算数据集 x 每列的标准差

def calculate_std_dev(x):

std_dev = np.sqrt(calculate_variance(x))

return std_dev

# 计算相关系数矩阵

def calculate_correlation_matrix(x, y=np.empty([0])):

# 先计算协方差矩阵

covariance_matrix = calculate_covariance_matrix(x, y)

# 计算 x, y 的标准差

std_dev_x = np.expand_dims(calculate_std_dev(x), 1)

std_dev_y = np.expand_dims(calculate_std_dev(y), 1)

correlation_matrix = np.divide(covariance_matrix, std_dev_x.dot(std_dev_y.t))

return np.array(correlation_matrix, dtype=float)

class pca:

"""

主成份分析算法 pca,非监督学习算法.

"""

def __init__(self):

self.eigen_values = none

self.eigen_vectors = none

self.k = 2

def transform(self, x):

"""

将原始数据集 x 通过 pca 进行降维

"""

covariance = calculate_covariance_matrix(x)

# 求解特征值和特征向量

self.eigen_values, self.eigen_vectors = np.linalg.eig(covariance)

# 将特征值从大到小进行排序,注意特征向量是按列排的,即 self.eigen_vectors 第 k 列是 self.eigen_values 中第 k 个特征值对应的特征向量

idx = self.eigen_values.argsort[::-1]

eigenvalues = self.eigen_values[idx][:self.k]

eigenvectors = self.eigen_vectors[:, idx][:, :self.k]

# 将原始数据集 x 映射到低维空间

x_transformed = x.dot(eigenvectors)

return x_transformed

def main:

# load the dataset

data = datasets.load_iris

x = data.data

y = data.target

# 将数据集 x 映射到低维空间

x_trans = pca.transform(x)

x1 = x_trans[:, 0]

x2 = x_trans[:, 1]

cmap = plt.get_cmap('viridis')

colors = [cmap(i) for i in np.linspace(0, 1, len(np.unique(y)))]

class_distr =

# plot the different class distributions

for i, l in enumerate(np.unique(y)):

_x1 = x1[y == l]

_x2 = x2[y == l]

_y = y[y == l]

class_distr.append(plt.scatter(_x1, _x2, color=colors[i]))

# add a legend

plt.legend(class_distr, y, loc=1)

# axis labels

plt.xlabel('principal component 1')

plt.ylabel('principal component 2')

plt.show

if __name__ == "__main__":

main

最终,我们将得到降维结果如下。其中,如果得到当特征数 (d) 远大于样本数 (n) 时,可以使用一点小技巧实现 pca 算法的复杂度转换。

pca 降维算法展示

当然,这一算法虽然经典且较为常用,其不足之处也非常明显。它可以很好的解除线性相关,但是面对高阶相关性时,效果则较差;同时,pca 实现的前提是假设数据各主特征是分布在正交方向上,因此对于在非正交方向上存在几个方差较大的方向,pca 的效果也会大打折扣。

kpca(kernel pca)

kpca 是核技术与 pca 结合的产物,它与 pca 主要差别在于计算协方差矩阵时使用了核函数,即是经过核函数映射之后的协方差矩阵。

引入核函数可以很好的解决非线性数据映射问题。kpca 可以将非线性数据映射到高维空间,在高维空间下使用标准 pca 将其映射到另一个低维空间。

kpca 降维算法展示

详细内容可参见 《python 机器学习》之特征抽取——kpca:

https://blog.csdn.net/weixin_40604987/article/details/79632888

代码地址:

https://github.com/heucoder/dimensionality_reduction_alo_codes/blob/master/codes/pca/kpca.py

lda(linear discriminant analysis)

lda 是一种可作为特征抽取的技术,其目标是向最大化类间差异,最小化类内差异的方向投影,以利于分类等任务即将不同类的样本有效的分开。lda 可以提高数据分析过程中的计算效率,对于未能正则化的模型,可以降低维度灾难带来的过拟合。

lda 降维算法展示

详细内容可参见《数据降维—线性判别分析(lda)》:

https://blog.csdn.net/chenvast/article/details/79227945

代码地址:

https://github.com/heucoder/dimensionality_reduction_alo_codes/tree/master/codes/lda

mds(multidimensional scaling)

mds 即多维标度分析,它是一种通过直观空间图表示研究对象的感知和偏好的传统降维方法。该方法会计算任意两个样本点之间的距离,使得投影到低维空间之后能够保持这种相对距离从而实现投影。

由于 sklearn 中 mds 是采用迭代优化方式,下面实现了迭代和非迭代的两种。

mds 降维算法展示

详细内容可参见《mds 算法》

https://blog.csdn.net/zhangweiguo_717/article/details/69663452

代码地址:

https://github.com/heucoder/dimensionality_reduction_alo_codes/tree/master/codes/mds

isomap

isomap 即等度量映射算法,该算法可以很好地解决 mds 算法在非线性结构数据集上的弊端。

mds 算法是保持降维后的样本间距离不变,isomap 算法则引进了邻域图,样本只与其相邻的样本连接,计算出近邻点之间的距离,然后在此基础上进行降维保距。

isomap 降维算法展示

详细内容可参见《isomap》

https://blog.csdn.net/zhangweiguo_717/article/details/69802312

代码地址:

https://github.com/heucoder/dimensionality_reduction_alo_codes/tree/master/codes/isomap

lle(locally linear embedding)

lle 即局部线性嵌入算法,它是一种非线性降维算法。该算法核心思想为每个点可以由与它相邻的多个点的线性组合而近似重构,然后将高维数据投影到低维空间中,使其保持数据点之间的局部线性重构关系,即有相同的重构系数。在处理所谓的流形降维的时候,效果比 pca 要好很多。

lle 降维算法展示

详细内容可参见《lle 原理及推导过程》

https://blog.csdn.net/scott198510/article/details/76099630

代码地址:

https://github.com/heucoder/dimensionality_reduction_alo_codes/tree/master/codes/lle

t-sne

t-sne 也是一种非线性降维算法,非常适用于高维数据降维到 2 维或者 3 维进行可视化。它是一种以数据原有的趋势为基础,重建其在低纬度(二维或三维)下数据趋势的无监督机器学习算法。

下面的结果展示参考了源代码,同时也可用 tensorflow 实现(无需手动更新参数)。

t-sne 降维算法展示

详细内容可参见《t-sne 使用过程中的一些坑》:

http://bindog.github.io/blog/2018/07/31/t-sne-tips/

代码地址:

https://github.com/heucoder/dimensionality_reduction_alo_codes/tree/master/codes/t-sne

le(laplacian eigenmaps)

le 即拉普拉斯特征映射,它与 lle 算法有些相似,也是以局部的角度去构建数据之间的关系。它的直观思想是希望相互间有关系的点(在图中相连的点)在降维后的空间中尽可能的靠近;以这种方式,可以得到一个能反映流形的几何结构的解。

le 降维算法展示

详细内容可参见《拉普拉斯特征图降维及其 python 实现》:

https://blog.csdn.net/hustlx/article/details/50850342

代码地址:

https://github.com/heucoder/dimensionality_reduction_alo_codes/tree/master/codes/le

lpp(locality preserving projections)

lpp 即局部保留投影算法,其思路和拉普拉斯特征映射类似,核心思想为通过最好的保持一个数据集的邻居结构信息来构造投影映射,但 lpp 不同于 le 的直接得到投影结果,它需要求解投影矩阵。

lpp 降维算法展示

详情请参见《局部保留投影算法 (lpp) 详解》:

https://blog.csdn.net/qq_39187538/article/details/90402961

代码地址:

https://github.com/heucoder/dimensionality_reduction_alo_codes/tree/master/codes/lpp

*《dimensionality_reduction_alo_codes》项目作者简介

heucoder,目前是哈尔滨工业大学计算机技术在读硕士生,主要活跃于互联网领域,知乎昵称为「超爱学习」,其 github 主页地址为: https://github.com/heucoder。

github 项目地址:

https://github.com/heucoder/dimensionality_reduction_alo_codes

雷锋网 ai 开发者

雷锋网年度评选——寻找19大行业的最佳ai落地实践

创立于2017年的「ai最佳掘金案例年度榜单」,是业内首个人工智能商业案例评选活动。雷锋网从商用维度出发,寻找人工智能在各个行业的最佳落地实践。

第三届评选已正式启动,关注微信公众号“雷锋网”,回复关键词“榜单”参与报名。详情可咨询微信号:xqxq_xq

最新排行
摄影师为流浪猫拍摄肖像,捕捉了这种无畏猫科动物的独特美
2019羽联总决赛:中国最大赢家 陈雨菲逆转世界第一夺冠
日右翼踹“慰安妇”像 台民众激愤蔡当局不吭声
3700亿购房资金要买车!国际大行发声力挺车市 上调4家车企评级!复苏还是寒冬?业内分歧巨大
港汇转弱金管局料港息缓升 新股小米集资或推高拆息
上半年63家IPO公司完成上市 仅为去年同期四成
中外车企集体发力夯实2018分水岭
男子深圳地铁丢失一个袋子 牵出17年前惊人命案
名人开飞机游艇参加气候变化峰会被批:虚伪到极点
从沉寂厂房到网红航空小镇 建德祖孙三代亲历家乡巨变
热门文章
政府再次面临“停摆” 特朗普拨款暂时阻止
被困电梯怎么办?给你支几招
丰田新4Runner,搭4.0L引擎,越野性能提升,比Jeep牧马人都硬派
12岁的李嫣第一次上杂志封面,我们却被李亚鹏深深感动!
我东风21D打航母全靠核弹头?印度军工称自己更先进
每日星语|1228双子座失去平衡;巨蟹座显得吝啬
露华浓集团亏损收窄 中国市场是助力之一
终于要推新款了,首曝长城魏派VV7新款造型,细节感官提升明显
做好猪瘟疫情防控 静待来年一季度猪价回升
9.7国服中路五大天王 天使姐妹强势称霸
© Copyright 2018-2019 saogh.com 王角资讯 Inc. All Rights Reserved.