PCA步骤

本文是PCA(主成分分析)的详细教程,通过图文并茂的方式深入浅出地介绍了PCA的基本原理和应用。从数据降维的角度出发,解释了PCA如何寻找数据的主要特征,并展示了PCA在实际数据分析中的实施步骤和效果。适合初学者和进阶者学习备忘。
### 主成分分析(PCA)的详细步骤及实现过程 主成分分析(PCA)是一种常见的数据降维技术,其核心目标是通过线性变换将高维数据投影到低维空间,同时最大程度地保留原始数据中的方差信息。以下是PCA的具体使用步骤及其实现方法: #### 1. 数据预处理 在应用PCA之前,通常需要对数据进行标准化处理。这是因为不同特征可能具有不同的量纲或尺度,如果不加以调整,可能导致某些特征占据主导地位而掩盖其他重要特征的信息[^2]。 ```python from sklearn.preprocessing import StandardScaler scaler = StandardScaler() X_scaled = scaler.fit_transform(X) ``` #### 2. 计算协方差矩阵 为了捕捉变量间的相关关系,需构建协方差矩阵 \( \Sigma \),其中每个元素表示两个随机向量之间的协方差。对于经过标准化后的数据集 \( X_{scaled} \),协方差矩阵可以通过以下方式计算得到[^3]。 \[ \Sigma = \frac{1}{N-1} X_{scaled}^\top X_{scaled} \] #### 3. 特征值分解 通过对协方差矩阵进行特征值分解,可以获得一组特征值以及对应的特征向量。这些特征向量定义了新坐标系的方向,而特征值则反映了各个方向上的方差大小。按照特征值从大到小排列,选取前 k 个最大的特征值所对应的方向作为主成分[^1]。 ```python import numpy as np cov_matrix = np.cov(X_scaled.T) eigenvalues, eigenvectors = np.linalg.eig(cov_matrix) ``` #### 4. 确定主成分数目 依据累计贡献率原则来决定最终保留下来的主成分数目 k 。一般情况下,当累积解释比例达到某个阈值(如85%-95%之间)即可停止增加维度数[^2]。 ```python total_variance = sum(eigenvalues) explained_variances = [(i / total_variance) * 100 for i in sorted(eigenvalues, reverse=True)] cumulative_explained_variance = np.cumsum(explained_variances) k = next(i for i,v in enumerate(cumulative_explained_variance) if v >= 90) + 1 print(f"Number of components to retain: {k}") ``` #### 5. 投影至低维子空间 利用选定的 k 个最大特征值对应的特征向量构成投影矩阵 W ,再把原数据乘以此矩阵完成降维操作[^3]。 ```python W = eigenvectors[:, :k] X_pca = X_scaled.dot(W) ``` 或者直接采用 `sklearn` 提供的功能简化流程: ```python from sklearn.decomposition import PCA pca = PCA(n_components=k) X_pca_sklearn = pca.fit_transform(X_scaled) ``` 以上便是完整的PCA算法理论基础与实践指南。
评论
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符
 
红包 添加红包
表情包 插入表情
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值