主成分分析在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
