Как я могу построить кластеры в fviz_pca_biplot в R? - PullRequest
Новая
       30

Как я могу построить кластеры в fviz_pca_biplot в R?

0 голосов
/ 06 июня 2019

Я хотел бы отобразить результаты кластера в моем fviz_pca_biplot

Я создал кластерный анализ с функцией eclust и PCA (prcomp) с библиотекой Factorextra, я строю свой pca с помощью fviz_pca_biplot, но яхотел бы добавить участки кластера к моему графику fviz_pca_biplot, это возможно? 

library(devtools)
library(ggbiplot)
library(readxl)

vv.poll6<-read.csv("vv_poll.csv")
vv.veg6<-read.csv("vv_veg_free100.csv")
vv.land = read.csv("vv_landscape.csv")

vv.poll6.p<-vv.poll6[,-1]/rowSums(vv.poll6[,-1])
vv.poll6.p2 = (sqrt(vv.poll6.p))
vv.veg6.p<-vv.veg6[ ,-c(1,2)]/rowSums(vv.veg6[ ,-c(1,2)])
vv.veg6.p2 = (sqrt(vv.veg6.p))
vv.poll6.pca = prcomp(vv.poll6.p2, center = TRUE, scale. = FALSE)

library("factoextra")

fviz_pca_biplot(vv.poll6.pca, col.var = "black", col.ind = vv.land$Landscape, palette = "Dark2", label = "all", repel = TRUE, legend.title = "Landscape", title = "PCA Vila Velha - Pollen")

# Cluster - I want to add this results to the previous plot 

veg.norm6<-decostand(sqrt(vv.veg6.p), "normalize")
clust.vv = eclust(veg.norm6, FUNcluster = c("hclust"), hc_metric = "euclidean", hc_method = "ward.D2")
res.hc6.vv <- eclust(veg.norm6, "hclust", hc_metric = "euclidean", hc_method = "ward.D2", k = 2) # compute hclust
fviz_dend(res.hc6.vv, rect = TRUE) # dendrogam
fviz_cluster(res.hc6.vv)
...