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

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

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

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

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

使用方法:

$Rscriptscripts/cnv_dumbbell_plot.r-husage:scripts/cnv_dumbbell_plot.r[-h]-cCNV_DATA-gGENE_SYMBOL[-HHEIGHT][-WWIDTH][-oOUTDIR][-pPREFIX]CNV_Dumbbell_plot:https://www.恰卡编程网.com/article/1572optionalarguments:-h,--helpshowthishelpmessageandexit-cCNV_DATA,--cnv_dataCNV_DATAinputcnvdatafilepath[required]-gGENE_SYMBOL,--gene_symbolGENE_SYMBOLEnterthefilepaththatcontainsthegenesymbol[required]-HHEIGHT,--heightHEIGHTtheheightofdumbbellplot[optional,default6]-WWIDTH,--widthWIDTHthewidthofdumbbellplot[optional,default8]-oOUTDIR,--outdirOUTDIRoutputfiledirectory[optional,defaultcwd]-pPREFIX,--prefixPREFIXoutfilenameprefix[optional,defaultm6a_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突变频率”这篇文章就分享到这里了,希望以上内容可以对大家有一定的帮助,让大家可以学到更多知识,如果觉得文章不错,可以把它分享出去让更多的人看到吧!

发布于 2022-03-17 21:16:56
收藏
分享
海报
0 条评论
36
上一篇:如何使用ggraph包绘制展示相关性的网络图 下一篇:python如何设置休眠时间
目录

    0 条评论

    本站已关闭游客评论,请登录或者注册后再评论吧~

    忘记密码?

    图形验证码