Mercurial > repos > greg > ideas
changeset 19:50f64f662726 draft
Uploaded
| author | greg |
|---|---|
| date | Fri, 11 Aug 2017 08:39:56 -0400 |
| parents | 54ec2ba55b87 |
| children | a70690dcf9ff |
| files | build_matrix.R |
| diffstat | 1 files changed, 6 insertions(+), 5 deletions(-) [+] |
line wrap: on
line diff
--- 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)