Using the data set fgl in MASS the following code
layout(matrix(1:9,3,3))
for(i in 1:9){
boxplot(fgl[,i] ~ type, data = fgl,main=dimnames(fgl)[[2]][i])}
produces a 3 by 3 array of plots, each one of which consists of six 
boxplots.
Is it possible to do this in lattice?
Steve
  "R version 2.7.2 (2008-08-25)" on Ubuntu 6.06