现代统计图形 - 科学网—博客

5.9 条件分割图 85
1 > par(mar = rep(0, 4), mgp = c(2, 0.5, 0))
2 > library(maps)
3 > coplot(lat ~ long | depth, data = quakes, number = 4,
4 + ylim = c(-45, -10.72), panel = function(x, y,
5 + ...) {
6 + map("world2", regions = c("New Zealand",
7 + "Fiji"), add = TRUE, lwd = 0.1, fill = TRUE,
8 + col = "lightgray")
9 + text(180, -13, "Fiji", adj = 1)
10 + text(170, -35, "NZ")
11 + points(x, y, col = rgb(0.2, 0.2, 0.2, 0.5))
12 + })
lat
-45 -35 -25 -15
165 170 175 180 185
NZ
NZ
Fiji
Fiji
Given : depth
100 200 300 400 500 600
NZ
NZ
Fiji
Fiji
165 170 175 180 185
long
图 5.14: 给定震源深图的地震经纬度条件分割图:四幅散点图有相同的坐标
系,震源深度按左下、右下、左上、右上的顺序逐渐增加,可以看到地震
发生地点逐渐在向斐济岛靠近。
-45 -35 -25 -15