这里有一个可能的解决方案。

library(FactoMineR)
library(factoextra)
library(dplyr)

res.ACP <- PCA(iris[,c(1:3)], scale.unit = TRUE, ncp= 5)

df <- data.frame(res.ACP$ind$coord[,1:2], Species=iris$Species)
df <- df %>%
      group_by(Species) %>%
      summarize(Dim1=mean(Dim.1), Dim2=mean(Dim.2))

fviz_pca_ind(res.ACP, geom.ind=c("point"), point.size=3, pointshape=16,
             legend.title="Species", addlabel=TRUE) +
  geom_text(data=df, aes(x=Dim1, y=Dim2, label=Species, color=Species), 
            size=6, nudge_y=0.25) +
  geom_point(data=df, aes(x=Dim1, y=Dim2, color=Species), size=5)

字符串

的数据