`2009Wages` <- read.csv("~/2009Wages.csv") dd <- `2009Wages` library(ggplot2) dbox <- ggplot(dd, aes(Rname, Vac, fill = Population)) + geom_boxplot(outlier.colour = "blue") dbox + xlab("Region") + ylab("Unfilled Vacancies of more than 4 Months")