flip_data <- tibble( n = 1:10, probability = (1/2)^n ) flip_data %>% ggplot() + aes(x = n, y = probability) + geom_col() + scale_x_continuous(breaks = 1:10)