msi_combine: msi_combine.xml comparison

comparison msi_combine.xml @ 2:00b6c61f5054 draft

planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/msi_combine commit 1c808d60243bb1eeda0cd26cb4b0a17ab05de2c0

author	galaxyp
date	Mon, 28 May 2018 12:35:49 -0400
parents	f3f6c32ab690
children	91bba2486773

comparison

equal deleted inserted replaced

-:f3f6c32ab690
+:00b6c61f5054
 #if str( $combine_conditional.combine_method ) == 'xy_shifts':
 input_list = read.delim("$combine_conditional.coordinates_file", header = FALSE,
 stringsAsFactors = FALSE)
 #end if
+loadRData <- function(fileName){
+#loads an RData file, and returns it
+load(fileName)
+get(ls()[ls() != "fileName"])
+}
 pixel_vector = numeric()
 names_vector = character()
 x_shifts = 0
 y_shifts = 0
 max_y = numeric()
 #set $msidata = []
 #set $pixelcoords = []
 #set $num_infiles = len($infiles)
+all_files = $num_infiles
 #for $i, $infile in enumerate($infiles):
 #if $infile.ext == 'imzml'
 msidata_$i <- readImzML('infile_${i}')
+#elif $infile.ext == 'analyze75'
+msidata_$i <- readAnalyze('infile_${i}')
+#else
+msidata_$i = loadRData('infile_${i}.RData')
+#end if
+########## wenn ncol msidata > 0 and nrow msidata > 0, dann alles abklappern andernfalls nicht in msiliste mitreinnehmen (unten hinter python vor end for: filenr ausprinten)
 sampleNames(msidata_$i) = "msidata"
 pixelcoords_$i = cbind(coord(msidata_$i)[,1:2], rep($i+1,ncol(msidata_$i)))
-#elif $infile.ext == 'analyze75'
-msidata_$i <- readAnalyze('infile_${i}')
-sampleNames(msidata_$i) = "msidata"
-pixelcoords_$i = cbind(coord(msidata_$i)[,1:2], rep($i+1,ncol(msidata_$i)))
-#else
-loadRData <- function(fileName){
-#loads an RData file, and returns it
-load(fileName)
-get(ls()[ls() != "fileName"])
-}
-msidata_$i = loadRData('infile_${i}.RData')
-sampleNames(msidata_$i) = "msidata"
-pixelcoords_$i = cbind(coord(msidata_$i)[,1:2], rep($i+1,ncol(msidata_$i)))
-#end if
-colnames(pixelcoords_$i)[3] = "file_number"
 #if str( $combine_conditional.combine_method ) == 'xy_shifts':
 coord(msidata_$i)\$x = coord(msidata_$i)\$x + input_list[$i+1,$combine_conditional.column_x]
 coord(msidata_$i)\$y = coord(msidata_$i)\$y + input_list[$i+1,$combine_conditional.column_y]
 pixel_vector = append(pixel_vector, rep(input_list[$i+1,$combine_conditional.column_names],times=ncol(msidata_$i)))
 #end if
 #silent $msidata.append('msidata_'+str($i))
 #silent $pixelcoords.append('pixelcoords_'+str($i))
 #end for
+colnames(pixelcoords_$i)[3] = "file_number"
 ###################### automatic combination ###################################
 ################################################################################
 #if str( $combine_conditional.combine_method ) == 'automatic_combine':
 print("automatic_combine")
 msidata_combined = do.call(combine, list(#echo ','.join($msidata)#))
+print("combination successful")
 sample_names = as.factor(names_vector)
 pData(msidata_combined)\$sample = sample_names
+print("names successful")
 ## create PDF to show pixels of each file
 pdf("combining_qc.pdf", width=15, height=15)
 position_df = cbind(coord(msidata_combined)[,1:2], pData(msidata_combined)\$sample)
 colnames(position_df)[3] = "sample_name"
+print("position_df")
 combine_plot = ggplot(position_df, aes(x=x, y=y, fill=sample_name))+
 geom_tile() +
 coord_fixed()+
 ggtitle("Spatial orientation of combined data")+
 <output name="combining_qc" file="12_auto_combined_QC.pdf" compare="sim_size" delta="20000"/>
 </test>
 </tests>
 <help>
 <![CDATA[
-This tool can combine several mass-spectrometry imaging files.
-1) m/z values need to be the same across all datasets
+Cardinal is an R package that implements statistical & computational tools for analyzing mass spectrometry imaging datasets. `More information on Cardinal <http://cardinalmsi.org//>`_
-2) pixels (defined by x and y coordinates) must be unique
+This tool uses the Cardinal combine function to combine several mass-spectrometry imaging data.
-1) Same m/z values/axis can be achieved with the filtering tool to get all datasets to the same m/z range and afterwards binning in the preprocessing tool to obtain the same bins for all dataset.
-2) The pixels (defined by x and y coordinates) must be unique across all datasets, therefore the option "Select the way you want to combine multiple files" is helpful:
-- "automatic combination": files are arranged in a grid with a distance in x and y direction which can be given by the user
-- "no coordinates shift": this option can only be used if all pixels are unique across datasets
-- "xy shifts by hand": each file can be moved in x and y direction according to the users need (define one tabular file in the order in which the files are loaded in the history (bottom to top) and define for each file the x and y coordinates shifts in separate columns
-- "check pixels before combination": no combination takes place. You will only get a pdf which shows the arrangement of the pixels (with or without additional xy shifts)
 Input data: 3 types of input data can be used:
-- imzml file (upload imzml and ibd file via the "composite" function) `Introduction to the imzml format <http://ms-imaging.org/wp/introduction/>`_
+- imzml file (upload imzml and ibd file via the "composite" function) `Introduction to the imzml format <https://ms-imaging.org/wp/imzml/>`_
 - Analyze7.5 (upload hdr, img and t2m file via the "composite" function)
 - Cardinal "MSImageSet" data (with variable name "msidata", saved as .RData)
-The output of this tool is a Cardinal "MSImageSet" with the variable name "msidata" saved as .RData.
+Prerequisite:
+- m/z values need to be the same across all datasets (before using this tool), this can be achieved with the filtering tool (use same m/z range) and the preprocessing tool (use same binning parameter)
+Options:
+- "automatic combination": files are arranged in a grid with a minimal distance in x and y direction which can be given by the user, subfiles are named according to input file name
+- "no coordinates shift": this option can only be used if all pixels are unique across datasets, no assignment of names to the subfiles
+- "xy shifts by hand": each file can be moved in x and y direction according to the users need (define one tabular file in the order in which the files are loaded in the history (bottom to top) and define for each file the x and y coordinates shifts in separate columns
+- "check pixels before combination": no combination takes place but a pdf shows the current arrangement of the pixels
+Output:
+- imzML file containing multiple subfiles
+- pdf that shows the pixel positions of the combined files
+- optional: intensity matrix as tabular file (intensities for masses in rows and pixel in columns)
 ]]>
 </help>
 <citations>
 <citation type="doi">10.1093/bioinformatics/btv146</citation>
 </citations>

Mercurial > repos > galaxyp > msi_combine

comparison msi_combine.xml @ 2:00b6c61f5054 draft