Объединить соседние регионы в R (объединить пространственные данные)?

Наверное, мне нужно было перефразировать свой ужасно сформулированный предыдущий вопрос (удалил). Вот еще одна попытка. Я хочу соединиться с соседними регионами так, чтобы исчезла их общая граница и была видна только их внешняя линия.

Вот воспроизводимый пример:

require(shapefiles)
require(sp)

xx <- readShapeSpatial(system.file("shapes/sids.shp", package="maptools")[1],
                   IDvar="FIPSNO", proj4string=CRS("+proj=longlat +ellps=clrk66"))

# show all the subregions
plot(xx)

Теперь давайте рассмотрим только регионы 3 и 5

plot(xx[c(3,5),])

Как я могу просто агрегировать эти регионы. На практике то, что я хочу сделать, похоже на то, чтобы иметь карту всего континента, показывающую все страны, и создать карту, показывающую Северную Америку и Южную Америку.

Мне кажется, это довольно распространенная задача, но я пока не могу найти подходящую функцию для ее выполнения. Я просто пропустил функцию или я могу просто вручную?