diff kegg_maps_visualization.R @ 6:eb140217a9f1 draft default tip

"planemo upload commit 1f9a98db2496834cc38a281cbe3d4b11de2bc3d2-dirty"

--- a/kegg_maps_visualization.R	Fri Jan 24 05:20:10 2020 -0500
+++ b/kegg_maps_visualization.R	Thu Mar 05 10:10:20 2020 -0500
@@ -172,11 +172,20 @@
 }
 
 get_limit <- function(mat) {
-  min = min(apply(mat,2,min))
-  max = max(apply(mat,2,max))
+  min = min(apply(mat,2,min,na.rm=TRUE))
+  max = max(apply(mat,2,max,na.rm=TRUE))
   return(c(min,max))
 }
 
+check_pathway_ids<- function(pathways_ids) {
+  problematic_pathways <- c("04215","04723")
+  to_remove <- intersect(pathways_ids,problematic_pathways)
+  if (length(to_remove) > 0){ print(paste("Pathway(s)",to_remove,"have been remove from input")) }
+  pathways_ids <- pathways_ids[which(!pathways_ids %in% problematic_pathways)]
+  if (length(pathways_ids) == 0){stop("No pathways ids to process")}
+  return (pathways_ids)
+}
+
 get_args <- function(){
   
   ## Collect arguments
@@ -238,7 +247,7 @@
     pathways_file = read_file(args$pathways_input,header2)
     ids <- sapply(rapply(strsplit(clean_bad_character(pathways_file[,pathway_col]),","),c), function(x) remove_kegg_prefix(x),USE.NAMES = FALSE)
   }
-  if (args$native_kegg) { ids <- ids[ids != "04215"] }
+  if (args$native_kegg) {ids = check_pathway_ids(ids)} #remove problematic pathways
   pathways_list <- read_file(args$pathways_list,F)
   if (!is.null(args$id_list)) {
     id_list <- get_list_from_cp(args$id_list)