factoextra是一个R软件包,可以轻松提取和可视化探索性多变量数据分析的输出,其中包括:
主成分分析(PCA),用于通过在不丢失重要信息的情况下减少数据的维度来总结连续(即定量)多变量数据中包含的信息。
对应分析(CA),它是适用于分析由两个定性变量(或分类数据)形成的大型列联表的主成分分析的扩展。
多重对应分析(MCA),它是将CA改编为包含两个以上分类变量的数据表格。
多因素分析(MFA)专用于数据集,其中变量按组(定性和/或定量变量)组织。
分层多因素分析(HMFA):在数据组织为分层结构的情况下,MFA的扩展。
混合数据因子分析(FAMD)是MFA的一个特例,致力于分析包含定量和定性变量的数据集。
有许多R包实现主要组件方法。这些软件包包括:FactoMineR,ade4,stats,ca,MASS和ExPosition。
- 为什么使用
- R package factoextra具有灵活且易于使用的方法,可以用人类可读的标准数据格式快速提取上述不同软件包的分析结果。
- 它以较少的输入产生了基于ggplot2的优雅数据可视化。
- 它还包含许多便于聚类分析和可视化的功能。
- factoextra R软件包可以处理来自多个软件包的PCA,CA,MCA,MFA,FAMD和HMFA的结果,用于提取和可视化数据中包含的最重要信息。
- 在PCA,CA,MCA,MFA,FAMD和HMFA之后,可以使用以下高亮显示最重要的行/列元素:
它们的cos2值对应于它们在因子图上的表示质量
他们对主要维度定义的贡献。
如果你想这样做,factoextra包提供了一个方便的解决方案。
- 如果您想使用PCA / MCA进行预测并使用ggplot2可视化补充变量/个体在因子图上的位置:那么factoextra可以为您提供帮助。
- 提取并可视化特征值/方差:
library("factoextra")
data("decathlon2")
df <- decathlon2[1:23, 1:10]
library("FactoMineR")
res.pca <- PCA(df, graph = T)
get_eig(res.pca)
fviz_screeplot(res.pca, addlabels = TRUE, ylim = c(0, 50))
# Extract the results for variables
var <- get_pca_var(res.pca)
var
#> Principal Component Analysis Results for variables
#> ===================================================
#> Name Description
#> 1 "$coord" "Coordinates for the variables"
#> 2 "$cor" "Correlations between variables and dimensions"
#> 3 "$cos2" "Cos2 for the variables"
#> 4 "$contrib" "contributions of the variables"
# Coordinates of variables
head(var$coord)
#> Dim.1 Dim.2 Dim.3 Dim.4 Dim.5
#> X100m -0.8506257 -0.17939806 0.3015564 0.03357320 -0.1944440
#> Long.jump 0.7941806 0.28085695 -0.1905465 -0.11538956 0.2331567
#> Shot.put 0.7339127 0.08540412 0.5175978 0.12846837 -0.2488129
#> High.jump 0.6100840 -0.46521415 0.3300852 0.14455012 0.4027002
#> X400m -0.7016034 0.29017826 0.2835329 0.43082552 0.1039085
#> X110m.hurdle -0.7641252 -0.02474081 0.4488873 -0.01689589 0.2242200
# Contribution of variables
head(var$contrib)
#> Dim.1 Dim.2 Dim.3 Dim.4 Dim.5
#> X100m 17.544293 1.7505098 7.338659 0.13755240 5.389252
#> Long.jump 15.293168 4.2904162 2.930094 1.62485936 7.748815
#> Shot.put 13.060137 0.3967224 21.620432 2.01407269 8.824401
#> High.jump 9.024811 11.7715838 8.792888 2.54987951 23.115504
#> X400m 11.935544 4.5799296 6.487636 22.65090599 1.539012
#> X110m.hurdle 14.157544 0.0332933 16.261261 0.03483735 7.166193
# Graph of variables: default plot
fviz_pca_var(res.pca, col.var = "black")
-
提取和可视化变量的结果:
可以使用它们对主轴的贡献(“contrib”)来控制变量颜色:
# Contributions of variables to PC1
fviz_contrib(res.pca, choice = "var", axes = 1, top = 10)
# Contributions of variables to PC2
fviz_contrib(res.pca, choice = "var", axes = 2, top = 10)
- 为每一个提取和可视化结果:
# Extract the results for individuals
ind <- get_pca_ind(res.pca)
ind
#> Principal Component Analysis Results for individuals
#> ===================================================
#> Name Description
#> 1 "$coord" "Coordinates for the individuals"
#> 2 "$cos2" "Cos2 for the individuals"
#> 3 "$contrib" "contributions of the individuals"
# Coordinates of individuals
head(ind$coord)
#> Dim.1 Dim.2 Dim.3 Dim.4 Dim.5
#> SEBRLE 0.1955047 1.5890567 0.6424912 0.08389652 1.16829387
#> CLAY 0.8078795 2.4748137 -1.3873827 1.29838232 -0.82498206
#> BERNARD -1.3591340 1.6480950 0.2005584 -1.96409420 0.08419345
#> YURKOV -0.8889532 -0.4426067 2.5295843 0.71290837 0.40782264
#> ZSIVOCZKY -0.1081216 -2.0688377 -1.3342591 -0.10152796 -0.20145217
#> McMULLEN 0.1212195 -1.0139102 -0.8625170 1.34164291 1.62151286
# Graph of individuals
# 1\. Use repel = TRUE to avoid overplotting
# 2\. Control automatically the color of individuals using the cos2
# cos2 = the quality of the individuals on the factor map
# Use points only
# 3\. Use gradient color
fviz_pca_ind(res.pca, col.ind = "cos2",
gradient.cols = c("#00AFBB", "#E7B800", "#FC4E07"),
repel = TRUE # Avoid text overlapping (slow if many points)
)
- Color individuals by groups:
# Compute PCA on the iris data set
# The variable Species (index = 5) is removed
# before PCA analysis
iris.pca <- PCA(iris[,-5], graph = FALSE)
# Visualize
# Use habillage to specify groups for coloring
fviz_pca_ind(iris.pca,
label = "none", # hide individual labels
habillage = iris$Species, # color by groups
palette = c("#00AFBB", "#E7B800", "#FC4E07"),
addEllipses = TRUE # Concentration ellipses
)
相关资料:无监督学习 聚类分析① - 简书
https://www.jianshu.com/p/921aa71afdf8
无监督学习 聚类分析② - 简书
https://www.jianshu.com/p/514d771c243e
无监督学习 聚类分析③ - 简书
https://www.jianshu.com/p/b8dcfc7c1d4c
无监督学习 聚类分析④ - 简书
https://www.jianshu.com/p/e9ba1f5f5ae7
factoextra包 聚类分析(2) - 简书
https://www.jianshu.com/p/18f14b6774f9