library(tidyverse)
penguins_wu <- palmerpenguins::penguins %>% 
  select(species, island, sex, body_mass_g, flipper_length_mm) %>% 
  filter(species %in% c('Gentoo', 'Chinstrap')) %>% 
  rename(penguin_species =species) %>%
  group_by(penguin_species, sex) %>% 
  summarize(mean_body = mean(body_mass_g, na.rm = TRUE),
            mean_flipper = mean(flipper_length_mm),
            n = n())
#> `summarise()` has grouped output by 'penguin_species'. You can override using
#> the `.groups` argument.

ggplot(penguins, aes(x = species, y = body_mass_g))+
  geom_jitter(width = 0.2, alpha = 0.5, color = 'navy')
#> Error in ggplot(penguins, aes(x = species, y = body_mass_g)): could not find function "ggplot"