R語言如何繪制GO弦圖

這篇文章主要介紹“R語言如何繪制GO弦圖”的相關知識，小編通過實際案例向大家展示操作過程，操作方法簡單快捷，實用性強，希望這篇“R語言如何繪制GO弦圖”文章能幫助大家解決問題。

R包GOplot將GO富集分析結果可視化，繪制弦圖

使用方法:

$Rscript $scriptdir/GOChord_plot.r -h
usage: /share/nas1/fangs/test/Ferroptosis-Related/TCGA-ESCC/scripts/GOChord_plot.r
       [-h] -g GO_DATA -d DEG_DATA [--no_logFC] [--gene GENE] [--GO GO]
       [-z GENE_SIZE] [-S GENE_SPACE] [-s GO_SPACE] [-t TITLE] [-o OUTDIR]
       [-p PREFIX] [-H HEIGHT] [-W WIDTH]

GOChord plot:https://www.億速云.com/article/1532

optional arguments:
  -h, --help            show this help message and exit
  -g GO_DATA, --GO_data GO_DATA
                        input a GO enrichment analysis data,the data with
                        columns for 'category','ID','term',adjusted
                        p-value('adj_pval') and 'genes'[required]
  -d DEG_DATA, --deg_data DEG_DATA
                        input the result of differential expression analysis
                        corresponding to GO enrichment analysis,the data with
                        columns for 'ID', 'logFC'[required]
  --no_logFC            There is no logFC value[optional,default False]
  --gene GENE           The number of GOterms associated with a
                        gene[optional,default 3]
  --GO GO               The number of genes associated with a
                        GOterms[optional,default 5]
  -z GENE_SIZE, --gene_size GENE_SIZE
                        The size of the gene name[optional,default 5]
  -S GENE_SPACE, --gene_space GENE_SPACE
                        The spacing between genes[optional,default 0.25]
  -s GO_SPACE, --GO_space GO_SPACE
                        The spacing between GOterms[optional,default 0.02]
  -t TITLE, --title TITLE
                        the title of the String figure[optional,default
                        GOChord_plot]
  -o OUTDIR, --outdir OUTDIR
                        output file directory[optional,default cwd]
  -p PREFIX, --prefix PREFIX
                        out file name prefix[optional,default GOChord_plot]
  -H HEIGHT, --height HEIGHT
                        the height of the String figure[optional,default 18]
  -W WIDTH, --width WIDTH
                        the width of the String figure[optional,default 20]

參數說明：

-g 輸入GO富集分析結果文件，文件的列中必須含有"category"(BP/CC/MF)、''ID"(GO ID)、"term"(功能描述)、"adj_pval"(校正后的p值)和"genes"(每個GO term中含有的gene ID)：

-d 輸入用于GO富集分析的基因文件，文件列中必須含有""ID"(gene ID)和"logFC"(差異倍數），若沒有差異倍數，會添加"logFC"列并將差異倍數設置為0：

--no_logFC

輸入文件是否有差異倍數，有則無需設置，沒有則需要設置該參數

--gene 5 --GO 4

--gene是設置每個gene至少對應的GO term數，如基因對應的GO term數沒有達到設置的值，那么該gene將不會在圖中顯示

--GO是設置每個GO term至少含有的gene數，如GO term含有的gene數沒有達到設置的值，那么該GO term也將不會在圖中顯示

-z 5 -S 0.25 -s 0.02

-z是設置圖中gene ID的字體大小，-S是設置圖中gene間的間隔大小，-s是設置圖中GO term間的間隔大小

使用舉例：

$Rscript $scriptdir/GOChord_plot.r -g GO/TP_vs_NT_GO_enrichment_stat.tsv \
    -d ../03.deg/fer_deg_gene.tsv --gene 5 --GO 4

關于“R語言如何繪制GO弦圖”的內容就介紹到這里了，感謝大家的閱讀。如果想了解更多行業(yè)相關的知識，可以關注億速云行業(yè)資訊頻道，小編每天都會為大家更新不同的知識點。