flip_data %>% mutate(cumulative_probability = cumsum(probability)) %>% ggplot() + aes(x = n, y = cumulative_probability) + geom_col() + geom_hline(yintercept = 0.9) + scale_x_continuous(breaks = 1:10)
