# HG changeset patch # User proteore # Date 1521544900 14400 # Node ID 8a229e9a0d4c0e0d643d1538a223320a40960a64 # Parent 69cf9e6283f8295bc9509c12f3d7aa4ce0ac5b2b planemo upload commit 8dc1650e0adc5a095a0d6457fdce849fe96c6b3b-dirty diff -r 69cf9e6283f8 -r 8a229e9a0d4c sel_ann_hpa.R --- a/sel_ann_hpa.R Tue Mar 20 07:01:46 2018 -0400 +++ b/sel_ann_hpa.R Tue Mar 20 07:21:40 2018 -0400 @@ -38,7 +38,6 @@ res.Level<-subset(res.Tissue, Level==level) } if (length(level)>1) { - print(level) res.Level<-subset(res.Tissue, Level %in% level) } @@ -46,7 +45,6 @@ res.Rel<-subset(res.Level, Reliability==reliability) } if (length(reliability)>1) { - print(reliability) res.Rel<-subset(res.Level, Reliability %in% reliability) } @@ -135,7 +133,7 @@ # Extract input input_type = args$input_type if (input_type == "list") { - list_id = strsplit(args$input, " +")[[1]] + list_id = strsplit(args$input, "[ \t\n]+")[[1]] } else if (input_type == "file") { filename = args$input @@ -143,14 +141,12 @@ header = args$header file = readfile(filename, header) list_id = c() - print(file) list_id = sapply(strsplit(file[,column_number], ";"), "[", 1) } input = list_id # Read reference file reference_file = read.table(args$ref_file, header = TRUE, sep = "\t", stringsAsFactors = FALSE, fill = TRUE) - print(colnames(reference_file)) # Extract other options atlas = args$atlas