make feature lists unique

epigen · May 27, 2024 · bd35288 · bd35288
1 parent fa938d9
commit bd35288
Showing 1 changed file with 2 additions and 0 deletions.
diff --git a/workflow/scripts/aggregate.R b/workflow/scripts/aggregate.R
@@ -87,6 +87,7 @@ for (group in unique(dea_filtered_results$group)){
     for (direction in unique(dea_filtered_results$direction)){
 
         tmp_features <- dea_filtered_results[(dea_filtered_results$group==group) & (dea_filtered_results$direction==direction),"feature"]
+        tmp_features <- unique(tmp_features)
         write(tmp_features, file.path(results_path,paste0(group,"_",direction,"_features.txt")))
     }
 }
@@ -97,6 +98,7 @@ if("feature_name" %in% colnames(dea_results)){
 
             tmp_features <- dea_filtered_results[(dea_filtered_results$group==group) & (dea_filtered_results$direction==direction),"feature_name"]
             tmp_features <- tmp_features[tmp_features != ""]
+            tmp_features <- unique(tmp_features)
             write(tmp_features, file.path(results_path,paste0(group,"_",direction,"_features_annot.txt")))
         }
     }