ggplotで複数プロット（スケール一緒）を描く

library(ggplot2)
library(reshape2)
tmp <- shp1@data[,c(13,19,24)]　##描きたい値（複数）
tmp.melt <- melt(tmp, id="id")　　　##idを使ってmelt. idがない場合はcbindでくっつけておく
tmp.ggmap <- fortify(shp1, region="id") ##shp1（描きたいベクターデータをggplot用に変換する）
tmp.ggmap <- join( tmp.ggmap, tmp.melt, by ="id")　##くっつける
PC1 <- ggplot( tmp.ggmap, aes(x = long, y = lat, group=group, fill=value)) +
geom_polygon() +
coord_map() +
facet_wrap(~variable) +　　##こいつが複数描く
scale_fill_gradient2(low="blue", high="red", midpoint=0) +　##色
theme(legend.position="right")

参考: http://www.cartopedia.co.uk/blog/2013/06/02/facet-wrapping-multivariate-data-reshape-and-ggplot/