From 2611e4e909adeb0e15c686891ca8c235dda046b9 Mon Sep 17 00:00:00 2001
From: aduvermy <arnaud.duvermy@ens-lyon.fr>
Date: Wed, 9 Feb 2022 16:21:19 +0000
Subject: [PATCH] add fun for visualizing res
---
src/visualization_fun.R | 18 ++++++++++++++++++
1 file changed, 18 insertions(+)
create mode 100644 src/visualization_fun.R
diff --git a/src/visualization_fun.R b/src/visualization_fun.R
new file mode 100644
index 0000000..fe11135
--- /dev/null
+++ b/src/visualization_fun.R
@@ -0,0 +1,18 @@
+library(ggplot2)
+
+# visualization functions
+mu_effect_visualization <- function(mu_effect_res){
+ label_wrap <- c("mu observed", "N gene DE", "min(|logFC|)", "var observed")
+ names(label_wrap) <- c("mu_observ", "res_DEA", "statistical_power", "var_observ")
+ figure = mu_effect_res %>% ggplot(., aes(x=vec_of_mu, y = value, col=factor(N_rep))) +
+ geom_point() + facet_wrap(~variable, scales = "free_y", labeller = labeller(variable = label_wrap)) + labs(color = "N replicates")
+ return(figure)
+}
+
+size_effect_visualization <- function(alpha_effect_res){
+ label_wrap <- c("mu observed", "N gene DE", "min(|logFC|)", "var observed")
+ names(label_wrap) <- c("mu_observ", "res_DEA", "statistical_power", "var_observ")
+ figure = alpha_effect_res %>% ggplot(., aes(x=vec_of_alpha, y = value, col=factor(N_rep))) +
+ geom_point() + facet_wrap(~variable, scales = "free_y", labeller = labeller(variable = label_wrap)) + labs(color = "N replicates")
+ return(figure)
+}
\ No newline at end of file
--
GitLab