je suis tout béotien pour les statistiques spatiales. J'essaye de reproduire la figure page 9 de ce rapport à propos des bébés sans main dans le département de l'Ain. Voici où j'en suis avec le code reproductible suivant :
Code : Tout sélectionner
load(url("http://pbil.univ-lyon1.fr/R/donnees/qrbb/AinCol.Rda"))
par(mar = c(0, 0, 2, 0) + 0.1)
plot(ain, col = ain$col)
grid(lty = 1)
main <- "Localisation des 7 cas d'agénésie des membres supérieurs\ndans le département de l'Ain entre 2009 et 2014, inclus"
title(main = main, line = -1)
par(lend = "butt")
legend("bottomleft", inset = 0, lwd = 10,
legend = c("Cas d'agénésie", "Cluster d'agénésie de la main"),
col = c("yellow", "lightblue"))
Je ne suis pas satisfait du rendu, j'aimerais avoir l'équivalent de asp = 1 pour les plots standards: mon grid() devrait me faire des carrés en première approximation. Je me doute bien que c'est compliqué dès lors que l'on projette un ellipsoïde sur un plan. Est-ce que quelqu'un aurait une solution simple pour moi ?
Amicalement,
jean