## Load districts of Germany
load(system.file("shapes", "districtsD.RData", package="surveillance"))
## indicate districts at the border
districtAtBorder <- polyAtBorder(districtsD)
## plot to check
plot(districtsD, col=districtAtBorder)Run the code above in your browser using DataLab