johrpan/ubigen
1
0
Fork 0
mirror of https://github.com/johrpan/ubigen.git synced 2025-10-26 19:57:24 +01:00
Code Issues Projects Releases Packages Wiki Activity Actions

scripts: Fix typo gsea/gtex

Browse source
This commit is contained in:
Elias Projahn 2024-12-01 17:02:46 +01:00
parent 7c8cc99d34
commit 54b3315041
1 changed files with 3 additions and 3 deletions
Download patch file Download diff file Expand all files Collapse all files

6
scripts/gsea_vs_cmap.R → scripts/gtex_vs_cmap.R
View file

@ -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