温馨提示×

温馨提示×

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

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

如何使用R包ggplot2绘制哑铃图展示基因的CNV突变频率

发布时间:2022-03-17 14:45:24 来源:亿速云 阅读:688 作者:小新 栏目:开发技术

这篇文章给大家分享的是有关如何使用R包ggplot2绘制哑铃图展示基因的CNV突变频率的内容。小编觉得挺实用的,因此分享给大家做个参考,一起跟随小编过来看看吧。

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 输入拷贝数变异数据经过GISTIC软件分析过后的结果文件:

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 指定生成图片的高和宽,默认高为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

感谢各位的阅读!关于“如何使用R包ggplot2绘制哑铃图展示基因的CNV突变频率”这篇文章就分享到这里了,希望以上内容可以对大家有一定的帮助,让大家可以学到更多知识,如果觉得文章不错,可以把它分享出去让更多的人看到吧!

向AI问一下细节

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

AI