R语言的enrichKEGG_pip.r如何使用

发布时间:2022-03-21 10:53:20 作者:iii
来源:亿速云 阅读:266

本文小编为大家详细介绍“R语言的enrichKEGG_pip.r如何使用”,内容详细,步骤清晰,细节处理妥当,希望这篇“R语言的enrichKEGG_pip.r如何使用”文章能帮助大家解决疑惑,下面跟着小编的思路慢慢深入,一起来学习新知识吧。

enrichKEGG_pip.r  KEGG 富集分析

使用说明:

$Rscript $scriptdir/enrichKEGG_pip.r -h
usage: /work/my_stad_immu/scripts/enrichKEGG_pip.r [-h] -g gene.list
                                                   [-d ann.db] [-t idtype]
                                                   [-s show.type]
                                                   [-O organism]
                                                   [-p pvalueCutoff]
                                                   [-q qvalueCutoff]
                                                   [-c showCategory]
                                                   [-n prefix] [-o outdir]
                                                   [-H height] [-W width]
KEGG enrich analysis : https://www.亿速云.com/article/1503
optional arguments:
  -h, --help            show this help message and exit
  -g gene.list, --gene.list gene.list
                        diff expressed gene list file, required
  -d ann.db, --ann.db ann.db
                        org.Hs.eg.db or org.Mm.eg.db ,for more info visit:
                        https://www.亿速云.com/article/1244 [optional,
                        default: org.Hs.eg.db ]
  -t idtype, --idtype idtype
                        deg gene id type: ENSEMBL SYMBOL ENTREZID GENENAME
                        [optional, default: SYMBOL ]
  -s show.type, --show.type show.type
                        set example ID type name [optional, default: ENTREZID
                        ]
  -O organism, --organism organism
                        organism ,for more info visit:
                        https://www.亿速云.com/article/787 [optional,
                        default: hsa ]
  -p pvalueCutoff, --pvalueCutoff pvalueCutoff
                        pvalue cutoff on enrichment tests to report,
                        [optional, default: 0.05 ]
  -q qvalueCutoff, --qvalueCutoff qvalueCutoff
                        qvalue cutoff on enrichment tests to report as
                        significant. Tests must pass i) pvalueCutoff on
                        unadjusted pvalues, ii) pvalueCutoff on adjusted
                        pvalues and iii) qvalueCutoff on qvalues to be
                        reported., [optional, default: 0.2 ]
  -c showCategory, --showCategory showCategory
                        how many KEGG Category to show, [optional, default: 10
                        ]
  -n prefix, --prefix prefix
                        the output file prefix [optional, default: KEGG ]
  -o outdir, --outdir outdir
                        output file directory [default cwd]
  -H height, --height height
                        the height of pic inches [default 5]
  -W width, --width width
                        the width of pic inches [default 7]

参数说明:

-g 输入基因列表文件:  脚本会读取第一列基因ID作为基因集:

-d 物种注释数据库: 人的:org.Hs.eg.db 

-t 指定输入基因列表的ID类型

使用举例:

Rscript $scriptdir/enrichKEGG_pip.r --gene.list $workdir/04.deg/S1_vs_S2.DEG.final.tsv \
  -o KEGG -n S1_vs_S2 --pvalueCutoff 0.01 --ann.db org.Hs.eg.db --organism hsa \
  --idtype SYMBOL

读到这里,这篇“R语言的enrichKEGG_pip.r如何使用”文章已经介绍完毕,想要掌握这篇文章的知识点还需要大家自己动手实践使用过才能领会,如果想了解更多相关内容的文章,欢迎关注亿速云行业资讯频道。

推荐阅读:
  1. R语言的Dataframe怎么使用
  2. 使用R语言怎么删除指定的数据

免责声明:本站发布的内容(图片、视频和文字)以原创、转载和分享为主,文章观点不代表本网站立场,如果涉及侵权请联系站长邮箱:is@yisu.com进行举报,并提供相关证据,一经查实,将立刻删除涉嫌侵权内容。

r语言

上一篇:R语言的enrichGSEA_pip.r怎么用

下一篇:R语言Cox与KM生存结果分析

相关阅读

您好,登录后才能下订单哦!

密码登录
登录注册
其他方式登录
点击 登录注册 即表示同意《亿速云用户服务条款》