R/venn_4.R

In cfbeuchel/CarlHelpR: A Package With Functions

#' 4-Way Venn Diagram
#'
#' Plot a non-proportional 4-Way Venn Diagram.
#' This function is from the package sytemPipeR found here:
#' http://faculty.ucr.edu/~tgirke/Documents/R_BioCond/My_R_Scripts/vennDia.R
#' or
#' https://bioconductor.org/packages/release/bioc/html/systemPipeR.html
#'
#' Please check the systemPipeR documentation for further details
#'
#' @param x1 first vector to be matched
#'
#' @param y1 second vector to be matched
#'
#' @param z1 third vector to be matched
#'
#' @param w1 third vector to be matched
#'
#' @param mytitle Title of the plot
#'
#' @param mylabels Labels for plotting
#'
#' @param plotte ...
#'
#' @export

venn4 <- function(x1, y1, z1, w1, mytitle = "4-Way Venn Diagram", mylabels = NA, plotte = T) {
    # 13/07/03 150119 vector check
    if (all(is.vector(x1) | is.factor(x1), is.vector(y1) | is.factor(y1), is.vector(z1) | is.factor(z1), is.vector(w1) | is.factor(w1)) == F)
        stop("All input data must be vectors...")

    if (is.na(mylabels[1]))
        mylabels = c(deparse(substitute(x1)), deparse(substitute(y1)), deparse(substitute(z1)), deparse(substitute(w1)))
    qlist <- venndiagram(x = x1, y = y1, z = z1, w = w1, unique = T, title = mytitle, labels = mylabels, plot = plotte, lines = c(2, 3, 4, 6), lcol = c(2, 3, 4, 6), tcol = 1,
        lwd = 3, cex = 1, printsub = T, type = "4el")
    qlist
}