# HG changeset patch
# User greg
# Date 1502390495 14400
# Node ID 79acab9c6524c96994c72a1021834a1fb71a9d61
# Parent  65739ddd00bebc04db951298707d5370b19eeb38
Uploaded

diff -r 65739ddd00be -r 79acab9c6524 build_matrix.R
--- a/build_matrix.R	Thu Aug 10 14:23:32 2017 -0400
+++ b/build_matrix.R	Thu Aug 10 14:41:35 2017 -0400
@@ -5,15 +5,15 @@
 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"),
+    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
 
-r_matrix = as.matrix(read.table(opt$input));
-status = match(r_matrix[,3], missing); 
-r_matrix[,3] = paste(opt$work_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);
+r_matrix <- as.matrix(read.table(opt$input));
+status <- match(r_matrix[,3], missing); 
+r_matrix[,3] <- paste(opt$work_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);