Merge pull request friendsofstrandseq#41 from p-smirnov/plot-sv-calls…

…_dev_p-smirnov1 Fixing small indexing errors in plot-sv-calls_dev
grimbough · Aug 9, 2023 · 05c7959 · 05c7959
2 parents 4463ca3 + 8d23cf5
commit 05c7959
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)