带有内部标签 r 的优化 VennDiagram

Question

我正在尝试以优化的方式（见下文）绘制维恩图，并将案例作为内部标签（而不是每个交叉点的案例数）。我知道他们每个人都有 post，但没有一个解决方案允许我同时做。

我有这个：

x <- list()
x$A <- as.character(c("Per_36","Cent","CeM","vDG","LAVL","RSGd"))
x$B <- as.character(c("vCA1","DLE","Per_36","vDG","DIE","Per_35"))
x$C <- as.character(c("vCA1","Cg1","LAVL", "RSGc", "RSGd","Per_35","Per_36"))
x$D <- as.character(c("Por","Cg1","RSGc","LAVL","Per_35","RSGd","Per_36"))


require(VennDiagram)
v0 <-venn.diagram(x, lwd = 3, col = c("red", "green", "orange", "blue"),
                fill = c("red", "blue", "green", "orange"), apha = 0.5, filename = NULL)
grid.draw(v0)
overlaps <- calculate.overlap(x)
overlaps <- rev(overlaps)
for (i in 1:length(overlaps)){
  v0[[i+8]]$label <- paste(overlaps[[i]], collapse = "\n")
}

grid.newpage()
grid.draw(v0)

我得到以下输出：

关于维恩图的组织方式我想这样做：

c <- venn(x, simplify = TRUE, small = 0.5, intersections = TRUE)

我使用 venn 函数和 simplify = TRUE 从包 gplots() 中获得。但是，在 venn 函数中，我似乎无法用标签名称替换计数。我使用了 intersections = TRUE，根据参数的描述，它应该起作用，但它不起作用（尽管如果我查看变量 c，信息就在那里）。

    Logical flag indicating if the returned object should have the attribute 
"individuals.in.intersections" featuring for every set a list of individuals
 that are assigned to it.

问题：使用 VennDiagram 包，有没有办法与 gplots 包的 venn 函数中的 simplify 参数完全相同？

问题 2：使用 gplots 包中的 venn 函数，有没有办法显示每个元素的名称而不是元素计数？就像我在 'venn.diagram' 函数中所做的那样？

提前致谢，

Answer 1

这是我的方法，到目前为止还没有解决方案，而是一种 hack。

# Print a venn and save it to an object
a <- venn(list(letters[1:5], letters[3:8]))
# save the intersections
b <- attr(a, "intersections")

# find the coordinates
s <- seq(0,500,100); abline(h=s); text(s, y=s, x=0)
s <- seq(0,500,50);  abline(v=s); text(s, y=0, x=s)

# the hack, destroy the venn to avoid the plotting of the internal numbers 
rownames(a) <- letters[1:nrow(a)]
a
plot.venn(a)
>Error in data[n, 1] : subscript out of bounds

# include the internal labels
text(200,300,paste(b$`01`,collapse = "\n"))
text(200,200,paste(b$`11`,collapse = "\n"))
text(200,100,paste(b$`10`,collapse = "\n"))

多个 venns 很烦人。否则，您可以将维恩保存为 .svg 并使用 inkscape 或类似软件对其进行编辑，或者通过电子邮件询问开发人员。

编辑：如果你的图看起来总是一样，你可以检查 venn 函数的源代码（在 RStudio 中点击 F2）并复制粘贴 4 和 5 圆 venns 的位置并将标签函数 lab("1000", data) 替换为你想要的标签。

4 个圈子：

      text(35, 250, lab("1000", data))
      text(140, 315, lab("0100", data))
      text(260, 315, lab("0010", data))
      text(365, 250, lab("0001", data))
      text(90, 280, lab("1100", data), cex = small)
      text(95, 110, lab("1010", data))
      text(200, 50, lab("1001", data), cex = small)
      text(200, 290, lab("0110", data))
      text(300, 110, lab("0101", data))
      text(310, 280, lab("0011", data), cex = small)
      text(130, 230, lab("1110", data))
      text(245, 75, lab("1101", data), cex = small)
      text(155, 75, lab("1011", data), cex = small)
      text(270, 230, lab("0111", data))
      text(200, 150, lab("1111", data))

编辑

现在我会改用 ggplot 解决方案

ggVennDiagram::ggVennDiagram(x)

带有内部标签 r 的优化 VennDiagram

optimized VennDiagram with internal labels r

r

venn-diagram

编辑