Skip to content
Snippets Groups Projects
Commit ebe9dbf8 authored by CAETANO Tomas's avatar CAETANO Tomas
Browse files

fix bug by_region FFT

parent 6b474206
No related branches found
No related tags found
No related merge requests found
......@@ -1367,8 +1367,8 @@ EMOTE_FFT = function(count_table = NULL,
summarise(value = median(count_norm_sum), .groups = "drop")
# Add 0 to missing positions
metaC_tot = tibble(reg = "start", start = feature_spaned_tb$start, strand = feature_spaned_tb$strand, width = q25+1, value = 0) %>%
rbind(tibble(reg = "stop", start = feature_spaned_tb$start, strand = feature_spaned_tb$strand, width = q25+1, value = 0)) %>%
metaC_tot = tibble(reg = "start", start = feature_spaned_tb$start, strand = feature_spaned_tb$strand, width = as.integer(q25+1), value = 0) %>%
rbind(tibble(reg = "stop", start = feature_spaned_tb$start, strand = feature_spaned_tb$strand, width = as.integer(q25+1), value = 0)) %>%
uncount(width) %>%
mutate(dist_to = c(rep(0:q25, nrow(feature_spaned_tb)),rep(-q25:0, nrow(feature_spaned_tb)) )) %>%
anti_join(metaC, by = c("start","strand","reg","dist_to")) %>% rbind(metaC) %>%
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment