Estimating effect sizes (part 2) using the emmeans package in R