diff egsea.R @ 1:73281fbdf6c1 draft

planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/egsea commit 225518a08941e7ef8e5c402e3696ec5fa6e592a0
author iuc
date Thu, 15 Feb 2018 02:34:59 -0500
parents a8a083193440
children ba2111ae6eb4
line wrap: on
line diff
--- a/egsea.R	Thu Jan 25 02:23:23 2018 -0500
+++ b/egsea.R	Thu Feb 15 02:34:59 2018 -0500
@@ -50,13 +50,14 @@
     make_option(c("-base_methods", "--base_methods"), type="character", help="Gene set testing methods"),
     make_option(c("-msigdb", "--msigdb"), type="character", help="MSigDB Gene Set Collections"),
     make_option(c("-keggdb", "--keggdb"), type="character", help="KEGG Pathways"),
+    make_option(c("-keggupdated", "--keggupdated"), type="logical", help="Use updated KEGG"),
     make_option(c("-gsdb", "--gsdb"), type="character", help = "GeneSetDB Gene Sets"),
     make_option(c("-display_top", "--display_top"), type="integer", help = "Number of top Gene Sets to display"),
     make_option(c("-min_size", "--min_size"), type="integer", help = "Minimum Size of Gene Set"),
     make_option(c("-fdr_cutoff", "--fdr_cutoff"), type="double", help = "FDR cutoff"),
     make_option(c("-combine_method", "--combine_method"), type="character", help="Method to use to combine the p-values"),
     make_option(c("-sort_method", "--sort_method"), type="character", help="Method to sort the results"),
-    make_option(c("-rdata", "--rdaOpt"), type="character", help="Output RData file")
+    make_option(c("-rdaOpt", "--rdaOpt"), type="character", help="Output RData file")
     )
 
 parser <- OptionParser(usage = "%prog [options] file", option_list=option_list)
@@ -188,10 +189,9 @@
     gsdb <- "none"
 }
 
-
 ## Index gene sets
 
-gs.annots <- buildIdx(entrezIDs=rownames(counts), species=args$species, msigdb.gsets=msigdb, gsdb.gsets=gsdb, kegg.exclude=kegg_exclude)
+gs.annots <- buildIdx(entrezIDs=rownames(counts), species=args$species, msigdb.gsets=msigdb, gsdb.gsets=gsdb, kegg.exclude=kegg_exclude, kegg.updated=args$keggupdated)
 
 
 ## Run egsea.cnt
@@ -201,6 +201,6 @@
 
 ## Output RData file
 
-if (!is.null(args$rdata)) {
+if (!is.null(args$rdaOpt)) {
   save.image(file = "EGSEA_analysis.RData")
 }
\ No newline at end of file