Picking out peaks that fit a pattern-CodePudding

The

However, these are the expected DNA fragment sizes (in bp):

ladder <- c(50, 75, 100, 125, 150, 200, 250, 300, 350, 400, 450, 475, 500, 550, 600, 650, 700, 750, 800, 850, 900, 950, 1000)

And here is a plot with the correct peaks picked out:

Is there any way to get a peak finder such as stat_peak or ggpmisc:::find_peaks to take an pattern of peaks into account when selecting peaks to achieve the second plot??

CodePudding user response：

Before plotting, doing some data manipulation to pull out the maximum value for each of the 23 DNA fragment groups with base R max function, and adding the max plot with additional geom_ layer for the max values.

Here is small reprex example that plots the max value for each group with "red".

data.frame(x = 1:9) %>% mutate(y = rep(LETTERS[1:3], c(3, 3, 3) )) %>% group_by(y) %>% mutate(max = max(x), .groups = "drop") %>% ggplot(aes(x = y, y = x)) geom_point() geom_point(aes(x = y, y = max), col = "red")

CodePudding user response：

I just added an example variable named frag_type that labels two out of 23.