monarchs-exploration.R

library(tidyverse)
library(lubridate)

source("https://raw.githubusercontent.com/mjfrigaard/RDataSets/master/DotR/monarchs.R")

monarchs <- 
  monarchs %>%
  mutate(reign_length = lubridate::interval(from, to) / years(1))

monarchs %>% 
  ggplot() +
  geom_segment(aes(x = from, y = 0,
                   xend = to, yend = reign_length)) +
  geom_point(aes(x = to, y = reign_length))

monarchs %>%
  ggplot(aes(x = reign_length)) +
  geom_histogram() +
  stat_function(fun = function(x, lambda, total) {
    lambda * exp(-lambda * x) * total
  }, args = c(lambda = 1 / mean(monarchs$reign_length, na.rm = T),
              total  = 200))