R语言画棒棒糖图展示snp在基因上的位置是怎样的

这篇文章将为大家详细讲解有关R语言画棒棒糖图展示snp在基因上的位置是怎样的,文章内容质量较高,因此小编分享给大家做个参考,希望大家阅读完这篇文章后对相关知识有一定的了解。

创新互联服务项目包括揭阳网站建设、揭阳网站制作、揭阳网页制作以及揭阳网络营销策划等。多年来,我们专注于互联网行业,利用自身积累的技术优势、行业经验、深度合作伙伴关系等,向广大中小型企业、政府机构等提供互联网行业的解决方案,揭阳网站推广取得了明显的社会效益与经济效益。目前,我们服务的客户以成都为中心已经辐射到揭阳省份的部分城市,未来相信会继续扩大服务区域并继续获得客户的支持与信任!

今天在一个微信群里看见有人问下面两幅图用R语言如何实现

R语言画棒棒糖图展示snp在基因上的位置是怎样的  
image.png
R语言画棒棒糖图展示snp在基因上的位置是怎样的  
image.png
 第一幅就是散点图

两列数据x,y用来控制散点的位置,一列变量用来映射颜色,Frequency用来映射大小。这个实现起来相对比较容易。还有一列用来添加文字标签

下面我们构造一份数据集用来画图

df<-data.frame(pvalue=sample(seq(0,1,by=0.001),20,replace=F),
              Effect_size=sample(seq(-2,2,by=0.05),20,replace = F),
              frequency=sample(seq(0,1,by=0.01),20,replace=F),
              color=sample(c("A","B"),20,replace = T),
              text=LETTERS[1:20])
df              
library(ggplot2)
library(ggrepel)
ggplot(df,aes(x=Effect_size,y=pvalue,size=frequency))+
 geom_point(aes(color=color),alpha=0.5)+
 theme_bw()+
 scale_size_continuous(range = c(5,10))+
 geom_text_repel(aes(label=text),size=5)
 
R语言画棒棒糖图展示snp在基因上的位置是怎样的  
image.png
 第二幅图如何实现我自己还真没有思路。

印象里好像是看见过有公众号推文介绍这个图。我翻了翻之前转发过的朋友圈还真找到了。有一个专门的包 trackViewer

这个包对应的论文是

trackViewer: a Bioconductor package for interactive and integrative visualization of multi-omics data. Nature Methods https://www.nature.com/articles/s41592-019-0430-y

github主页 https://github.com/jianhong/trackViewer

一份很详细的帮助文档 https://bioconductor.org/packages/release/bioc/vignettes/trackViewer/inst/doc/trackViewer.html

文档中实现的图片有

R语言画棒棒糖图展示snp在基因上的位置是怎样的  
image.png
R语言画棒棒糖图展示snp在基因上的位置是怎样的  
image.png

等等。

今天就不重复文档中的每个例子了。知道有这个包可以实现这个功能就可以了。重复一个最基本的例子

library(trackViewer)
SNP <- c(10, 12, 1400, 1402)
sample.gr <- GRanges("chr1", IRanges(SNP, width=1, names=paste0("snp", SNP)))
features <- GRanges("chr1", IRanges(c(1, 501, 1001),
                                   width=c(120, 400, 405),
                                   names=paste0("block", 1:3)))
features$fill <- c("#FF8833", "#51C6E6", "#DFA32D")
sample.gr$color <- sample.int(6, length(SNP), replace=TRUE)
sample.gr$border <- sample(c("gray80", "gray30"), length(SNP), replace=TRUE)
lolliplot(sample.gr, features)
 
R语言画棒棒糖图展示snp在基因上的位置是怎样的

关于R语言画棒棒糖图展示snp在基因上的位置是怎样的就分享到这里了,希望以上内容可以对大家有一定的帮助,可以学到更多知识。如果觉得文章不错,可以把它分享出去让更多的人看到。

分享文章:R语言画棒棒糖图展示snp在基因上的位置是怎样的
网站路径:https://www.cdcxhl.com/article26/jocccg.html

成都网站建设公司_创新互联,为您提供App开发标签优化云服务器网站收录网站策划建站公司

广告

声明:本网站发布的内容(图片、视频和文字)以用户投稿、用户转载内容为主,如果涉及侵权请尽快告知,我们将会在第一时间删除。文章观点不代表本网站立场,如需处理请联系客服。电话:028-86922220;邮箱:631063699@qq.com。内容未经允许不得转载,或转载时需注明来源: 创新互联

搜索引擎优化