如何在PCA双向图中显示观测结果?双向、图中、结果、如何在

2023-09-03 13:54:38 作者:陪你度过漫长岁月

我希望使用FactoExtra程序包在主成分分析双曲线图上显示观测行名。

library(factoextra)
## organize dataset to obtain unique rownames
iris$id <- rownames(iris)
iris$id <- paste(iris$Species, iris$id, sep="_")
rownames(iris) <- iris$id
iris <- iris[,1:4]

## Run PCA
pca <- prcomp(iris, center=T, scale. = T)

## On this plot observation names are displayed
fviz_pca_ind(pca,
         col.ind = "cos2", # Color by the quality of representation
         gradient.cols = c("#00AFBB", "#E7B800", "#FC4E07"),
         repel = TRUE
)

)

## Is there any way to display the same row names in this plot? (instead of points)
fviz_pca_biplot(pca,
            col.var = "#2E9FDF", # Variables color
            col.ind = "#696969",  # Individuals color
            label = 'var',
            repel = TRUE
)
因子分析︱使用Stata做主成分分析

提前谢谢

推荐答案

您可以指定geom.ind=";txt";

fviz_pca_biplot(pca,geom.ind="text",labelsize=2,
col.var = "#2E9FDF",col.ind = "#696969")

并尝试添加ggrepel:

options(ggrepel.max.overlaps = 20)
fviz_pca_biplot(pca,geom.ind="text",labelsize=2,
col.var = "#2E9FDF",col.ind = "#696969",repel=TRUE)