diff goseq.xml @ 3:783e8b70b047 draft

planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/goseq commit 46c4278d292ab4d76dc5f3f74c3109c3179be7ef
author iuc
date Mon, 24 Sep 2018 06:29:03 -0400
parents ab492df30cdf
children ae39895af5fe
line wrap: on
line diff
--- a/goseq.xml	Mon Oct 23 11:19:12 2017 -0400
+++ b/goseq.xml	Mon Sep 24 06:29:03 2018 -0400
@@ -1,12 +1,14 @@
-<tool id="goseq" name="goseq" version="1.26.0">
+<tool id="goseq" name="goseq" version="1.30.0">
     <description>tests for overrepresented gene categories</description>
     <requirements>
-        <requirement type="package" version="1.3.2">r-optparse</requirement>
-        <requirement type="package" version="1.26.0">bioconductor-goseq</requirement>
-        <requirement type="package" version="3.3.0">bioconductor-org.hs.eg.db</requirement>
-        <requirement type="package" version="3.4.0">bioconductor-org.dm.eg.db</requirement>
-        <requirement type="package" version="3.4.1">bioconductor-org.dr.eg.db</requirement>
-        <requirement type="package" version="3.4.0">bioconductor-org.mm.eg.db</requirement>
+        <requirement type="package" version="1.6.0">r-optparse</requirement>
+        <requirement type="package" version="1.30.0">bioconductor-goseq</requirement>
+        <requirement type="package" version="3.5.0">bioconductor-org.hs.eg.db</requirement>
+        <requirement type="package" version="3.5.0">bioconductor-org.dm.eg.db</requirement>
+        <requirement type="package" version="3.5.0">bioconductor-org.dr.eg.db</requirement>
+        <requirement type="package" version="3.5.0">bioconductor-org.mm.eg.db</requirement>
+        <requirement type="package" version="0.7.6">r-dplyr</requirement>
+        <requirement type="package" version="3.0.0">r-ggplot2</requirement>
     </requirements>
     <stdio>
         <regex match="Execution halted"
@@ -23,7 +25,7 @@
                description="An undefined error occured, please check your input carefully and contact your administrator." />
     </stdio>
     <version_command><![CDATA[
-echo $(R --version | grep version | grep -v GNU)", goseq version" $(R --vanilla --slave -e "library(goseq); cat(sessionInfo()\$otherPkgs\$goseq\$Version)" 2> /dev/null | grep -v -i "WARNING: ")", optparse version" $(R --vanilla --slave -e "library(optparse); cat(sessionInfo()\$otherPkgs\$optparse\$Version)" 2> /dev/null | grep -v -i "WARNING: ")", org.Hs.eg.db version" $(R --vanilla --slave -e "library(org.Hs.eg.db); cat(sessionInfo()\$otherPkgs\$org.Hs.eg.db\$Version)" 2> /dev/null | grep -v -i "WARNING: ")", org.Dr.eg.db version" $(R --vanilla --slave -e "library(org.Dr.eg.db); cat(sessionInfo()\$otherPkgs\$org.Dr.eg.db\$Version)" 2> /dev/null | grep -v -i "WARNING: ")", org.Dm.eg.db version" $(R --vanilla --slave -e "library(org.Dm.eg.db); cat(sessionInfo()\$otherPkgs\$org.Dm.eg.db\$Version)" 2> /dev/null | grep -v -i "WARNING: ")", org.Mm.eg.db version" $(R --vanilla --slave -e "library(org.Mm.eg.db); cat(sessionInfo()\$otherPkgs\$org.Mm.eg.db\$Version)" 2> /dev/null | grep -v -i "WARNING: ")
+echo $(R --version | grep version | grep -v GNU)", goseq version" $(R --vanilla --slave -e "library(goseq); cat(sessionInfo()\$otherPkgs\$goseq\$Version)" 2> /dev/null | grep -v -i "WARNING: ")", optparse version" $(R --vanilla --slave -e "library(optparse); cat(sessionInfo()\$otherPkgs\$optparse\$Version)" 2> /dev/null | grep -v -i "WARNING: ")", org.Hs.eg.db version" $(R --vanilla --slave -e "library(org.Hs.eg.db); cat(sessionInfo()\$otherPkgs\$org.Hs.eg.db\$Version)" 2> /dev/null | grep -v -i "WARNING: ")", org.Dr.eg.db version" $(R --vanilla --slave -e "library(org.Dr.eg.db); cat(sessionInfo()\$otherPkgs\$org.Dr.eg.db\$Version)" 2> /dev/null | grep -v -i "WARNING: ")", org.Dm.eg.db version" $(R --vanilla --slave -e "library(org.Dm.eg.db); cat(sessionInfo()\$otherPkgs\$org.Dm.eg.db\$Version)" 2> /dev/null | grep -v -i "WARNING: ")", org.Mm.eg.db version" $(R --vanilla --slave -e "library(org.Mm.eg.db); cat(sessionInfo()\$otherPkgs\$org.Mm.eg.db\$Version)" 2> /dev/null | grep -v -i "WARNING: ")", dplyr version" $(R --vanilla --slave -e "library(dplyr); cat(sessionInfo()\$otherPkgs\$dplyr\$Version)" 2> /dev/null | grep -v -i "WARNING: ")", ggplot2 version" $(R --vanilla --slave -e "library(ggplot2); cat(sessionInfo()\$otherPkgs\$ggplot2\$Version)" 2> /dev/null | grep -v -i "WARNING: ")
     ]]></version_command>
     <command><![CDATA[
 Rscript '$__tool_directory__/goseq.r'
@@ -51,11 +53,15 @@
 --make_plots '$out.make_plots'
 --length_bias_plot '$length_bias_plot'
 --sample_vs_wallenius_plot '$sample_vs_wallenius_plot'
+
 --rdata '$out.rdata_out'
-
 --p_adj_method '$adv.p_adj_method'
 --use_genes_without_cat '$adv.use_genes_without_cat'
 
+#if $out.topgo_plot:
+    --top_plot '$out.topgo_plot'
+#end if
+
     ]]></command>
 
     <!-- Input Files-->
