tidyplot package

devtools::install_github("jbengler/tidyplots")
library(tidyverse)
library(tidyplots)

study %>%
  tidyplot(x = treatment, y = score, color = treatment) %>%
  add_mean_bar(alpha = 0.4) %>%
  add_sem_errorbar() %>%
  add_data_points_beeswarm()

Error in Ops.data.frame(guide_loc, panel_loc) :
  ‘==’ only defined for equally-sized data frames