使用 ggsubplot 指定绘图顺序

Question

我正在使用 ggplot 和 ggsubplot 向地图添加条形图，但无法弄清楚如何指定首先绘制哪个。我想先绘制北方的，这样它们就位于任何重叠的地块后面。使用低填充阿尔法，这些应该仍然是可见的。这是工作流程：

library(ggsubplot)
library(ggplot2)
library(maps)
library(plyr)

world_map = map_data("world")
(p = ggplot() + geom_polygon(data = world_map, aes(x=long, y=lat,group=group)))

d = ddply(world_map,.(region),summarize,long=mean(long),lat=mean(lat))
d = d[sample(1:nrow(d), 50),]
d = rbind(d,d)
d$cat = rep(c('A','B'), each=nrow(d)/2)
d$value = sample(1:10, nrow(d), rep=T)
head(d)

p + geom_subplot(data=d, aes(long, lat, group=region, subplot = geom_bar(aes(cat, value, fill=cat), 
                 col='black', alpha=0.9, stat="identity")), width = 30, height=30)

如您所见，情节顺序似乎很随机。所以我尝试将地区（国家）更改为有序因素：

d$region = factor(d$region, ordered=T)
(ord = count(d[,c('region','lat')], vars=c('region','lat')))
ordered_levels = order(ord$lat, decreasing=T)
print(ord[ordered_levels,])
levels(d$region) = levels(d$region)[ordered_levels]
levels(d$region)

p + geom_subplot(data=d, aes(long, lat, group=region, subplot = geom_bar(aes(cat, value, fill=cat), 
                 col='black', alpha=0.9, stat="identity")), width = 30, height=30)

但这似乎并不能解决问题。非常感谢任何建议。

Answer 1

这是你的想法吗？

正如您所指出的，您需要按纬度订购d，并在对geom_subplot(...) 的调用中使用group=lat。

set.seed(1)
d = ddply(world_map,.(region),summarize,long=mean(long),lat=mean(lat))
d = d[sample(1:nrow(d), 50),]
d = rbind(d,d)
d$cat = rep(c('A','B'), each=nrow(d)/2)
d$value = sample(1:10, nrow(d), rep=T)
d <- d[order(d$lat),]
head(d)
p + geom_subplot(data=d, aes(long, lat, group=lat, 
      subplot = geom_bar(aes(cat, value, fill=cat), 
        col='black', alpha=0.9, stat="identity")), width = 30, height=30)

Answer 2

在生成的图中，底层地图几乎是一团糟，但这种预购似乎将高纬度项目带到了最前面：

world_map = world_map[order(world_map$lat),]

不清楚您是否希望在靠近赤道的纬度下绘制负纬度，因此您还可以选择使用 abs(world_map$lat) 作为顺序。