Repository 'msi_qualitycontrol'
hg clone https://toolshed.g2.bx.psu.edu/repos/galaxyp/msi_qualitycontrol

Changeset 8:52ef77866de8 (2018-05-28)
Previous changeset 7:b86a66dd1a16 (2018-05-08) Next changeset 9:963c7ec00141 (2018-06-11)
Commit message:
planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/msi_qualitycontrol commit 1c808d60243bb1eeda0cd26cb4b0a17ab05de2c0
modified:
msi_qualitycontrol.xml
test-data/Analyze75.hdr
test-data/Analyze75.img
test-data/Analyze75.t2m
test-data/inputcalibrantfile2.txt
test-data/inputpeptides.txt
added:
test-data/QC_analyze75.pdf
test-data/QC_empty_spectra.pdf
test-data/QC_imzml.pdf
test-data/QC_rdata.pdf
test-data/empty_spectra.rdata
test-data/preprocessed.RData
removed:
test-data/LM8_file16.rdata
test-data/LM8_file16output.pdf
test-data/Testfile_qualitycontrol_analyze75.pdf
test-data/Testfile_qualitycontrol_imzml.pdf
test-data/Testfile_qualitycontrol_rdata.pdf
test-data/example_continousS042.RData
test-data/inputpeptides.csv
test-data/preprocessing_results1.RData
b
diff -r b86a66dd1a16 -r 52ef77866de8 msi_qualitycontrol.xml
--- a/msi_qualitycontrol.xml Tue May 08 02:36:43 2018 -0400
+++ b/msi_qualitycontrol.xml Mon May 28 12:38:50 2018 -0400
[
b'@@ -1,9 +1,9 @@\n-<tool id="mass_spectrometry_imaging_qc" name="MSI Qualitycontrol" version="1.7.0.5">\n+<tool id="mass_spectrometry_imaging_qc" name="MSI Qualitycontrol" version="1.10.0.0">\n     <description>\n         mass spectrometry imaging QC\n     </description>\n     <requirements>\n-        <requirement type="package" version="1.7.0">bioconductor-cardinal</requirement>\n+        <requirement type="package" version="1.10.0">bioconductor-cardinal</requirement>\n         <requirement type="package" version="2.2.1">r-ggplot2</requirement>\n         <requirement type="package" version="1.1_2">r-rcolorbrewer</requirement>\n         <requirement type="package" version="2.2.1">r-gridextra</requirement>\n@@ -12,12 +12,12 @@\n     <command detect_errors="exit_code">\n     <![CDATA[\n         #if $infile.ext == \'imzml\'\n-            cp \'${infile.extra_files_path}/imzml\' infile.imzML &&\n-            cp \'${infile.extra_files_path}/ibd\' infile.ibd &&\n+            ln -s \'${infile.extra_files_path}/imzml\' infile.imzML &&\n+            ln -s \'${infile.extra_files_path}/ibd\' infile.ibd &&\n         #elif $infile.ext == \'analyze75\'\n-            cp \'${infile.extra_files_path}/hdr\' infile.hdr &&\n-            cp \'${infile.extra_files_path}/img\' infile.img &&\n-            cp \'${infile.extra_files_path}/t2m\' infile.t2m &&\n+            ln -s \'${infile.extra_files_path}/hdr\' infile.hdr &&\n+            ln -s \'${infile.extra_files_path}/img\' infile.img &&\n+            ln -s \'${infile.extra_files_path}/t2m\' infile.t2m &&\n         #else\n             ln -s \'$infile\' infile.RData &&\n         #end if\n@@ -37,10 +37,9 @@\n ## Read MALDI Imaging dataset\n \n #if $infile.ext == \'imzml\'\n-    msidata = readMSIData(\'infile.imzML\')\n+    msidata = readImzML(\'infile\')\n #elif $infile.ext == \'analyze75\'\n-    msidata = readMSIData(\'infile.hdr\')\n-\n+    msidata = readAnalyze(\'infile\')\n #else\n     load(\'infile.RData\')\n #end if\n@@ -107,7 +106,7 @@\n   peakpickinginfo=processinginfo@peakPicking\n }\n \n-### Read tabular file with peptide masses for plots and heatmap images: \n+### Read tabular file with masses for plots and heatmap images: \n \n #if $peptide_file:\n \n@@ -143,10 +142,8 @@\n         number_calibrants_in = length(calibrant_list[,1])\n         number_calibrants_valid = length(inputcalibrants[,1])\n #else\n-    ###inputcalibrants = data.frame(0,0)\n \n     inputcalibrants = as.data.frame(matrix(, nrow = 0, ncol = 2))\n-\n     number_calibrants_in = 0\n     number_calibrants_valid = 0\n #end if\n@@ -372,7 +369,7 @@\n         {\n           image(msidata, mz=inputmasses[mass], plusminus=$plusminus_dalton, \n           main= paste0(inputnames[mass], " (", round(inputmasses[mass], digits = 2)," \xc2\xb1 ", $plusminus_dalton, " Da)"),\n-                contrast.enhance = "histogram", ylim=c(maximumy+2, 0))\n+                contrast.enhance = "histogram", ylim= c(maximumy+0.2*maximumy,minimumy-0.2*minimumy))\n         }\n     } else {print("3) The inputpeptide masses were not provided or outside the mass range")}\n \n@@ -428,7 +425,7 @@\n     pca = PCA(msidata, ncomp=2) \n     par(mfrow = c(2,1))\n     plot(pca, col=c("black", "darkgrey"), main="PCA for two components")\n-    image(pca, col=c("black", "white"),ylim=c(maximumy+2, 0),  strip=FALSE)\n+    image(pca, col=c("black", "white"),ylim= c(maximumy+0.2*maximumy,minimumy-0.2*minimumy),  strip=FALSE)\n \n \n     ############################# III) properties over acquisition (spectra index)##########\n@@ -600,12 +597,12 @@\n         <param name="infile" type="data" format="imzml,rdata,analyze75" label="Inputfile as imzML, Analyze7.5 or Cardinal MSImageSet saved as RData"\n             help="Upload composite datatype imzml (ibd+imzML) or analyze75 (hdr+img+t2m) or regular upload .RData (Cardinal MSImageSet)"/>\n         <param name="filename" type="text" value="" optional="true" label="Title" help="will appear in the quality report. If nothing given it will take the dataset name."/>\n-        <param name="peptide_file" type="data" optional="true" format="tabular" label="Text '..b'75.img"/>\n                 <composite_data value="Analyze75.t2m"/>\n             </param>\n-            <param name="peptide_file" value="inputpeptides.txt" ftype="txt"/>\n-            <param name="calibrant_file" ftype="txt" value="inputcalibrantfile2.txt"/>\n+            <param name="peptide_file" value="inputpeptides.txt"/>\n+            <param name="calibrant_file" value="inputcalibrantfile2.txt"/>\n             <param name="plusminus_dalton" value="0.5"/>\n             <param name="filename" value="Testfile_analyze75"/>\n-            <output name="plots" file="Testfile_qualitycontrol_analyze75.pdf" compare="sim_size" delta="20000"/>\n+            <output name="plots" file="QC_analyze75.pdf" compare="sim_size" delta="20000"/>\n         </test>\n-\n         <test>\n-            <param name="infile" value="preprocessing_results1.RData" ftype="rdata"/>\n+            <param name="infile" value="preprocessed.RData" ftype="rdata"/>\n+            <param name="plusminus_dalton" value="0"/>\n+            <param name="filename" value="Testfile_rdata"/>\n+            <output name="plots" file="QC_rdata.pdf" compare="sim_size" delta="20000"/>\n+        </test>\n+        <test>\n+            <param name="infile" value="empty_spectra.rdata" ftype="rdata"/>\n+            <param name="peptide_file" value="inputpeptides.txt"/>\n+            <param name="calibrant_file" value="inputcalibrantfile2.txt"/>\n             <param name="plusminus_dalton" value="0.1"/>\n             <param name="filename" value="Testfile_rdata"/>\n-            <output name="plots" file="Testfile_qualitycontrol_rdata.pdf" compare="sim_size" delta="20000"/>\n-        </test>\n-        <test>\n-            <param name="infile" value="LM8_file16.rdata" ftype="rdata"/>\n-            <param name="peptide_file" value="inputpeptides.txt" ftype="txt"/>\n-            <param name="calibrant_file" ftype="txt" value="inputcalibrantfile2.txt"/>\n-            <param name="plusminus_dalton" value="0.1"/>\n-            <param name="filename" value="Testfile_rdata"/>\n-            <output name="plots" file="LM8_file16output.pdf" compare="sim_size" delta="20000"/>\n+            <output name="plots" file="QC_empty_spectra.pdf" compare="sim_size" delta="20000"/>\n         </test>\n     </tests>\n     <help>\n         <![CDATA[\n-Quality control for maldi imaging mass spectrometry data. The output of this tool contains key values and plots of the imaging data as pdf. \n-For additional beautiful heatmap images use the MSI ion images tool and to plot more mass spectra use the MSI massspectra tool. \n+Cardinal is an R package that implements statistical & computational tools for analyzing mass spectrometry imaging datasets. `More information on Cardinal <http://cardinalmsi.org//>`_\n+\n+This tool uses some Cardinal functions to create a quality control report with descriptive plots for mass-spectrometry imaging data. \n \n Input data: 3 types of input data can be used:\n \n-- imzml file (upload imzml and ibd file via the "composite" function) `Introduction to the imzml format <http://ms-imaging.org/wp/introduction/>`_\n+- imzml file (upload imzml and ibd file via the "composite" function) `Introduction to the imzml format <https://ms-imaging.org/wp/imzml/>`_\n - Analyze7.5 (upload hdr, img and t2m file via the "composite" function)\n - Cardinal "MSImageSet" data (with variable name "msidata", saved as .RData)\n \n+Options: \n \n+- masses of interest as tabular file, used to generate heatmap images\n+- internal calibrants as tabular file, used for the following plots: Number of calibrant per spectrum, heatmap images, mass-spectrum plot zoomed in for calibrant region, ppm accuracy\n+- fold change plot: draws a heatmap of the fold change of two masses (log2(intensity ratio))\n+\n+Output: \n+\n+- pdf with numbers and descriptive plots to check the quality of the mass-spectrometry imaging data\n+\n+Tip: \n+\n+- For additional heatmap images use the MSI ion images tool and to plot more mass spectra use the MSI massspectra tool. \n \n         ]]>\n     </help>\n'
b
diff -r b86a66dd1a16 -r 52ef77866de8 test-data/Analyze75.hdr
b
Binary file test-data/Analyze75.hdr has changed
b
diff -r b86a66dd1a16 -r 52ef77866de8 test-data/Analyze75.img
b
Binary file test-data/Analyze75.img has changed
b
diff -r b86a66dd1a16 -r 52ef77866de8 test-data/Analyze75.t2m
b
Binary file test-data/Analyze75.t2m has changed
b
diff -r b86a66dd1a16 -r 52ef77866de8 test-data/LM8_file16.rdata
b
Binary file test-data/LM8_file16.rdata has changed
b
diff -r b86a66dd1a16 -r 52ef77866de8 test-data/LM8_file16output.pdf
b
Binary file test-data/LM8_file16output.pdf has changed
b
diff -r b86a66dd1a16 -r 52ef77866de8 test-data/QC_analyze75.pdf
b
Binary file test-data/QC_analyze75.pdf has changed
b
diff -r b86a66dd1a16 -r 52ef77866de8 test-data/QC_empty_spectra.pdf
b
Binary file test-data/QC_empty_spectra.pdf has changed
b
diff -r b86a66dd1a16 -r 52ef77866de8 test-data/QC_imzml.pdf
b
Binary file test-data/QC_imzml.pdf has changed
b
diff -r b86a66dd1a16 -r 52ef77866de8 test-data/QC_rdata.pdf
b
Binary file test-data/QC_rdata.pdf has changed
b
diff -r b86a66dd1a16 -r 52ef77866de8 test-data/Testfile_qualitycontrol_analyze75.pdf
b
Binary file test-data/Testfile_qualitycontrol_analyze75.pdf has changed
b
diff -r b86a66dd1a16 -r 52ef77866de8 test-data/Testfile_qualitycontrol_imzml.pdf
b
Binary file test-data/Testfile_qualitycontrol_imzml.pdf has changed
b
diff -r b86a66dd1a16 -r 52ef77866de8 test-data/Testfile_qualitycontrol_rdata.pdf
b
Binary file test-data/Testfile_qualitycontrol_rdata.pdf has changed
b
diff -r b86a66dd1a16 -r 52ef77866de8 test-data/empty_spectra.rdata
b
Binary file test-data/empty_spectra.rdata has changed
b
diff -r b86a66dd1a16 -r 52ef77866de8 test-data/example_continousS042.RData
b
Binary file test-data/example_continousS042.RData has changed
b
diff -r b86a66dd1a16 -r 52ef77866de8 test-data/inputcalibrantfile2.txt
--- a/test-data/inputcalibrantfile2.txt Tue May 08 02:36:43 2018 -0400
+++ b/test-data/inputcalibrantfile2.txt Mon May 28 12:38:50 2018 -0400
b
@@ -1,3 +1,3 @@
 869.51 mass1
-1001.62 mass2
-1023.6 mass3
+1111.1 mass2
+1296.7 mass3
b
diff -r b86a66dd1a16 -r 52ef77866de8 test-data/inputpeptides.csv
--- a/test-data/inputpeptides.csv Tue May 08 02:36:43 2018 -0400
+++ /dev/null Thu Jan 01 00:00:00 1970 +0000
b
@@ -1,3 +0,0 @@
-152 mass1
-328.9 mass2
-185.2 mass3
b
diff -r b86a66dd1a16 -r 52ef77866de8 test-data/inputpeptides.txt
--- a/test-data/inputpeptides.txt Tue May 08 02:36:43 2018 -0400
+++ b/test-data/inputpeptides.txt Mon May 28 12:38:50 2018 -0400
b
@@ -1,3 +1,6 @@
+152
+328.9
+185.9
 854.5
-1296.7
+1250.1
 2000.8
b
diff -r b86a66dd1a16 -r 52ef77866de8 test-data/preprocessed.RData
b
Binary file test-data/preprocessed.RData has changed
b
diff -r b86a66dd1a16 -r 52ef77866de8 test-data/preprocessing_results1.RData
b
Binary file test-data/preprocessing_results1.RData has changed