Merge pull request #41 from p-smirnov/plot-sv-calls_dev_p-smirnov1

Fixing small indexing errors in plot-sv-calls_dev
friendsofstrandseq · Aug 9, 2023 · 74419f6 · 74419f6
2 parents 30e7535 + d54ec74
commit 74419f6
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/workflow/scripts/plotting/plot-sv-calls_dev.R b/workflow/scripts/plotting/plot-sv-calls_dev.R
@@ -363,7 +363,7 @@ if (!is.null(f_segments)) {
 
 
   seg <-
-    merge(seg, bins[, .N, by = chrom][, .(chrom, N = c(0, cumsum(N))[1:(.N - 1)])], by = "chrom")
+    merge(seg, bins[, .N, by = chrom][, .(chrom, N = c(0, cumsum(N)[- .N]))], by = "chrom")
 
   # print(c(1, bps[1:(.N - 1)] + 1))
   # print(bps)
@@ -372,7 +372,7 @@ if (!is.null(f_segments)) {
   # stop()
 
 
-  seg[, `:=`(from = c(1, bps[1:(.N - 1)] + 1), to = bps), by = chrom]
+  seg[, `:=`(from = c(1, bps[-length(bps)] + 1), to = bps), by = chrom]
 
   # print(seg)