diff --git a/src/visualization_fun.R b/src/visualization_fun.R
new file mode 100644
index 0000000000000000000000000000000000000000..fe111357b1a6d97984fc906c7198760f61521ebe
--- /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