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 :

Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
library(maps)
france<-map("france")
france$names
je ne sais pas du tout comment faire le lien entre mon jeu de données et cette carte de France...

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 :

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)
j'ai une ligne d'erreur qui me dit : Erreur dans map.poly(database, regions, exact, xlim, ylim, boundary, interior, :
no recognized region names

alors que DPT1 a exactement les mêmes noms que les france$names

Merci pour votre aide