作品分享
问答交流
发现
任务
客服工单
library("FactoMineR") library("factoextra") library("dplyr") data = read.csv("G:/TPTNW.csv",row.names = "gene_id") 这是分组信息 datagroup = read.csv("G:/TPTNWgroup.csv") t()用于转换矩阵行列,用起来方便省事 verdata = t(data) res.pca <- PCA(verdata, graph = FALSE) 将同一组的individuals圈在一起 fviz_pca_ind(res.pca, geom.ind = c("point","text"), # 表现点与笔墨 col.ind = datagroup$group, #以组分颜色,这里直接上面导入的group信息即可 palette = c("#00AFBB", "#E7B800", "#FC4E07"), #颜色 addEllipses = TRUE, #画椭圆 legend.title = "Groups" #分组名)
举报
本版积分规则 回帖后跳转到最后一页
Powered by CangBaoKu v1.0 小黑屋藏宝库It社区( 冀ICP备14008649号 )
GMT+8, 2025-8-9 20:48, Processed in 0.161855 second(s), 35 queries.© 2003-2025 cbk Team.