#
# Script to demonstrate what I am trying to do.
#

# Simulate some data:
Year <- factor(rep(4:8,each=50,times=2))
Type <- rep(c("National","Local"),each=250)
M0   <- 1300+50*(0:4)
set.seed(42)
M1   <- M0 + runif(5,-100,-50)
X0   <- rnorm(250,rep(M0,each=50),150)
X1   <- rnorm(250,rep(M1,each=50),100)
DemoDat <- data.frame(Year=Year,Score=c(X0,X1),Type=Type)

# Grouped boxplots:
library(ggplot2)
print(ggplot(data=DemoDat) +
    geom_boxplot(aes(x=Year, y=Score, color=Type,fill=Type),
                 position=position_dodge(1),alpha=0) +
    theme_minimal() +
    scale_fill_discrete(name="National v. Local") +
    ylim(700,2100))