分享是一种态度
1引言
有小伙伴在对 clusterData 聚类后的亚群,使用 enrichCluster 来做富集分析时,想用自己的基因集来分析,但是不支持。另外一种办法就是你自己拿亚群去使用 clusterProfiler::enricher 函数自己来做富集分析,然后整理一下数据,做成和 enrichCluster 输出结果一样的格式就行。不过还是需要周转一番。我直接加点代码到 enrichCluster 函数里也不是特别麻烦。
2安装
重新安装获取新功能:
3使用介绍
聚类:
library(ClusterGVis)
library(org.Mm.eg.db)
然后自己去 GSEA 官网下载一下基因集数据:
然后准备一下 enricher 函数需要的 TERM2GENE 和 TERM2NAME 数据即可。设置一下富集类型 type = “ownSet” 。
4结尾
如果有任何疑问或者建议,欢迎在 github 留言!
欢迎加入生信交流群。加我微信我也拉你进 微信群聊 老俊俊生信交流群 (微信交流群需收取 20 元入群费用,一旦交费,拒不退还!(防止骗子和便于管理)) 。QQ 群可免费加入, 记得进群按格式修改备注哦。
老俊俊微信:
知识星球:
往期回顾目录
R 语言做竞争风险回归分析 听说你还想添加 KEGG 注释? BioSeqUtils 延伸序列及提取启动子序列 ClusterGVis 对接单细胞啦 About BioSeqUtils 基因特征序列提取 R 包 BioSeqUtils 关于 chatGPT 无法连接的问题 问问 chatGPT 来提取最长转录本及基因序列 GseaVis 一键对接 GSEA 软件结果并可视化 同学你又在画 GSEA?
声明:文中观点不代表本站立场。本文传送门:https://eyangzhen.com/190726.html