depth.col <-
with(quakes,
level.colors(depth, at = do.breaks(range(depth), 30),
col.regions = terrain.colors))
xyplot(lat ~ long | equal.count(stations), quakes,
strip = strip.custom(var.name = "Stations"),
colours = depth.col,
panel = function(x, y, colours, subscripts, ...) {
panel.xyplot(x, y, pch = 21, col = "transparent",
fill = colours[subscripts], ...)
})
Run the code above in your browser using DataLab