Estoy tratando de averiguar cómo realizar las operaciones (tan fáciles en GIS) en R.

Tomemos algunos ejemplos de datos de polígonos despdep paquete

library("spdep")
c <- readShapePoly(system.file("etc/shapes/columbus.shp", package="spdep")[1])
plot(c)

Me las arreglé para averiguar que puedo elegir polígonos con declaraciones lógicas usandosubset. Por ejemplo:

cc <- subconjunto (c, c @ data $ POLYID <5) gráfico (cc)

Ahora, supongamos que tengo otro marco de datos que me gustaría unir a mis datos espaciales:

POLYID=1:9
TO.LINK =101:109
link.data <- data.frame(POLYID=POLYID, TO.LINK=TO.LINK)

Usando estos dos conjuntos de datos, ¿cómo puedo obtener dos marcos de datos espaciales:

Primero, consiste en polígonos que tienen su ID en el segundo marco de datosSegundo, consiste en el conjunto opuesto - polígonos que no existen en el segundo marco de datos.

¿Cómo podría llegar a este punto?