R语言画棒棒糖图展示snp在基因上的位置是怎样的
发表于:2024-09-24 作者:千家信息网编辑
千家信息网最后更新 2024年09月24日,这篇文章将为大家详细讲解有关R语言画棒棒糖图展示snp在基因上的位置是怎样的,文章内容质量较高,因此小编分享给大家做个参考,希望大家阅读完这篇文章后对相关知识有一定的了解。今天在一个微信群里看见有人问
千家信息网最后更新 2024年09月24日R语言画棒棒糖图展示snp在基因上的位置是怎样的
这篇文章将为大家详细讲解有关R语言画棒棒糖图展示snp在基因上的位置是怎样的,文章内容质量较高,因此小编分享给大家做个参考,希望大家阅读完这篇文章后对相关知识有一定的了解。
第一幅就是散点图今天在一个微信群里看见有人问下面两幅图用R语言如何实现
两列数据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)
第二幅图如何实现我自己还真没有思路。印象里好像是看见过有公众号推文介绍这个图。我翻了翻之前转发过的朋友圈还真找到了。有一个专门的包
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
文档中实现的图片有
等等。
今天就不重复文档中的每个例子了。知道有这个包可以实现这个功能就可以了。重复一个最基本的例子
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在基因上的位置是怎样的就分享到这里了,希望以上内容可以对大家有一定的帮助,可以学到更多知识。如果觉得文章不错,可以把它分享出去让更多的人看到。
位置
语言
文档
基因
棒棒糖
例子
内容
数据
文章
更多
知识
篇文章
帮助
不错
主页
公众
功能
印象
变量
图片
数据库的安全要保护哪些东西
数据库安全各自的含义是什么
生产安全数据库录入
数据库的安全性及管理
数据库安全策略包含哪些
海淀数据库安全审计系统
建立农村房屋安全信息数据库
易用的数据库客户端支持安全管理
连接数据库失败ssl安全错误
数据库的锁怎样保障安全
论文版权和数据库有关系吗
网络安全设备销售合同
苹果电脑怎么登admin服务器
上海浦东曙光服务器回收
天津服务器回收价格
厦门中尧软件开发有限公司
网络安全为人民的实施方案
常州大学嵌入式软件开发
t310服务器怎么做raid
软件开发项目怎么分割
软件开发价格计算公式
上海商标互联网科技有限公司
数据库取千万数据的速度
数据库自增长10w条数据
企业网络安全论文怎么选题
数据库的设计会遇到什么问题
哈尔滨电子病例软件开发
化解网络安全风险隐患
信元网络安全等保
数据库组成的特征
自动操作服务器
山西服务器虚拟化系统云主机
公用网是网络安全产品吗
access数据库管理
软件开发的生命周期是指什么
继续教育 网络安全法题库
招联金融软件开发工资
sql有批量附加数据库吗
欧阳裕尧 广州腾科网络技术有限公司
数据库组成的特征