Bonjour à tous,
j'ai un jeu de données qui fournit pour différentes régions françaises des teneurs en toxines trouvées au champs... j'ai classé ces teneurs en 3 classes et j'aimerai créer une carte où j'ai la répartition géographique de ces trois classes.
Je me suis donc lancée dans la fonction maps :
je ne sais pas du tout comment faire le lien entre mon jeu de données et cette carte de France...
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3 library(maps) france<-map("france") france$names
voici mon jeu de données :
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
8
9
10 > summary(data) ESPECE ANNEE REGION DPT DPT1 CLASSE_PLUIE CLASSES_DON BLE TENDRE:2349 Min. :2001 BRETAGNE :553 COTES D ARMOR : 279 Cotes-Darmor : 279 <10 : 452 <1250 :1878 1st Qu.:2003 CENTRE :337 SEINE ET MARNE: 198 Seine-et-Marne : 198 >40 : 445 >1750 : 113 Median :2004 ILE DE France:282 ILE ET VILAINE: 162 Ille-et-Vilaine: 162 10-40:1282 1250-1750: 44 Mean :2005 NORMANDIE :191 CALVADOS : 137 Calvados : 137 NA's : 170 NA's : 314 3rd Qu.:2007 PICARDIE :167 EURE ET LOIR : 114 Eure-et-Loir : 114 Max. :2012 (Other) :783 (Other) :1423 (Other) :1423 NA's : 36 NA's : 36 NA's : 36 >
j'ai essayé ceci, mais en vain :
j'ai une ligne d'erreur qui me dit : Erreur dans map.poly(database, regions, exact, xlim, ylim, boundary, interior, :
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3 col <- c(1,2,3) map ('france', regions=data$DPT1, col=col[data$CLASSES_DON], fill=TRUE, add=TRUE) legend ('bottomleft', legend=attr (df$CLASSES_DON, 'levels'), col=col, lty=1, lwd=10)
no recognized region names
alors que DPT1 a exactement les mêmes noms que les france$names
Merci pour votre aide
Partager