R/get_venn.R

In frhl/our: genetic data analysis

# load packages
#' @title get_venn
#' @description get's venn statistics
#' @param lst a names list
#' @export

get_venn <- function(lst){

  #haplo = fread(IN_FILE1)
  #gnomad = fread(IN_FILE2)
  #mydat = fread(IN_FILE3)

  require(gplots)

  print(names(lst))
  write('\nCombining venn..',stdout())
  d <- venn(lst, show.plot = F)
  myattr = lengths(attributes(d)$intersections)
  write('\nVENN DIAGRAM OPTIONS [start]',stdout())
  print(myattr)
  write('\nVENN DIAGRAM OPTIONS [end]',stdout())
  res = stack(attributes(d)$intersections)
  return(res)

}