r - 页面中的多个 spplots 和单个索引

Question

我有一个名为的多边形形状文件CROPS.shp和以下键 -

1 - Rice
2 - Wheat 
3 - Corn 
4 - Millets 

我使用包将形状文件读rgdal入一个名为my_crops. 现在my_crops@data有六个不同的字段，即Y1、Y2、Y3、和。这些字段中的每一个都具有从 1 到 4 的值。Y4Y5Y6

例如，前几行crops@data就像

Y1   Y2   Y3   Y4   Y5    Y6 
1    2    4    1    1    2
3    4    1    1    1    2
4    2    2    2    1    3

我想做什么？

我想要一个包含 6 个图（2 行 3 列）的页面，其中第一个图将基于 field Y1，最后一个图将基于 field Y6。在每个 shapefile 图中，绘图（多边形形状）的颜色应基于值 1 到 4。

这些图应在页面底部附有一个图例。

Answer 1

你的意思是这样的吗？

library(maptools)
nc.sids <- readShapePoly(system.file("etc/shapes/sids.shp", package="spdep")[1],
  ID="FIPSNO", proj4string=CRS("+proj=longlat +ellps=clrk66"))

set.seed(1)
nc.sids$a = factor(sample(letters[1:4], 100, replace=T))
nc.sids$b = factor(sample(letters[1:4], 100, replace=T))
nc.sids$c = factor(sample(letters[1:4], 100, replace=T))
nc.sids$d = factor(sample(letters[1:4], 100, replace=T))
nc.sids$e = factor(sample(letters[1:4], 100, replace=T))
nc.sids$f = factor(sample(letters[1:4], 100, replace=T))
spplot(nc.sids[c("a", "b", "c", "d", "e", "f")], 
    colorkey = list(space = "bottom"), layout = c(3, 2),
    as.table = TRUE)