# HG changeset patch
# User greg
# Date 1502455196 14400
# Node ID 50f64f662726a0d3bf3981ff7a540b7b2586fc6b
# Parent  54ec2ba55b8746fcd5a54d65fbc307f942f98f76
Uploaded

diff -r 54ec2ba55b87 -r 50f64f662726 build_matrix.R
--- a/build_matrix.R	Thu Aug 10 14:53:46 2017 -0400
+++ b/build_matrix.R	Fri Aug 11 08:39:56 2017 -0400
@@ -12,8 +12,9 @@
 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)==TRUE,];
-write.table(array(r_matrix, dim=c(length(r_matrix)/3, 3)), file=opt$output, quote=FALSE, row.names=FALSE, col.names=FALSE);
+data_table <- read.table(opt$input)
+as.matrix(data_table)
+status <- match(data_table[,3], missing)
+data_table[,3] <- paste(opt$work_dir, data_table[,1], ".", data_table[,2], ".bed.gz", sep="")
+data_table <- data_table[is.na(status)==TRUE,]
+write.table(array(data_table, dim=c(length(data_table)/3, 3)), file=opt$output, quote=FALSE, row.names=FALSE, col.names=FALSE)