主成分分析在R语言里面的实现（PCA学习笔记）

简介

主成分分析（Principal Component Analysis，PCA），顾名思义，即是拿来分析’‘主成分“的。通常和 PCA 联系在一起的是降维，当手里的数据集有成千上万个特征时，PCA可以减少数据集的维数，同时保留对数据集贡献最大的特征。

本质上主成分分析是找到一个欧式空间的线性变换，把原始数据从“一组旧的标准正交基下的表示”转化成“另一组新的标准正交基下的表示”，降维发生在新的标准正交基下的表示，直接去掉了后面几个维度的坐标值。简单来说就是利用线性变换，将分析数据的方差投影到二维的坐标上。

在学生信学习过程中，PCA是我们经常用到的分析方法，目的是为了找到有共同特征的不同聚类，在处理RNA-seq数据中发挥作用，可用于判断批次效应或者离群点。

PCA用到的R包

在pca常用的R包就俩个，一个是FactoMineR包，此包常用于分析；另外一个是factoextra包，是用来做可视化的，factoextra包内含了基于ggplot2的数据可视化的函数，是一个非常实用的包。

以iris数据集为例，提取并可视化特征值

1589253676500.png
1589253950787.png

代码示例

library("FactoMineR")
library(factoextra)
iris.pca <- PCA(iris[,-5],  graph = T)
fviz_screeplot(iris.pca, addlabels = TRUE, 
               ylim = c(0, 75)
               )

提取可视化变量的结果（coord，cor，cos2，contribution）

var <- get_pca_var(iris.pca)
View(var)
head(var$coord)
#                  Dim.1      Dim.2       Dim.3       Dim.4
#Sepal.Length  0.8901688 0.36082989 -0.27565767 -0.03760602
#Sepal.Width  -0.4601427 0.88271627  0.09361987  0.01777631
#Petal.Length  0.9915552 0.02341519  0.05444699  0.11534978
#Petal.Width   0.9649790 0.06399985  0.24298265 -0.07535950
head(var$cos2)
#                 Dim.1       Dim.2       Dim.3        Dim.4
#Sepal.Length 0.7924004 0.130198208 0.075987149 0.0014142127
#Sepal.Width  0.2117313 0.779188012 0.008764681 0.0003159971
#Petal.Length 0.9831817 0.000548271 0.002964475 0.0133055723
#Petal.Width  0.9311844 0.004095980 0.059040571 0.0056790544
head(var$contrib)
#                 Dim.1       Dim.2     Dim.3     Dim.4
#Sepal.Length 27.150969 14.24440565 51.777574  6.827052
#Sepal.Width   7.254804 85.24748749  5.972245  1.525463
#Petal.Length 33.687936  0.05998389  2.019990 64.232089
#Petal.Width  31.906291  0.44812296 40.230191 27.415396

利用fviz_pca_ind函数进行可视化

fviz_pca_ind(iris.pca,
             geom.ind = 'point',
             habillage = iris$Species, # color by groups
             palette = c("#00AFBB", "#E7B800", "#FC4E07"),
             addEllipses = T # Concentration ellipses)```

1589254966726.png