pg_adelie <- pg %>% filter(species == "Adelie")
cor(pg_adelie$bill_length_mm, pg_adelie$bill_depth_mm, use = "complete.obs") # [1] 0.3914917
.EOD.