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