diff --git a/src/emote-tk.R b/src/emote-tk.R index b2b297d107925d852b76fc20e8b70e46b4303e85..df8bac6a6a7486a0278fe996b1c3bf32704e23fd 100644 --- a/src/emote-tk.R +++ b/src/emote-tk.R @@ -1413,7 +1413,7 @@ EMOTE_frame_preference = function(count_table = NULL, regions_count_table = lapply(1:nrow(feature_spaned_tb), function(i) EMOTE_parse_region(region = feature_spaned_tb[i,], count_table = count_table)) %>% bind_rows() - + # Retrieve percent of read mapped on each 3 translation frames if(!by_region){ frame = regions_count_table %>% @@ -1424,8 +1424,8 @@ EMOTE_frame_preference = function(count_table = NULL, else{ frame = regions_count_table %>% mutate(frame = abs(dist_to%%3) + 1) %>% - group_by(start, strand,frame) %>% - summarise(percent = sum(count_norm_sum)/sum(regions_count_table$count_norm_sum), .groups = "drop") %>% + group_by(start, strand, frame) %>% + summarise(percent = sum(count_norm_sum), .groups = "drop") %>% inner_join(annot, by = c("start","strand")) } return(frame)