Диаграмма Венна с метками предметов

Предположим, у меня есть два вектора

foo <- c('a','b','c','d')
baa <- c('a','e','f','g')

Кто-нибудь знает способ создания диаграммы Венна, но визуализирует элементы вектора внутри диаграммы?

Вот так? (сделано в powerpoint)

3 ответа

Решение

С использованием 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)

введите описание изображения здесь

Быстрое решение с использованием 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)

Который производит

Очевидно, что этот метод быстро выйдет из-под контроля с большим количеством категорий и пересечений.

Другой способ, который я считаю более простым, это использование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))

Другие вопросы по тегам