Mercurial > repos > galaxyp > cardinal_segmentations
changeset 9:c68bbd8414d8 draft
"planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/cardinal commit e499c9124d3fd85a7fc47b95c206ce91a5e3678c-dirty"
line wrap: on
line diff
--- a/macros.xml Wed May 13 17:52:58 2020 +0000 +++ b/macros.xml Tue Nov 03 22:30:47 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@"><