library(maptools)
cook <- readShapePoly(system.file("maps/CookCensusTracts.shp",
package="McSpatial"))
sampvar <- cook$CHICAGO==1|(!is.na(cook$CAREA)&cook$CAREA=="O'Hare")
cook <- cook[sampvar==TRUE,]
Run the code above in your browser using DataLab