It is easy to control the main title with plot, but you will get much better looking plots if you use xyplot. library(lattice) tmp <- data.frame(x=rnorm(100), y=rnorm(100,4,2), ID=paste("Patient", rep(1:5, rep(20,5)))) xyplot(y ~ x | ID, data=tmp) Rich