Minor correction of the male / female GD ranking.

16-Data_integration_all/refine_GD_rankings.R

@@ -64,6 +64,7 @@ for (i in seq_len(length(config$integrations))) {
   # We then merge across the two genders.
   GD_FM <- merge(x = GD_females, y = GD_males, by = "Gene", all = TRUE, suffixes = c("",".M")) %>%
     mutate(ranking_value = ranking_value + ranking_value.M) %>%
+    mutate(adj_P_value = pmin(adj_P_value, adj_P_value.M)) %>%
     select(Gene, log_fold_change, P_value, adj_P_value, pi_value, ranking_value) %>%
     arrange(desc(ranking_value))
   k              <- 5
@@ -77,6 +78,7 @@ for (i in seq_len(length(config$integrations))) {
   rm(res_females_fn, k, limma_analysis)
   GD_MF <- merge(x = GD_males, y = GD_females, by = "Gene", all = TRUE, suffixes = c("",".F")) %>%
     mutate(ranking_value = ranking_value + ranking_value.F) %>%
+    mutate(adj_P_value = pmin(adj_P_value, adj_P_value.F)) %>%
     select(Gene, log_fold_change, P_value, adj_P_value, pi_value, ranking_value) %>%
     arrange(desc(ranking_value))
   k              <- 6