主成分分析介绍可以看此视频
library(FactoMineR)
data("iris")
iris_pca <- PCA(iris[,-5], graph =F)
summary(iris_pca)
library(factoextra)
eig_val <- get_eigenvalue(iris_pca) # 获得特征值,其实summary里也有
eig_val
scree_plot <-fviz_eig(iris_pca)
scree_plot
fviz_pca_var(iris_pca, col.var = "black")
# 个体PCA可视化
iris_pca
ind_plot <- fviz_pca_ind(iris_pca,
geom.ind = "point", # 点图
col.ind = iris$Species, # 颜色分类
palette = c("#00AFBB", "#E7B800", "#FC4E07"),#颜色设置
addEllipses = TRUE, # 添加椭圆
legend.title = "Groups") #添加标题
ind_plot
代码比较简单,主要还是理解原理