Mercurial > repos > proteore > proteore_expression_levels_by_tissue

changeset 6:8a229e9a0d4c draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload commit 8dc1650e0adc5a095a0d6457fdce849fe96c6b3b-dirty
author proteore
date Tue, 20 Mar 2018 07:21:40 -0400
parents 69cf9e6283f8
children f7e93829327e
files sel_ann_hpa.R
diffstat 1 files changed, 1 insertions(+), 5 deletions(-) [+]
line wrap: on
line diff
--- 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