It doesn't seem to work. This what I have though.
I have a STATE_NAME variable in the data that I could potentially use,
assuming I know the write syntax.
Any ideas?
Thanks!
Pete
nc_file <- system.file("/LynchingHotSpots.shp", package="maptools")[1]
llCRS <- CRS("+proj=longlat +datum=NAD27")
nc <- readShapePoly("LynchingHotSpots.shp", proj4string=llCRS)
statesWanted <- c("MN","WI")
nc <- nc[nc$STATE_NAME %in% statesWanted,]
print(spplot(nc, c("Lynchings")))
_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo