df |> 
  dplyr::filter(!is.na(Species)) |> 
  ggplot() +
  aes(x = Species, y = Petal.Length, fill = Species) +
  geom_violin()