# HG changeset patch
# User mvdbeek
# Date 1456504882 18000
# Node ID 81283f3d65c77f1aa77677a1fa653591a7761928
# Parent  0283381bdafdd30854c3f9422af0b0ea3e166bfc
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/goseq_1_22_0 commit fdd0811efc61c31f88ff17096fbe8ee8cfacd766-dirty

diff -r 0283381bdafd -r 81283f3d65c7 goseq.r
--- a/goseq.r	Fri Feb 26 08:49:05 2016 -0500
+++ b/goseq.r	Fri Feb 26 11:41:22 2016 -0500
@@ -1,6 +1,6 @@
 sink(stdout(), type = "message")
-library(goseq)
-library(optparse)
+suppressWarnings(suppressMessages(library(goseq)))
+suppressWarnings(suppressMessages(library(optparse)))
 
 option_list <- list(
     make_option(c("-d", "--dge_file"), type="character", help="Path to file with differential gene expression result"),
@@ -34,15 +34,13 @@
 sample_vs_wallenius_plot = args$sample_vs_wallenius_plot
 repcnt = args$repcnt
 
-
 # format DE genes into vector suitable for use with goseq
 dge_table = read.delim(dge_file, header = TRUE, sep="\t", check.names = FALSE)
 genes = as.integer(dge_table[,p_adj_column]<p_adj_cutoff)
 names(genes) = dge_table[,1] # Assuming first row contains gene names
 
 # Get gene lengths
-
-if (length_file) {
+if (length_file !=FALSE) {
   length_table = read.delim(length_file, header=TRUE, sep="\t", check.names=FALSE, row.names=TRUE)
   gene_lengths = length_table[names(genes),]$length
   } else {
@@ -53,7 +51,7 @@
 
 pdf(length_bias_plot)
 pwf=nullp(genes, genome, gene_id, gene_lengths)
-dev.off()
+message = dev.off()
 # wallenius approximation of p-values
 GO.wall=goseq(pwf, genome, gene_id)
 
@@ -68,7 +66,7 @@
      xlab="log10(Wallenius p-values)",ylab="log10(Sampling p-values)",
      xlim=c(-3,0))
      abline(0,1,col=3,lty=2)
-  dev.off()
+  message = dev.off()
   write.table(GO.samp, sampling_tab, sep="\t", row.names = FALSE, quote = FALSE)
 }
 
diff -r 0283381bdafd -r 81283f3d65c7 goseq.xml
--- a/goseq.xml	Fri Feb 26 08:49:05 2016 -0500
+++ b/goseq.xml	Fri Feb 26 11:41:22 2016 -0500
@@ -48,11 +48,11 @@
         <data format="pdf" label="Plot P-value from sampling against wallenius distribution" name="sample_vs_wallenius_plot">
             <filter>repcnt != 0</filter>
         </data>
-        <data format="tabular" label="Ranked category list - wallenius" name="wallenius_tab" />
+        <data format="tabular" label="Ranked category list - no length bias correction" name="nobias_tab" />
         <data format="tabular" label="Ranked category list - sampling" name="sampling_tab">
             <filter>repcnt != 0</filter>
         </data>
-        <data format="tabular" label="Ranked category list - no length bias correction" name="nobias_tab" />
+        <data format="tabular" label="Ranked category list - wallenius approx. of p-values" name="wallenius_tab" />
     </outputs>
     <tests>
         <test>