Le forum du Master ESA économétrie et statistique appliquée - Université d'Orléans

Vous n'êtes pas identifié.

Annonce

Vous êtes sur le forum du master ESA !

Le site du master ESA - description de la formation, notes de cours, contacts... est ici :

http://www.univ-orleans.fr/deg/masters/ESA/

#1 23-07-2007 15:22:24

alaa-eddine
Member
Date d'inscription: 07-03-2007
Messages: 398

Proc GMAP : ça peut toujours servir...

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 :

Code:

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 :

Code:

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... :

Code:

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 » :


Code:

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 !

Hors ligne

 

#2 23-07-2007 15:24:43

alaa-eddine
Member
Date d'inscription: 07-03-2007
Messages: 398

Re: Proc GMAP : ça peut toujours servir...

arf !! désolé je me suis tromper de rubrique ! je voulais mettre ça dans trucs et astuces !

Hors ligne

 

Pied de page des forums

Powered by PunBB
© Copyright 2002–2005 Rickard Andersson

[ Generated in 0.007 seconds, 8 queries executed ]