QQ登录

只需要一步,快速开始

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

python 解决主成分分析 pca

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

1189

主题

4

听众

2934

积分

该用户从未签到

跳转到指定楼层
1#
发表于 2024-3-21 11:03 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
主成分分析(Principal Component Analysis,PCA)是一种常用的降维技术和数据预处理方法,它通过线性变换将高维数据映射到低维空间中,以找到数据中的主要特征。3 G$ g" ~6 Q2 u4 i; O& i
主成分分析的基本思想是将原始数据投影到一个新的坐标系中,使得投影后的数据具有最大的方差。这些新的坐标轴被称为主成分,而每个主成分都是原始特征的线性组合。主成分按照其所解释的方差贡献程度进行排序,最重要的主成分排在前面。
1 ~2 z, ^4 J% g- U& a! z. Q6 B主成分分析的步骤如下:& Y6 z3 R! r; Q$ [7 w

& V4 S! D5 m$ _* l2 V1.标准化数据:将原始数据进行标准化处理,使得每个特征的均值为0,方差为1。这样可以避免某些特征由于量级差异造成的影响。2 P+ e  A- i! F- S1 J
2.计算协方差矩阵:根据标准化后的数据计算协方差矩阵。协方差矩阵描述了不同特征之间的相关性。
% c! K3 U! ]- q+ r- M% a+ i3.计算特征值和特征向量:对协方差矩阵进行特征值分解,得到特征值和对应的特征向量。特征值表示每个主成分所解释的方差,特征向量表示各个主成分的方向。
' H* c2 S' M( o  h% c. V4.选择主成分:根据特征值的大小选择要保留的主成分的数量。通常选择保留累计贡献率较高的主成分。% z  h  E' m1 T
5.数据转换:将原始数据投影到选定的主成分上,得到降维后的数据。) N6 Q. S) W" X! N

+ |; ]3 J5 O# l, O, K主成分分析的主要应用包括降维、可视化、特征提取和去除数据中的噪音。通过降低数据的维度,主成分分析可以简化数据集并去除冗余信息,从而提高后续分析的效率和准确性。
- B4 o# X2 B+ R; }. P0 p6 g7 w9 J! ?! J& @% ^  q! D  U
逐行解释代码的含义:3 e: z/ N3 I; Z) f  D. m
import numpy as np# S( q3 o9 N# u" B* j- f& W/ F. l
import pandas as pd1 F3 X" J9 h7 \+ u. W
from sklearn.decomposition import PCA0 T7 F; z7 U% n1 G5 P

2 B6 P2 ?7 X" A, v0 E这些是导入所需的库。numpy用于数值计算,pandas用于数据处理,sklearn.decomposition中的PCA用于主成分分析。
6 Q& K' D6 Z" o0 m6 r6 T7 V0 b% N8 W' p' Kdf = pd.DataFrame({
! |: S5 ]  f8 L+ {9 r6 b, s    'x1': [149.5, 162.5, 162.7, 162.2, 156.5],
  p% ]; R- l6 ]" t9 G* e    'x2': [69.5, 77, 78.5, 87.5, 74.5],
! ]/ ~3 Z1 j: c6 g) M/ r+ S% g    'x3': [38.5, 55.5, 50.8, 65.5, 49]! U" ]0 G  W% b& x4 ]" N! d/ c; R
})
6 r  j- n+ \8 x' w2 g% B' S) }$ u+ c( a' v# J
这里创建了一个数据帧df,包含了3个变量 x1、x2、x3 的观测值。数据集中每一列代表一个变量,每一行代表一个观测值。  c. d( |) s4 J$ F0 a) |6 O
model = PCA().fit(np.array(df))! T/ |8 r9 k) R* Y3 V% Q8 C
- Z# J' ?+ y' i3 n( P5 N
这行代码创建了一个PCA对象,并使用fit方法拟合数据。fit方法将数据df作为输入,并根据数据计算主成分分析模型。1 K. {/ q+ [* _& J# i' d6 |/ Q  c+ |
print('特征值:', model.explained_variance_)
8 Q. G5 F4 J$ h0 [7 J& Mprint('贡献率:', model.explained_variance_ratio_)
/ U: B! k* S3 n4 p7 Y3 I( i3 q6 hprint('各主成分的系数:', model.components_)
+ p  K0 X9 V% v: N$ g: [2 ^) H" x
这几行代码分别打印了主成分分析模型的三个重要属性:/ A5 r" M0 B, A( N" |: r3 s+ e$ f

0 j0 A9 o( X. d$ T% z2 F: P1.explained_variance_:特征值,表示每个主成分的方差。4 t& R; d; H$ ]$ s8 O7 j
2.explained_variance_ratio_:贡献率,表示每个主成分的方差占总方差的比例。. i: t4 }. N" o  I
3.components_:各主成分的系数,表示每个主成分在原始变量空间中的权重。
6 t' Q* r/ p! b/ s
$ t" V7 p/ @6 Q4 r4 c9 l! H7 Wpca_df = pd.DataFrame(model.transform(np.array(df)))
% z! ^: ~$ D) G; W2 Lpca_df.columns = ['F1', 'F2', 'F3']
) t; ?6 `, V% ]8 u& X: n; Fpca_df
% j1 n5 l: O, |5 B1 t" {
6 e6 C/ z/ m' H# r0 T) W这几行代码使用model.transform方法将原始数据进行主成分转换,并将结果存储到一个新的数据帧pca_df中。pca_df包含三个列,分别命名为'F1'、'F2'、'F3',分别表示三个主成分的值。
' g; I0 f+ D" ~; C& I$ I; X希望这个逐行解释对你有帮助!如果你还有其他问题,请随时提问。1 V9 K% K+ c3 O' M( a2 D

9 R1 G$ R( K6 }" V
6 M4 R0 Z2 N) P" l  D4 ^3 h

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-7-10 06:58 , Processed in 0.445554 second(s), 54 queries .

回顶部