Я хочу отдельно отметить значения NA для SpatialPointsDataFrame в карте spplot и критически включить эти точки в мою легенду.
В настоящее время я отмечаю значения NA красным X, отображая весь SpatialPointsDataFrame как набор красных X с помощью опции sp.layout, а затем снова отображаю шейп-файл поверх красных X. Но я не могу понять, как включить элементы sp.layout в легенду. И я также не знаю, как обозначать значения NA красным крестиком, а также отображать значения не-NA по цвету.
tzaea_latlon <- read.dta13("PATH/TZMapFlourPrices.dta")
tzdata <- tzaea_latlon[ , c("flourEA", "lflourEA","rnd")]
tzshp = SpatialPointsDataFrame(cbind(tzaea_latlon$lon_modr,tzaea_latlon$lat_modr),
tzdata, proj4string=CRS("+proj=longlat +datum=WGS84"))
tzshape1 <- tzshp[tzshp@data$rnd==1, ]
spplot(tzshape1,zcol="lflourEA",cuts=8,cex = .75,
key.space="right", digits=1,col.regions=brewer.pal(8, "YlGnBu"),
par.settings = list(axis.line = list(col = 'transparent')),
sp.layout = list(tz1))
Так что я просто хочу найти способ добавить красные буквы X под текущей легендой и пометить их как «отсутствующие данные». Данные для воспроизведения текущей карты находятся здесь: https://www.dropbox.com/sh/mysj7zqrfckb05h/AAAr0VX9Rd0NnCpzhz58MzbPa?dl=0