Объединяйте полигоны и рисуйте с помощью spplot()

Я не уверен в том, что вы пытаетесь сделать здесь.

Ошибка связана с тем, что spplot используется для отрисовки пространственных объектов с атрибутами, т.е. с ассоциированными данными. Ваш объект gadm относится к классу SpatialPolygonsDataFrame, поэтому он определяет полигоны и связанные с ними данные, доступ к которым можно получить через слот gadm@data. Когда вы используете UnionSpatialPolygons, вы получаете только объект класса SpatialPolygons, который можно построить с помощью plot, но не с spplot:

IDs <- gadm$ID_2
IDs[IDs %in% c(11500:11521)] <- "11500"
gadm_new <- unionSpatialPolygons(gadm, IDs)
plot(gadm_new)

Если вы хотите использовать spplot, вам нужно вручную объединить связанные данные так же, как вы объединили полигоны, а затем построить обратно SpatialPolygonsDataFrame. Один из способов сделать это заключается в следующем:

gadm_new <- gadm
## Change IDs
gadm_new$ID_2[gadm_new$ID_2 %in% c(11500:11521)] <- "11500"
## Merge Polygons
gadm_new.sp <- unionSpatialPolygons(gadm_new, gadm_new$ID_2)
## Merge data
gadm_new.data <- unique(gadm_new@data[,c("ID_2", "ENGTYPE_2")])
## Rownames of the associated data frame must be the same as polygons IDs
rownames(gadm_new.data) <- gadm_new.data$ID_2
## Build the new SpatialPolygonsDataFrame
gadm_new <- SpatialPolygonsDataFrame(gadm_new.sp, gadm_new.data) 

Затем вы можете использовать spplot для построения карты со связанным атрибутом:

spplot(gadm_new, "ENGTYPE_2", main="Test", lwd=.4, col="white")

Обратите внимание, что здесь я использовал только переменную ENGTYPE_2 ваших данных, а не переменную NAME_2, так как я не вижу смысла представлять переменную, где каждое значение кажется уникальным для каждого полигона.