Mercurial > repos > iuc > heinz
changeset 3:2b80a2596064 draft default tip
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/heinz commit 8007f71281553ddfa45e6f8e1172952d956bb000"
author | iuc |
---|---|
date | Thu, 11 Jun 2020 07:18:17 -0400 |
parents | 704c401e0afb |
children | |
files | bum.R |
diffstat | 1 files changed, 14 insertions(+), 14 deletions(-) [+] |
line wrap: on
line diff
--- a/bum.R Thu Jan 23 14:07:54 2020 -0500 +++ b/bum.R Thu Jun 11 07:18:17 2020 -0400 @@ -2,15 +2,16 @@ # Usage: Rscript bum.R --input p-values.txt --output result.txt --verbose TRUE # Set up R error handling to go to stderr -options(show.error.messages=F, error=function(){cat(geterrmessage(),file=stderr());q("no",1,F)}) - -# Avoid crashing Galaxy with an UTF8 error on German LC settings -#loc <- Sys.setlocale("LC_MESSAGES", "en_US.UTF-8") +err_foo <- function() { + cat(geterrmessage(), file = stderr()); + q("no", 1, F) +} +options(show.error.messages = F, error = err_foo) # Import required libraries suppressPackageStartupMessages({ - library('getopt') - library('BioNet') + library("getopt") + library("BioNet") }) # Take in trailing command line arguments @@ -18,17 +19,16 @@ # Get options using the spec as defined by the enclosed list # Read the options from the default: commandArgs(TRUE) -option_specification <- matrix(c( - 'input', 'i', 2, 'character', - 'output', 'o', 2, 'character' -), byrow=TRUE, ncol=4); +option_specification <- matrix(c("input", "i", 2, "character", + "output", "o", 2, "character" + ), byrow = TRUE, ncol = 4); # Parse options options <- getopt(option_specification); pvals <- read.table(options$input) -bum <- fitBumModel(pvals,plot=FALSE) +bum <- fitBumModel(pvals, plot = FALSE) mat <- c(bum$lambda, bum$a) -#bumtablename <- paste(options$output,sep="\t") -write.table(x=mat, file=options$output,quote=FALSE, row.names=FALSE, col.names=FALSE) -message ("Parameters have been estimated successfully!") +write.table(x = mat, file = options$output, quote = FALSE, + row.names = FALSE, col.names = FALSE) +message("Parameters have been estimated successfully!")