From 789a4bf91ffd6b4b4ccbc739c7ec6a20895d26ec Mon Sep 17 00:00:00 2001
From: AlexWeinreb <alexis.weinreb@yale.edu>
Date: Wed, 12 Oct 2022 11:52:59 -0400
Subject: [PATCH] add some test code

---
 Functions.R | 28 ++++++++++++++++++++++++++--
 1 file changed, 26 insertions(+), 2 deletions(-)

diff --git a/Functions.R b/Functions.R
index d9a854d..f083d98 100644
--- a/Functions.R
+++ b/Functions.R
@@ -224,10 +224,34 @@ perform_de <- function(ident.1, ident.2, method, ...){
 # head(res2)
 # head(res3)
 # 
-# list(wilcox = res2$gene[res2$p_val_adj < .05],
-#      edgeR = res1$gene[res1$p_val_adj < .05],
+# list(wilcox = res2$gene[res2$FDR < .05],
+#      edgeR = res1$gene[res1$FDR < .05],
 #      sc = res3$gene[res3$p_val_adj < .05]) |>
 #   eulerr::euler() |>
 #   plot()
+# 
+# left_join(res3 |>
+#                   select(gene, sc_wilcox = avg_logFC),
+#                 res1 |>
+#                   select(gene, pb_edger = logFC),
+#                 by = "gene") |>
+#   ggplot() +
+#   geom_point(aes(x = sc_wilcox, y = pb_edger))
+# 
+# left_join(res3 |>
+#             select(gene, sc_wilcox = avg_logFC),
+#           res2 |>
+#             select(gene, pb_wilcox = log2FC),
+#           by = "gene") |>
+#   ggplot() +
+#   geom_point(aes(x = sc_wilcox, y = pb_wilcox))
+# left_join(res1 |>
+#             select(gene, pb_edger = logFC),
+#           res2 |>
+#             select(gene, pb_wilcox = log2FC),
+#           by = "gene") |>
+#   ggplot() +
+#   geom_point(aes(x = pb_edger, y = pb_wilcox))
+