airport_delay.md

library(tidyverse)
library(nycflights13)

delay <- flights |> 
    group_by(dest) |> 
    summarize(means = mean(dep_delay, na.rm = TRUE))

joined <- inner_join(delay, flights, by = "dest")

airports |> 
    inner_join(joined, c("faa" = "dest")) |> 
    group_by(means) |> 
    ggplot(aes(lon, lat, col = means)) + 
    borders("state") +
    geom_point() + 
    coord_quickmap()

^{Created on 2022-01-21 by the reprex package (v2.0.0)}