假设我有两个向量
foo <- c('a','b','c','d')
baa <- c('a','e','f','g')
Run Code Online (Sandbox Code Playgroud)
有没有人知道产生维恩图的方法,但在图中可视化的矢量项.
像这样?(用powerpoint制作)
use*_*650 15
使用包中的venn.diagram
功能的快速解决方案VennDiagram
.标签(计数)在函数中是硬编码的,因此不能使用函数参数进行更改.但是对于这样一个简单的例子,你可以改变grobs
自己.
library(VennDiagram)
# your data
foo <- c('a','b','c','d')
baa <- c('a','e','f','g')
# Generate plot
v <- venn.diagram(list(foo=foo, baa=baa),
fill = c("orange", "blue"),
alpha = c(0.5, 0.5), cat.cex = 1.5, cex=1.5,
filename=NULL)
# have a look at the default plot
grid.newpage()
grid.draw(v)
# have a look at the names in the plot object v
lapply(v, names)
# We are interested in the labels
lapply(v, function(i) i$label)
# Over-write labels (5 to 7 chosen by manual check of labels)
# in foo only
v[[5]]$label <- paste(setdiff(foo, baa), collapse="\n")
# in baa only
v[[6]]$label <- paste(setdiff(baa, foo) , collapse="\n")
# intesection
v[[7]]$label <- paste(intersect(foo, baa), collapse="\n")
# plot
grid.newpage()
grid.draw(v)
Run Code Online (Sandbox Code Playgroud)
哪个产生
显然,这种方法很快就会失控,因为更多的类别和交叉点.
使用RAM
包:
library(RAM)
foo <- c('a','b','c','d')
baa <- c('a','e','f','g')
group.venn(list(foo=foo, baa=baa), label=TRUE,
fill = c("orange", "blue"),
cat.pos = c(0, 0),
lab.cex=1.1)
Run Code Online (Sandbox Code Playgroud)
小智 5
我发现更简单的另一种方法是将ggvenn
包与show_elements = TRUE
.
library(ggvenn)
library(RColorBrewer)
AA <- c("hi","foo", "bar","yep","woo","hoo")
BB <- c("baa","yep", "woo","yes")
CC <- c("yes","foo","hi","woo", "huh")
x <- list(AA=AA , BB=BB , CC=CC)
ggvenn(x, show_elements = T, label_sep = "\n", fill_color = brewer.pal(name="Set2",n=3))
Run Code Online (Sandbox Code Playgroud)