Я не думаю, что вы можете изменить вывод с fviz_pca_ind()
, поэтому вам нужно будет извлечь данные из результатов и построить их снова, используя ggplot2:
library(factoextra)
library(ggplot2)
data <- iris
colnames(data)[5] <- "A"
data$B <- sample(letters[1:2],nrow(data),replace=TRUE)
res.pca <- prcomp(data[,1:4], scale = TRUE)
basic_plot <- fviz_pca_ind(res.pca, label="none")
ggplot(cbind(basic_plot$data,data[,c("A","B")]),
aes(x=x,y=y,col=A,shape=B)) + geom_point() + theme_bw()