R语言

关注公众号 jb51net

关闭
首页 > 软件编程 > R语言 > R语言绘制维恩图ggvenn

R语言绘制维恩图ggvenn示例详解

作者:Wei_Sun

这篇文章主要为大家介绍了R语言绘制维恩图ggvenn示例详解,有需要的朋友可以借鉴参考下,希望能够有所帮助,祝大家多多进步,早日升职加薪

引言

韦恩图,Venn diagram,常用图的一种,用来展示集合之间的特异性和共同性。现在有很多在线的网站都可以绘制,但是R来画也方便,其中ggvenn是基于ggplot2的专门绘制韦恩图的R包。

官方网站:https://github.com/yanlinlin82/ggvenn

1.安装

ggvenn在CRAN上,直接用Install.packages就可以完成安装:

> install.packages("ggvenn")
> library(ggvenn)

2.基础用法

ggvenn支持list和data.frame两种数据格式。这里以三个基因文件为例:

读取三个基因文件:

> set1<-read.csv("gene_a.csv")
> set2<-read.csv("gene_b.csv")
> set3<-read.csv("gene_c.csv")

提取每个文件的基因id,创建list:

> dat <- list( A = set1$gene_id,  B = set2$gene_id, C = set3$gene_id)

绘图:

> ggvenn(dat)

绘制部分list:

> ggvenn(dat, c("A", "B")) 

3.图形美化

填充

边框

集合名

集合内文本

百分比

> ggvenn(dat,show_percentage = T,
  stroke_color = "white",
  stroke_size = 0.5,
  fill_color = c("#E41A1C","#1E90FF","#FF8C00"),
  set_name_color =c("#E41A1C","#1E90FF","#FF8C00"), 
  set_name_size = 15,text_size=6)

4.提取交集部分并输出

> A_B <- as.data.frame(intersect(set1$gene_id, set2$gene_id))
> write.csv(A_B,"A_B_common_gene.csv",row.names = F)

以上就是R语言绘制维恩图ggvenn示例详解的详细内容,更多关于R绘制维恩图ggvenn的资料请关注脚本之家其它相关文章!

您可能感兴趣的文章:
阅读全文