navinlabcode/copykit

runMetrics

Closed this issue · 1 comments

Breakpoint count is doing on a chromosome by chromosome basis, but it needs to be on chr arm basis, otherwise it is counting centromere as a breakpoint.

Turns out I had already considered that when I wrote the code

copykit/R/runMetrics.R

Lines 80 to 82 in eff169f

purrr::map_dfc(x, function(i) {
length(rle(i)$values) - 1
}) %>%