# Scatter plot matrix
data("iris")
pairs(iris[1:4], main = "Edgar Anderson's Iris Data", font.main = 4, 
        pch = 19)
pairs(iris[1:4], main = "Edgar Anderson's Iris Data", pch = 21, 
        bg = c("red", "green3", "blue")[codes(iris$Species)])

#Coplots
data(quakes)
coplot(long ~ lat | depth, data = quakes, pch = 21, bg = "green3")