Mercurial > repos > galaxyp > cardinal_spectra_plots
changeset 12:10566783d718 draft
"planemo upload for repository https://github.com/galaxyproteomics/tools-galaxyp/tree/master/tools/cardinal commit e499c9124d3fd85a7fc47b95c206ce91a5e3678c-dirty"
| author | galaxyp |
|---|---|
| date | Tue, 03 Nov 2020 23:11:11 +0000 |
| parents | d2596e1be091 |
| children | c5aa9fca18c6 |
| files | spectra_plots.xml |
| diffstat | 1 files changed, 8 insertions(+), 8 deletions(-) [+] |
line wrap: on
line diff
--- a/spectra_plots.xml Tue Oct 06 08:09:08 2020 +0000 +++ b/spectra_plots.xml Tue Nov 03 23:11:11 2020 +0000 @@ -1,4 +1,4 @@ -<tool id="cardinal_spectra_plots" name="MSI plot spectra" version="@VERSION@.0"> +<tool id="cardinal_spectra_plots" name="MSI plot spectra" version="@VERSION@.1"> <description> mass spectrometry imaging mass spectra plots </description> @@ -82,7 +82,7 @@ merged_annotation = merge(msidata_coordinates, spectra_input, by=c("x", "y"), all.x=TRUE) merged_annotation[is.na(merged_annotation)] = "NA" merged_annotation = merged_annotation[order(merged_annotation\$pixel_index),] - msidata\$annotation = as.factor(merged_annotation[,4]) + msidata\$annotation = factor(merged_annotation[,4], levels = unique(as.character(merged_annotation[,4]))) ## keep the right order ## overview plot over annotated samples number_combined = length(levels(msidata\$annotation)) @@ -214,7 +214,7 @@ #if $fullmz: ## plot full mz for single tabular mz - print(plot(msidata, coord=list(x=x_coord, y=y_coord), key=TRUE, grid = $grid_variable)) + print(plot(msidata, coord=list(x=x_coord, y=y_coord), key=TRUE, grid = $grid_variable, col="black")) #end if pixelname = paste0("x = ", x_coord,", ", "y = ", y_coord) @@ -261,7 +261,7 @@ } ## print single tabular mz; manual zoom - print(plot(msidata[minmasspixel:maxmasspixel,], grid = $grid_variable, coord=list(x=x_coord, y=y_coord), key=TRUE, xlim= c($token.xlimmin,$token.xlimmax))) + print(plot(msidata[minmasspixel:maxmasspixel,], grid = $grid_variable, coord=list(x=x_coord, y=y_coord), key=TRUE, xlim= c($token.xlimmin,$token.xlimmax), col="black")) #end for @@ -306,7 +306,7 @@ } ## print single tabular mz; tabular zoom - print(plot(msidata[minmasspixel:maxmasspixel,], grid = $grid_variable, coord=list(x=x_coord, y=y_coord), key=TRUE, xlim= c(xlimmin,xlimmax))) + print(plot(msidata[minmasspixel:maxmasspixel,], grid = $grid_variable, coord=list(x=x_coord, y=y_coord), key=TRUE, xlim= c(xlimmin,xlimmax), col="black")) } } @@ -329,7 +329,7 @@ #if $fullmz: - print(plot(msidata, run="infile", key=TRUE, strip=FALSE, main="Average spectrum", grid = $grid_variable)) + print(plot(msidata, run="infile", key=TRUE, strip=FALSE, main="Average spectrum", grid = $grid_variable, col="black")) #end if ##################### II) Sample: plot zoom-in mass spectrum ########## @@ -366,7 +366,7 @@ } } - print(plot(msidata[minmasspixel:maxmasspixel,], grid = $grid_variable, main="Average spectrum", run = "infile", strip=FALSE, key=TRUE, xlim= c($token.xlimmin,$token.xlimmax))) + print(plot(msidata[minmasspixel:maxmasspixel,], grid = $grid_variable, main="Average spectrum", run = "infile", strip=FALSE, key=TRUE, xlim= c($token.xlimmin,$token.xlimmax), col="black")) #end for @@ -401,7 +401,7 @@ } ) - print(plot(msidata[minmasspixel:maxmasspixel,], grid = $grid_variable, run="infile", key=TRUE, strip=FALSE, main="Average spectrum", xlim= c(xlimmin,xlimmax))) + print(plot(msidata[minmasspixel:maxmasspixel,], grid = $grid_variable, run="infile", key=TRUE, strip=FALSE, main="Average spectrum", xlim= c(xlimmin,xlimmax), col="black")) } }