¿Cómo crear un atractivo mapa de coropletas en R?

He estado experimentando con técnicas de mapeo de coropletas en R, habiendo obtenido los beneficios de ggplot2 para crear gráficos hermosos dentro de un poderoso paquete de análisis de datos.

Cuando se trata de mapear, no he podido producir resultados comparables a los que obtengo de forma habitual de QGIS. Usando datos de flujo en Sheffield (datos completamente replicables y archivo .qgs aquí ), QGIS se produjo fácilmente lo siguiente:

LomejorquepudeproducirusandoR(usandoelcódigodescrito aquí ) fueron estos dos mapas, el primero con GISTools, el segundo utilizando ggplot2:

Para mí, la opción ggplot2 parece mucho más atractiva, si tan solo pudiera resolver el problema de las líneas defectuosas (probablemente un problema con el comando fortify (), o no leer en shapefiles usando readOGR () descrito aquí .)

Entonces, la pregunta es doble: ¿es la opción ggplot2 la mejor solución de mapeo de coropletas en R y, si es así, cómo puedo resolver el problema de las líneas blancas defectuosas?

El código replicable para descubrir lo que he hecho es aquí .

Editar: desde que la opción de salida de choropleth () es más atractiva:

Lafealdaddeestaopciónsepuedereducirexportandoaunaresoluciónmásaltayeliminandolaleyenda( otras leyendas de R están disponibles). Todavía no conseguí cambiar el grosor de la línea pero está mejorando! Las líneas rojas representan flujos a zonas que emplean > 5000 personas.