#create grouped boxplots with custom legend labels p <- ggplot(data, aes(x=team, y=values, fill=program)) + geom_boxplot() + scale_fill_discrete(labels=c('High Program', 'Low Program')) #display grouped boxplots p