annotate maldi_quant_preprocessing.xml @ 4:82b1de5c142b draft default tip

planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f1e1cd260ef2884d0ba12e2b614df3c72d0934dc
author galaxyp
date Sat, 04 Mar 2023 19:13:31 +0000
parents aa81463e21ea
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
4
82b1de5c142b planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f1e1cd260ef2884d0ba12e2b614df3c72d0934dc
galaxyp
parents: 3
diff changeset
1 <tool id="maldi_quant_preprocessing" name="MALDIquant preprocessing" version="@VERSION@.0">
3
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
2 <description>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
3 Preprocessing of mass-spectrometry imaging data
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
4 </description>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
5 <macros>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
6 <import>maldi_macros.xml</import>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
7 </macros>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
8 <expand macro="requirements"/>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
9 <command detect_errors="exit_code">
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
10 <![CDATA[
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
11 cat '${maldi_quant_preprocessing}' &&
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
12 #if $infile.ext == 'imzml'
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
13 cp '${infile.extra_files_path}/imzml' infile.imzML &&
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
14 cp '${infile.extra_files_path}/ibd' infile.ibd &&
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
15 #elif $infile.ext == 'analyze75'
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
16 cp '${infile.extra_files_path}/hdr' infile.hdr &&
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
17 cp '${infile.extra_files_path}/img' infile.img &&
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
18 cp '${infile.extra_files_path}/t2m' infile.t2m &&
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
19 du infile.hdr &&
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
20 du infile.img &&
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
21 du -s -B1 infile.hdr &&
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
22 #end if
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
23 Rscript "${maldi_quant_preprocessing}" &&
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
24
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
25 mkdir $outfile_imzml.files_path &&
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
26 mv ./out.imzMl "${os.path.join($outfile_imzml.files_path, 'imzml')}" | true &&
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
27 mv ./out.ibd "${os.path.join($outfile_imzml.files_path, 'ibd')}" | true &&
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
28 echo "imzML file:" > $outfile_imzml &&
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
29 ls -l "$outfile_imzml.files_path" >> $outfile_imzml
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
30 ]]>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
31 </command>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
32 <configfiles>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
33 <configfile name="maldi_quant_preprocessing"><![CDATA[
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
34
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
35 @R_IMPORTS@
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
36
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
37 #if $restriction_conditional.restriction == 'restrict':
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
38
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
39 print('Reading mask region')
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
40 ## Import imzML file
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
41
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
42 coordinate_matrix = as.matrix(read.delim("$restriction_conditional.coordinates_file", header = $restriction_conditional.coordinates_header, stringsAsFactors = FALSE))[,1:2]
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
43 coordinate_matrix = coordinate_matrix[,c($restriction_conditional.column_x, $restriction_conditional.column_y)]
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
44
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
45 maldi_data = importImzMl('infile.imzML',
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
46 coordinates = coordinate_matrix)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
47 pixelnames = paste0("x = ", coordinates(maldi_data)[,1],", y = ", coordinates(maldi_data)[,2])
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
48
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
49 #else:
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
50
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
51 print('Reading entire file')
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
52 #if $infile.ext == 'imzml'
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
53 ## Import imzML file
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
54 maldi_data = import( 'infile.imzML', type="imzML" )
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
55 coordinates_info = cbind(coordinates(maldi_data)[,1:2], c(1:length(maldi_data)))
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
56 #elif $infile.ext == 'analyze75'
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
57 ## Import analyze7.5 file
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
58 maldi_data = importAnalyze( 'infile.hdr' )
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
59 coordinates_info = cbind(coordinates(maldi_data)[,1:2], c(1:length(maldi_data)))
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
60 #end if
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
61
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
62 #end if
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
63
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
64 ## Quality control plots during preprocessing
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
65
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
66 pdf("prepro_qc_plot.pdf", fonts = "Times", pointsize = 12)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
67 plot(0,type='n',axes=FALSE,ann=FALSE)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
68
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
69 ## if no filename is given, name of file in Galaxy history is used
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
70 #set $filename = $infile.display_name
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
71 title(main=paste("$filename"))
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
72
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
73
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
74 #################### Preprocessing methods #####################################
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
75
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
76 ## QC plot on input file
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
77 avgSpectra = averageMassSpectra(maldi_data,method="mean")
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
78 plot(avgSpectra, main="Average spectrum for input file")
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
79
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
80 pixel_number = length(maldi_data)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
81 minmz = round(min(unlist(lapply(maldi_data,mass))), digits=4)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
82 maxmz = round(max(unlist(lapply(maldi_data,mass))), digits=4)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
83 mean_features = round(length(unlist(lapply(maldi_data,mass)))/length(maldi_data), digits=2)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
84 number_features = length(unique(unlist(lapply(maldi_data,mass))))
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
85 medint = round(median(unlist(lapply(maldi_data,intensity))), digits=2)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
86 inputdata = c(minmz, maxmz,number_features,mean_features,medint, pixel_number)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
87 QC_numbers= data.frame(inputdata = c(minmz, maxmz,number_features, mean_features, medint, pixel_number))
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
88 vectorofactions = "inputdata"
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
89
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
90
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
91 #for $method in $methods:
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
92
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
93 #if str( $method.methods_conditional.method ) == 'Transformation':
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
94
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
95 print('transforming')
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
96 ##transformation
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
97 maldi_data = transformIntensity(maldi_data, method="$method.methods_conditional.transform_method")
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
98 ## QC plot and numbers
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
99 avgSpectra = averageMassSpectra(maldi_data,method="mean")
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
100 plot(avgSpectra, main="Average spectrum after transformation")
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
101 pixel_number = length(maldi_data)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
102 minmz = round(min(unlist(lapply(maldi_data,mass))), digits=4)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
103 maxmz = round(max(unlist(lapply(maldi_data,mass))), digits=4)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
104 mean_features = round(length(unlist(lapply(maldi_data,mass)))/length(maldi_data), digits=2)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
105 medint = round(median(unlist(lapply(maldi_data,intensity))), digits=2)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
106 number_features = length(unique(unlist(lapply(maldi_data,mass))))
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
107 transformed = c(minmz, maxmz,number_features,mean_features,medint,pixel_number)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
108 QC_numbers= cbind(QC_numbers, transformed)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
109 vectorofactions = append(vectorofactions, "transformed")
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
110
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
111
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
112 #elif str( $method.methods_conditional.method ) == 'Smoothing':
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
113
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
114 print('smoothing')
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
115 ##smoothing
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
116
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
117 #if str($method.methods_conditional.methods_for_smoothing.smooth_method ) == 'SavitzkyGolay':
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
118 print('SavitzkyGolay')
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
119
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
120 maldi_data = smoothIntensity(maldi_data,
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
121 method="SavitzkyGolay", polynomialOrder=$method.methods_conditional.methods_for_smoothing.polynomial,
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
122 halfWindowSize=$method.methods_conditional.halfWindowSize)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
123
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
124 #elif str($method.methods_conditional.methods_for_smoothing.smooth_method ) == 'MovingAverage':
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
125 print('MovingAverage')
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
126
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
127 maldi_data = smoothIntensity(maldi_data,
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
128 method="MovingAverage", weighted=$method.methods_conditional.methods_for_smoothing.weighted,
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
129 halfWindowSize=$method.methods_conditional.halfWindowSize)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
130
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
131 #end if
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
132
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
133 ## QC plot and numbers
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
134 avgSpectra = averageMassSpectra(maldi_data,method="mean")
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
135 plot(avgSpectra, main="Average spectrum after smoothing", sub="")
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
136 pixel_number = length(maldi_data)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
137 minmz = round(min(unlist(lapply(maldi_data,mass))), digits=4)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
138 maxmz = round(max(unlist(lapply(maldi_data,mass))), digits=4)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
139 mean_features = round(length(unlist(lapply(maldi_data,mass)))/length(maldi_data), digits=2)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
140 medint = round(median(unlist(lapply(maldi_data,intensity))), digits=2)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
141 number_features = length(unique(unlist(lapply(maldi_data,mass))))
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
142 smoothed = c(minmz, maxmz,number_features,mean_features,medint,pixel_number)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
143 QC_numbers= cbind(QC_numbers, smoothed)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
144 vectorofactions = append(vectorofactions, "smoothed")
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
145
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
146
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
147 #elif str( $method.methods_conditional.method ) == 'Baseline':
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
148
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
149 print('baseline removing')
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
150 ## Remove baseline
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
151
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
152 ## Choose random spectra for QC plots
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
153 random_spectra = sample(1:length(maldi_data), 4, replace=FALSE)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
154
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
155 #if str($method.methods_conditional.methods_for_baseline.baseline_method ) == 'SNIP':
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
156 print('SNIP')
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
157
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
158 par(mfrow = c(2,2))
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
159 for (random_sample in random_spectra){
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
160 maldi_data_baseline = estimateBaseline(maldi_data[[random_sample]],
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
161 method="SNIP", iterations=$method.methods_conditional.methods_for_baseline.iterations)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
162 plot(maldi_data[[random_sample]], sub="", main=paste0("Estimated baseline for spectrum ", random_sample))
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
163 lines(maldi_data_baseline, col="blue", lwd=2)}
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
164
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
165 maldi_data = removeBaseline(maldi_data,
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
166 method="SNIP",
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
167 iterations=$method.methods_conditional.methods_for_baseline.iterations)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
168
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
169
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
170 #elif str($method.methods_conditional.methods_for_baseline.baseline_method ) == 'TopHat':
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
171 print('TopHat')
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
172
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
173 par(mfrow = c(2,2))
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
174 for (random_sample in random_spectra){
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
175 maldi_data_baseline = estimateBaseline(maldi_data[[random_sample]],
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
176 method="TopHat", halfWindowSize=$method.methods_conditional.methods_for_baseline.tophat_halfWindowSize)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
177 plot(maldi_data[[random_sample]], sub="", main=paste0("Estimated baseline for spectrum ", random_sample))
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
178 lines(maldi_data_baseline, col="blue", lwd=2)}
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
179
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
180 maldi_data = removeBaseline(maldi_data,
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
181 method="TopHat",
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
182 halfWindowSize=$method.methods_conditional.methods_for_baseline.tophat_halfWindowSize)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
183
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
184
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
185 #elif str($method.methods_conditional.methods_for_baseline.baseline_method ) == 'ConvexHull':
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
186 print('ConvexHull')
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
187
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
188 par(mfrow = c(2,2))
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
189 for (random_sample in random_spectra){
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
190 maldi_data_baseline = estimateBaseline(maldi_data[[random_sample]],
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
191 method="ConvexHull")
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
192 plot(maldi_data[[random_sample]], sub="", main=paste0("Estimated baseline for spectrum ", random_sample))
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
193 lines(maldi_data_baseline, col="blue", lwd=2)}
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
194
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
195 maldi_data = removeBaseline(maldi_data,
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
196 method="ConvexHull")
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
197
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
198
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
199 #elif str($method.methods_conditional.methods_for_baseline.baseline_method ) == 'median':
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
200 print('median')
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
201
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
202 par(mfrow = c(2,2))
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
203 for (random_sample in random_spectra){
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
204 maldi_data_baseline = estimateBaseline(maldi_data[[random_sample]],
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
205 method="median", halfWindowSize=$method.methods_conditional.methods_for_baseline.median_halfWindowSize)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
206 plot(maldi_data[[random_sample]], sub="", main=paste0("Estimated baseline for spectrum ", random_sample))
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
207 lines(maldi_data_baseline, col="blue", lwd=2)}
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
208
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
209 maldi_data = removeBaseline(maldi_data,
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
210 method="median",
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
211 halfWindowSize=$method.methods_conditional.methods_for_baseline.median_halfWindowSize)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
212
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
213 #end if
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
214
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
215 ## QC plot and numbers
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
216 par(mfrow = c(1,1))
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
217 avgSpectra = averageMassSpectra(maldi_data,method="mean")
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
218 plot(avgSpectra, main="Average spectrum after baseline removal")
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
219 pixel_number = length(maldi_data)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
220 minmz = round(min(unlist(lapply(maldi_data,mass))), digits=4)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
221 maxmz = round(max(unlist(lapply(maldi_data,mass))), digits=4)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
222 mean_features = round(length(unlist(lapply(maldi_data,mass)))/length(maldi_data), digits=2)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
223 medint = round(median(unlist(lapply(maldi_data,intensity))), digits=2)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
224 number_features = length(unique(unlist(lapply(maldi_data,mass))))
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
225 baseline_removed = c(minmz, maxmz,number_features,mean_features,medint,pixel_number)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
226 QC_numbers= cbind(QC_numbers, baseline_removed)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
227 vectorofactions = append(vectorofactions, "bl_removed")
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
228
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
229
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
230 #elif str( $method.methods_conditional.method ) == 'Calibrate':
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
231
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
232 print('calibrate')
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
233 ##calibrate
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
234
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
235 #if str($method.methods_conditional.cond_calibration_range) == "yes":
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
236 ## calibrate only given m/z range
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
237 maldi_data = calibrateIntensity(maldi_data,
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
238 method="$method.methods_conditional.calibrate_method",
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
239 range=c($method.methods_conditional.cond_calibration_range.mass_start, $method.methods_conditional.cond_calibration_range.mass_end))
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
240 #else:
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
241 maldi_data = calibrateIntensity(maldi_data,
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
242 method="$method.methods_conditional.calibrate_method")
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
243 #end if
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
244 ## QC plot and numbers
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
245 avgSpectra = averageMassSpectra(maldi_data,method="mean")
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
246 plot(avgSpectra, main="Average spectrum after normalization")
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
247 pixel_number = length(maldi_data)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
248 minmz = round(min(unlist(lapply(maldi_data,mass))), digits=4)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
249 maxmz = round(max(unlist(lapply(maldi_data,mass))), digits=4)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
250 mean_features = round(length(unlist(lapply(maldi_data,mass)))/length(maldi_data), digits=2)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
251 medint = round(median(unlist(lapply(maldi_data,intensity))), digits=2)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
252 number_features = length(unique(unlist(lapply(maldi_data,mass))))
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
253 intensity_calibrated = c(minmz, maxmz,number_features,mean_features,medint,pixel_number)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
254 QC_numbers= cbind(QC_numbers, intensity_calibrated)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
255 vectorofactions = append(vectorofactions, "calibrated")
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
256
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
257
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
258 #elif str( $method.methods_conditional.method ) == 'Align':
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
259
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
260 print('align')
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
261 ##align spectra with 3 separate functions
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
262
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
263 ## create reference if needed
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
264
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
265 ## 1) detect peaks:
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
266 peaks <- detectPeaks(maldi_data, halfWindowSize=$method.methods_conditional.halfWindowSize,
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
267 method="$method.methods_conditional.peak_method", SNR=$method.methods_conditional.snr)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
268
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
269
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
270 #if str($method.methods_conditional.reference_for_alignment.align_ref) == 'no_reference':
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
271
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
272 ## 2) calculate warping:
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
273 warping_function <- determineWarpingFunctions(peaks,
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
274 tolerance=$method.methods_conditional.tolerance, method="$method.methods_conditional.warping_method",
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
275 allowNoMatches=$method.methods_conditional.allow_nomatch, minFrequency = $method.methods_conditional.reference_for_alignment.min_frequency)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
276
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
277 ## 3) warp spectra:
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
278 maldi_data = warpMassSpectra(maldi_data, warping_function, emptyNoMatches=$method.methods_conditional.empty_nomatch)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
279
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
280
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
281 #elif str($method.methods_conditional.reference_for_alignment.align_ref) == 'yes_reference':
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
282
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
283 ## create reference mass_vector from tabular file
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
284 mass_vector = read.delim("$method.methods_conditional.reference_for_alignment.reference_file", header = $method.methods_conditional.reference_for_alignment.reference_header, stringsAsFactors = FALSE)[,$method.methods_conditional.reference_for_alignment.mz_column]
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
285 int_vector = rep(1,length(mass_vector))
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
286 mass_list = createMassPeaks(mass_vector, int_vector)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
287
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
288 #if str($method.methods_conditional.reference_for_alignment.separate_alignment) == "FALSE"
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
289 print('default alignment')
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
290
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
291 ## 2) calculate warping:
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
292 warping_function <- determineWarpingFunctions(peaks,
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
293 tolerance=$method.methods_conditional.tolerance, method="$method.methods_conditional.warping_method",
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
294 allowNoMatches=$method.methods_conditional.allow_nomatch, reference = mass_list)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
295
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
296 ## 3) warp spectra:
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
297 maldi_data = warpMassSpectra(maldi_data, warping_function, emptyNoMatches=$method.methods_conditional.empty_nomatch)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
298
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
299 #elif str($method.methods_conditional.reference_for_alignment.separate_alignment) == "TRUE"
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
300 print('spectra wise alignment')
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
301
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
302 maldi_data_new_list =list()
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
303
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
304 for (pixelnb in 1:length(peaks))
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
305 {
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
306 ## 2) calculate warping:
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
307 warping_function <- determineWarpingFunctions(peaks[[pixelnb]],
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
308 tolerance=$method.methods_conditional.tolerance, method="$method.methods_conditional.warping_method",
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
309 allowNoMatches=$method.methods_conditional.allow_nomatch, reference = mass_list)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
310
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
311 ## 3) warp spectra:
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
312 maldi_data_new = warpMassSpectra(list(maldi_data[[pixelnb]]), warping_function, emptyNoMatches=$method.methods_conditional.empty_nomatch)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
313 maldi_data_new_list = c(maldi_data_new_list, maldi_data_new)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
314
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
315 }
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
316 maldi_data = maldi_data_new_list
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
317 #end if
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
318 #end if
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
319
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
320
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
321 #if $method.methods_conditional.remove_empty:
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
322 print(paste(length(findEmptyMassObjects(maldi_data)), " empty spectra were removed", sep=" "))
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
323
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
324 ## only if there are empty spectra to remove
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
325
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
326 if (length(findEmptyMassObjects(maldi_data))>0)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
327
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
328 {
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
329 maldi_data = removeEmptyMassObjects(maldi_data)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
330 }
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
331 #end if
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
332
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
333 ## QC plot
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
334
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
335 if (length(maldi_data)>0){
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
336 avgSpectra = averageMassSpectra(maldi_data,method="mean")
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
337 plot(avgSpectra, main="Average spectrum after alignment")
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
338 medint = round(median(unlist(lapply(maldi_data,intensity))), digits=2)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
339 }else{print("All spectra are empty")}
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
340
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
341 pixel_number = length(maldi_data)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
342 minmz = round(min(unlist(lapply(maldi_data,mass))), digits=4)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
343 maxmz = round(max(unlist(lapply(maldi_data,mass))), digits=4)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
344 mean_features = round(length(unlist(lapply(maldi_data,mass)))/length(maldi_data), digits=2)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
345
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
346 number_features = length(unique(unlist(lapply(maldi_data,mass))))
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
347 spectra_aligned = c(minmz, maxmz,number_features,mean_features, medint,pixel_number)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
348 QC_numbers= cbind(QC_numbers, spectra_aligned)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
349 vectorofactions = append(vectorofactions, "aligned")
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
350
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
351 #elif str( $method.methods_conditional.method ) == 'skip_preprocessing':
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
352 ##for now as option to filter large files
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
353
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
354 #end if
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
355
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
356 #end for
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
357
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
358 rownames(QC_numbers) = c("min m/z", "max mz", "# features", "median \n# peaks (int.>0)", "median\nintensity", "pixel\nnumber")
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
359 colnames(QC_numbers) = vectorofactions
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
360 plot(0,type='n',axes=FALSE,ann=FALSE)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
361 grid.table(t(QC_numbers))
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
362
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
363 dev.off()
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
364
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
365 ## export imzML file
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
366 if (length(maldi_data)>0){
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
367 MALDIquantForeign::exportImzMl(maldi_data, file="out.imzMl", processed=$export_processed)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
368
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
369 }else{"All spectra are empty, outputfiles will be empty,too."}
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
370
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
371 ]]>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
372 </configfile>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
373 </configfiles>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
374 <inputs>
4
82b1de5c142b planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f1e1cd260ef2884d0ba12e2b614df3c72d0934dc
galaxyp
parents: 3
diff changeset
375 <param name="infile" type="data" format="imzml,analyze75" label="MSI data" help="Input file as imzML (composite upload). The file must be in profile mode, not centroided."/>
3
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
376 <conditional name="restriction_conditional">
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
377 <param name="restriction" type="select" label="Use only spectra of interest" help="This option only works for imzML files">
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
378 <option value="no_restriction" selected="True">No, calculate on entire file</option>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
379 <option value="restrict">Yes, restrict to spectra of interest</option>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
380 </param>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
381 <when value="restrict">
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
382 <param name="coordinates_file" type="data" format="tabular" label="Tabular file with coordinates"/>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
383 <param name="column_x" data_ref="coordinates_file" label="Column with x values" type="data_column"/>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
384 <param name="column_y" data_ref="coordinates_file" label="Column with y values" type="data_column"/>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
385 <param name="coordinates_header" type="boolean" label="File contains a header line" truevalue="TRUE" falsevalue="FALSE"/>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
386 </when>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
387 <when value="no_restriction"/>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
388 </conditional>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
389 <repeat name="methods" title="Method" min="1">
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
390 <conditional name="methods_conditional">
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
391 <param name="method" type="select" label="Select a method">
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
392 <option value="Transformation" selected="True">Transformation</option>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
393 <option value="Smoothing">Smoothing</option>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
394 <option value="Baseline">Baseline removal</option>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
395 <option value="Calibrate">Intensity calibration (normalization)</option>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
396 <option value="Align">Align spectra (warping/phase correction)</option>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
397 <option value="skip_preprocessing">Skip preprocessing</option>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
398 <validator type="empty_field" />
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
399 </param>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
400 <when value="Transformation">
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
401 <param name="transform_method" type="select" label="Transformation method">
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
402 <option value="sqrt" selected="True">sqrt</option>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
403 <option value="log">log</option>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
404 <option value="log2">log2</option>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
405 <option value="log10">log10</option>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
406 <validator type="empty_field" />
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
407 </param>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
408 </when>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
409 <when value="Smoothing">
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
410 <conditional name="methods_for_smoothing">
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
411 <param name="smooth_method" type="select" label="Smoothing method" help="This method smoothes the intensity values of a MassSpectrum object.">
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
412 <option value="SavitzkyGolay" selected="True">SavitzkyGolay</option>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
413 <option value="MovingAverage">MovingAverage</option>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
414 </param>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
415 <when value="SavitzkyGolay">
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
416 <param name="polynomial" value="3" type="text" label="Polynomial order"
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
417 help="Controls the order of the filter, should be smaller than the resulting window."/>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
418 </when>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
419 <when value="MovingAverage">
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
420 <param name="weighted" type="boolean" label="Weighted average" help = "Indicates if the average should be equal weight or if it should have weights depending on the distance from the center as calculated as 1/2^abs(-halfWindowSize:halfWindowSize) with the sum of all weigths normalized to 1" truevalue="TRUE" falsevalue="FALSE"/>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
421 </when>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
422 </conditional>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
423 <param name="halfWindowSize" type="integer" value="10"
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
424 label="Half window size"
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
425 help="Number of data points, the resulting window reaches from
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
426 mass[currentIndex-halfWindowSize] to mass[currentIndex+halfWindowSize]
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
427 (window size is 2*halfWindowSize+1)."/>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
428 </when>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
429 <when value="Baseline">
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
430 <conditional name="methods_for_baseline">
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
431 <param name="baseline_method" type="select" label="Baseline removal method">
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
432 <option value="SNIP" selected="True">SNIP</option>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
433 <option value="TopHat">TopHat</option>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
434 <option value="ConvexHull">ConvexHull</option>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
435 <option value="median">median</option>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
436 <validator type="empty_field" />
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
437 </param>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
438 <when value="SNIP">
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
439 <param name="iterations" type="integer" value="100"
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
440 label="Number of iterations" help="Corresponds to half window size: The resulting window reaches from mass[cur_index-iterations] to mass[cur_index+iterations]"/>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
441 </when>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
442 <when value="TopHat">
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
443 <param name="tophat_halfWindowSize" type="integer" value="10"
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
444 label="Half window size"
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
445 help="Number of data points, the resulting window reaches from
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
446 mass[currentIndex-halfWindowSize] to mass[currentIndex+halfWindowSize]"/>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
447 </when>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
448 <when value="ConvexHull"/>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
449 <when value="median">
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
450 <param name="median_halfWindowSize" type="integer" value="10"
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
451 label="Half window size"
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
452 help="Number of data points, the resulting window reaches from
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
453 mass[currentIndex-halfWindowSize] to mass[currentIndex+halfWindowSize]"/>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
454 </when>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
455 </conditional>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
456 </when>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
457 <when value="Calibrate">
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
458 <param name="calibrate_method" type="select" label="Intensity calibration method" help="Intensity normalization">
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
459 <option value="TIC" selected="True">TIC</option>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
460 <option value="PQN">PQN</option>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
461 <option value="median">median</option>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
462 <validator type="empty_field" />
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
463 </param>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
464 <conditional name="cond_calibration_range">
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
465 <param name="calibration_range" type="select" label="m/z range" help="Instead of the whole m/z range, a specified m/z range can be used to calculate the scaling factor.">
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
466 <option value="no" selected="True">complete m/z range</option>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
467 <option value="yes">specify a m/z range</option>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
468 </param>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
469 <when value="no"/>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
470 <when value="yes">
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
471 <param name="mass_start" type="integer" value="800"
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
472 label="Start of m/z range, has to be inside m/z range"
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
473 help="Scaling factor is calculated on the mass range and applied to the whole spectrum."/>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
474 <param name="mass_end" type="integer" value="3000"
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
475 label="End of m/z range, has to be inside m/z range"/>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
476 </when>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
477 </conditional>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
478 </when>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
479 <when value="Align">
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
480 <param name="warping_method" type="select" label="Alignment method">
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
481 <option value="lowess" selected="True">Lowess</option>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
482 <option value="linear">Linear</option>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
483 <option value="quadratic">Quadratic</option>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
484 <option value="cubic">Cubic</option>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
485 </param>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
486
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
487 <param name="tolerance" type="float" value="0.00005"
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
488 label="Tolerance"
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
489 help="abs(mz1 - mz2)/mz2, maximal relative deviation of a peak position (m/z) to be considered as identical. For 50ppm use 0.00005 or 5e-5" />
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
490
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
491 <param name="halfWindowSize" type="integer" value="20"
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
492 label="Half window size"
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
493 help="Number of data points, the resulting window reaches from
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
494 mass[currentIndex-halfWindowSize] to mass[currentIndex+halfWindowSize]
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
495 (window size is 2*halfWindowSize+1).
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
496 The best size differs depending on the selected smoothing method."/>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
497
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
498 <param name="peak_method" type="select" label="Noise estimation function">
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
499 <option value="MAD" selected="True">MAD</option>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
500 <option value="SuperSmoother">SuperSmoother</option>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
501 </param>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
502
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
503 <param name="snr" type="integer" value="2" label="Signal-to-noise-ratio"/>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
504 <param name="allow_nomatch" type="boolean" label="Allow no matches" help="Don't throw an error when less than 2 reference m/z were found in a spectrum" truevalue="TRUE" falsevalue="FALSE"/>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
505 <param name="empty_nomatch" type="boolean" label="Empty no matches" help="If TRUE the intensity values of MassSpectrum or MassPeaks objects with missing (NA) warping functions are set to zero" truevalue="TRUE" falsevalue="FALSE"/>
4
82b1de5c142b planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f1e1cd260ef2884d0ba12e2b614df3c72d0934dc
galaxyp
parents: 3
diff changeset
506 <param name="remove_empty" type="boolean" label="Remove empty spectra" truevalue="TRUE" falsevalue="FALSE"/>
3
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
507
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
508 <conditional name="reference_for_alignment">
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
509 <param name="align_ref" type="select" label="Reference" help="If given, samples will be aligned to reference, use internal calibrants to perform m/z calibration">
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
510 <option value="no_reference" selected="True">no reference</option>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
511 <option value="yes_reference">reference from tabular file</option>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
512 </param>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
513 <when value="no_reference">
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
514 <param name="min_frequency" type="float" value="0.9" label = "minFrequency" help="Removal of all peaks which occur in less than minFrequency spectra to generate the reference m/z"/>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
515 </when>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
516 <when value="yes_reference">
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
517 <param name="reference_file" type="data" format="tabular"
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
518 label="Reference m/z values"
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
519 help="Tabular file"/>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
520 <param name="mz_column" data_ref="reference_file" label="Column with m/z values" type="data_column"/>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
521 <param name="reference_header" type="boolean" label="File contains a header line" truevalue="TRUE" falsevalue="FALSE"/>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
522 <param name="separate_alignment" type="boolean" label="Spectrum wise alignment" help="Internal binning is omitted to avoid interaction between spectra" truevalue="TRUE" falsevalue="FALSE"/>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
523 </when>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
524 </conditional>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
525 </when>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
526 <when value="skip_preprocessing"/>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
527 </conditional>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
528 </repeat>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
529 <param name="export_processed" type="boolean" label="Export processed imzML" help="otherwise continuous imzML will be exported" truevalue="TRUE" falsevalue="FALSE"/>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
530 </inputs>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
531 <outputs>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
532 <data format="imzml" name="outfile_imzml" label="${tool.name} on ${on_string}" />
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
533 <data format="pdf" name="plots" from_work_dir="prepro_qc_plot.pdf" label="${tool.name} on ${on_string}: QC"/>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
534 </outputs>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
535 <tests>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
536 <test>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
537 <param name="infile" value="" ftype="imzml">
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
538 <composite_data value="Example_Continuous.imzML"/>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
539 <composite_data value="Example_Continuous.ibd"/>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
540 </param>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
541 <conditional name="restriction_conditional">
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
542 <param name="restriction" value="restrict"/>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
543 <param name="coordinates_file" value="restricted_pixels.tabular"/>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
544 <param name="column_x" value="1"/>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
545 <param name="column_y" value="2"/>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
546 </conditional>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
547 <conditional name="methods_conditional">
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
548 <param name="method" value="Transformation"/>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
549 <param name="transform_method" value="log2"/>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
550 <param name="method" value="Smoothing"/>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
551 <param name="smooth_method" value="SavitzkyGolay"/>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
552 <param name="method" value="Basline"/>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
553 <param name="baseline_method" value ="TopHat"/>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
554 </conditional>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
555 <output name="outfile_imzml" ftype="imzml" file="preprocessing1.imzml.txt" lines_diff="4">
4
82b1de5c142b planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f1e1cd260ef2884d0ba12e2b614df3c72d0934dc
galaxyp
parents: 3
diff changeset
556 <extra_files type="file" file="outfile1.imzml" name="imzml" lines_diff="8"/>
3
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
557 <extra_files type="file" file="outfile1.ibd" name="ibd" compare="sim_size"/>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
558 </output>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
559 <output name="plots" file="Preprocessing1_QC.pdf" compare="sim_size"/>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
560 </test>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
561 <test>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
562 <param name="infile" value="" ftype="imzml">
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
563 <composite_data value="Example_Continuous.imzML"/>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
564 <composite_data value="Example_Continuous.ibd"/>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
565 </param>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
566 <conditional name="methods_conditional">
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
567 <param name="method" value="Calibrate"/>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
568 <param name="calibrate_method" value="median"/>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
569 </conditional>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
570 <output name="outfile_imzml" ftype="imzml" file="preprocessing3.imzml.txt" lines_diff="4">
4
82b1de5c142b planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f1e1cd260ef2884d0ba12e2b614df3c72d0934dc
galaxyp
parents: 3
diff changeset
571 <extra_files type="file" file="outfile3.imzml" name="imzml" lines_diff="8"/>
3
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
572 <extra_files type="file" file="outfile3.ibd" name="ibd" compare="sim_size"/>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
573 </output>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
574 <output name="plots" file="Preprocessing3_QC.pdf" compare="sim_size"/>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
575 </test>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
576 </tests>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
577 <help><![CDATA[
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
578
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
579 @MADLI_QUANT_DESCRIPTION@
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
580
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
581 -----
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
582
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
583 **Input data**
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
584
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
585 - MSI data: 2 types of input data can be used:
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
586
4
82b1de5c142b planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f1e1cd260ef2884d0ba12e2b614df3c72d0934dc
galaxyp
parents: 3
diff changeset
587 - imzml file (upload imzml and ibd file via the "composite" function) `Introduction to the imzml format <https://ms-imaging.org/imzml/>`_
3
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
588
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
589 - Optional: Tabular file with pixel coordinates to restrict reading of imzML files to coordinates of interest. Tabular files with any header name or no header at all are supported.
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
590
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
591 ::
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
592
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
593 x_coord y_coord
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
594 1 1
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
595 2 1
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
596 3 1
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
597 ...
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
598 ...
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
599
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
600 - Optional: Tabular file with reference m/z for the spectra align function. At least 2 m/z values of the input list must be present in every spectrum to peform the alignment. First column must contain m/z values, without empty fields or letters. Tabular files with any header name or no header at all are supported.
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
601
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
602 ::
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
603
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
604
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
605 m/z
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
606 100.0
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
607 100.01
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
608 100.02
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
609 ...
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
610 ...
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
611
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
612
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
613
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
614 **Options**
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
615
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
616 - Transformation: Variance stabilization through intensity transformation:'log', 'log2', 'log10' and 'squareroot' (sqrt) are available
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
617 - Smoothing: Smoothing of the peaks reduces noise and improves peak detection. Available smoothing methods are 'SavitzkyGolay' and 'Moving Average'
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
618
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
619 - For all smoothing methods: The larger the 'Half window size', the stronger the smoothing. The resulting window should be smaller than the FWHM (full width at half maximum) of the typical peaks. Moving average needs smaller window size than SavitzkyGolay.
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
620 - Moving average: Recommended for broader peaks/high m/z range spectra. Weighted moving average: Points in the center get larger weight factors than points away from the center.
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
621 - SavitzkyGolay: Recommended for sharp peaks/low m/z range, preserves the shape of the local maxima. The PolynomialOrder should be smaller than the resulting window. Negative values will be replaced with 0.
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
622
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
623 - Baseline reduction: Baseline reduction removes background intensity generated by chemical noise (common in MALDI datasets).
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
624
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
625 - Available methods are SNIP, TopHat,ConvexHull and median:
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
626 - SNIP is the default baseline reduction method in MALDIquant.
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
627 - ConvexHull is not appropriate for MALDI-TOF baseline removal.
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
628 - The moving median may generate negative intensities.
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
629 - Except for the ConvexHull all methods have a parameter for the 'Half window size' (in SNIP it is called 'iterations'). The smaller the window the more baseline will be removed but also parts of the peaks. Wider windows preserve the peak height better and produce a smoother baseline, but some local background variation will remain.
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
630
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
631 - Intensity calibration (normalization): Normalization of intensities to Total Ion Current (TIC), median spectrum, Probabilistic Quotient Normalization (PQN)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
632
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
633 - TIC and median are local calibration methods: each spectrum is normalized on its own (each peak is divided by the TIC or median of the spectrum)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
634 - PQN is a global calibration method: In PQN all spectra are calibrated using the TIC calibration first. Subsequently, a median reference spectrum is created and the intensities in all spectra are standardized using the reference spectrum and a spectrum-specific median is calculated for each spectrum. Finally, each spectrum is rescaled by the median of the ratios of its intensity values and that of the reference spectrum
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
635
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
636 - Spectra alignment (warping): alignment for (re)calibration of m/z values.
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
637
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
638 - peak detection is performed, the reference peaks will be matched to those detected peaks
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
639 - without external reference m/z: internal reference is obtained by filtering and binning the picked peaks to find landmark peaks and their average m/z
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
640 - with external reference m/z: the given m/z are used as a reference, at least 10 reference values are recommended
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
641 - non linear warping (parametric time warping plus binning) to match the reference peaks (internal or external) to the present (picked) peaks with the given tolerance. At least two m/z per spectrum are needed for the alignment. To prevent an error when this criterium is not fullfilled, "Don't throw an error when less than 2 reference m/z were found in a spectrum" should be set to yes. If the not aligned spectra should be set to zero select yes in "If TRUE the intensity values of MassSpectrum or MassPeaks objects with missing (NA) warping functions are set to zero". In order to remove such empty spectra set "Should empty spectra be removed" to yes.
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
642
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
643
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
644 **Output**
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
645
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
646 - imzML file (imzML format can be continuous or processed)
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
647 - PDF with average mass spectra after each preprocessing step
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
648
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
649 .. _MALDIquant: http://strimmerlab.org/software/maldiquant/
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
650
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
651 ]]>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
652 </help>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
653 <expand macro="citation"/>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
654 </tool>
aa81463e21ea planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
655