数据降维之线性判别分析(LDA)——基本原理与基于python sklearn库的LDA实现

目录简介算法流程基于python sklearn库的LDA例程简介线性判别分析(Linear Discriminate Analysis, LDA)通过正交变换将一组可能存在相关性的变量降维变量,目标是将高维数据投影至低维后,同类的数据之间距离尽可能近、不同类数据之间距离尽可能远。应用场景:对于拟合、分类算法,可以利用主成分分析对输入数据实现降维,去除冗余数据可以提高计算效率并提高计算精度...
摘要由CSDN通过智能技术生成

简介

线性判别分析(Linear Discriminate Analysis, LDA)通过正交变换将一组可能存在相关性的变量降维变量,目标是将高维数据投影至低维后,同类的数据之间距离尽可能近、不同类数据之间距离尽可能远。

应用场景:
对于拟合、分类算法,可以利用主成分分析对输入数据实现降维,去除冗余数据可以提高计算效率并提高计算精度;
对于难以进行可视化的高维数据可以利用主成分分析映射至二维进行可视化,方便进行展示,如果需要展示的数据从属于多个类别,LDA的降维效果优于PCA。
由于LDA本身降维时就以分离不同类数据为目的,因此对于二分类问题可以直接降至一维进行分类。

缺点:
利用LDA降维时有最大维数限制,若输入数据的总类别数为 C C C则最大维数是 C − 1 C-1 C1,这意味着在对高维但类别数较少的问题进行降维时LDA可能会导致大量信息的丢失。

算法流程

符号约定:
记总样本个数为 n n n,每个初始样本所含特征数为 N N N,最终选择的主成分所含特征数为 K K K
记第 i i i个样本的第 j j j个特征为: x i j x_{ij} xij
n n n个初始样本为: X 1 , X 2 , . . . , X n X_1,X_2,...,X_n X1,X2,...,Xn
每个初始样本为一个 N N N维向量,第 i i i个样本为 X i = [ x i 1 x i 1 . . . x i N ] X_i=\begin{bmatrix}x_{i1} &x_{i1} & ... &x_{iN} \end{bmatrix} Xi=[xi1xi1...xiN]
记降维后的样本为: Y 1 , Y 2 , . . . , Y

  • 1
    点赞
  • 26
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
这里是一个简单的Python实现线性判别分析LDA的代码示例: ```python import numpy as np class LDA: def __init__(self): self.mean_1 = None self.mean_2 = None self.covariance_1 = None self.covariance_2 = None def fit(self, X, y): X_1 = X[y == 0] X_2 = X[y == 1] self.mean_1 = np.mean(X_1, axis=0) self.mean_2 = np.mean(X_2, axis=0) self.covariance_1 = np.cov(X_1.T) self.covariance_2 = np.cov(X_2.T) # calculate the within-class scatter matrix Sw = self.covariance_1 + self.covariance_2 # calculate the between-class scatter matrix Sb = np.outer((self.mean_1 - self.mean_2), (self.mean_1 - self.mean_2)) # calculate the eigenvalues and eigenvectors of (Sw^-1)*Sb eig_vals, eig_vecs = np.linalg.eig(np.linalg.inv(Sw).dot(Sb)) # sort the eigenvalues in descending order idx = eig_vals.argsort()[::-1] eig_vals = eig_vals[idx] eig_vecs = eig_vecs[:, idx] # select the first eigenvector as the projection axis self.w = eig_vecs[:, 0] def predict(self, X): return np.dot(X, self.w) ``` 这个类包含了`fit`和`predict`方法。在`fit`方法中,我们首先将数据集划分成两个类别,并计算每个类别的均值和协方差矩阵。然后,我们计算类内散布矩阵和类间散布矩阵,并使用它们计算(LDA)的投影向量。在`predict`方法中,我们通过将数据集投影到投影向量上来预测类别。 请注意,这是一个相对简单的实现,可能不适用于所有情况。在实际应用中,可能需要更复杂的实现来处理更大的数据集或更高的维数。

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值