The most obvious methods that come to my mind are to use either Lattice or ggplot2. Here's an example using lattice:
library(lattice)
depthgroup<-equal.count(quakes$depth, number=3, overlap=0)
magnitude<-equal.count(quakes$mag, number=2, overlap=0)
xyplot(lat ~ long | depthgroup*magnitude,
data=quakes,
main="Fiji Earthquakes",
ylab="latitude", xlab="longitude",
pch=".",
scales=list(x=list(alternating=c(1,1,1))),
between=list(y=1),
par.strip.text=list(cex=0.7),
par.settings=list(axis.text=list(cex=0.7)))
In lattice you would change the main= parameter.
The above example was lifted from here.
I don't have a good ggplot2 example, but there are a metricasston of examples with ggpolot2 over at the learn r blog.
One option might be this example where they use ggplot2 and
opts (title = "RSS and NINO3.4 Temperature Anomalies \nand SATO Index Trends Since 1980")
But you would have to have all three graphs created in gg2plot, naturally.
I think you should be fine with either lattice or ggplot2.