Mercurial > repos > galaxyp > cardinal_filtering
changeset 9:c3d2ca5f9dd0 draft
"planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/cardinal commit e499c9124d3fd85a7fc47b95c206ce91a5e3678c-dirty"
line wrap: on
line diff
--- a/filtering.xml Wed May 13 17:54:30 2020 +0000 +++ b/filtering.xml Tue Nov 03 22:26:34 2020 +0000 @@ -5,7 +5,7 @@ </macros> <expand macro="requirements"> <requirement type="package" version="2.3">r-gridextra</requirement> - <requirement type="package" version="3.2.1">r-ggplot2</requirement> + <requirement type="package" version="3.3.2">r-ggplot2</requirement> </expand> <expand macro="print_version"/> <command detect_errors="exit_code"> @@ -365,7 +365,7 @@ theme(plot.title = element_text(hjust = 0.5))+ theme(legend.position="bottom",legend.direction="vertical")+ guides(fill=guide_legend(ncol=5,byrow=TRUE)) - print(pixel_image) + print(pixel_image + scale_fill_manual(values=c("#00BFC4", "#F8766D"))) ### visual mz feature control @@ -384,10 +384,9 @@ theme(legend.position="bottom",legend.direction="vertical")+ labs(title="Overview of filtered m/z", x="m/z", y = "count") + guides(fill=guide_legend(ncol=5,byrow=TRUE)) - print(histogram_mz) + print(histogram_mz + scale_fill_manual(values=c("#00BFC4", "#F8766D"))) - - dev.off() + dev.off() ## save msidata as imzML file, will only work if there is at least 1 m/z left @@ -395,6 +394,7 @@ ## make sure that coordinates are integers coord(msidata)\$y = as.integer(coord(msidata)\$y) coord(msidata)\$x = as.integer(coord(msidata)\$x) + msidata = as(msidata, "MSContinuousImagingExperiment") writeImzML(msidata, "out")} @@ -525,7 +525,7 @@ <expand macro="processed_infile_imzml"/> <conditional name="processed_cond"> <param name="processed_file" value="processed"/> - <param name="accuracy" value="100"/> + <param name="accuracy" value="200"/> <param name="units" value="ppm"/> </conditional> <param name="pixel_filtering" value="two_columns"/>
--- a/macros.xml Wed May 13 17:54:30 2020 +0000 +++ b/macros.xml Tue Nov 03 22:26:34 2020 +0000 @@ -1,10 +1,10 @@ <macros> - <token name="@VERSION@">2.4.0</token> + <token name="@VERSION@">2.6.0</token> <xml name="requirements"> <requirements> <requirement type="package" version="@VERSION@">bioconductor-cardinal</requirement> - <requirement type="package" version="3.6.1">r-base</requirement> + <!--requirement type="package" version="3.6.1">r-base</requirement--> <yield/> </requirements> </xml> @@ -117,6 +117,13 @@ <token name="@DATA_PROPERTIES_INRAM@"><![CDATA[ ########################### QC numbers ######################## ## including intensity calculations which need data in RAM + + int_matrix = as.matrix(spectra(msidata)) ## only load once into RAM, then re-use + ## Number of NA in spectra matrix + NAcount = sum(is.na(int_matrix)) + ## replace NA with zero to calculate data properties based on intensity matrix, no change in msidata + int_matrix[is.na(int_matrix)] <- 0 + ## Number of features (mz) maxfeatures = length(features(msidata)) ## Range mz @@ -131,14 +138,12 @@ minimumy = min(coord(msidata)[,2]) maximumy = max(coord(msidata)[,2]) ## Range of intensities - minint = round(min(as.matrix(spectra(msidata)), na.rm=TRUE), digits=2) - maxint = round(max(as.matrix(spectra(msidata)), na.rm=TRUE), digits=2) + minint = round(min(int_matrix), digits=2) + maxint = round(max(int_matrix), digits=2) ## Number of intensities > 0, for if conditions - npeaks= sum(as.matrix(spectra(msidata))>0, na.rm=TRUE) + npeaks= sum(int_matrix>0) ## Number of NA in spectra matrix - NAcount = sum(is.na(spectra(msidata))) - ## Number of NA in spectra matrix - infcount = sum(is.infinite(as.matrix(spectra(msidata)))) + infcount = sum(is.infinite(int_matrix)) ## Number of duplicated coordinates dupl_coord = sum(duplicated(coord(msidata))) properties = c("Number of m/z features", @@ -175,7 +180,7 @@ - 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 saved as .RData + - Cardinal "MSImageSet" or "MSImagingExperiment" saved as .RData ]]></token> <token name="@MZ_TABULAR_INPUT_DESCRIPTION@"><