Mercurial > repos > galaxyp > maldi_quant_peak_detection

annotate maldi_quant_peakdetection.xml @ 4:e9300ef37403 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit ecdc3a64aa245d80dbc5487b2bf10a85a43adc6d
author galaxyp
date Fri, 22 Mar 2019 08:26:23 -0400
parents 36d38d2cf88c
children e66f552a3c47
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
4
e9300ef37403 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit ecdc3a64aa245d80dbc5487b2bf10a85a43adc6d
galaxyp
parents: 3
diff changeset
1 <tool id="maldi_quant_peak_detection" name="MALDIquant peak detection" version="@VERSION@.4">
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
2 <description>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
3 Peak detection, binning and filtering for mass-spectrometry imaging data
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
4 </description>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
5 <macros>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
6 <import>maldi_macros.xml</import>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
7 </macros>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
8 <expand macro="requirements"/>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
9 <command detect_errors="exit_code">
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
10 <![CDATA[
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
11 #if $infile.ext == 'imzml'
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
12 cp '${infile.extra_files_path}/imzml' infile.imzML &&
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
13 cp '${infile.extra_files_path}/ibd' infile.ibd &&
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
14 #elif $infile.ext == 'analyze75'
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
15 cp '${infile.extra_files_path}/hdr' infile.hdr &&
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
16 cp '${infile.extra_files_path}/img' infile.img &&
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
17 cp '${infile.extra_files_path}/t2m' infile.t2m &&
1
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
18 #else
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
19 ln -s $infile infile.RData &&
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
20 #end if
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
21 Rscript '${maldi_quant_peak_detection}'&&
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
22 mkdir $outfile_imzml.files_path &&
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
23 mv ./out.imzMl "${os.path.join($outfile_imzml.files_path, 'imzml')}" | true &&
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
24 mv ./out.ibd "${os.path.join($outfile_imzml.files_path, 'ibd')}" | true &&
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
25 echo "imzML file:" > $outfile_imzml &&
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
26 ls -l "$outfile_imzml.files_path" >> $outfile_imzml
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
27 ]]>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
28 </command>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
29 <configfiles>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
30 <configfile name="maldi_quant_peak_detection"><![CDATA[
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
31
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
32 @R_IMPORTS@
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
33
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
34 #if $restriction_conditional.restriction == 'restrict':
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
35
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
36 print('Reading mask region')
1
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
37
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
38 ## Import imzML file
1
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
39 coordinate_matrix = as.matrix(read.delim("$restriction_conditional.coordinates_file", header = $restriction_conditional.coordinates_header, stringsAsFactors = FALSE))[,1:2]
3
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
40 coordinate_matrix = coordinate_matrix[,c($restriction_conditional.column_x, $restriction_conditional.column_y)]
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
41
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
42 maldi_data <- importImzMl('infile.imzML',
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
43 coordinates = coordinate_matrix, centroided = $centroids)
1
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
44 pixelnames = paste("xy", coordinates(maldi_data)[,1],coordinates(maldi_data)[,2], sep="_")
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
45
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
46
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
47 #else:
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
48
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
49 print('Reading entire file')
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
50 ## Import imzML file
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
51
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
52 #if $infile.ext == 'imzml'
1
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
53 print('imzML file')
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
54 #if str($centroids) == "TRUE"
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
55 peaks <- importImzMl('infile.imzML', centroided = $centroids)
2
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
56 pixelnames = paste("xy", coordinates(peaks)[,1],coordinates(peaks)[,2], sep="_")
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
57 coordinates_info = cbind(coordinates(peaks)[,1:2], c(1:length(peaks)))
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
58 #else
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
59 maldi_data <- importImzMl('infile.imzML', centroided = $centroids)
1
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
60 pixelnames = paste("xy", coordinates(maldi_data)[,1],coordinates(maldi_data)[,2], sep="_")
2
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
61 coordinates_info = cbind(coordinates(maldi_data)[,1:2], c(1:length(maldi_data)))
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
62 #end if
2
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
63
1
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
64
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
65 #elif $infile.ext == 'tabular'
1
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
66 print('tabular file')
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
67 #set $centroids = "TRUE" ## will be used in some if conditions
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
68 peak_tabular = read.delim("$infile", header = TRUE, stringsAsFactors = FALSE)
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
69 peak_list = split(peak_tabular, f = peak_tabular\$spectrum) ## will be ordered according to spectrum
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
70 pixelnames = unique(peak_tabular\$spectrum)
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
71
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
72 peaks = list()
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
73 for (spectra in 1:length(peak_list))
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
74 {
1
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
75 single_peaks = createMassPeaks(peak_list[[spectra]]\$mass, peak_list[[spectra]]\$intensity, snr=peak_list[[spectra]]\$snr)
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
76 peaks[[spectra]] = single_peaks
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
77 }
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
78
1
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
79 #else
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
80 print('rdata file')
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
81 loadRData <- function(fileName){
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
82 #loads an RData file, and returns it
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
83 load(fileName)
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
84 get(ls()[ls() != "fileName"])
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
85 }
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
86 msidata = loadRData('infile.RData')
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
87 centroided(msidata) = $centroids
2
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
88 ## change to correct pixelnames
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
89
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
90 x_coords = unlist(lapply(strsplit(names(Cardinal::pixels(msidata)), ","), `[[`, 1))
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
91 y_coords = unlist(lapply(strsplit(names(Cardinal::pixels(msidata)), ","), `[[`, 2))
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
92 x_coordinates = gsub("x = ","",x_coords)
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
93 y_coordinates = gsub(" y = ","",y_coords)
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
94 pixelnames = paste0("xy_", x_coordinates, "_", y_coordinates)
1
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
95
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
96 cardinal_coordinates = as.matrix(Cardinal::coord(msidata)[,1:2])
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
97
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
98 if (centroided(msidata) == FALSE){
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
99 ## create mass spectrum object
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
100 cardinal_mzs = Cardinal::mz(msidata)
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
101 maldi_data = list()
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
102 for(number_spectra in 1:ncol(msidata)){
3
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
103 maldi_data[[number_spectra]] = createMassSpectrum(mass = cardinal_mzs, intensity = iData(msidata)[,number_spectra])}
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
104 coordinates_info = cbind(cardinal_coordinates, c(1:length(maldi_data)))
1
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
105 }else{
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
106 peaks = list()
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
107 for (spectra in 1:ncol(msidata))
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
108 {
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
109 single_peaks = createMassPeaks(Cardinal::mz(msidata), Cardinal::spectra(msidata)[,spectra], snr=as.numeric(rep("NA", nrow(msidata))))
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
110 peaks[[spectra]] = single_peaks
3
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
111 }
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
112 coordinates_info = cbind(cardinal_coordinates, c(1:length(peaks)))}
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
113 #end if
1
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
114 #end if
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
115
1
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
116 ## default summarized = FALSE
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
117 summarized_spectra = FALSE
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
118
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
119 ## Quality control plots during peak detection
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
120 pdf("peaks_qc_plot.pdf", fonts = "Times", pointsize = 12)
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
121 plot(0,type='n',axes=FALSE,ann=FALSE)
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
122
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
123 ## if no filename is given, name of file in Galaxy history is used
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
124 #set $filename = $infile.display_name
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
125
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
126 title(main=paste("$filename"))
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
127
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
128 ## plot input file spectrum:
2
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
129 #if $centroids:
3
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
130 ## Choose random spectra for QC plots
4
e9300ef37403 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit ecdc3a64aa245d80dbc5487b2bf10a85a43adc6d
galaxyp
parents: 3
diff changeset
131 print(length(peaks))
3
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
132 random_spectra = sample(1:length(peaks), 4, replace=FALSE)
4
e9300ef37403 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit ecdc3a64aa245d80dbc5487b2bf10a85a43adc6d
galaxyp
parents: 3
diff changeset
133 random_spectra_name = pixelnames[random_spectra]
3
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
134 par(mfrow = c(2, 2), oma=c(0,0,2,0))
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
135 for (random_sample in random_spectra){
4
e9300ef37403 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit ecdc3a64aa245d80dbc5487b2bf10a85a43adc6d
galaxyp
parents: 3
diff changeset
136 plot(peaks[[random_sample]],sub="", main=paste0("spectrum ", pixelnames[random_sample]))}
3
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
137 title("Input spectra", outer=TRUE, line=0)
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
138
1
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
139 #else
3
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
140 ## Choose random spectra for QC plots
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
141 random_spectra = sample(1:length(maldi_data), 4, replace=FALSE)
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
142 par(mfrow = c(2, 2), oma=c(0,0,2,0))
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
143 for (random_sample in random_spectra){
4
e9300ef37403 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit ecdc3a64aa245d80dbc5487b2bf10a85a43adc6d
galaxyp
parents: 3
diff changeset
144 plot(maldi_data[[random_sample]],sub="", main=paste0("spectrum ", pixelnames[random_sample]))
e9300ef37403 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit ecdc3a64aa245d80dbc5487b2bf10a85a43adc6d
galaxyp
parents: 3
diff changeset
145 }
3
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
146 title("Input spectra", outer=TRUE, line=0)
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
147 #end if
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
148
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
149
1
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
150 ## QC numbers for input file
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
151 #if str($centroids) == "TRUE"
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
152 pixel_number = length(peaks)
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
153 minmz = round(min(unlist(lapply(peaks,mass))), digits=4)
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
154 maxmz = round(max(unlist(lapply(peaks,mass))), digits=4)
2
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
155 mean_features = round(length(unlist(lapply(peaks,mass)))/length(peaks), digits=2)
1
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
156 medint = round(median(unlist(lapply(peaks,intensity))), digits=2)
2
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
157 number_features = length(unique(unlist(lapply(peaks,mass))))
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
158 inputdata = c(minmz, maxmz,number_features,mean_features, medint)
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
159 QC_numbers= data.frame(inputdata = c(minmz, maxmz,number_features, mean_features, medint))
1
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
160 vectorofactions = "inputdata"
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
161 #else
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
162 pixel_number = length(maldi_data)
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
163 minmz = round(min(unlist(lapply(maldi_data,mass))), digits=4)
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
164 maxmz = round(max(unlist(lapply(maldi_data,mass))), digits=4)
2
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
165 mean_features = round(length(unlist(lapply(maldi_data,mass)))/length(maldi_data), digits=2)
1
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
166 medint = round(median(unlist(lapply(maldi_data,intensity))), digits=2)
2
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
167 number_features = length(unique(unlist(lapply(maldi_data,mass))))
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
168 inputdata = c(minmz, maxmz,number_features,mean_features, medint)
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
169 QC_numbers= data.frame(inputdata = c(minmz, maxmz,number_features, mean_features, medint))
1
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
170 vectorofactions = "inputdata"
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
171 #end if
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
172
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
173 #if str($tabular_annotation.load_annotation) == 'yes_annotation':
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
174
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
175 ## read and extract x,y,annotation information
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
176 input_tabular = read.delim("$tabular_annotation.annotation_file", header = $tabular_annotation.tabular_header, stringsAsFactors = FALSE)
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
177 annotation_input = input_tabular[,c($tabular_annotation.column_x, $tabular_annotation.column_y, $tabular_annotation.column_names)]
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
178 colnames(annotation_input) = c("x", "y", "annotation") ## rename annotations header to default name "annotation"
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
179
3
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
180 ## merge provided annotation with coordinate information of MSI data
1
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
181 colnames(coordinates_info)[3] = "pixel_index"
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
182 merged_annotation = merge(coordinates_info, annotation_input, by=c("x", "y"), all.x=TRUE)
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
183 merged_annotation[is.na(merged_annotation)] = "NA"
3
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
184 ## order coordinate information according to pixel index to make sure that the order stays the same
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
185 merged_annotation = merged_annotation[order(merged_annotation\$pixel_index),]
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
186 samples = as.factor(merged_annotation\$annotation)
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
187
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
188 ## print annotation overview into PDF output
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
189
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
190 combine_plot = ggplot(merged_annotation, aes(x=x, y=y, fill=annotation))+
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
191 geom_tile() +
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
192 coord_fixed()+
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
193 ggtitle("Spatial orientation of annotated data")+
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
194 theme_bw()+
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
195 theme(plot.title = element_text(hjust = 0.5))+
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
196 theme(text=element_text(family="ArialMT", face="bold", size=12))+
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
197 theme(legend.position="bottom",legend.direction="vertical")+
2
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
198 theme(legend.key.size = unit(0.2, "line"), legend.text = element_text(size = 8))+
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
199 theme(legend.position="bottom",legend.direction="vertical")+
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
200 guides(fill=guide_legend(ncol=4,byrow=TRUE))
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
201
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
202 print(combine_plot)
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
203
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
204 #end if
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
205
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
206 #################### Preprocessing methods #####################################
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
207
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
208 #for $method in $methods:
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
209
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
210 #if str( $method.methods_conditional.method ) == 'Peak_detection':
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
211 print('peak detection')
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
212 ##peak detection
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
213
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
214 #if $method.methods_conditional.use_annotations:
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
215 maldi_data <- averageMassSpectra(maldi_data, labels=samples,method="mean") ## use average spectra for peak picking
1
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
216 pixelnames = levels(samples)
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
217 summarized_spectra = TRUE
3
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
218 random_spectra = sample(1:length(maldi_data), 4, replace=TRUE)
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
219 #end if
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
220
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
221 peaks <- detectPeaks(maldi_data, method="$method.methods_conditional.peak_method",
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
222 halfWindowSize=$method.methods_conditional.halfWindowSize,SNR=$method.methods_conditional.snr)
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
223
1
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
224 ## QC plot and numbers
3
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
225 ## plot old spectra with baseline in blue and picked peaks in green
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
226 par(mfrow = c(2, 2), oma=c(0,0,2,0))
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
227 for (random_sample in random_spectra){
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
228 noise = estimateNoise(maldi_data[[random_sample]], method= "$method.methods_conditional.peak_method")
4
e9300ef37403 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit ecdc3a64aa245d80dbc5487b2bf10a85a43adc6d
galaxyp
parents: 3
diff changeset
229 plot(maldi_data[[random_sample]], sub="", main=paste0("spectrum ", pixelnames[random_sample]))
3
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
230 lines(noise[,1], noise[,2]*$method.methods_conditional.snr, col="blue")
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
231 points(peaks[[random_sample]], col="green", pch=20)}
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
232 title("S/N in blue and picked peaks in green", outer=TRUE, line=0)
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
233
2
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
234 ## plot new spectrum
3
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
235 par(mfrow = c(2, 2), oma=c(0,0,2,0))
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
236 for (random_sample in random_spectra){
4
e9300ef37403 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit ecdc3a64aa245d80dbc5487b2bf10a85a43adc6d
galaxyp
parents: 3
diff changeset
237 plot(peaks[[random_sample]], sub="", main=paste0("spectrum ", pixelnames[random_sample]))}
3
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
238 title("Picked peaks", outer=TRUE, line=0)
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
239
1
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
240 pixel_number = length(peaks)
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
241 minmz = round(min(unlist(lapply(peaks,mass))), digits=4)
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
242 maxmz = round(max(unlist(lapply(peaks,mass))), digits=4)
2
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
243 mean_features = round(length(unlist(lapply(peaks,mass)))/length(peaks), digits=2)
1
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
244 medint = round(median(unlist(lapply(peaks,intensity))), digits=2)
2
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
245 number_features = length(unique(unlist(lapply(peaks,mass))))
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
246 peaks_picked = c(minmz, maxmz,number_features,mean_features, medint)
1
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
247 QC_numbers= cbind(QC_numbers, peaks_picked)
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
248 vectorofactions = append(vectorofactions, "peaks_picked")
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
249
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
250 if (length(peaks[!sapply(peaks, isEmpty)])>0){
2
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
251 featureMatrix <- intensityMatrix(peaks)
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
252 #if $infile.ext == 'imzml'
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
253 #if str($centroids) == "FALSE"
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
254 featureMatrix <- intensityMatrix(peaks, maldi_data)
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
255 #end if
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
256 #end if
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
257 featureMatrix2 =cbind(pixelnames, featureMatrix)
1
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
258 colnames(featureMatrix2)[1] = c("mz")
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
259 featureMatrix2 = t(featureMatrix2)
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
260 write.table(featureMatrix2, file="$intensity_matrix", quote = FALSE, row.names = TRUE, col.names=FALSE, sep = "\t")
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
261 }else{print("There are no spectra with peaks left")}
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
262
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
263
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
264 #elif str( $method.methods_conditional.method ) == 'monoisotopic_peaks':
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
265
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
266 print('monoisotopic peaks')
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
267 ##monoisotopic peaks
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
268
3
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
269 ## keep peaks to plot them with monoisotopic peaks
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
270 picked_peaks = peaks
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
271
3
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
272 peaks = monoisotopicPeaks(peaks, minCor=$method.methods_conditional.minCor,
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
273 tolerance=$method.methods_conditional.tolerance,
4
e9300ef37403 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit ecdc3a64aa245d80dbc5487b2bf10a85a43adc6d
galaxyp
parents: 3
diff changeset
274 distance=c($method.methods_conditional.distance),
3
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
275 size=$method.methods_conditional.size)
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
276
2
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
277 ## plot old spectrum with picked isotopes as green dots
3
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
278 par(mfrow = c(2, 2), oma=c(0,0,2,0))
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
279 for (random_sample in random_spectra){
4
e9300ef37403 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit ecdc3a64aa245d80dbc5487b2bf10a85a43adc6d
galaxyp
parents: 3
diff changeset
280 plot(picked_peaks[[random_sample]], sub="", main=paste0("spectrum ", pixelnames[random_sample]))
3
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
281 points(peaks[[random_sample]], col="green", pch=20)}
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
282 title(paste0("Monoisotopic peaks in green"), outer=TRUE, line=0)
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
283
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
284
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
285 par(mfrow = c(2, 2), oma=c(0,0,2,0))
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
286 for (random_sample in random_spectra){
4
e9300ef37403 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit ecdc3a64aa245d80dbc5487b2bf10a85a43adc6d
galaxyp
parents: 3
diff changeset
287 plot(peaks[[random_sample]], sub="", main=paste0("spectrum ", pixelnames[random_sample]))}
3
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
288 title("Monoisotopic peaks", outer=TRUE, line=0)
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
289
1
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
290 minmz = round(min(unlist(lapply(peaks,mass))), digits=4)
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
291 maxmz = round(max(unlist(lapply(peaks,mass))), digits=4)
2
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
292 mean_features = round(length(unlist(lapply(peaks,mass)))/length(peaks), digits=2)
1
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
293 medint = round(median(unlist(lapply(peaks,intensity))), digits=2)
2
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
294 number_features = length(unique(unlist(lapply(peaks,mass))))
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
295 monoisotopes = c(minmz, maxmz,number_features,mean_features, medint)
1
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
296 QC_numbers= cbind(QC_numbers, monoisotopes)
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
297 vectorofactions = append(vectorofactions, "monoisotopes")
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
298
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
299 if (length(peaks[!sapply(peaks, isEmpty)])>0){
2
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
300 featureMatrix <- intensityMatrix(peaks)
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
301 ## only for profile imzML file: featurematrix is overwritten:
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
302 #if $infile.ext == 'imzml'
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
303 #if str($centroids) == "FALSE"
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
304 featureMatrix <- intensityMatrix(peaks, maldi_data)
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
305 #end if
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
306 #end if
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
307 featureMatrix2 =cbind(pixelnames, featureMatrix)
1
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
308 colnames(featureMatrix2)[1] = c("mz")
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
309 featureMatrix2 = t(featureMatrix2)
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
310 write.table(featureMatrix2, file="$intensity_matrix", quote = FALSE, row.names = TRUE, col.names=FALSE, sep = "\t")
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
311 }else{print("There are no spectra with peaks left")}
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
312
3
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
313
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
314
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
315 #elif str( $method.methods_conditional.method ) == 'Align':
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
316
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
317 print('align')
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
318 ##align spectra with 2 separate functions
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
319
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
320
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
321 #if str($method.methods_conditional.reference_for_alignment.align_ref) == 'no_reference':
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
322 ## 1) calculate warping:
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
323 warping_function <- determineWarpingFunctions(peaks,
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
324 tolerance=$method.methods_conditional.tolerance, method="$method.methods_conditional.warping_method",
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
325 allowNoMatches=$method.methods_conditional.allow_nomatch, minFrequency = $method.methods_conditional.reference_for_alignment.min_frequency)
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
326 ## 2) warp spectra:
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
327 peaks = warpMassPeaks(peaks, warping_function, emptyNoMatches=$method.methods_conditional.empty_nomatch)
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
328
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
329
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
330 #elif str($method.methods_conditional.reference_for_alignment.align_ref) == 'yes_reference':
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
331
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
332 ## create reference mass_vector from tabular file
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
333 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]
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
334 int_vector = rep(1,length(mass_vector))
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
335 mass_list = createMassPeaks(mass_vector, int_vector)
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
336
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
337 #if str($method.methods_conditional.reference_for_alignment.separate_alignment) == "FALSE"
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
338 print('default alignment')
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
339
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
340 ## 1) calculate warping:
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
341 warping_function <- determineWarpingFunctions(peaks,
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
342 tolerance=$method.methods_conditional.tolerance, method="$method.methods_conditional.warping_method",
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
343 allowNoMatches=$method.methods_conditional.allow_nomatch, reference = mass_list)
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
344 ## 2) warp spectra:
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
345 peaks = warpMassPeaks(peaks, warping_function, emptyNoMatches=$method.methods_conditional.empty_nomatch)
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
346
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
347 #elif str($method.methods_conditional.reference_for_alignment.separate_alignment) == "TRUE"
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
348 print('spectra wise alignment')
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
349
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
350 peaks_new_list =list()
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
351
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
352 for (pixelnb in 1:length(peaks))
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
353 {
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
354
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
355 ## 1) calculate warping:
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
356 warping_function <- determineWarpingFunctions(peaks[[pixelnb]],
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
357 tolerance=$method.methods_conditional.tolerance, method="$method.methods_conditional.warping_method",
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
358 allowNoMatches=$method.methods_conditional.allow_nomatch, reference = mass_list)
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
359 ## 2) warp spectra:
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
360 peaks_new_list = warpMassPeaks(list(peaks[[pixelnb]]), warping_function, emptyNoMatches=$method.methods_conditional.empty_nomatch)
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
361 }
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
362 peaks_new = peaks_new_list
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
363
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
364 #end if
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
365
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
366 #end if
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
367
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
368 ## QC plot and numbers
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
369 par(mfrow = c(2, 2), oma=c(0,0,2,0))
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
370 for (random_sample in random_spectra){
4
e9300ef37403 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit ecdc3a64aa245d80dbc5487b2bf10a85a43adc6d
galaxyp
parents: 3
diff changeset
371 plot(peaks[[random_sample]], sub="", main=paste0("spectrum ", pixelnames[random_sample]))}
3
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
372 title("Aligned spectra", outer=TRUE, line=0)
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
373 minmz = round(min(unlist(lapply(peaks,mass))), digits=4)
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
374 maxmz = round(max(unlist(lapply(peaks,mass))), digits=4)
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
375 mean_features = round(length(unlist(lapply(peaks,mass)))/length(peaks), digits=2)
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
376 medint = round(median(unlist(lapply(peaks,intensity))), digits=2)
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
377 number_features = length(unique(unlist(lapply(peaks,mass))))
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
378 aligned = c(minmz, maxmz,number_features,mean_features, medint)
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
379 QC_numbers= cbind(QC_numbers, aligned)
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
380 vectorofactions = append(vectorofactions, "aligned")
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
381
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
382 if (length(peaks[!sapply(peaks, isEmpty)])>0){
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
383 featureMatrix <- intensityMatrix(peaks)
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
384 ## only for profile imzML file: featurematrix is overwritten:
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
385 #if $infile.ext == 'imzml'
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
386 #if str($centroids) == "FALSE"
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
387 featureMatrix <- intensityMatrix(peaks, maldi_data)
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
388 #end if
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
389 #end if
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
390 featureMatrix2 =cbind(pixelnames, featureMatrix)
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
391 colnames(featureMatrix2)[1] = c("mz")
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
392 featureMatrix2 = t(featureMatrix2)
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
393 write.table(featureMatrix2, file="$intensity_matrix", quote = FALSE, row.names = TRUE, col.names=FALSE, sep = "\t")
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
394 }else{print("There are no spectra with peaks left")}
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
395
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
396 #elif str( $method.methods_conditional.method ) == 'Binning':
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
397
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
398 print('binning')
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
399 ##m/z binning
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
400
2
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
401 peaks <- binPeaks(peaks, tolerance=$method.methods_conditional.bin_tolerance, method="$method.methods_conditional.bin_method")
1
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
402
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
403 ## QC plot and numbers
3
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
404 par(mfrow = c(2, 2), oma=c(0,0,2,0))
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
405 for (random_sample in random_spectra){
4
e9300ef37403 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit ecdc3a64aa245d80dbc5487b2bf10a85a43adc6d
galaxyp
parents: 3
diff changeset
406 plot(peaks[[random_sample]], sub="", main=paste0("spectrum ", pixelnames[random_sample]))}
3
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
407 title("Binned spectra", outer=TRUE, line=0)
1
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
408 minmz = round(min(unlist(lapply(peaks,mass))), digits=4)
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
409 maxmz = round(max(unlist(lapply(peaks,mass))), digits=4)
2
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
410 mean_features = round(length(unlist(lapply(peaks,mass)))/length(peaks), digits=2)
1
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
411 medint =round( median(unlist(lapply(peaks,intensity))), digits=2)
2
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
412 number_features = length(unique(unlist(lapply(peaks,mass))))
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
413 binned = c(minmz, maxmz,number_features,mean_features, medint)
1
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
414 QC_numbers= cbind(QC_numbers, binned)
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
415 vectorofactions = append(vectorofactions, "binned")
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
416
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
417 if (length(peaks[!sapply(peaks, isEmpty)])>0){
2
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
418 featureMatrix <- intensityMatrix(peaks)
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
419 #if $infile.ext == 'imzml'
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
420 #if str($centroids) == "FALSE"
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
421 featureMatrix <- intensityMatrix(peaks, maldi_data)
2
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
422 #elif str($centroids) == "TRUE"
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
423 featureMatrix <- intensityMatrix(peaks)
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
424 #end if
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
425 #end if
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
426 featureMatrix2 =cbind(pixelnames, featureMatrix)
1
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
427 colnames(featureMatrix2)[1] = c("mz")
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
428 featureMatrix2 = t(featureMatrix2)
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
429 write.table(featureMatrix2, file="$intensity_matrix", quote = FALSE, row.names = TRUE, col.names=FALSE, sep = "\t")
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
430 }else{print("There are no spectra with peaks left")}
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
431
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
432
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
433 #elif str( $method.methods_conditional.method ) == 'Filtering':
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
434
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
435 print('filtering')
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
436 ##m/z filtering
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
437
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
438 ## filtering on all pixels or on pixel groups:
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
439 #if str($method.methods_conditional.filter_annot_groups ) == 'FALSE':
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
440
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
441 peaks <- filterPeaks(peaks,
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
442 minFrequency=$method.methods_conditional.minFrequency,
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
443 minNumber=$method.methods_conditional.minNumber,
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
444 mergeWhitelists=$method.methods_conditional.mergeWhitelists)
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
445
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
446 #elif str( $method.methods_conditional.filter_annot_groups ) == 'TRUE':
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
447
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
448 peaks <- filterPeaks(peaks,
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
449 minFrequency=$method.methods_conditional.minFrequency,
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
450 minNumber=$method.methods_conditional.minNumber,
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
451 mergeWhitelists=$method.methods_conditional.mergeWhitelists, label = samples)
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
452 #end if
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
453
1
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
454 ##QC plot and numbers
3
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
455 par(mfrow = c(2, 2), oma=c(0,0,2,0))
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
456 for (random_sample in random_spectra){
4
e9300ef37403 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit ecdc3a64aa245d80dbc5487b2bf10a85a43adc6d
galaxyp
parents: 3
diff changeset
457 plot(peaks[[random_sample]], sub="", main=paste0("spectrum ", pixelnames[random_sample]))}
3
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
458 title("Filtered spectra", outer=TRUE, line=0)
1
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
459 minmz = round(min(unlist(lapply(peaks,mass))), digits=4)
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
460 maxmz = round(max(unlist(lapply(peaks,mass))), digits=4)
2
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
461 mean_features = round(length(unlist(lapply(peaks,mass)))/length(peaks), digits=2)
1
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
462 medint = round(median(unlist(lapply(peaks,intensity))), digits=2)
2
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
463 number_features = length(unique(unlist(lapply(peaks,mass))))
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
464 filtered = c(minmz, maxmz,number_features,mean_features, medint)
1
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
465 QC_numbers= cbind(QC_numbers, filtered)
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
466 vectorofactions = append(vectorofactions, "filtered")
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
467
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
468 if (length(peaks[!sapply(peaks, isEmpty)])>0){
2
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
469 featureMatrix <- intensityMatrix(peaks)
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
470 #if $infile.ext == 'imzml'
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
471 #if str($centroids) == "FALSE"
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
472 featureMatrix <- intensityMatrix(peaks, maldi_data)
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
473 #end if
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
474 #end if
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
475 featureMatrix2 =cbind(pixelnames, featureMatrix)
1
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
476 colnames(featureMatrix2)[1] = c("mz")
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
477 featureMatrix2 = t(featureMatrix2)
1
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
478 }else{print("There are no spectra with peaks left")
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
479 featureMatrix2 = matrix(0, ncol=1, nrow=1)}
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
480 write.table(featureMatrix2, file="$intensity_matrix", quote = FALSE, row.names = TRUE, col.names=FALSE, sep = "\t")
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
481 #end if
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
482 #end for
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
483
1
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
484 if (length(peaks[!sapply(peaks, isEmpty)])>0){
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
485 ## mass peaks output
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
486 mass_peaks = data.frame(matrix(,ncol=3, nrow=0))
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
487 for (spectrum in 1:length(peaks)){
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
488 spectrum_df = data.frame(peaks[[spectrum]]@snr, peaks[[spectrum]]@mass, peaks[[spectrum]]@intensity)
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
489 spectrum_df\$spectrum_id = rep(pixelnames[[spectrum]], length(peaks[[spectrum]]@mass))
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
490 mass_peaks = rbind(mass_peaks,spectrum_df)
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
491 }
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
492 colnames(mass_peaks) = c("snr", "mass", "intensity", "spectrum")
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
493 write.table(mass_peaks, file="$masspeaks", quote = FALSE, row.names = FALSE, col.names=TRUE, sep = "\t")
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
494 }else{print("There are no spectra with peaks left")}
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
495
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
496 ## print table with QC values
3
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
497 rownames(QC_numbers) = c("min m/z", "max mz", "# features", "median \n# peaks (int.>0)", "median\nintensity")
1
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
498 plot(0,type='n',axes=FALSE,ann=FALSE)
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
499 grid.table(t(QC_numbers))
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
500
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
501 dev.off()
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
502
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
503 if (summarized_spectra == FALSE){
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
504 #if $infile.ext == 'imzml'
3
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
505 MALDIquantForeign::exportImzMl(peaks, file="out.imzMl", processed=TRUE)
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
506 #elif $infile.ext == 'tabular'
1
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
507 masspeaks_coordinates = matrix(unlist(strsplit(as.character(pixelnames), "\\_")), ncol=3, byrow=TRUE)
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
508 ## extract x and y values and create the coordinate matrix in case tabular was input
1
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
509 peaklist_coordinates = unique(cbind(as.numeric(masspeaks_coordinates[,2]), as.numeric(masspeaks_coordinates[,3])))
3
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
510 exportImzMl(peaks, file="out.imzMl", processed=TRUE, coordinates=peaklist_coordinates)
1
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
511 #elif $infile.ext == 'rdata'
3
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
512 MALDIquantForeign::exportImzMl(peaks, file="out.imzMl", processed=TRUE, coordinates=cardinal_coordinates)
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
513 #end if
1
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
514
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
515 }
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
516
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
517 ]]>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
518 </configfile>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
519 </configfiles>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
520 <inputs>
3
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
521 <param name="infile" type="data" format="imzml,tabular,rdata" label="MSI data" help="Input file as imzML (composite upload), tabular peaklist or Cardinal MSImageSet saved as RData (regular upload)"/>
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
522 <param name="centroids" type="boolean" label="Centroided input" help="Choose Yes if peak detection has already been done. Peak detection cannot be run again on centroided data" truevalue="TRUE" falsevalue="FALSE"/>
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
523 <conditional name="restriction_conditional">
3
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
524 <param name="restriction" type="select" label="Use only spectra of interest" help="This option only works for imzML files">
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
525 <option value="no_restriction" selected="True">No, calculate on entire file</option>
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
526 <option value="restrict">Yes, restrict to spectra of interest</option>
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
527 </param>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
528 <when value="restrict">
3
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
529 <param name="coordinates_file" type="data" format="tabular" label="Tabular file with coordinates"/>
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
530 <param name="column_x" data_ref="coordinates_file" label="Column with x values" type="data_column"/>
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
531 <param name="column_y" data_ref="coordinates_file" label="Column with y values" type="data_column"/>
2
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
532 <param name="coordinates_header" type="boolean" label="File contains a header line" truevalue="TRUE" falsevalue="FALSE"/>
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
533 </when>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
534 <when value="no_restriction"/>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
535 </conditional>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
536
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
537 <conditional name="tabular_annotation">
3
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
538 <param name="load_annotation" type="select" label="Spectra annotations" help="Annotations can be used for group wise peak detection or filtering">
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
539 <option value="no_annotation" selected="True">No</option>
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
540 <option value="yes_annotation">Yes</option>
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
541 </param>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
542 <when value="yes_annotation">
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
543 <param name="annotation_file" type="data" format="tabular" label="Use annotations from tabular file"
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
544 help="Tabular file with three columns: x values, y values and pixel annotations"/>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
545 <param name="column_x" data_ref="annotation_file" label="Column with x values" type="data_column"/>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
546 <param name="column_y" data_ref="annotation_file" label="Column with y values" type="data_column"/>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
547 <param name="column_names" data_ref="annotation_file" label="Column with pixel annotations" type="data_column"/>
2
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
548 <param name="tabular_header" type="boolean" label="File contains a header line" truevalue="TRUE" falsevalue="FALSE"/>
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
549 </when>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
550 <when value="no_annotation"/>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
551 </conditional>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
552 <repeat name="methods" title="Method" min="1">
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
553 <conditional name="methods_conditional">
1
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
554 <param name="method" type="select" label="Select a method">
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
555 <option value="Peak_detection">Peak detection</option>
4
e9300ef37403 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit ecdc3a64aa245d80dbc5487b2bf10a85a43adc6d
galaxyp
parents: 3
diff changeset
556 <option value="monoisotopic_peaks">Monoisotopic peaks</option>
3
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
557 <option value="Align">Align Spectra (warping/phase correction)</option>
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
558 <option value="Binning">Binning</option>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
559 <option value="Filtering">Filtering</option>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
560 </param>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
561 <when value="Peak_detection">
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
562 <param name="peak_method" type="select" label="Noise estimation function">
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
563 <option value="MAD" selected="True">MAD</option>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
564 <option value="SuperSmoother">SuperSmoother</option>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
565 </param>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
566 <param name="halfWindowSize" type="integer" value="20"
2
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
567 label="Half window size (number of data points)"
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
568 help="The resulting window reaches from
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
569 mass[currentIndex-halfWindowSize] to mass[currentIndex+halfWindowSize]
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
570 (window size is 2*halfWindowSize+1).
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
571 The best size differs depending on the selected smoothing method."/>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
572 <param name="snr" type="integer" value="2" label="Signal-to-noise-ratio" help=""/>
2
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
573 <param name="use_annotations" type="boolean" label="Detect peaks on average mass spectra" help="Spectra with same annotation are averaged before peak detection, no imzML output" truevalue="TRUE" falsevalue="FALSE"/>
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
574 </when>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
575 <when value="monoisotopic_peaks">
2
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
576 <param name="minCor" type="float" value="0.95" label="Minimal correlation"
4
e9300ef37403 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit ecdc3a64aa245d80dbc5487b2bf10a85a43adc6d
galaxyp
parents: 3
diff changeset
577 help="Minimal correlation between the peak pattern generated by the model and the experimental peaks in the MassPeaks object to be recognized as isotopic pattern."/>
e9300ef37403 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit ecdc3a64aa245d80dbc5487b2bf10a85a43adc6d
galaxyp
parents: 3
diff changeset
578 <param name="tolerance" type="float" label="Tolerance" value="0.0001"
e9300ef37403 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit ecdc3a64aa245d80dbc5487b2bf10a85a43adc6d
galaxyp
parents: 3
diff changeset
579 help="Maximal relative deviation of a peak position (m/z) to be considered as identical: abs(((mass[i]+distance)-mass[i+1])/mass[i]) smaller than 'tolerance'. For 100ppm use 0.0001" />
e9300ef37403 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit ecdc3a64aa245d80dbc5487b2bf10a85a43adc6d
galaxyp
parents: 3
diff changeset
580 <param name="distance" type="text" label="Distance" value="1.00235" help="Distance between two consecutive peaks in an isotopic pattern. 1.00235 is average distance for polypeptides. Multiple values can be used to find multiple charged pattern e.g. 1, 0.5 ,0.33">
e9300ef37403 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit ecdc3a64aa245d80dbc5487b2bf10a85a43adc6d
galaxyp
parents: 3
diff changeset
581 <sanitizer invalid_char="">
e9300ef37403 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit ecdc3a64aa245d80dbc5487b2bf10a85a43adc6d
galaxyp
parents: 3
diff changeset
582 <valid initial="string.digits">
e9300ef37403 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit ecdc3a64aa245d80dbc5487b2bf10a85a43adc6d
galaxyp
parents: 3
diff changeset
583 <add value="," />
e9300ef37403 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit ecdc3a64aa245d80dbc5487b2bf10a85a43adc6d
galaxyp
parents: 3
diff changeset
584 <add value=":" />
e9300ef37403 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit ecdc3a64aa245d80dbc5487b2bf10a85a43adc6d
galaxyp
parents: 3
diff changeset
585 <add value="." />
e9300ef37403 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit ecdc3a64aa245d80dbc5487b2bf10a85a43adc6d
galaxyp
parents: 3
diff changeset
586 </valid>
e9300ef37403 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit ecdc3a64aa245d80dbc5487b2bf10a85a43adc6d
galaxyp
parents: 3
diff changeset
587 </sanitizer>
e9300ef37403 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit ecdc3a64aa245d80dbc5487b2bf10a85a43adc6d
galaxyp
parents: 3
diff changeset
588 </param>
e9300ef37403 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit ecdc3a64aa245d80dbc5487b2bf10a85a43adc6d
galaxyp
parents: 3
diff changeset
589 <param name="size" type="text" label="Size" value="3:10" help="Size (length) of isotopic pattern, longer patterns are prefered over shorter ones, min size is 2, a range can be used.">
e9300ef37403 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit ecdc3a64aa245d80dbc5487b2bf10a85a43adc6d
galaxyp
parents: 3
diff changeset
590 <sanitizer invalid_char="">
e9300ef37403 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit ecdc3a64aa245d80dbc5487b2bf10a85a43adc6d
galaxyp
parents: 3
diff changeset
591 <valid initial="string.digits">
e9300ef37403 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit ecdc3a64aa245d80dbc5487b2bf10a85a43adc6d
galaxyp
parents: 3
diff changeset
592 <add value=":" />
e9300ef37403 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit ecdc3a64aa245d80dbc5487b2bf10a85a43adc6d
galaxyp
parents: 3
diff changeset
593 </valid>
e9300ef37403 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit ecdc3a64aa245d80dbc5487b2bf10a85a43adc6d
galaxyp
parents: 3
diff changeset
594 </sanitizer>
e9300ef37403 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit ecdc3a64aa245d80dbc5487b2bf10a85a43adc6d
galaxyp
parents: 3
diff changeset
595 </param>
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
596 </when>
3
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
597
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
598 <when value="Align">
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
599 <param name="warping_method" type="select" label="Warping methods">
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
600 <option value="lowess" selected="True">Lowess</option>
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
601 <option value="linear">Linear</option>
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
602 <option value="quadratic">Quadratic</option>
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
603 <option value="cubic">Cubic</option>
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
604 </param>
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
605 <param name="tolerance" type="float" value="0.00005"
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
606 label="Tolerance = abs(mz1 - mz2)/mz2"
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
607 help="Maximal relative deviation of a peak position (m/z) to be considered as identical. For 50ppm use 0.00005 or 5e-5" />
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
608 <param name="allow_nomatch" type="boolean" label="Don't throw an error when less than 2 reference m/z were found in a spectrum" truevalue="TRUE" falsevalue="FALSE"/>
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
609 <param name="empty_nomatch" type="boolean" label="If TRUE the intensity values of MassSpectrum or MassPeaks objects with missing (NA) warping functions are set to zero" truevalue="TRUE" falsevalue="FALSE"/>
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
610 <param name="remove_empty" type="boolean" label="Should empty spectra be removed" truevalue="TRUE" falsevalue="FALSE" help="For Cardinal RData files this step can only be performed if pixel annotations were provided"/>
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
611
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
612 <conditional name="reference_for_alignment">
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
613 <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">
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
614 <option value="no_reference" selected="True">no reference</option>
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
615 <option value="yes_reference">reference from tabular file</option>
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
616 </param>
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
617 <when value="no_reference">
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
618 <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"/>
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
619 </when>
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
620 <when value="yes_reference">
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
621 <param name="reference_file" type="data" format="tabular"
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
622 label="Reference m/z values"
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
623 help="Tabular file"/>
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
624 <param name="mz_column" data_ref="reference_file" label="Column with m/z values" type="data_column"/>
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
625 <param name="reference_header" type="boolean" label="File contains a header line" truevalue="TRUE" falsevalue="FALSE"/>
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
626 <param name="separate_alignment" type="boolean" label="Spectrum wise alignment" help="Internal binning is omitted to avoid interaction between spectra" truevalue="TRUE" falsevalue="FALSE"/>
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
627 </when>
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
628 </conditional>
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
629 </when>
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
630
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
631 <when value="Binning">
3
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
632 <param name="bin_tolerance" type="float" value="0.002" label="Tolerance"
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
633 help="After the alignment the peak positions (mass) are very similar but not identical. The binning is needed to make similar peak mass values identical."/>
4
e9300ef37403 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit ecdc3a64aa245d80dbc5487b2bf10a85a43adc6d
galaxyp
parents: 3
diff changeset
634 <param name="bin_method" display="radio" type="select" label="Bin creation rule" help="strict - creates bins never containing two or more peaks of the same sample. relaxed - allows multiple peaks of the same sample in one bin.">
2
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
635 <option value="strict" selected="True" >strict</option>
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
636 <option value="relaxed" >relaxed</option>
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
637 </param>
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
638 </when>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
639 <when value="Filtering">
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
640 <param name="minFrequency" type="float" value="0.25"
2
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
641 label="Removal of all peaks which occur in less than minFrequency spectra" help="Relative threshold. The higher value from relative and absolute threshold is taken. Set one value to zero to be sure it will not be used."/>
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
642 <param name="minNumber" type="float" value="1.0"
2
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
643 label="Removal of all peaks which occur in less than minNumber spectra" help="Absolute threshold. The higher value from relative and absolute threshold is taken. Set one value to zero to be sure it will not be used."/>
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
644 <param name="filter_annot_groups" type="boolean" label="Group wise filtering with pixel annotations." help="If not specified a single group is assumed or when filtering has been done group wise it will automatically be group wise when selecting filtering on all pixel" truevalue="TRUE" falsevalue="FALSE"/>
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
645 <param name="mergeWhitelists" type="boolean" truevalue="TRUE" falsevalue="FALSE"
2
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
646 label="mergeWhitelists" help="Yes means that peaks that survive the filtering in one annotation group are also kept in other groups regardless if the filtering criteria are met in these groups"/>
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
647 </when>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
648 </conditional>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
649 </repeat>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
650 </inputs>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
651 <outputs>
3
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
652 <data format="imzml" name="outfile_imzml" label="${tool.name} on ${on_string}: imzml">
2
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
653 <!--filter>methods_conditional['method'] == 'Peak_detection'</filter-->
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
654 </data>
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
655 <data format="pdf" name="plots" from_work_dir="peaks_qc_plot.pdf" label = "${tool.name} on ${on_string}: QC"/>
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
656 <data format="tabular" name="masspeaks" label="${tool.name} on ${on_string}: peaklist"/>
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
657 <data format="tabular" name="intensity_matrix" label="${tool.name} on ${on_string}: intensity_matrix"/>
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
658 </outputs>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
659 <tests>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
660 <test>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
661 <param name="infile" value="" ftype="imzml">
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
662 <composite_data value="Example_Continuous.imzML"/>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
663 <composite_data value="Example_Continuous.ibd"/>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
664 </param>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
665 <conditional name="tabular_annotation">
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
666 <param name="load_annotation" value="yes_annotation"/>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
667 <param name="annotation_file" value="pixel_annotations.tabular"/>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
668 <param name="column_x" value="1"/>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
669 <param name="column_y" value="2"/>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
670 <param name="column_names" value="3"/>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
671 <param name="tabular_header" value="TRUE"/>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
672 </conditional>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
673 <repeat name="methods">
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
674 <conditional name="methods_conditional">
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
675 <param name="method" value="Peak_detection"/>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
676 <param name="peak_method" value="SuperSmoother"/>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
677 <param name="halfWindowSize" value="1"/>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
678 <param name="snr" value="5"/>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
679 <param name="use_annotations" value="TRUE"/>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
680 </conditional>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
681 </repeat>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
682 <output name="plots" file="peakdetection1_QC.pdf" compare="sim_size"/>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
683 <output name="masspeaks" file="masspeaks1.tabular"/>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
684 <output name="intensity_matrix" file="int1.tabular"/>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
685 </test>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
686 <test>
1
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
687 <param name="infile" value="masspeaks3_forinput.tabular"/>
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
688 <param name="centroids" value="TRUE"/>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
689 <repeat name="methods">
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
690 <conditional name="methods_conditional">
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
691 <param name="method" value="monoisotopic_peaks"/>
2
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
692 <param name="tolerance" value="0.0004"/>
4
e9300ef37403 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit ecdc3a64aa245d80dbc5487b2bf10a85a43adc6d
galaxyp
parents: 3
diff changeset
693 <param name="size" value="3"/>
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
694 </conditional>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
695 </repeat>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
696 <output name="plots" file="peakdetection2_QC.pdf" compare="sim_size"/>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
697 <output name="masspeaks" file="masspeaks2.tabular"/>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
698 <output name="intensity_matrix" file="int2.tabular"/>
3
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
699 <output name="outfile_imzml" ftype="imzml" file="peak_detection2.imzml.txt" lines_diff="4">
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
700 <extra_files type="file" file="peak_detection2.imzml" name="imzml" lines_diff="6"/>
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
701 <extra_files type="file" file="peak_detection2.ibd" name="ibd" compare="sim_size"/>
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
702 </output>
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
703 </test>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
704 <test>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
705 <param name="infile" value="" ftype="imzml">
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
706 <composite_data value="Example_Continuous.imzML"/>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
707 <composite_data value="Example_Continuous.ibd"/>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
708 </param>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
709 <conditional name="tabular_annotation">
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
710 <param name="load_annotation" value="yes_annotation"/>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
711 <param name="annotation_file" value="pixel_annotations.tabular"/>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
712 <param name="column_x" value="1"/>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
713 <param name="column_y" value="2"/>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
714 <param name="column_names" value="3"/>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
715 <param name="tabular_header" value="TRUE"/>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
716 </conditional>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
717 <repeat name="methods">
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
718 <conditional name="methods_conditional">
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
719 <param name="method" value="Peak_detection"/>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
720 <param name="peak_method" value="MAD"/>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
721 <param name="halfWindowSize" value="1"/>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
722 <param name="snr" value="2"/>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
723 </conditional>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
724 </repeat>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
725 <repeat name="methods">
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
726 <conditional name="methods_conditional">
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
727 <param name="method" value="Binning"/>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
728 <param name="bin_tolerance" value="0.01"/>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
729 </conditional>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
730 </repeat>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
731 <repeat name="methods">
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
732 <conditional name="methods_conditional">
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
733 <param name="method" value="Filtering"/>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
734 <param name="minFrequency" value="0.5"/>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
735 <param name="minNumber" value="3"/>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
736 <param name="filter_annot_groups" value="TRUE"/>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
737 <param name="mergeWhitelists" value="FALSE"/>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
738 </conditional>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
739 </repeat>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
740 <output name="plots" file="peakdetection3_QC.pdf" compare="sim_size"/>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
741 <output name="intensity_matrix" file="intensity_matrix3.tabular"/>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
742 <output name="masspeaks" file="masspeaks3.tabular"/>
3
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
743 <output name="outfile_imzml" ftype="imzml" file="peak_detection3.imzml.txt" lines_diff="4">
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
744 <extra_files type="file" file="peak_detection3.imzml" name="imzml" lines_diff="6"/>
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
745 <extra_files type="file" file="peak_detection3.ibd" name="ibd" compare="sim_size"/>
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
746 </output>
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
747 </test>
2
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
748 <test>
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
749 <param name="infile" value="testfile_squares.rdata" ftype="rdata"/>
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
750 <param name="method" value="Peak_detection"/>
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
751 <param name="peak_method" value="MAD"/>
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
752 <param name="halfWindowSize" value="20"/>
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
753 <param name="snr" value="2"/>
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
754 <output name="plots" file="peakdetection4_QC.pdf" compare="sim_size"/>
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
755 <output name="intensity_matrix" file="intensity_matrix4.tabular"/>
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
756 <output name="masspeaks" file="masspeaks4.tabular"/>
3
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
757 <output name="outfile_imzml" ftype="imzml" file="peak_detection4.imzml.txt" lines_diff="4">
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
758 <extra_files type="file" file="peak_detection4.imzml" name="imzml" lines_diff="6"/>
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
759 <extra_files type="file" file="peak_detection4.ibd" name="ibd" compare="sim_size"/>
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
760 </output>
2
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
761 </test>
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
762 </tests>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
763 <help>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
764 <![CDATA[
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
765
1
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
766 @MADLI_QUANT_DESCRIPTION@
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
767
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
768 -----
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
769
1
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
770 **Input data**
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
771
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
772 - MSI data: 3 types of input data can be used:
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
773
1
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
774 - imzml file (upload imzml and ibd file via the "composite" function) `Introduction to the imzml format <https://ms-imaging.org/wp/imzml/>`_
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
775 - Cardinal "MSImageSet" data saved as .RData
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
776 - MSI data as peak list (tabular file) with the columns named "snr", "mass", "intensity" and "spectrum". The spectrum has to be in the following format: xy_1_1 (for pixel coordinates x1y1). The header must have exactly the four column names.
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
777
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
778 ::
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
779
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
780 snr mass intensity spectrum
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
781 5.34 304.16 0.10 xy_1_1
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
782 12.09 305 0.2 xy_1_1
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
783 6.80 306.25 0.133 xy_1_1
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
784 ...
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
785 ...
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
786
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
787
3
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
788 - 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.
1
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
789
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
790 ::
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
791
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
792 x_coord y_coord
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
793 1 1
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
794 2 1
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
795 3 1
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
796 ...
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
797 ...
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
798
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
799
2
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
800 - Optional: Tabular file(s) containing pixel coordinates and annotation. X and y values in separate columns and the corresponding annotation in a third column. Tabular files with any header name or no header at all are supported. The annotations can be used to either detect peak on the average mass spectrum of each annotation group ('Detect peaks on average mass spectra') or to filter for peaks within pixel groups ('Group wise filtering'). This option is not compatible with tabular peak list inputs.
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
801
1
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
802 ::
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
803
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
804 x_coord y_coord annotation
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
805 1 1 healthy
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
806 2 1 healthy
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
807 3 1 disease
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
808 ...
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
809 ...
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
810
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
811 **Options**
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
812
2
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
813 - Peak detection: detection of peaks, only possible with profile mode input. A peak is a local maximum above a user defined noise threshold.
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
814
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
815 - Two functions exist to estimate the noise: MAD and supersmoother.
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
816 - SNR is an abbreviation for signal-to-noise-ratio. A local maximum has to be higher than SNR*noise to be recognize as peak.
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
817 - 'halfWindowSize': The resulting window reaches from mass[currentIndex-halfWindowSize] to mass[currentIndex+halfWindowSize]. A local maximum has to be the highest one in the given window to be recognized as peak.
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
818 - Detect peaks on average mass spectra: Spectra with the same annotation (taken from the annotation tabular input) are averaged and peak picking is performed on the average spectrum of each annotation group. The exported imzML is empty and cannot be used for further analysis steps. The peaklist and intensity matrix outputs contain the annotation group names with their averaged intensity values. Filtering steps have to be done in the same run as the peak picking.
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
819
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
820
3
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
821 - Monoisotopic peaks: Keeps only the monoisotopic peaks
2
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
822
4
e9300ef37403 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit ecdc3a64aa245d80dbc5487b2bf10a85a43adc6d
galaxyp
parents: 3
diff changeset
823 - Based on poisson model for isotopic patterns as decribed in (`Breen et al. <https://doi.org/10.1002/1522-2683(20000601)21:11%3C2243::AID-ELPS2243%3E3.0.CO;2-K>`_)
e9300ef37403 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit ecdc3a64aa245d80dbc5487b2bf10a85a43adc6d
galaxyp
parents: 3
diff changeset
824 - Isotopic pattern can be characterized and recognized by
e9300ef37403 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit ecdc3a64aa245d80dbc5487b2bf10a85a43adc6d
galaxyp
parents: 3
diff changeset
825
e9300ef37403 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit ecdc3a64aa245d80dbc5487b2bf10a85a43adc6d
galaxyp
parents: 3
diff changeset
826 - the similarity of the experimental pattern with the modelled pattern
e9300ef37403 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit ecdc3a64aa245d80dbc5487b2bf10a85a43adc6d
galaxyp
parents: 3
diff changeset
827 - the distance between consecutive isotopic peaks. For polypeptides the average distance is 1.00235 (`Park et al. <https://pubs.acs.org/doi/abs/10.1021/ac800913b>`_). Multiply charged analytes have smaller distances between the peaks (e.g. z = 1 distance = ~1; z = 2: distance = ~0.5; z = 3: distance = ~0.3333) To search for differently charged isotopic pattern multiple distances can be applied - the order matters because the first distance that matches is reported (1, 0.5, 0.3333).
e9300ef37403 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit ecdc3a64aa245d80dbc5487b2bf10a85a43adc6d
galaxyp
parents: 3
diff changeset
828 - the size (length) of the pattern, multiple values can be applied, longer patterns are prefered over shorter ones.
e9300ef37403 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit ecdc3a64aa245d80dbc5487b2bf10a85a43adc6d
galaxyp
parents: 3
diff changeset
829
2
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
830
3
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
831 - Spectra alignment (warping): alignment for (re)calibration of m/z values.
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
832
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
833 - without external reference m/z: internal reference is obtained by filtering (default 90%) and binning the peaks to find landmark peaks and their average m/z
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
834 - with external reference m/z: the m/z provided in a tabular file are used as a reference, at least 10 reference values are recommended
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
835 - non linear warping (parametric time warping plus binning) to match the reference peaks (internal or external) to the peaks with the given tolerance. At least two m/z per spectrum are needed for the alignment. This requirement can be skipped by setting "Don't throw an error when less than 2 reference m/z were found in a spectrum" 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.
2
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
836
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
837
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
838 - Peak binning: After the alignment the peak positions (m/z) are very similar but not identical. The binning is needed to make similar peak m/z values identical. The algorithm is based on the following workflow:
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
839
3
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
840 1. Put all m/z in a sorted vector.
2
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
841 2. Calculate differences between each neighbor.
3
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
842 3. Divide the m/z vector at the largest gap (largest difference) and form a left and a right bin.
2
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
843 4. Rerun step 3 for the left and/or the right bin if they don't fulfill the following criteria:
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
844
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
845 - All peaks in a bin are near to the mean (abs(mass-meanMass)/meanMass < tolerance).
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
846 - method == "strict": The bin doesn't contain two or more peaks of the same sample.
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
847
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
848
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
849 - Peak filtering: Removal of less frequent m/z features:
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
850
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
851 - minFrequency : between 0 and 1: m/z has to occur in 0 - 100% of all spectra; minNumber: m/z has to occur in at least this amount of spectra --> out of those two criteria the stricter value will be used
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
852 - Group wise filtering with pixel annotations: 'Yes' means that the filtering criteria are applied for each annotation group separately.
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
853 - mergeWhitelists: 'Yes' means that peaks that survive the filtering in one annotation group are also kept in other groups regardless if the filtering criteria are met in these groups
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
854 - To filter data that was averaged before peak detection: Filtering has to be done in the same tool run as the peak detection. The filtering criteria are automatically applied per annotation group (Group wise filtering can be 'No') and not per pixel. Example: to keep only m/z that were detected in at least half of all annotation groups set minFrequency to 0.5.
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
855
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
856
1
eaaa73b043e6 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 0825a4ccd3ebf4ca8a298326d14f3e7b25ae8415
galaxyp
parents: 0
diff changeset
857 **Output**
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
858
3
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
859 - centroided, processed imzML file, imzML file is empty when 'Detect peaks on average mass spectra' is chosen.
36d38d2cf88c planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit f127be2141cf22e269c85282d226eb16fe14a9c1
galaxyp
parents: 2
diff changeset
860 - pdf with mass spectra plots of four random spectra and a table with key values after each preprocessing step
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
861 - peak list (tabular file) with the columns "snr", "mass", "intensity" and "spectrum"
2
17c54820f3be planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit d2f311f7fff24e54c565127c40414de708e31b3c
galaxyp
parents: 1
diff changeset
862 - tabular file with intensity matrix (m/z in rows and spectra in columns). If the input file was imzML in profile mode the intensities before peak picking are also stored in the matrix . For all other inputs not picked values are set to NA. For peak picking on the average of multiple spectra, each spectra group is a column with mean intensities for each m/z.
0
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
863
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
864 .. _MALDIquant: http://strimmerlab.org/software/maldiquant/
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
865
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
866 ]]>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
867 </help>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
868 <expand macro="citation"/>
01212bf66f61 planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/MALDIquant commit 5feaf3d0e0da8cef1241fecc1f4d6f81324594e6
galaxyp
parents:
diff changeset
869 </tool>