Zach, try: pdf("Name.pdf") par(mfrow=c(2,2)) par(mai=c(0.5,0.3, 0.6, 0.4)) plot(x, y, ..........) plot(x, z ...........) plot(z, y,...........) plot(z+x, x-y, .....) dev.off() Andrew