@@ -100,6 +106,7 @@
 
         <!-- Output Options -->
         <section name="out" title="Output Options">
+            <param name="topgo_plot" type="boolean" checked="false" label="Output Top GO terms plot?" help="Output a PDF plot of the Top 10 over-represented GO terms. Default: No" />
             <param name="make_plots" type="boolean" checked="false" label="Produce diagnostic plots?" help="This will produce the length bias (PWF) plot. If both sampling and wallenius methods are selected, it will also produce a plot comparing their p-values. These plots may help you compare the different p-value estimation methods that goseq can use. Default: No" />
             <param name="rdata_out" type="boolean" checked="false" label="Output RData file?" help="Output all the data used by R to construct the tables and plots, can be loaded into R. Default: No" />
         </section>
@@ -120,7 +127,7 @@
 
     <outputs>
         <data name="wallenius_tab" format="tabular" label="${tool.name} on ${on_string}: Ranked category list - Wallenius method">
-            <filter>methods['wallenius]'</filter>
+            <filter>methods['wallenius']</filter>
         </data>
         <data name="sampling_tab" format="tabular" label="${tool.name} on ${on_string}: Ranked category list - Sampling method">
             <filter>methods['repcnt'] != 0</filter>
@@ -139,9 +146,23 @@
         <data name="rdata" format="rdata" from_work_dir="goseq_analysis.RData" label="${tool.name} on ${on_string}: RData file">
             <filter>out['rdata_out']</filter>
         </data>
+        <data name="top_plot" format="pdf" from_work_dir="top10.pdf" label="${tool.name} on ${on_string}: Top over-represented GO terms plot">
+            <filter>methods['wallenius']</filter>
+            <filter>out['topgo_plot']</filter>
+        </data>
     </outputs>
 
     <tests>
+        <!-- Ensure top plot is output -->
+        <test expect_num_outputs="2">
+            <param name="dge_file" value="dge_list.tab" ftype="tabular" />
+            <param name="length_file" value="gene_length.tab" ftype="tabular" />
+            <param name="catSource" value="history" />
+            <param name="category_file" value="category.tab" ftype="tabular" />
+            <param name="use_genes_without_cat" value="true" />
+            <param name="topgo_plot" value="true" />
+            <output name="top_plot" ftype="pdf" file="topgo.pdf" compare="sim_size"/>
+        </test>
         <!-- Ensure Wallenius table is output -->
         <test expect_num_outputs="1">
             <param name="dge_file" value="dge_list.tab" ftype="tabular" />
@@ -149,7 +170,12 @@
             <param name="catSource" value="history" />
             <param name="category_file" value="category.tab" ftype="tabular" />
             <param name="use_genes_without_cat" value="true" />
-            <output name="wallenius_tab" file="wal.tab" compare="contains" />
+            <output name="wallenius_tab">
+                <assert_contents>
+                    <has_text_matching expression="category.*over_represented_pvalue.*under_represented_pvalue.*numDEInCat.*numInCat.*term.*ontology.*p.adjust.over_represented.*p.adjust.under_represented" />
+                    <has_text_matching expression="GO:0000278.*0.01" />
+                </assert_contents>
+            </output>
         </test>
         <!-- Ensure getting GO categories works -->
         <test expect_num_outputs="1">
@@ -159,7 +185,12 @@
             <param name="genome" value="hg38" />
             <param name="gene_id" value="ensGene" />
             <param name="use_genes_without_cat" value="true" />
-            <output name="wallenius_tab" ftype="tabular" file="getgo.hg38.tab" compare="contains"/>
+            <output name="wallenius_tab">
+                <assert_contents>
+                    <has_text_matching expression="category.*over_represented_pvalue.*under_represented_pvalue.*numDEInCat.*numInCat.*term.*ontology.*p.adjust.over_represented.*p.adjust.under_represented" />
+                    <has_text_matching expression="GO:0005576.*5.2" />
+                </assert_contents>
+            </output>
         </test>
         <!-- Ensure getting GO categories for another genome (zebrafish) works -->
         <test expect_num_outputs="1">
@@ -169,7 +200,12 @@
             <param name="genome" value="danRer10"/>
             <param name="gene_id" value="ensGene" />
             <param name="use_genes_without_cat" value="true" />
-            <output name="wallenius_tab" ftype="tabular" file="getgo.danRer10.tab" compare="contains"/>
+            <output name="wallenius_tab">
+                <assert_contents>
+                    <has_text_matching expression="category.*over_represented_pvalue.*under_represented_pvalue.*numDEInCat.*numInCat.*term.*ontology.*p.adjust.over_represented.*p.adjust.under_represented" />
+                    <has_text_matching expression="GO:0031324.*0.08" />
+                </assert_contents>
+            </output>
         </test>
         <!-- Ensure length bias plot works -->
         <test expect_num_outputs="2">
@@ -277,7 +313,7 @@
 **Outputs**
 
 * This tool outputs a tabular file containing a ranked list of gene categories, similar to below. The default output is the Wallenius method table. If the Sampling and/or Hypergeometric methods are also selected, additional tables are produced.
-* Optionally, this tool can also output some diagnostic plots and an RData file, see **Output Options** above.
+* Optionally, this tool can also output a plot of the top 10 over-represented categories, some diagnostic plots and an RData file, see **Output Options** above.
 
 Example: