ogive

R code of the ggplot2 example

library(ggplot2)
library(dplyr)
library(HistData)

## Stigler's corrected table. Each cumulative share goes at the upper limit
## of its class, half an inch above the recorded size; the polygon starts at 0
## at the lower limit of the first class.
ogive <- ChestStigler |>
  arrange(chest) |>
  mutate(upper = chest + 0.5, share = cumsum(count) / sum(count))
ogive <- bind_rows(tibble(upper = min(ogive$chest) - 0.5, share = 0), ogive)

at_40 <- filter(ogive, upper == 40.5)

p <- ggplot(ogive, aes(x = upper, y = share)) +
  geom_line(linewidth = 0.7, colour = "#9b2367") +
  geom_point(size = 1.6, colour = "#9b2367") +
  annotate("segment", x = -Inf, xend = at_40$upper, y = at_40$share, yend = at_40$share,
           colour = "#7d7d7d", linetype = "dashed", linewidth = 0.3) +
  annotate("segment", x = at_40$upper, xend = at_40$upper, y = -Inf, yend = at_40$share,
           colour = "#7d7d7d", linetype = "dashed", linewidth = 0.3) +
  scale_x_continuous(labels = lex_number, breaks = lex_break_breaks, expand = lex_break_expand) +
  scale_y_continuous(labels = scales::label_percent(), expand = expansion(mult = c(0, 0.05))) +
  labs(title = paste0(scales::percent(at_40$share, accuracy = 1),
                      " of the chests measured 40.5 inches or less"),
       subtitle = paste0(format(sum(ChestStigler$count), big.mark = "\u2009"),
                         " Scottish militiamen; each point at the upper limit of its class"),
       x = "Chest girth (inches)", y = "Militiamen",
       caption = "Data: Quetelet 1846, corrected by Stigler 1986, via the HistData R package") +
  theme_lex_viz(break_not_zero = "x")