df |> 
  tidyr::drop_na(Species) |> 
  ggplot() +
  aes(x = Species, y = Petal.Length, fill = Species) +
  geom_violin()