Поскольку не предоставлен набор образцов данных, его немного сложно настроить, поэтому я предоставляю выборку с 4 значениями, и каждый из этих цветов имеет 1000 значений.
library(tidyverse)
colorrs <- c("red","green","blue","orange")#load the 4 colors
colordata <- tibble::tibble(color = rep(colorrs, 1000),#make a tibble
values = rnorm(4000))
#plot normally distributed values by color and adjust the colors
colordata %>%
#Don't forget to set the aes values. X is first, y is second
ggplot(aes(color,values,fill = color))+
#geom_boxplot to create boxplots
geom_boxplot()+
scale_fill_manual(values = c("blue","green","orange","red"))