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

ddbox2 <- ggplot(dd, aes(Rname, Vac)) + geom_boxplot(aes(fill = factor(Population)))
ddbox2 + xlab("Region") + ylab("Unfilled Vacancies of more than 4 Months")