R语言的enrichGSEA_pip.r怎么用

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

这篇文章主要介绍“R语言的enrichGSEA_pip.r怎么用”,在日常操作中,相信很多人在R语言的enrichGSEA_pip.r怎么用问题上存在疑惑,小编查阅了各式资料,整理出简单好用的操作方法,希望对大家解答”R语言的enrichGSEA_pip.r怎么用”的疑惑有所帮助!接下来,请跟着小编一起来学习吧!

enrichGSEA_pip.R  GSEA富集分析

使用说明:

$Rscript $scriptdir/enrichGSEA_pip.r -h
usage: /work/my_stad_immu/scripts/enrichGSEA_pip.r [-h] -a all.deg.file -g
                                                   gmtfile [-p pvalueCutoff]
                                                   [-t pvalueCutoff]
                                                   [-n prefix] [-o outdir]
                                                   [-H height] [-W width]
GSEA enrich analysis :https://www.亿速云.com/article/1504
optional arguments:
  -h, --help            show this help message and exit
  -a all.deg.file, --all.deg.file all.deg.file
                        all diff express gene list file,must include log2FC
                        column, required
  -g gmtfile, --gmtfile gmtfile
                        GSEA gmtfile function class file, required
  -p pvalueCutoff, --pvalueCutoff pvalueCutoff
                        pvalue cutoff on enrichment tests to report,
                        [optional, default: 0.1 ]
  -t top, --top top
                        top NES for barplot [optional, default:10 ]
  -n prefix, --prefix prefix
                        the output file prefix [optional, default: GSEA ]
  -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 5]

参数说明:

-a 输入差异基因分析所以的结果; 必须含有log2FC这列差异倍数信息,用于GSEA排序;

-g 指定 gmt文件  基因集: 更多GSEA功能富集数据下载:http://software.broadinstitute.org/gsea/downloads.jsp#msigdb

使用举例:

#更多GSEA功能富集数据下载:http://software.broadinstitute.org/gsea/downloads.jsp#msigdb
wget -c https://data.broadinstitute.org/gsea-msigdb/msigdb/release/7.4/c2.cp.kegg.v7.4.symbols.gmt

Rscript $scriptdir/enrichGSEA_pip.r   --all.deg.file $workdir/04.deg/S1_vs_S2.all.tsv \
  --gmtfile  c2.cp.kegg.v7.4.symbols.gmt -o GSEA -n S1_vs_S2_KEGG -p 0.05

到此,关于“R语言的enrichGSEA_pip.r怎么用”的学习就结束了,希望能够解决大家的疑惑。理论与实践的搭配能更好的帮助大家学习,快去试试吧!若想继续学习更多相关知识,请继续关注亿速云网站,小编会继续努力为大家带来更多实用的文章!

推荐阅读:
  1. R语言的glm()函数怎么用
  2. R语言的Dataframe怎么用

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

r语言

上一篇:OTUtable标准化方法是什么

下一篇:R语言的enrichKEGG_pip.r如何使用

相关阅读

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

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