Mercurial > repos > greg > ideas
changeset 2:7a77ec6342c1 draft
Uploaded
author | greg |
---|---|
date | Thu, 10 Aug 2017 12:34:42 -0400 |
parents | d75cbb2db2c4 |
children | 9874800487e6 |
files | build_matrix.R |
diffstat | 1 files changed, 20 insertions(+), 0 deletions(-) [+] |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/build_matrix.R Thu Aug 10 12:34:42 2017 -0400 @@ -0,0 +1,20 @@ +#!/usr/bin/env Rscript + +suppressPackageStartupMessages(library("optparse")) + +option_list <- list( + make_option(c("-i", "--input"), action="store", dest="input", help="Input .bed.gz file produced by prepMat"), + make_option(c("-o", "--output"), action="store", dest="output", help="Output file"), + make_option(c("-w", "--work_dir"), action="store", dest="work_dir", help="Working directory"), +) + +parser <- OptionParser(usage="%prog [options] file", option_list=option_list) +args <- parse_args(parser, positional_arguments=TRUE) +opt <- args$options + +input_dir = paste(opt$work_dir, "Input/", sep=""); +r_matrix = as.matrix(read.table(opt$input)); +status = match(r_matrix[,3], missing); +r_matrix[,3] = paste(input_dir, r_matrix[,1], ".", r_matrix[,2], ".bed.gz", sep=""); +r_matrix = r_matrix[is.na(status)==T,]; +write.table(array(r_matrix, dim = c(length(r_matrix) / 3, 3)), opt$output, quote=F, row.names=F, col.names=F);