From 54b331504175cfe022fb65c91e2745f4aba4df88 Mon Sep 17 00:00:00 2001 From: Elias Projahn Date: Sun, 1 Dec 2024 17:02:46 +0100 Subject: [PATCH] scripts: Fix typo gsea/gtex --- scripts/{gsea_vs_cmap.R => gtex_vs_cmap.R} | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) rename scripts/{gsea_vs_cmap.R => gtex_vs_cmap.R} (96%) diff --git a/scripts/gsea_vs_cmap.R b/scripts/gtex_vs_cmap.R similarity index 96% rename from scripts/gsea_vs_cmap.R rename to scripts/gtex_vs_cmap.R index 8fccfef..11ead5f 100644 --- a/scripts/gsea_vs_cmap.R +++ b/scripts/gtex_vs_cmap.R @@ -1,7 +1,7 @@ library(data.table) library(here) -i_am("scripts/gsea_vs_cmap.R") +i_am("scripts/gtex_vs_cmap.R") ranking_gtex <- ubigen::rank_genes(ubigen::gtex_all) ranking_cmap <- ubigen::rank_genes(ubigen::cmap) @@ -48,7 +48,7 @@ data[percentile_gtex < 0.95 & percentile_cmap >= 0.95, group := "0_1"] data[percentile_gtex >= 0.95 & percentile_cmap < 0.95, group := "1_0"] data[percentile_gtex >= 0.95 & percentile_cmap >= 0.95, group := "1_1"] -fwrite(data, file = here("scripts/output/gsea_vs_cmap_groups.csv")) +fwrite(data, file = here("scripts/output/gtex_vs_cmap_groups.csv")) write(data[group == "0_0", gene], here("scripts/output/genes_0_0.txt")) write(data[group == "0_1", gene], here("scripts/output/genes_0_1.txt")) @@ -151,7 +151,7 @@ fig <- plotly::plot_ly() |> plotly::save_image( fig, - file = here("scripts/output/gsea_vs_cmap.svg"), + file = here("scripts/output/gtex_vs_cmap.svg"), width = 6.27 * 72, height = 6.27 * 72, scale = 96 / 72