Bonjour à tous !
J'ai fait une petite découverte en fouinant dans la doc SAS.. du coté SAS/GRAPH, et donc je suis tombé sur la PROC GMAP qui permet de s'amuser avec les .shp ( cartes ).

Tout d'abord il faut se munir d'un .shp ( exemple : le « south.shp » dans les samples fournis avec GeoDa )

On commence par une PROC MAPIMPORT pour charger la carte dans SAS :

proc mapimport out=south datafile='C:\Program Files\GeoDa\Sample Data\south.shp' ;
run;

NOTE:  41163 observations copiées de la table WORK.SOUTH.
NOTE: La table WORK.SOUTH a 41163 observations et 72 variables.
NOTE: La procédure SORT a utilisé (Durée totale du processus) :
      temps réel                   3.67 secondes
      temps processeur   0.26 secondes

Pas bon !! ça le fait avec toutes les cartes je ne sais pas pourquoi...

astuce :

proc sort data=south ;by fipsno;run;

data south1 ; set south;
by fipsno ;
if first.fipsno then output;run;

NOTE:  41163 observations copiées de la table WORK.SOUTH.
NOTE: La table WORK.SOUTH1 a 1412 observations et 72 variables.
NOTE: L'étape DATA a utilisé (Durée totale du processus) :
      temps réel                   1.15 secondes
      temps processeur   0.04 secondes

On dispose d'un SOUTH1.WORK avec 1412 obs ( plus logique car il n'y a que 1412 contés dans le sud des usa )

Et pour afficher... :

proc gmap data=south1 map=south;
id fipsno;
choro HR80 / levels=5 coutline=black woutline=3 ;
run;
quit;

id : clé primaire de la table (unique)

choro : pour afficher la carte en couleur selon les niveaux de la variable HR80

on peut choisir l'option « PRISM » au lieu de choro pour afficher la carte en... prism lol ( ça fait une carte en 3D )

« BLOCK » pour avoir des barres verticales..

ex :

Block HR80 / levels=5 shape=cylinder


Dans le cas des variables discrètes on rajoute l'option « discrete » :

title1 'elections 2000';
footnote1 'usa';

pattern1 value=solid color=lime;
pattern2 value=solid color=cyan; 

proc gmap data=elections map=usa;
id fipsno;
choro winner / discrete;
run;
quit;

...etc ! ça marche à peu prêt comme une PROC GPLOT.

voili voualou ! Plus d'infos « The GMAP Procedure » dans l'aide SAS.

oulaaa 15h25... sauve-qui-peut !