{"word_count": 307, "character_count": 1950}
Suite à ma note sur le megacombo “R + Google Map Charts” et aux commentaires avisés d'Éric Verdeil, un autre blogueur du coin, j'essaie de fabriquer des cartes avec maptools sous R et spmap sous Stata. J'essaie, plus exactement, de reproduire ce type de cartes, qui ont un lien avec ma thèse (message subliminal à quelques lecteurs et à moi-même).Après quelques heures de code et une bonne dizaine d'emails échangés avec Baptiste Coulmont, dont les conseils avisés m'ont permis de régler plusieurs problèmes que je n'aurais jamais réussis à élucider tout seul, je pense avoir simplement réussi à me prouver que ce n'était pas simple.Mais après quelques minutes de code supplémentaire, toujours avec l'aide inestimable de Baptiste Coulmont, je peux enfin mettre cette note à jour et visualiser correctement, voire même élégamment (si le terme est approprié), le taux de mortalité par cancer du poumon chez les femmes en 2008 dans les 27 États-membres de l'Union européenne, en suivant un modèle de l'European Cancer Observatory :J'ai utilisé un fonds de carte EuroGeographics, des données de l'International Agency for Cancer Research, et la patience de Baptiste, plus quelques autres correspondants qui recevront bientôt tous un email de remerciements avec un lien vers une page permettant de répliquer l'intégralité de ma démarche.En attendant, le code est disponible, avec la carte et des explications je ne sais pas vraiment si j'ai le droit de distribuer les fichiers SHP (shapefile) utilisés pour la carte, même avec une mention complète de la source.Comme l'écrit Baptiste, le code final n'est pas très impressionnant, mais il a fallu de nombreuses étapes avant d'arriver à intégrer exactement les bonnes commandes, et à résoudre les situations illogiques. Le même résultat sous Stata (non reproduit ici) est moins élégant, mais s'obtient en une demi-douzaine de lignes à partir d'un modèle disponible dans la documentation du module approprié.