df <- iris
df$Species[c(5,23,44,45,12,125)] <- NA 
#
df |>  
  ggplot() +
  aes(x = Species, y = Petal.Length, fill = Species) +
  geom_violin()