survival curve

R code of the ggplot2 example

library(ggplot2)
library(HistData)

half <- HalleyLifeTable[HalleyLifeTable$age == 33, ]

p <- ggplot(HalleyLifeTable, aes(x = age, y = number)) +
  geom_step(linewidth = 0.7, colour = "#9b2367") +
  geom_point(data = half, size = 2, colour = "#9b2367") +
  annotate("text", x = half$age + 3, y = half$number + 40, label = "507 alive at 33",
           hjust = 0, vjust = 0, size = 3.2, colour = "#535353") +
  scale_x_continuous(labels = lex_number, expand = expansion(mult = c(0, 0.05))) +
  scale_y_continuous(labels = lex_number, expand = expansion(mult = c(0, 0.05)), limits = c(0, 1000)) +
  labs(title = "Half of Breslau's one-year-olds lived to 33",
       subtitle = "Halley's life table (1693): people alive at each age, out of 1\u2009000 aged one",
       x = "Age (years)", y = "People alive",
       caption = "Data: Halley 1693, via the HistData R package") +
  theme_lex_viz()