Lab Hw11

Create a airport map similar to the one above with the color representing the number of flights per day from NYC to each airport.

library(nycflights13) 
library(tidyverse)  

## ── Attaching core tidyverse packages ──────────────────────── tidyverse 2.0.0 ──
## ✔ dplyr     1.1.4     ✔ readr     2.1.6
## ✔ forcats   1.0.1     ✔ stringr   1.6.0
## ✔ ggplot2   4.0.1     ✔ tibble    3.3.1
## ✔ lubridate 1.9.4     ✔ tidyr     1.3.2
## ✔ purrr     1.2.1     
## ── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
## ✖ dplyr::filter() masks stats::filter()
## ✖ dplyr::lag()    masks stats::lag()
## ℹ Use the conflicted package (<http://conflicted.r-lib.org/>) to force all conflicts to become errors

flights %>%
  group_by(dest) %>%
  summarise(flights_per_day = n() / 365) -> flights_per_day_data


airports %>%
  semi_join(flights_per_day_data, by = c("faa" = "dest")) %>%
  left_join(flights_per_day_data, by = c("faa" = "dest")) %>%
  ggplot(aes(lon, lat)) +
  borders("state") +
  geom_point(aes(color = flights_per_day)) +
  scale_color_gradient(low = "green", high = "red") +
  coord_quickmap() +
  labs(color = "Flights/day")

## Warning: `borders()` was deprecated in ggplot2 4.0.0.
## ℹ Please use `annotation_borders()` instead.
## This warning is displayed once per session.
## Call `lifecycle::last_lifecycle_warnings()` to see where this warning was
## generated.