优化的带有内部标签的VennDiagram r

Question

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

我有这个：

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)

我使用带有gplots()的{​​{1}}函数从包venn获得的。但是，在venn函数中，我似乎无法用标签的名称替换计数。我使用了十字形= TRUE，它通过参数的描述应该有效，但它没有（尽管如果我查看变量c，信息就在那里）。

simplify = TRUE

问题：使用 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包中的simplify函数中的venn参数完全相同？

问题2：使用gplots包中的venn函数，有没有办法显示每个元素的名称而不是元素计数？就像我在'venn.diagram＆＃39;功能

提前致谢，

Answer 1

这是我的方法，到目前为止还没有解决方案而是黑客攻击。

# 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"))

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

编辑： 如果你的情节看起来是一样的，你可以检查一下venn函数的源代码（通过点击F2在RStudio中）并复制粘贴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))