Load packages library(tidyverse) library(reshape2) # for data library(mosaic) library(sjmisc) library(skimr) Data setup data(tips) Aggregate data per group tips_aggr <- tips %>% group_by(smoker) %>% summarise(tip_avg = mean(tip), tip_md = median(tip), tip_sd = sd(tip), tip_iqr = IQR(tip)) tips_aggr #> # A tibble: 2 x 5 #> smoker tip_avg tip_md tip_sd tip_iqr #> <fct> <dbl> <dbl> <dbl> <dbl> #> 1 No 2.99 2.74 1.38 1.50 #> 2 Yes 3.01 3 …