Hay una forma, si bien trabajosa, al menos sí, dinámica. En primer lugar, vamos a necesitar los conjuntos en una lista, y configurar un tamaño de letra para las etiquetas un numero más pequeño que en tu ejemplo:
library(VennDiagram)
BM1 <- c('Cistaceae', 'Oleaceae','Fagaceae','Rhamnaceae','Asteraceae','Asteraceae')
BM2 <- c('Boraginaceae','Poaceae','Poaceae')
BM3 <- c('Poaceae','Asteraceae','Chenopodiaceae','Asteraceae','Poaceae','Solanaceae','Solanaceae','Asteraceae','Rubiaceae')
BM4<- c('Amaranthaceae','Chenopodiaceae','Convolvulaceae','Convolvulaceae','Cucurbitaceae','Cucurbitaceae','Cucurbitaceae','Cucurbitaceae','Boraginaceae','Boraginaceae','Solanaceae')
lista_conjuntos <- list(Oak=BM1, Wasteland=BM2,Edge=BM3, Crop=BM4)
v2 <- venn.diagram(lista_conjuntos,
fill = c("red", "blue","green","purple"),
alpha = c(0.5, 0.5, 0.5,0.5), cat.cex = 1.5, cex=.8,
filename=NULL)
Ahora, teniendo lista_conjuntos
, vamos a construir las etiquetas de cada área de intersección, para esto podemos usar calculate.overlap()
del propio paquete VennDiagram
:
labels <- lapply(calculate.overlap(lista_conjuntos), function(x) {paste(unique(x), collapse='\n')})
tail(labels, 3)
$`a14`
[1] ""
$a1
[1] "Rubiaceae"
$a3
[1] "Amaranthaceae\nConvolvulaceae\nCucurbitaceae"
Lo que obtenemos es una lista con los nombres de los elementos de cada intersección, dónde además en el nombre del elemento (a1..an
) tenemos un índice que no permitira luego acceder a los elementos del diagrama de venn, es decir el objeto v2
. Este objeto es un Grob List
, una lista de objetos gráficos, los primeros objetos de la lista son los polígonos que se dibujan y los siguientes son las etiquetas (que por defecto tienen las cantidad de elementos de cada área). Estos elementos don del tipo text
, y nos interesará ubicar el primero de estos:
min(which(unlist(lapply(v2, function(x) {class(x)[[1]]}=="text"))))
[1] 9
El elemento 9 será la intersección a1
, el 10 es a2
y así. Obviamente esto variará en función a la cantidad de conjuntos. Entonces podremos configurar la etiqueta deseada y dibujar ahora sí el diagrama:
for (x in names(labels)) {
v2[[ start + as.integer(gsub("a", "",x)) ]]$label = labels[[x]]
}
grid.newpage()
grid.draw(v2)
Salida: