Mercurial > repos > hackdna > pathprint
diff consensus.r @ 0:0cebe436a553 draft default tip
Initial upload.
author | hackdna |
---|---|
date | Fri, 17 May 2013 14:29:33 -0400 |
parents | |
children |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/consensus.r Fri May 17 14:29:33 2013 -0400 @@ -0,0 +1,32 @@ +#!/usr/bin/env Rscript + +#sink(file("/dev/null", "w"), type = "message"); + +library(pathprint); + +# Figure out the relative path to the galaxy-pathprint.r library. +script.args <- commandArgs(trailingOnly = FALSE); +script.name <- sub("--file=", "", script.args[grep("--file=", script.args)]) +script.base <- dirname(script.name) +library.path <- file.path(script.base, "galaxy-pathprint.r"); +source(library.path) + +data(GEO.metadata.matrix); + +usage <- function() { + stop("Usage: consensus.r <input> <threshold> <output>", call. = FALSE) +} + +## Get the command line arguments. +args <- commandArgs(trailingOnly = TRUE) + +input <- ifelse(! is.na(args[1]), args[1], usage()) +threshold <- ifelse(! is.na(args[2]), as.numeric(args[2]), usage()) +output <- ifelse(! is.na(args[3]), args[3], usage()) + +fingerprint <- loadFingerprint(input); +consensus <- consensusFingerprint(fingerprint, threshold); + +saveConsensus(consensus, output); + +quit("no", 0)