【原】R语言GEO数据处理（五）

# 4. PCA分析 ----------------------------------------------------------------

library(FactoMineR)

library(factoextra)

dat=as.data.frame(t(exp2))

dat.pca <- PCA(dat, graph = FALSE)

pca_plot <- fviz_pca_ind(dat.pca,

                         geom.ind = "point",

                         col.ind = group_list,

                         palette = c("#00AFBB", "#E7B800"),

                         addEllipses = TRUE,

                         legend.title = "Groups"

)

pca_plot2 <- pca_plot+

  ggtitle("GSE13850_PCA") +

  theme(plot.title = element_text(hjust = 0.5)) #设置标题居中

pca_plot2