Repository 'mspurity_puritya'
hg clone https://toolshed.g2.bx.psu.edu/repos/computational-metabolomics/mspurity_puritya

Changeset 1:3d6356bcf998 (2019-12-05)
Previous changeset 0:56cce1a90b73 (2019-11-27) Next changeset 2:9b067f4c5f7f (2020-04-01)
Commit message:
"planemo upload for repository https://github.com/computational-metabolomics/mspurity-galaxy commit 20d432cdc6326114d05c26fd51889b872513b57d"
modified:
flagRemove.R
macros.xml
purityA.R
b
diff -r 56cce1a90b73 -r 3d6356bcf998 flagRemove.R
--- a/flagRemove.R Wed Nov 27 14:26:04 2019 -0500
+++ b/flagRemove.R Thu Dec 05 07:42:36 2019 -0500
[
@@ -138,7 +138,7 @@
     chosen_blank <- samplelist_blank[samplelist_blank %in% xset@phenoData$class]
     if (length(chosen_blank)>1){
         print('ERROR: only 1 blank is currently allowed to be used with this tool')
-        exit()
+        quit()
     }
     blank_class <- as.character(chosen_blank)
     print(blank_class)
@@ -171,7 +171,7 @@
     grp_peaklist <- ffrm_out[[2]]
     removed_peaks <- ffrm_out[[3]]
 
-    save.image(file=file.path(opt$out_dir, 'xset_filtered.RData'))
+    save.image(file=file.path(opt$out_dir, 'xset_filtered.RData'), version=2)
 
     # grpid needed for mspurity ID needed for deconrank... (will clean up at some up)
     peak_pth <- file.path(opt$out_dir, 'peaklist_filtered.tsv')
b
diff -r 56cce1a90b73 -r 3d6356bcf998 macros.xml
--- a/macros.xml Wed Nov 27 14:26:04 2019 -0500
+++ b/macros.xml Thu Dec 05 07:42:36 2019 -0500
b
@@ -1,7 +1,7 @@
 <?xml version="1.0"?>
 <macros>
     <token name="@TOOL_VERSION@">1.12.0</token>
-    <token name="@GALAXY_TOOL_VERSION@">0</token>
+    <token name="@GALAXY_TOOL_VERSION@">1</token>
 
     <xml name="requirements">
         <requirements>
b
diff -r 56cce1a90b73 -r 3d6356bcf998 purityA.R
--- a/purityA.R Wed Nov 27 14:26:04 2019 -0500
+++ b/purityA.R Thu Dec 05 07:42:36 2019 -0500
b
@@ -6,9 +6,9 @@
   make_option(c("-o", "--out_dir"), type="character"),
   make_option("--mzML_files", type="character"),
   make_option("--galaxy_names", type="character"),
-  make_option("--minOffset", default=0.5),
-  make_option("--maxOffset", default=0.5),
-  make_option("--ilim", default=0.05),
+  make_option("--minOffset", type="numeric"),
+  make_option("--maxOffset", type="numeric"),
+  make_option("--ilim", type="numeric"),
   make_option("--iwNorm", default="none", type="character"),
   make_option("--exclude_isotopes", action="store_true"),
   make_option("--isotope_matrix", type="character"),
@@ -22,18 +22,26 @@
 opt <- parse_args(OptionParser(option_list=option_list))
 print(opt)
 
-minOffset = as.numeric(opt$minOffset)
-maxOffset = as.numeric(opt$maxOffset)
 
 if (opt$iwNorm=='none'){
     iwNorm = FALSE
     iwNormFun = NULL
 }else if (opt$iwNorm=='gauss'){
     iwNorm = TRUE
-    iwNormFun = msPurity::iwNormGauss(minOff=-minOffset, maxOff=maxOffset)
+    if (is.null(opt$minOffset) || is.null(opt$maxOffset)){
+      print('User has to define offsets if using Gaussian normalisation')
+    }else{
+      iwNormFun = msPurity::iwNormGauss(minOff=-as.numeric(opt$minOffset), 
+                                      maxOff=as.numeric(opt$maxOffset))
+    }
 }else if (opt$iwNorm=='rcosine'){
     iwNorm = TRUE
-    iwNormFun = msPurity::iwNormRcosine(minOff=-minOffset, maxOff=maxOffset)
+    if (is.null(opt$minOffset) || is.null(opt$maxOffset)){
+      print('User has to define offsets if using R-cosine normalisation')
+    }else{
+      iwNormFun = msPurity::iwNormRcosine(minOff=-as.numeric(opt$minOffset), 
+                                          maxOff=as.numeric(opt$maxOffset))
+    }
 }else if (opt$iwNorm=='QE5'){
     iwNorm = TRUE
     iwNormFun = msPurity::iwNormQE.5()