R code used to produce the graph shown below.

data(airquality) 
par(mar=c(4,4.5,1.5,0)+.1) 
plot(airquality$Temp, airquality$Ozone, 
     xlab="Temperature", ylab="Ozone", 
     main="Airquality relationship?") 
pairs(airquality)