R语言读取hallmarks的gmt文档的不同姿势整理
不同格式各有所用
1.读取数据框格式的
hallmarks <- clusterProfiler::read.gmt("~/genelist/h.all.v7.4.symbols.gmt") #返回的是表格
hallmarks$term<- gsub('HALLMARK_',"",hallmarks$term)
适配Y叔的clusterProfiler的后续分析,比如整理后genelist后做gsea
gseahallmarks <- GSEA(geneList,
TERM2GENE =hallmarks)
2.读取列表格式的
hallmarks <- fgsea::gmtPathways("~/genelist/h.all.v7.4.symbols.gmt")
names(hallmarks)<- gsub('HALLMARK_',"",names(hallmarks))
原文地址:https://blog.csdn.net/m0_58549466/article/details/144569169
本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.kler.cn/a/448139.html 如若内容造成侵权/违法违规/事实不符,请联系邮箱:809451989@qq.com进行投诉反馈,一经查实,立即删除!
本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.kler.cn/a/448139.html 如若内容造成侵权/违法违规/事实不符,请联系邮箱:809451989@qq.com进行投诉反馈,一经查实,立即删除!