Mercurial > repos > artbio > deseq2_normalization

changeset 2:8faa01102fbc draft default tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/main/tools/deseq2_normalization commit 770f615298b9314baec54a539adcb2c3d0e729a7
author artbio
date Sat, 14 Oct 2023 18:00:34 +0000
parents ca845e480c89
children
files deseq2_normalization.xml
diffstat 1 files changed, 16 insertions(+), 10 deletions(-) [+]
line wrap: on
line diff
--- a/deseq2_normalization.xml	Mon Jan 08 04:32:10 2018 -0500
+++ b/deseq2_normalization.xml	Sat Oct 14 18:00:34 2023 +0000
@@ -1,7 +1,7 @@
-<tool id="deseq2_normalization" name="DESeq2 normalization" version="1.0.1">
+<tool id="deseq2_normalization" name="DESeq2 normalization" version="1.40.2+galaxy0" profile="22.01">
 	<description>of hit lists (geometric method)</description>
     <requirements>
-      <requirement type="package" version="1.8.2gx=r3.2.1_0">bioconductor-deseq2</requirement>
+      <requirement type="package" version="1.40.2">bioconductor-deseq2</requirement>
     </requirements>
     <stdio>
         <exit_code range="1:" level="fatal" description="Tool exception" />
@@ -10,19 +10,25 @@
     <configfiles>
         <configfile name="normalizationDESeq">
       ## Setup R error handling to go to stderr
-      options( show.error.messages=F,  error = function () { cat( geterrmessage(), file=stderr() ); q( "no", 1, F ) } )
+      options(show.error.messages =FALSE, error = function(){
+          cat(geterrmessage(), file=stderr())
+          q("no", 1, FALSE)
+          }
+      )
       suppressMessages(require(DESeq2))
-      countsTable = read.delim("${input}", header=TRUE, check.names=FALSE)
-      rownames( countsTable )= countsTable[,1]
-      countsTable= countsTable[ , -1 ]
-      conds = factor(c(rep ("C", length(countsTable[1,])-1 ), "T"))
+      countsTable = read.delim("${input}", header = TRUE, check.names = FALSE)
+      rownames(countsTable ) = countsTable[,1]
+      countsTable = countsTable[, -1]
+      conds = factor(c(rep("C", length(countsTable[1, ]) - 1), "T"))
       cds = DESeqDataSetFromMatrix(countsTable, DataFrame(conds), ~ conds)
       cds = estimateSizeFactors( cds )
-      write.table (as.data.frame(sizeFactors(cds)), file="${output_factors}", row.names=TRUE, col.names=TRUE, quote= FALSE, dec = ".", sep = "\t", eol = "\n")
+      write.table(as.data.frame(sizeFactors(cds)), file = "${output_factors}", row.names = TRUE,
+                  col.names = TRUE, quote = FALSE, dec = ".", sep = "\t", eol = "\n")
       norm_table = as.data.frame(counts( cds, normalized=TRUE ))
       norm_table = cbind(rownames(norm_table), norm_table)
-      colnames(norm_table) = c("gene", colnames(countsTable) )
-      write.table (norm_table, file="${norm_hit_table}", row.names=FALSE, col.names=TRUE, quote= FALSE, dec = ".", sep = "\t", eol = "\n")
+      colnames(norm_table) = c("gene", colnames(countsTable))
+      write.table(norm_table, file = "${norm_hit_table}", row.names = FALSE, col.names = TRUE,
+                  quote = FALSE, dec = ".", sep = "\t", eol = "\n")
         </configfile>
     </configfiles>
 	<inputs>