The grouped boxplot is one of the features included in the HH package. An example is posted on my site http://astro.ocis.temple.edu/~rmh/HH/bwplot-color.pdf See also the R-help email with an example vaguely similar to yours: https://stat.ethz.ch/pipermail/r-help/2009-March/190541.html