Mercurial > repos > greg > ideas
changeset 21:7e22d880e551 draft
Uploaded
| author | greg |
|---|---|
| date | Fri, 11 Aug 2017 09:29:18 -0400 |
| parents | a70690dcf9ff |
| children | f8e36f4b71c7 |
| files | build_matrix.R |
| diffstat | 1 files changed, 12 insertions(+), 6 deletions(-) [+] |
line wrap: on
line diff
--- a/build_matrix.R Fri Aug 11 09:11:40 2017 -0400 +++ b/build_matrix.R Fri Aug 11 09:29:18 2017 -0400 @@ -12,9 +12,15 @@ 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) +#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)