# HG changeset patch # User greg # Date 1502388558 14400 # Node ID c59cbd44dcff69728b3a6de709083fac4a0025ae # Parent 5b5733e0934436888571ba1ffa4a421486a6bcdd Uploaded diff -r 5b5733e09344 -r c59cbd44dcff build_matrix.R --- a/build_matrix.R Thu Aug 10 14:09:11 2017 -0400 +++ b/build_matrix.R Thu Aug 10 14:09:18 2017 -0400 @@ -5,16 +5,14 @@ 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[,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);