Commit fd8d9310 authored by Lena Hersemann's avatar Lena Hersemann

adjusted geom_text position of volcano plots

parent 9f08a922
......@@ -40,7 +40,7 @@ devtools::source_url("https://raw.githubusercontent.com/holgerbrandl/datautils/v
# devtools::source_url("https://raw.githubusercontent.com/holgerbrandl/datautils/v1.40/R/bio/diffex_commons.R")
# source("/lustre/projects/bioinfo/herseman/scripts/ngs_tools/dge_workflow/diffex_commons.R", chdir = T)
devtools::source_url("https://git.mpi-cbg.de/bioinfo/ngs_tools/raw/v4/dge_workflow/diffex_commons.R", chdir = T)
devtools::source_url("https://git.mpi-cbg.de/bioinfo/ngs_tools/raw/v5/dge_workflow/diffex_commons.R", chdir = T)
## process input arguments
......@@ -594,9 +594,9 @@ hitCounts = filter(deResults, is_hit) %>%
rename(hits = n) %>%
## complement hit counts with directed contrasts without any diffex genes
right_join(merge(mutate(contrasts),data_frame(c1_overex=c(T,F)), by=NULL)) %>% mutate(hits=if_else(is.na(hits), as.integer(0), hits)) %>%
merge(data.frame(c1_overex = c(T, F), x_pos = c(maxX + 2, -maxX - 2)))
merge(data.frame(c1_overex = c(T, F), x_pos = c(maxX - maxX*0.1, -maxX + maxX*0.1)))
hitCounts %<>% mutate(y_pos = quantile(- log10(deResults$pvalue), 0.99, na.rm = TRUE))
hitCounts %<>% mutate(y_pos = quantile(- log10(deResults$pvalue), 0.987, na.rm = TRUE))
#+ fig.width=16, fig.height=14
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment