df %>% ggplot(aes(x=age_group, y=height, fill=age_group)) + geom_boxplot(width=0.5,lwd=1)+ labs(subtitle="Filling Boxplot with Colors by a Variable")