par(mfrow=c(1,4), mar=c (5, 4, 4, .5))
with(anscombe, plot(x1, y1, xlab = expression(Y[1]), ylab = expression(Y[2]), main = paste("r = ",round(cor(x1, y1),2)))); abline(3,0.5)
with(anscombe, plot(x2, y2, xlab = expression(Y[1]), ylab = expression(Y[2]), main = paste("r = ",round(cor(x2, y2),2)))); abline(3,0.5)
with(anscombe, plot(x3, y3, xlab = expression(Y[1]), ylab = expression(Y[2]), main = paste("r = ",round(cor(x3, y3),2)))); abline(3,0.5)
with(anscombe, plot(x4, y4, xlab = expression(Y[1]), ylab = expression(Y[2]), main = paste("r = ",round(cor(x4, y4),2)))); abline(3,0.5)
mtext("(a)",side = 3, at = -50, line = 1); mtext("(b)",side = 3, at = -32, line = 1)
mtext("(c)",side = 3, at = -14, line = 1); mtext("(d)",side = 3, at = 5, line = 1)Run the code above in your browser using DataLab