Mercurial > repos > proteore > proteore_maps_visualization
changeset 10:a4db9efe1e43 draft
"planemo upload commit 78a50fa16dcea6afd93a2f4de17ed193a2aa7dd3-dirty"
| author | proteore |
|---|---|
| date | Mon, 02 Mar 2020 09:23:05 +0000 |
| parents | 9ac0ade20b3f |
| children | 7d04efde2526 |
| files | kegg_maps_visualization.R kegg_maps_visualization.xml |
| diffstat | 2 files changed, 18 insertions(+), 3 deletions(-) [+] |
line wrap: on
line diff
--- a/kegg_maps_visualization.R Wed Jan 08 09:03:41 2020 +0000 +++ b/kegg_maps_visualization.R Mon Mar 02 09:23:05 2020 +0000 @@ -177,6 +177,15 @@ 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)
--- a/kegg_maps_visualization.xml Wed Jan 08 09:03:41 2020 +0000 +++ b/kegg_maps_visualization.xml Mon Mar 02 09:23:05 2020 +0000 @@ -1,9 +1,10 @@ -<tool id="kegg_maps_visualization" name="Pathways visualization" version="2020.01.08"> - <description>of (differentially expressed) genes/proteins [KEGG]</description> +<tool id="kegg_maps_visualization" name="KEGG pathways mapping and visualization" version="2020.03.02"> + <description>[PathView]</description> <requirements> <requirement type="package" version="3.6.2" >R</requirement> <requirement type="package" version="1.26.0">bioconductor-pathview</requirement> <requirement type="package" version="1.26.0">bioconductor-keggrest</requirement> + <requirement type="package" version="1.30.10">r-biocmanager</requirement> </requirements> <command detect_errors="exit_code"><