Here’s what I tried. I can plot a selection of regions, but I can’t seem to remove an arbitrary list of region numbers, unless I’ve done something wrong by selecting the regions I want to plot with departements[-exclude]. I also get an error when I try to use map.text to label a map with only the regions I’m selecting. >, departements <,- map(‘france’,namesonly=T, plot=FALSE) >, # returns a vector of names of regions >, >, exclude <,- c(47, #Alpes-Maritimes + 66, # Haute-Savoie + 76, # Savoie + 95, # Territore-de-Belfort + 109, 110, 111, # Var: Iles d’Hyeres + 49, 53, 54, 55, # Moribhan: Isles + 62, 64, # Vendee: Isles + 72, 75 # Charente-Maritime: Isles + ) >, >, depts <,- departements[-exclude] >, gfrance <,-map(‘france’, regions=depts) >, labels <,- (as.character(1:length(departements)))[-exclude] >, gfrance <,-map.text(‘france’, regions=depts, add=FALSE, labels=labels) Error in map.text(‘france’, regions = depts, add = FALSE, labels = labels) : map object must have polygons (fill=TRUE) Another problem, potentially more difficult for mapping data on the map of France is that the ‘departements’ are actually just the polygons in the map, arbitrarily numbered from east to west, and from north to south — they don’t correspond to the ‘official’ administrative region numbers. As well, the departement names don’t always match exactly (ignoring accents, e.g., Val-d’Oise vs. Val-Doise) so it would be another challenge to plot my historical data on the map of France. -Michael Stephane DRAY wrote: >, Source.