From 765760df33e855d2ce74c46230f73b0376fa4961 Mon Sep 17 00:00:00 2001
From: Elias Projahn <elias@johrpan.de>
Date: Wed, 2 Oct 2024 11:01:58 +0200
Subject: [PATCH] scripts: Add Venn diagram

---
 scripts/comparison.R | 7 +++++++
 1 file changed, 7 insertions(+)

diff --git a/scripts/comparison.R b/scripts/comparison.R
index 44710a2..ece759c 100644
--- a/scripts/comparison.R
+++ b/scripts/comparison.R
@@ -17,6 +17,13 @@ datasets <- list(
   "Joshi 2022" = j2022
 )
 
+VennDiagram::venn.diagram(datasets, filename = NULL, disable.logging = TRUE) |>
+  ggplot2::ggsave(file = here("scripts/output/venn.svg"), device = "svg")
+
+partitions <- VennDiagram::get.venn.partitions(datasets) |> data.table()
+genes_venn <- partitions[1]$..values..[[1]]
+write(genes_venn, file = here("scripts/output/genes_venn.txt"))
+
 ranking_gtex <- ubigen::rank_genes(ubigen::gtex_all)
 ranking_cmap <- ubigen::rank_genes(ubigen::cmap)