Mercurial > repos > greg > ideas
changeset 22:f8e36f4b71c7 draft
Uploaded
| author | greg |
|---|---|
| date | Fri, 11 Aug 2017 09:39:53 -0400 |
| parents | 7e22d880e551 |
| children | 64c372086d41 |
| files | build_matrix.R |
| diffstat | 1 files changed, 6 insertions(+), 12 deletions(-) [+] |
line wrap: on
line diff
--- a/build_matrix.R Fri Aug 11 09:29:18 2017 -0400 +++ b/build_matrix.R Fri Aug 11 09:39:53 2017 -0400 @@ -12,15 +12,9 @@ args <- parse_args(parser, positional_arguments=TRUE) opt <- args$options -#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) - -x=as.matrix(read.table(opt$input)) -t=match(x[,3],missing) -x[,3]=paste(opt$work_dir, x[,1],".",x[,2],".bed.gz", sep="") -x=x[is.na(t)==T,] -write.table(array(x,dim=c(length(x)/3,3)), opt$output, quote=F,row.names=F,col.names=F) +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)