Pessoal como destacar o contorno de um conjunto de municípios de um mapa .shp importado pelo R? Criei uma função chamada mapa.coropletico definida abaixo só que os municípios apresentam divididos por exemplo em macro regiões. Gostaria de contornar as macros desses municípios por linhas vermelhas por exemplo.<div>
<br></div><div><br></div><div><div>mapa.coropletico <- function(dados){</div><div>      library(maptools)</div><div>      library(RColorBrewer)</div><div>      library(classInt)</div><div>      mapa <- readShapePoly("C:\\Users\\GORR-PEDRO\\INDICADORES\\MALHA MUNICIPAL (MAPA)\\PB.shp")  </div>
<div>      colours <- brewer.pal(5, "Blues")</div><div>      brks<-classIntervals(dados, n=5, style="quantile")</div><div>      plot(brks, pal=colours)</div><div>      brks<- brks$brk</div><div>
      plot(mapa, col=colours[findInterval(dados, brks,all.inside=TRUE)], axes=TRUE)  </div><div>      title(paste (main = "MAPA DA PARAÍBA"), xlab = "Longitude", ylab = "Latitude")</div><div>
      SpatialPolygonsRescale(layout.north.arrow(1), offset= c(-38.5,-8.5),plot.grid=FALSE, scale = 0.4) # Plota a direção do Norte.</div><div>      legend(locator(1), legend=leglabs(brks), fill=colours, bty="o", cex = 0.7) # Plota a legenda ao lado do gráfico.</div>
<div>      SpatialPolygonsRescale(layout.scale.bar(height = 0.01), offset= c(-38.9,-8.0), fill= c("transparent", "black"), plot.grid=F)</div><div>      text(-36.84889, -8.70114, "10KM", cex= 1)</div>
<div>      text(534000,152000, "Boundary Data Crown Copyright Ordnance Survey 2009.", cex= 1)</div><div>      text(556500, 166000, "Coberturas", cex= 1)          </div><div>}</div><div><br></div>-- <br>
<font color="#000099">Saudações,</font><div><font color="#000099">Pedro Rafael Diniz Marinho.</font></div><div><font color="#000099">Currículo Lattes: <a href="http://buscatextual.cnpq.br/buscatextual/visualizacv.do?id=K4250792T6" target="_blank">http://buscatextual.cnpq.br/buscatextual/visualizacv.do?id=K4250792T6</a></font></div>
<font color="#000099"><br></font><br>
</div>