R語言的enrichGSEA_pip.r怎么用

這篇文章主要介紹“R語言的enrichGSEA_pip.r怎么用”，在日常操作中，相信很多人在R語言的enrichGSEA_pip.r怎么用問題上存在疑惑，小編查閱了各式資料，整理出簡單好用的操作方法，希望對大家解答”R語言的enrichGSEA_pip.r怎么用”的疑惑有所幫助！接下來，請跟著小編一起來學(xué)習(xí)吧！

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]

參數(shù)說明：

-a 輸入差異基因分析所以的結(jié)果； 必須含有l(wèi)og2FC這列差異倍數(shù)信息，用于GSEA排序；

-g 指定 gmt文件  基因集： 更多GSEA功能富集數(shù)據(jù)下載：http://software.broadinstitute.org/gsea/downloads.jsp#msigdb

使用舉例：

#更多GSEA功能富集數(shù)據(jù)下載：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

到此，關(guān)于“R語言的enrichGSEA_pip.r怎么用”的學(xué)習(xí)就結(jié)束了，希望能夠解決大家的疑惑。理論與實(shí)踐的搭配能更好的幫助大家學(xué)習(xí)，快去試試吧！若想繼續(xù)學(xué)習(xí)更多相關(guān)知識，請繼續(xù)關(guān)注億速云網(wǎng)站，小編會(huì)繼續(xù)努力為大家?guī)砀鄬?shí)用的文章！