Solução alternativa [e mais bonitinha] pra questão 5.5

[diogoprov]

penguins %>%
       nest_by(species) %>%
       mutate(corr_especies = cor(data$bill_length_mm, data$bill_depth_mm, use = "na.or.complete"))%>%
       summarise(corr_especies, .groups = "keep")

vai me retornar uma tabela com os nomes das espécies e as respectivas correlações:

# A tibble: 3 × 2
# Groups:   species [3]
  species   corr_especies
  <fct>             <dbl>
1 Adelie            0.391
2 Chinstrap         0.654
3 Gentoo            0.643