setwd("F://张秀秀//过程性文件//8//8//OX1931B10")
new_gremline<-read.table("germline_mutations_hg38_new.txt",sep = "\t",header = T)
old_gremline<-read.table("germline_mutations_hg38_old.txt",sep = "\t",header = T)
dim(new_gremline) #2572
dim(old_gremline) #2367
#为了确保准确性,这次不进行位点的比较(因为不同染色体会有不同的位点。这也没办法说服我自己),
#而是将我们关注的几个信息合并后进行比较
new_gremline<-within(new_gremline, mutation<- paste(Chr,Start,Ref,Alt,sep=" "))
old_gremline<-within(old_gremline, mutation<- paste(Chr,Start,Ref,Alt,sep=" "))
new_mutation<-new_gremline$mutation #需要进行画图的对象1
old_mutation<-old_gremline$mutation #需要进行画图的对象2
#然后就是一个简单的venn plot
library(VennDiagram)
# 二维韦恩图
venn.plot <- venn.diagram(x = list(Removal=new_mutation,Original=old_mutation),filename = NULL,lwd = 4,fill = c("cornflowerblue", "darkorchid1"),alpha = 0.75,label.col = "black",cex = 2,fontfamily = "serif",fontface = "bold",cat.col = c("cornflowerblue", "darkorchid1"),cat.cex = 1,cat.fontfamily = "serif",cat.fontface = "bold",cat.dist = c(0.03, 0.03),cat.pos = c(-20, 14),main = "OX1931B10 patient(9.6%)"
)
grid.draw(venn.plot)
grid.newpage()

最后的结果:

R语言绘图 | Venn图相关推荐

  1. R语言绘图—热图绘制

    R语言绘图-热图绘制 原创 R语言与医学生  2022-06-24 20:59 发表于广东 热图使用颜色来表示二维图中第三个变量的变化和大小.热图能清楚直观的看到数据的变化.今天就简单介绍下R语言中热 ...

  2. R语言绘图——小提琴图

    在过去的某一天,许久没有灵感的我,忽然想起来,这个世界上还有一个网站叫做R Gallery,于是就有了这个系列! 系列1:小提琴图--所谓小提琴图,其实是最简单的箱线图的一种变种,在箱线图显示的中位数 ...

  3. R语言绘图 | 折线图画法,如何画出你满意的图?

    作者:丁点helper 来源:丁点帮你 上篇文章中,我们用世界银行的例子给大家介绍了散点图怎么画,那折线图呢? 绘制折线图 如果将散点图上的点从左往右连接起来,就会得到一个折线图.今天我们以R中自带的 ...

  4. R语言绘图—“金字塔图”

    这里我主要模拟一些数据,大家可以根据自己数据的情况来绘制. #模拟数据 set.seed(1234) mydata<-data.frame(age=rnorm(n=6000,mean=40,sd ...

  5. (生物信息学)R语言绘图初-中-高级——3-10分文章必备——Venn图(韦恩图)(初级)

    生物信息学文章的发表要求除了思路和热点以外,图片绘制是否精美也是十分重要的,本专栏为(生物信息学)R语言绘图初-中-高级--3-10分文章必备,主要通过大量文献,总结3-10分文章中高频出现的各种图片 ...

  6. R语言可视化韦恩图实战:VENN DIAGRAM

    R语言可视化韦恩图实战:VENN DIAGRAM #韦恩图 # load Venn diagram package library("VennDiagram")# move to ...

  7. R语言绘图patchwork拼图详解快速实现组合图拼接

    pdf教程下载 此文内容来自微信公众号:R语言搬运工,扫码关注公众号浏览更多精彩内容** 我们在实际科研绘图或者写报告的时候,往往会同时出多幅统计图,如何将他们拼在一起是即将或者正在面临的问题.R语言 ...

  8. R语言绘图——组间差异箱线图图,小提琴图、聚类热力图、相关系数热力图

    生物医学工程常用图 箱线图:使用ggboxplot()函数来展示数据在多组间的差异 小提琴图:使用ggviolin()函数来展示数据在多组间的比较 热力图:pheatmap()展示聚类热图反映数据在多 ...

  9. R语言学习 - 热图美化 (数值标准化和调整坐标轴顺序)

    生物信息学习的正确姿势 NGS系列文章包括NGS基础.在线绘图.转录组分析 (Nature重磅综述|关于RNA-seq你想知道的全在这).ChIP-seq分析 (ChIP-seq基本分析流程).单细胞 ...

  10. R语言绘图:28个实用程序包

    本文转载自" 阿虎定量笔记",作者郑连虎,转载己获授权.其他人转载请联系原作者,点我跳转原文. 全文注释 # 示例参考帮助文档 # 我用Rstudio重现了所有示例 # 初识R语言 ...

最新文章

  1. SQLite的sqlite_master表
  2. Cisco1242胖AP转瘦AP并加入Cisco无线控制器
  3. linux命令uname
  4. Reflection
  5. 性能测试工具Curl-Loader
  6. 判断一段文件是UTF-8编码还是GB2312的编码方式
  7. gradle依赖冲突的解决方式
  8. 计算机主机箱外面有什么区别吗,电脑主机箱后面有三个插头空 都是干什么用的啊?...
  9. 分析一个文本文件中各个单词出现的频率,把频率最高的10个词打印出来
  10. 120 Python程序中的线程操作-队列
  11. 使用apache benchmark(ab) 测试报错: apr_socket_recv: Connection timed out (110)
  12. Android Studio下载Gradle一直失败,使用正确的下载地址
  13. Java 通过银行卡号获取银行名称
  14. linux红帽子镜像下载,红帽Red Hat Linux相关产品iso镜像下载
  15. 小米手机证书信任设置在哪里_通过此设置帮助更快地启动小米手机
  16. 银河系召唤“四有青年”,宇宙原力即将觉醒
  17. 微信公众号与微信小程序相辅相成
  18. idm 无法下载此文件 备注: SSL_connect error 5
  19. MWC 2018火热开幕,AI取代手机成主角,华为超越苹果称霸!
  20. 用sed命令shell扩大王林快码词库

热门文章

  1. noip2011 公交观光
  2. DSPE-PEG3-Mal,带有DSPE的基团小分子PEG试剂
  3. 基于Python+kociemba+opencv的3阶魔方自动还原
  4. 用v-modle做一个简单的计算器
  5. 怎么用wifi进e站_科普微e站 | 注意!路由器千万别再这样用,当心信号越来越差!...
  6. 语音识别sdk_语音识别 sdk_离线语音识别sdk - 云+社区 - 腾讯云
  7. iOS - 手机摇一摇
  8. uni-app以h5网页的方式发布企业微信应用
  9. 生产力工具:功能强大又好用的浏览器网页截屏工具FireShot
  10. jquery 校验中国身份证号码