热点新闻
R语言用列数据画PCA带椭圆图
2023-07-09 05:05  浏览:1499  搜索引擎搜索“手机晒展网”
温馨提示:信息一旦丢失不一定找得到,请务必收藏信息以备急用!本站所有信息均是注册会员发布如遇到侵权请联系文章中的联系方式或客服删除!
联系我时,请说明是在手机晒展网看到的信息,谢谢。
展会发布 展会网站大全 报名观展合作 软文发布

之前找了很多信息,见到的都是用行画PCA,用列的大多都是使用fviz_pca_var这个函数,但这个函数椭圆我又死活加不上,再加上画椭圆时的分组信息又没法加在原数据里,纠结了很久。但最后还是被我试出来了。

先放个原数据格式:




image




image

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" #分组名

)

结果图:




image

生信

作者:清珺
链接:https://www.jianshu.com/p/be963a3b0e81

发布人:6dc9****    IP:117.173.23.***     举报/删稿
展会推荐
让朕来说2句
评论
收藏
点赞
转发