Mercurial > repos > iuc > egsea
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