1、为什么要画韦恩图?
在对两个或多个样品之间基因集合比较,查看多少基因是几个样品共有,多少基因是几个样品特有,基因之间的交集是怎样的,我们会采用韦恩图来表示。
2、实现过程
这是我们需要复现的图,主要是复现上面部分
1、先用非生物学的例子来演示
#install.packages("ggVennDiagram") ##下载并加载包
library(ggVennDiagram)
mydata<-list(Agroup=1:31,Bgroup=20:80)
ggVennDiagram(mydata)
这是直接画出来的图
#改变颜色
ggVennDiagram(mydata)+scale_fill_gradient(low = "#F4C4B6", high = "#34C4B6")
还是很丑,好一点,但不多
#绘制3个集合的图
mydata<list(Agroup=1:31,Bgroup=20:80,Cgroup=26:55)
ggVennDiagram(mydata)+scale_fill_gradient(low = "#F4C4B6", high = "#34C4B6")
以上参考了别人,但是我找不到了,不是有意,侵删
2、拿需要复现的来练习练习
DOI: 10.1002/pmic.202000210(补充文件表格中的Tab1)
install.packages("ggVennDiagram") ##下载并加载包
library(ggVennDiagram)
library(readr)
test <- read_csv("R/test.csv")#导入数据
#按筛选条件分组数据
sEV <- filter(test,Mean_EV > 1 )
WCL <- filter(test,Mean_WCL > 1)
#合并两组数据
veen_data <- list(sEV_G=sEV$Gene.name,WCL_G=WCL$Gene.name)
#画图
ggVennDiagram(veen_data,
set_color = "black",#标签颜色
set_size =5,#标签大小
label = c("both")#选择标签形式"both", "count", "percent", "none"
)+
scale_fill_gradient(low = "#EDF8E9", high = "#74C476")