本文小编为大家详细介绍“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如何使用”文章已经介绍完毕,想要掌握这篇文章的知识点还需要大家自己动手实践使用过才能领会,如果想了解更多相关内容的文章,欢迎关注亿速云行业资讯频道。
亿速云「云服务器」,即开即用、新一代英特尔至强铂金CPU、三副本存储NVMe SSD云盘,价格低至29元/月。点击查看>>
免责声明:本站发布的内容(图片、视频和文字)以原创、转载和分享为主,文章观点不代表本网站立场,如果涉及侵权请联系站长邮箱:is@yisu.com进行举报,并提供相关证据,一经查实,将立刻删除涉嫌侵权内容。