# HG changeset patch # User iuc # Date 1548695016 18000 # Node ID 4f497a8a72e86280c764e66a4de0420632df0f9e # Parent 13f324f9cdb7cb8178b63dd4777a9be657126b88 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 05bc13637dc9f8e523bc72844ff6eb0794f98ad3 diff -r 13f324f9cdb7 -r 4f497a8a72e8 deseq2.R --- a/deseq2.R Thu Dec 06 15:48:59 2018 -0500 +++ b/deseq2.R Mon Jan 28 12:03:36 2019 -0500 @@ -46,7 +46,8 @@ spec <- matrix(c( "quiet", "q", 0, "logical", "help", "h", 0, "logical", - "batch_factors", "", 1, "character", + "cores", "s", 0, "integer", + "batch_factors", "w", 1, "character", "outfile", "o", 1, "character", "countsfile", "n", 1, "character", "rlogfile", "r", 1, "character", @@ -105,6 +106,14 @@ library("gplots") }) +if (opt$cores > 1) { + library("BiocParallel") + register(MulticoreParam(opt$cores)) + parallel = TRUE +} else { + parallel = FALSE +} + # build or read sample table trim <- function (x) gsub("^\\s+|\\s+$", "", x) @@ -273,7 +282,7 @@ if (verbose) cat(paste("using disperion fit type:",fitType,"\n")) # run the analysis -dds <- DESeq(dds, fitType=fitType, betaPrior=betaPrior, minReplicatesForReplace=minRep) +dds <- DESeq(dds, fitType=fitType, betaPrior=betaPrior, minReplicatesForReplace=minRep, parallel=parallel) # create the generic plots and leave the device open if (!is.null(opt$plots)) { diff -r 13f324f9cdb7 -r 4f497a8a72e8 deseq2.xml --- a/deseq2.xml Thu Dec 06 15:48:59 2018 -0500 +++ b/deseq2.xml Mon Jan 28 12:03:36 2019 -0500 @@ -1,4 +1,4 @@ - + Determines differentially expressed features from count tables bioconductor-deseq2 @@ -41,6 +41,7 @@ #import json #import os Rscript '${__tool_directory__}/deseq2.R' + --cores \${GALAXY_SLOTS:-1} -o '$deseq_out' #if $pdf: -p '$plots'