On Fri, 28 Sep 2012, David Winsemius wrote: > bwplot(quant ~ param , data=mg.d, main='Magnesium', ylab='Concentration > (mg/L)', subset= era=='Pre-mining' ) David, Don: Thank you. I tried subset= and era== separately, not together. Now I know. Much appreciated, Rich