# HG changeset patch # User dktanwar # Date 1513616705 18000 # Node ID aa29b20bbb45def8bb016df6fe95eb043f06aea0 # Parent 609b5a2ae0e8b68f993be3b7dae6a06c587b6e0d Uploaded diff -r 609b5a2ae0e8 -r aa29b20bbb45 csaw.R --- a/csaw.R Mon Dec 18 11:20:20 2017 -0500 +++ b/csaw.R Mon Dec 18 12:05:05 2017 -0500 @@ -48,14 +48,19 @@ table <- data.frame(df, assay(windows), stringsAsFactors = F, check.names = F) colnames(table)[4:ncol(table)] <- file_names - # Remove spaces in the table ---- setDT(table) for (j in names(table)) set(table, j = j, value = table[[trimws(j)]]) -table_sp <- data.frame(table) # Save final table ---- -fwrite(x = table_sp, file = output, quote = F, row.names = F, sep = "\t") +# fwrite(x = table_sp, file = output, quote = F, row.names = F, sep = "\t") +dt <- table[,regions:=paste0(seqnames,"-", start, "-", end)] +table_sp <- data.frame(dt) + +for(i in 4:(ncol(table_sp)-1)){ + tmp <- table_sp[,c(ncol(table_sp), i)] + fwrite(x = tmp, file = output, quote = F, row.names = F, sep = "\t") +} # # Save individual files ---- # Sys.time()