R Диаграмма Венна Ошибка - тройная диаграмма Венна

Question

Я изо всех сил пытался нарисовать простую диаграмму Венна, и я получаю ошибки, которые я не понимаю.

Я не могу понять, что не так в моих кодах (ниже) по сравнению с написанным в Rdocumentation для функции диаграммы draw.triple.venn - по моему мнению, код в порядке ( https://www.rdocumentation.org/packages/VennDiagram/versions/1.6.20/topics/draw.triple.venn).

library(VennDiagram)
library(gridExtra)

body(draw.triple.venn)[[78]] <- substitute(cell.labels <- paste0(areas," : ", round( 100*areas/sum(areas), 1), "%"))
# this line of code adds %s to each of the cells in the Venn Diagram

g <- draw.triple.venn(
    area1 = 2951, 
    area2 = 2764, 
    area3 = 2764,
    n12 = 719, 
    n23 = 807, 
    n13 = 1034, 
    n123 = 325,
    category = c("Mental Health", "Community", "Social Care"),
    fill = c("blue", "red", "green"), 
    cat.col = c("blue", "red", "green"), 
    lty = "blank",
    euler.d = TRUE, 
    scaled = TRUE, 
    cex = 2, 
    cat.cex = 2);

grid.arrange(gTree(children = g), main = "Breakdown of WSIC Service Users", sub = "By table")

Сообщение об ошибке, которое я получаю и не могу понять, что это значит:

Error in if (max.x - min.x >= max.y - min.y) { : 
 missing value where TRUE/FALSE needed
In addition: Warning messages:
1: In sqrt(r1^2 - (l.x.cept.13 - x.centres[1])^2) : NaNs produced
2: In sqrt(r1^2 - (l.x.cept.13 - x.centres[1])^2) : NaNs produced

Спасибо!

Answer 1

Существует сообщение об ошибке с этой проблемой здесь. Я попытался воспроизвести пример, и я получаю ту же ошибку. Если вы хотите попробовать мою библиотеку nVennR, вы можете получить следующую диаграмму:

library(nVennR)
myV <- createVennObj(nSets = 3, sNames = c("Mental Health", "Community", "Social Care"), sSizes = c(0, 598, 913, 807, 873, 1034, 719, 325))
myV <- plotVenn(nVennObj = myV, borderWidth = 2, setColors = c('blue', 'red', 'green'))

В результате получается файл SVG:

Есть некоторые параметры визуализации, как вы можете видеть в виньетке. Там вы также можете найти обоснование чисел в sSizes, Есть и другие режимы ввода данных. Самый простой способ - предоставить список списков с элементами каждого набора (вы можете увидеть один пример в начале виньетки). Пакет рассчитает все регионы на диаграмме для вас, а затем вы можете запросить элементы, принадлежащие каждому региону.