scatter plot
plot(faithful$eruptions, faithful$waiting)
cor(faithful$eruptions, faithful$waiting)
## [1] 0.9008112