data(geyser, package="MASS")
library(ggplot2) ggplot(geyser, aes(duration, waiting)) + geom_point() + geom_density2d()