您好,登錄后才能下訂單哦!
這篇文章給大家分享的是有關(guān)如何使用R包ggplot2繪制啞鈴圖展示基因的CNV突變頻率的內(nèi)容。小編覺得挺實(shí)用的,因此分享給大家做個(gè)參考,一起跟隨小編過來看看吧。
R包ggplot2繪制啞鈴圖展示基因的CNV突變頻率
$Rscript scripts/cnv_dumbbell_plot.r -h usage: scripts/cnv_dumbbell_plot.r [-h] -c CNV_DATA -g GENE_SYMBOL [-H HEIGHT] [-W WIDTH] [-o OUTDIR] [-p PREFIX] CNV_Dumbbell_plot:https://www.億速云.com/article/1572 optional arguments: -h, --help show this help message and exit -c CNV_DATA, --cnv_data CNV_DATA input cnv data file path[required] -g GENE_SYMBOL, --gene_symbol GENE_SYMBOL Enter the file path that contains the gene symbol[required] -H HEIGHT, --height HEIGHT the height of dumbbell plot[optional,default 6] -W WIDTH, --width WIDTH the width of dumbbell plot[optional,default 8] -o OUTDIR, --outdir OUTDIR output file directory[optional,default cwd] -p PREFIX, --prefix PREFIX out file name prefix[optional,default m6a_cnv]
-c 輸入拷貝數(shù)變異數(shù)據(jù)經(jīng)過GISTIC軟件分析過后的結(jié)果文件:
Gene Symbol | Locus ID | Cytoband | TCGA-3Z-A93Z-01A-11D-A36W-01 | TCGA-6D-AA2E-01A-11D-A36W-01 | TCGA-A3-3306-01A-01D-0858-01 | TCGA-A3-3307-01A-01D-0858-01 | TCGA-A3-3308-01A-02D-1322-01 |
ACAP3 | 116983 | 1p36.33 | -0.001 | 0.002 | 0.009 | 0.015 | 0.005 |
ACTRT2 | 140625 | 1p36.32 | -0.001 | 0.002 | 0.009 | 0.015 | 0.005 |
-g 輸入用于繪制啞鈴圖的基因文件:
ID | group |
RBM15 | W |
RBM15B | W |
METTL14 | W |
-H和-W 指定生成圖片的高和寬,默認(rèn)高為6,寬為8
Rscript ../scripts/cnv_dumbbell_plot.r -g ../m6a.tsv \ -c ../01.TCGA_data/TP_CNV_Gistic2_Level_4/all_data_by_genes.txt
感謝各位的閱讀!關(guān)于“如何使用R包ggplot2繪制啞鈴圖展示基因的CNV突變頻率”這篇文章就分享到這里了,希望以上內(nèi)容可以對(duì)大家有一定的幫助,讓大家可以學(xué)到更多知識(shí),如果覺得文章不錯(cuò),可以把它分享出去讓更多的人看到吧!
免責(zé)聲明:本站發(fā)布的內(nèi)容(圖片、視頻和文字)以原創(chuàng)、轉(zhuǎn)載和分享為主,文章觀點(diǎn)不代表本網(wǎng)站立場,如果涉及侵權(quán)請(qǐng)聯(lián)系站長郵箱:is@yisu.com進(jìn)行舉報(bào),并提供相關(guān)證據(jù),一經(jīng)查實(shí),將立刻刪除涉嫌侵權(quán)內(nèi)容。