R 在一个页面中有多个SPPLOT，并且具有单个索引

我有一个名为

CROPS.shp

的多边形形状文件和以下键-

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

我使用

rgdal

包将shape文件读入名为

my_crops

的变量中。现在

my_crops@data

有六个不同的字段，即

Y1

，

Y2

，

Y3

，

Y4

，

Y5

和

Y6

。每个字段都有1到4的值

例如，

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列）的页面，其中第一个绘图基于字段

Y1

，最后一个绘图基于字段

Y6

。绘图（多边形形状）的颜色应基于每个shapefile绘图中的值1到4

---

绘图应在页面底部附有一个图例

你的意思是这样的吗

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)

库（maptools）
北卡罗来纳州小岛屿发展中国家