QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 2556|回复: 0
打印 上一主题 下一主题

python 解决主成分分析 pca

[复制链接]
字体大小: 正常 放大

1189

主题

4

听众

2934

积分

该用户从未签到

跳转到指定楼层
1#
发表于 2024-3-21 11:03 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
主成分分析(Principal Component Analysis,PCA)是一种常用的降维技术和数据预处理方法,它通过线性变换将高维数据映射到低维空间中,以找到数据中的主要特征。
. X' P8 t1 L) x主成分分析的基本思想是将原始数据投影到一个新的坐标系中,使得投影后的数据具有最大的方差。这些新的坐标轴被称为主成分,而每个主成分都是原始特征的线性组合。主成分按照其所解释的方差贡献程度进行排序,最重要的主成分排在前面。% [( h# }  U6 q6 Q2 J. Z2 \; w: L
主成分分析的步骤如下:0 x5 Y5 B: t; H8 D! S
0 w5 G) ?( F6 k- t9 s! m" e
1.标准化数据:将原始数据进行标准化处理,使得每个特征的均值为0,方差为1。这样可以避免某些特征由于量级差异造成的影响。
* d8 [( I! k. b4 p1 K# d2.计算协方差矩阵:根据标准化后的数据计算协方差矩阵。协方差矩阵描述了不同特征之间的相关性。
+ w* c8 y1 L' v" ~& l3 l+ T3.计算特征值和特征向量:对协方差矩阵进行特征值分解,得到特征值和对应的特征向量。特征值表示每个主成分所解释的方差,特征向量表示各个主成分的方向。$ V0 H& s3 }2 {2 X" d
4.选择主成分:根据特征值的大小选择要保留的主成分的数量。通常选择保留累计贡献率较高的主成分。
5 j/ ~) O; r# W" {/ F$ O4 e" n5.数据转换:将原始数据投影到选定的主成分上,得到降维后的数据。
+ i8 ?# C# |/ ~) |3 o$ k1 j$ @
1 ]& ?! W5 Z  }1 W; V/ u主成分分析的主要应用包括降维、可视化、特征提取和去除数据中的噪音。通过降低数据的维度,主成分分析可以简化数据集并去除冗余信息,从而提高后续分析的效率和准确性。
0 a% k# A. q$ @3 W4 I' L3 Q; Y4 O6 ^2 P4 r7 c
逐行解释代码的含义:
1 u* R3 E9 W; q; |3 ^& B- Eimport numpy as np
# A# A+ X8 |9 U2 ~import pandas as pd
# v+ `" B/ Y3 B: _from sklearn.decomposition import PCA7 K, R5 \4 o% d6 ]2 q
- K# y+ F+ Q, h: v
这些是导入所需的库。numpy用于数值计算,pandas用于数据处理,sklearn.decomposition中的PCA用于主成分分析。, _8 f! r) |; X9 |
df = pd.DataFrame({9 ^( K7 U0 v$ p0 x4 [* G+ q, s0 x
    'x1': [149.5, 162.5, 162.7, 162.2, 156.5],8 A4 \9 W1 D" w; Q+ _) M: C
    'x2': [69.5, 77, 78.5, 87.5, 74.5],4 J! Q& {7 m- ~  W4 \* b  `5 R
    'x3': [38.5, 55.5, 50.8, 65.5, 49]
* p) C; I( F, W$ ?  h  T}). ^0 |' u$ {$ w0 o3 p
2 ~. j+ n1 O- H& i. R8 G- e. E* p
这里创建了一个数据帧df,包含了3个变量 x1、x2、x3 的观测值。数据集中每一列代表一个变量,每一行代表一个观测值。/ o- H! z# v- r4 r3 I. t$ R
model = PCA().fit(np.array(df))2 j  R; ~! q( G  K3 T) R
) c4 d4 J/ x+ v; |& U4 }2 {+ F2 z
这行代码创建了一个PCA对象,并使用fit方法拟合数据。fit方法将数据df作为输入,并根据数据计算主成分分析模型。
6 ]; T0 E+ e7 S% A: Pprint('特征值:', model.explained_variance_): `5 }- w. p6 e3 ^/ {
print('贡献率:', model.explained_variance_ratio_)
( G) R7 h1 {; V4 c5 i) F  Lprint('各主成分的系数:', model.components_)
( _+ Z, U9 ^! }4 {' n/ S) H0 }
这几行代码分别打印了主成分分析模型的三个重要属性:
) K. [0 w+ u. Q* k: C' p" `# H# \8 R: ^# C3 c. `% K% y
1.explained_variance_:特征值,表示每个主成分的方差。1 B' u- \) J" U4 b3 u2 c* C. o# N
2.explained_variance_ratio_:贡献率,表示每个主成分的方差占总方差的比例。
, t; l: {0 W3 p$ g/ x% o3.components_:各主成分的系数,表示每个主成分在原始变量空间中的权重。. ?4 W7 c% i* o
' L2 `) D7 _0 e9 f* i. E0 e
pca_df = pd.DataFrame(model.transform(np.array(df)))4 S; M- I0 ^: N% N# F
pca_df.columns = ['F1', 'F2', 'F3']
( L4 J5 l" W  j2 i3 ~, ]7 _$ N0 Npca_df+ x5 ?' Z' @% \' y$ o3 e' P
7 q. {8 B, N8 \& A& O* W  r+ K2 v
这几行代码使用model.transform方法将原始数据进行主成分转换,并将结果存储到一个新的数据帧pca_df中。pca_df包含三个列,分别命名为'F1'、'F2'、'F3',分别表示三个主成分的值。
& I* x! w% [, i$ \( G希望这个逐行解释对你有帮助!如果你还有其他问题,请随时提问。$ p8 Y7 s, h+ m% o: n  j$ P

$ f! W1 Y6 n; [  x; U2 a% y3 K7 c$ ]+ i" b

29.pca.py

671 Bytes, 下载次数: 0, 下载积分: 体力 -2 点

售价: 2 点体力  [记录]  [购买]

zan
转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信
您需要登录后才可以回帖 登录 | 注册地址

qq
收缩
  • 电话咨询

  • 04714969085
fastpost

关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

手机版|Archiver| |繁體中文 手机客户端  

蒙公网安备 15010502000194号

Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

GMT+8, 2026-8-5 01:39 , Processed in 0.562842 second(s), 55 queries .

回顶部