library(palmerpenguins)
library(dplyr)
mean_sd <- function(x, na.rm = TRUE) {
c(mean = mean(x, na.rm = na.rm),
sd = sd(x, na.rm = na.rm))
}
penguins |> group_by(species, sex) |>
summarize(mean_sd(body_mass_g))
# A tibble: 16 × 3
# Groups: species, sex [8]
species sex `mean_sd(body_mass_g)`
<fct> <fct> <dbl>
1 Adelie female 3369.
2 Adelie female 269.
3 Adelie male 4043.
4 Adelie male 347.
5 Adelie <NA> 3540
6 Adelie <NA> 477.
7 Chinstrap female 3527.
8 Chinstrap female 285.
9 Chinstrap male 3939.
10 Chinstrap male 362.
11 Gentoo female 4680.
12 Gentoo female 282.
13 Gentoo male 5485.
14 Gentoo male 313.
15 Gentoo <NA> 4588.
16 Gentoo <NA> 338.