enrichCluster 关于非模式物种富集分析的使用

引言


其实 enrichCluster 既可以针对于模式物种进行富集分析,当然也是可以针对非模式物种的,不过使用者们针对非模式物种时,往往会报错。这里做一个示例说明,演示如何对非模式物种进行分析。
提供自己的 orgdb 文件


假如你已经有了自己构建好的 orgdb 文件,只需要设置一下 id.trans/toType/readable 几个参数即可:
enrich <- enrichCluster(object = cm,
OrgDb = org.myorg.eg.db,
type = “BP”,
id.trans = F,
toType = “GID”,
readable = F,
pvalueCutoff = 0.5,
topn = 5)
提供 gene-term-name 文件


你也可以用 基因-通路 id 和 基因-通路名称 这两个文件来进行富集分析,可以参考 enrichr 函数 需要准备的文件格式:
enrich <- enrichCluster(object = cm,
type = “ownSet”,
id.trans = F,
TERM2GENE = ko2gene,
TERM2NAME = ko2name,
pvalueCutoff = 0.5,
topn = 5)

上面对于非模式物种就可以顺利运行了。
结尾


路漫漫其修远兮,吾将上下而求索。
欢迎加入生信交流群。加我微信我也拉你进 微信群聊 老俊俊生信交流群 (微信交流群需收取 20 元入群费用,一旦交费,拒不退还!(防止骗子和便于管理)) 。QQ 群可免费加入, 记得进群按格式修改备注哦。

声明:文中观点不代表本站立场。本文传送门:https://eyangzhen.com/419356.html

(0)
联系我们
联系我们
分享本页
返回顶部