annotate deseq_normalization.xml @ 0:557c3cdf5b0a draft default tip

planemo upload for repository https://bitbucket.org/drosofff/gedtools/
author mvdbeek
date Mon, 29 Jun 2015 05:50:30 -0400
parents
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
557c3cdf5b0a planemo upload for repository https://bitbucket.org/drosofff/gedtools/
mvdbeek
parents:
diff changeset
1 <tool id="DESeq_normalization" name="DESeq Normalization" version="1.0.2">
557c3cdf5b0a planemo upload for repository https://bitbucket.org/drosofff/gedtools/
mvdbeek
parents:
diff changeset
2 <description>of hit lists (geometric method)</description>
557c3cdf5b0a planemo upload for repository https://bitbucket.org/drosofff/gedtools/
mvdbeek
parents:
diff changeset
3 <requirements>
557c3cdf5b0a planemo upload for repository https://bitbucket.org/drosofff/gedtools/
mvdbeek
parents:
diff changeset
4 <requirement type="package" version="3.1.2">R</requirement>
557c3cdf5b0a planemo upload for repository https://bitbucket.org/drosofff/gedtools/
mvdbeek
parents:
diff changeset
5 <requirement type="package" version="2.14">biocbasics</requirement>
557c3cdf5b0a planemo upload for repository https://bitbucket.org/drosofff/gedtools/
mvdbeek
parents:
diff changeset
6 </requirements>
557c3cdf5b0a planemo upload for repository https://bitbucket.org/drosofff/gedtools/
mvdbeek
parents:
diff changeset
7 <command interpreter="Rscript"> $normalizationDESeq </command>
557c3cdf5b0a planemo upload for repository https://bitbucket.org/drosofff/gedtools/
mvdbeek
parents:
diff changeset
8 <inputs>
557c3cdf5b0a planemo upload for repository https://bitbucket.org/drosofff/gedtools/
mvdbeek
parents:
diff changeset
9 <param name="input" type="data" format="tabular" label="Hit lists of items"/>
557c3cdf5b0a planemo upload for repository https://bitbucket.org/drosofff/gedtools/
mvdbeek
parents:
diff changeset
10 </inputs>
557c3cdf5b0a planemo upload for repository https://bitbucket.org/drosofff/gedtools/
mvdbeek
parents:
diff changeset
11 <outputs>
557c3cdf5b0a planemo upload for repository https://bitbucket.org/drosofff/gedtools/
mvdbeek
parents:
diff changeset
12 <data name="output_factors" format="tabular" label="Normalization Factors (Geometric mean method)" />
557c3cdf5b0a planemo upload for repository https://bitbucket.org/drosofff/gedtools/
mvdbeek
parents:
diff changeset
13 <data name="norm_hit_table" format="tabular" label="Normalized Hit Table (Geometric mean method)" />
557c3cdf5b0a planemo upload for repository https://bitbucket.org/drosofff/gedtools/
mvdbeek
parents:
diff changeset
14 </outputs>
557c3cdf5b0a planemo upload for repository https://bitbucket.org/drosofff/gedtools/
mvdbeek
parents:
diff changeset
15 <tests>
557c3cdf5b0a planemo upload for repository https://bitbucket.org/drosofff/gedtools/
mvdbeek
parents:
diff changeset
16 <test>
557c3cdf5b0a planemo upload for repository https://bitbucket.org/drosofff/gedtools/
mvdbeek
parents:
diff changeset
17 <param name="input" value="counts.tab" ftype="tabular"/>
557c3cdf5b0a planemo upload for repository https://bitbucket.org/drosofff/gedtools/
mvdbeek
parents:
diff changeset
18 <output name="output_factors" file="norm_factors.tab" ftype="tabular"/>
557c3cdf5b0a planemo upload for repository https://bitbucket.org/drosofff/gedtools/
mvdbeek
parents:
diff changeset
19 <output name="norm_hit_table" file="norm_table.tab" ftype="tabular"/>
557c3cdf5b0a planemo upload for repository https://bitbucket.org/drosofff/gedtools/
mvdbeek
parents:
diff changeset
20 </test>
557c3cdf5b0a planemo upload for repository https://bitbucket.org/drosofff/gedtools/
mvdbeek
parents:
diff changeset
21 </tests>
557c3cdf5b0a planemo upload for repository https://bitbucket.org/drosofff/gedtools/
mvdbeek
parents:
diff changeset
22 <configfiles>
557c3cdf5b0a planemo upload for repository https://bitbucket.org/drosofff/gedtools/
mvdbeek
parents:
diff changeset
23 <configfile name="normalizationDESeq">
557c3cdf5b0a planemo upload for repository https://bitbucket.org/drosofff/gedtools/
mvdbeek
parents:
diff changeset
24 ## Setup R error handling to go to stderr
557c3cdf5b0a planemo upload for repository https://bitbucket.org/drosofff/gedtools/
mvdbeek
parents:
diff changeset
25 options( show.error.messages=F, error = function () { cat( geterrmessage(), file=stderr() ); q( "no", 1, F ) } )
557c3cdf5b0a planemo upload for repository https://bitbucket.org/drosofff/gedtools/
mvdbeek
parents:
diff changeset
26 suppressMessages(require(DESeq))
557c3cdf5b0a planemo upload for repository https://bitbucket.org/drosofff/gedtools/
mvdbeek
parents:
diff changeset
27 countsTable = read.delim("${input}", header=TRUE, check.names=FALSE)
557c3cdf5b0a planemo upload for repository https://bitbucket.org/drosofff/gedtools/
mvdbeek
parents:
diff changeset
28 rownames( countsTable )= countsTable[,1]
557c3cdf5b0a planemo upload for repository https://bitbucket.org/drosofff/gedtools/
mvdbeek
parents:
diff changeset
29 countsTable= countsTable[ , -1 ]
557c3cdf5b0a planemo upload for repository https://bitbucket.org/drosofff/gedtools/
mvdbeek
parents:
diff changeset
30 conds = c(rep ("C", length(countsTable[1,])-1 ), "T")
557c3cdf5b0a planemo upload for repository https://bitbucket.org/drosofff/gedtools/
mvdbeek
parents:
diff changeset
31 cds = newCountDataSet( countsTable, conds )
557c3cdf5b0a planemo upload for repository https://bitbucket.org/drosofff/gedtools/
mvdbeek
parents:
diff changeset
32 cds = estimateSizeFactors( cds )
557c3cdf5b0a planemo upload for repository https://bitbucket.org/drosofff/gedtools/
mvdbeek
parents:
diff changeset
33 write.table (as.data.frame(sizeFactors(cds)), file="${output_factors}", row.names=TRUE, col.names=TRUE, quote= FALSE, dec = ".", sep = "\t", eol = "\n")
557c3cdf5b0a planemo upload for repository https://bitbucket.org/drosofff/gedtools/
mvdbeek
parents:
diff changeset
34 norm_table = as.data.frame(counts( cds, normalized=TRUE ))
557c3cdf5b0a planemo upload for repository https://bitbucket.org/drosofff/gedtools/
mvdbeek
parents:
diff changeset
35 norm_table = cbind(rownames(norm_table), norm_table)
557c3cdf5b0a planemo upload for repository https://bitbucket.org/drosofff/gedtools/
mvdbeek
parents:
diff changeset
36 colnames(norm_table) = c("gene", colnames(countsTable) )
557c3cdf5b0a planemo upload for repository https://bitbucket.org/drosofff/gedtools/
mvdbeek
parents:
diff changeset
37 write.table (norm_table, file="${norm_hit_table}", row.names=FALSE, col.names=TRUE, quote= FALSE, dec = ".", sep = "\t", eol = "\n")
557c3cdf5b0a planemo upload for repository https://bitbucket.org/drosofff/gedtools/
mvdbeek
parents:
diff changeset
38 </configfile>
557c3cdf5b0a planemo upload for repository https://bitbucket.org/drosofff/gedtools/
mvdbeek
parents:
diff changeset
39 </configfiles>
557c3cdf5b0a planemo upload for repository https://bitbucket.org/drosofff/gedtools/
mvdbeek
parents:
diff changeset
40 <help>
557c3cdf5b0a planemo upload for repository https://bitbucket.org/drosofff/gedtools/
mvdbeek
parents:
diff changeset
41
557c3cdf5b0a planemo upload for repository https://bitbucket.org/drosofff/gedtools/
mvdbeek
parents:
diff changeset
42 **What it does**
557c3cdf5b0a planemo upload for repository https://bitbucket.org/drosofff/gedtools/
mvdbeek
parents:
diff changeset
43
557c3cdf5b0a planemo upload for repository https://bitbucket.org/drosofff/gedtools/
mvdbeek
parents:
diff changeset
44 DESeq Normalization Factors (geometic mean method)
557c3cdf5b0a planemo upload for repository https://bitbucket.org/drosofff/gedtools/
mvdbeek
parents:
diff changeset
45 and
557c3cdf5b0a planemo upload for repository https://bitbucket.org/drosofff/gedtools/
mvdbeek
parents:
diff changeset
46 Normalized hit lists
557c3cdf5b0a planemo upload for repository https://bitbucket.org/drosofff/gedtools/
mvdbeek
parents:
diff changeset
47
557c3cdf5b0a planemo upload for repository https://bitbucket.org/drosofff/gedtools/
mvdbeek
parents:
diff changeset
48 </help>
557c3cdf5b0a planemo upload for repository https://bitbucket.org/drosofff/gedtools/
mvdbeek
parents:
diff changeset
49 </tool>