1 Tutoriel 3 : Grahiques avec SAS/GRAPH

faut des sorties graphiques à partir des principales procédures de Statistique

Tutoriel 3 : Grahiques avec SAS/GRAPH (module SAS/Stat) sans nécessité le module SAS/Graph.
Le choix est fait ici de laisser SAS produire les graphiques par défaut dans la
plupart des procédures, c’est tout à fait utile dans uns phase exploratoire, mais
d’apprendre à utiliser au moins certaines des possibilités offertes par le module
Résumé SAS/Graph. Celui-ci, très complet, permet en effet de produire la plupart des
graphes nécessaires à une illustration pertinente et efficace des résultats sta-
Le module SAS/GRAPH permet de tracer des graphes dont la réso- tistiques avec une très (trop) large variété d’options. Il serait en effet frustrant
lution est adaptée au périphérique utilisé (écran, imprimante, page de commencer à produire un graphique dans l’environnement ODS, afin de
web. . . ) et avec toutes les options possibles. Cette vignette décrit les profiter de l’apparente simplicité, pour finalement s’apercevoir que les bonnes
principales procédures (gchart, gplot, annotate) et options parmi un options nécessaires à la réalisation du graphique souhaité ne sont disponible
nombre disponible considérable. Bien que la plupart des procédures que dans SAS/Graph !
statistiques fournissent (open delevery system) des graphiques par
défaut, la construction spécifique du graphique adapté au problème 1.2 SAS/Graph
posé est souvent incontournable.
Les programmes SAS ci-dessous fournissent les exemples types de gra-
Plan des tuteuriels : phiques uni et bidimensionnels à partie de SAS/Graph. Certains (histo-
• Prise en main grammes, profils) ne sont pas une simple utilisation d’une procédure standard,
• Gestion des données ils nécessitent la construction de tables spécifiques. Il est important de savoir
• Graphiques contrôler les options de ses graphiques afin de s’affranchir de celles imposées
• Macros-commandes par la plupart des logiciels guidés par menu. La plupart de ces graphiques sont
• Bases de données ceux qui illustrent les vignettes de statistique descriptive élémentaire.
Un soin particulier doit être apporté à la génération de fichiers graphiques.
1 Objectifs En effet, l’apparente simplicité du copié/collé entre SAS et un traitement
de texte sous Windows masque une limitation ; il s’agit d’une simple re-
1.1 ODS vs. SAS/Graph copie d’écran qui impose la définition du graphique. En créant un fichier
postscript ou jpeg, on utilise, au moment de l’impression ou de la visua-
Deux systèmes de production de graphiques sont en "concurrence". Le
lisation, la définition optimale du support. Par ailleurs, la production "program-
premier historique est partie intégrante du module SAS/Graph tandis que le
mée " de fichiers graphiques permet d’inclure ceux-ci dans une page html ou
deuxième, très récent, est en production depuis la version 9.3 dans le module
xml pour la mise à jour automatique et en temps réel de sites web.
SAS de base à travers l’ODS (output delivery system) au moment de créa-
tion de fichiers de sortie au format HTML par défaut et paramétrable en jpeg, Des macros seront mises au point lors d’une prochaine séance afin de facili-
png, gif.... La plupart des procédures graphiques apparaissent donc en double : ter l’exécution de ces tâches.
i.e. gplot de SAS/Graph et sgplot de SAS de base, avec évidemment des Les petits jeux de données sont saisis au clavier tandis que les plus volumi-
syntaxes et possibilités différentes mais des principes similaires, comme l’ad- neux sont dans des fichiers accessibles dans le répertoire data du siteWikistat.
jonction d’un ensemble d’annotations.
L’intérêt majeur des sorties ODS est de fournir systématiquement et par dé-
 2 Tutoriel 3 : Grahiques avec SAS/GRAPH

2 Environnement ODS RTF CLOSE;

2.1 Taille des graphiques 3 Commandes globales

Le graphique demandé est tracé dans une zone dont les dimensions sont défi-
nies par les paramètres hsize et vsize de la commande globale goptions Elles définissent des objets (axes, symboles, trames, légendes) et les options
(les valeurs maximales sont prises par défaut) diminués de l’espace nécessaire utilisés pour les tracés ; elles demeurent valables jusqu’à une nouvelle défini-
à l’édition des titres, sous-titres, notes, légendes,. . . . Les dimensions peuvent tion ou la fin de la session sas.
être exprimées en trois unités : pouce, cm ou pct qui signifie “pourcen- 3.1 Axes
tage de la dimension totale”. Cette dernière unité est préférable pour exprimer
les tailles de caractères et symboles lorsque les dimensions globales, liés au Des types d’axes, numérotés de 1 à 99 sont définis avant de pouvoir être
périphérique de sortie, sont sujettes à modifications. utilisés dans les différents graphiques. Ils précisent l’échelle (liste de valeurs,
logarithmique), l’apparence (longueur, couleur, épaisseur, style de ligne, ori-
2.2 Sauvegarde des graphiques gine), les marques d’échelle (nombre, couleur, épaisseur, hauteur), les valeurs
En l’absence de commande explicite, le graphe apparaît sur le périphérique des échelles (format), le libellé (police,. . . ).
par défaut, l’écran de l’ordinateur, dans la fenêtre de visualisation des résultats.
Une fois les graphiques mis au point, ils peuvent être sauvés dans des fichiers axis1 order=(1973 to 1981 by 2)
ou un traitement de texte. label=(’annee’)
minor=(number=1)
Un fichier de format .jpg est systématiquement créé (sous windows) dans width=3;
un répertoire temporaire ainsi qu’un fichier [Link] contenant le gra- axis2 order=(0 to 10000 by 1000)
phique ; consulter la fenêtre du journal. Il devrait être possible de contrôler la label=(’Revenu en francs’)
destination de ce fichier. minor=none
2.3 Image écran width=3;

Un clic droit sur le graphique (windows mais unix ?) ouvre un menu qui 3.2 Légendes
permet de sauver l’image dans le format .png avec la définition de l’écran.
Comme pour les axes, différents types de légendes (de 1 à 99) sont définis-
2.4 ODS sables. Ils spécifient positions et textes des libellés qui identifient les différents
graphismes et symboles utilisés.
Comme vu en introduction, le graphique peut être automatiquement orienté
dans une fichier au bon format (.rtf, .html...) et sans doute avec une meilleure 3.3 Symboles
définition que celle de l’écran.
Les différents types de symboles (1 à 99) sont définis afin de décrire les
ODS RTF BODY=’[Link]’; modes de représentation recherchés. Sont concernés : le symbole (forme,
ODS GRAPHICS ON; taille, couleur) utilisé pour représenter un point, le type de lignes reliant les
/* Programme SAS */ points (couleur, continue, hachurée, pointillée,. . . ), la façon ou mode d’inter-
ODS GRAPHICS OFF; polation incluant barres, boîtes à moustaches, escaliers, splines, intervalles de
 3 Tutoriel 3 : Grahiques avec SAS/GRAPH

confiance, régression (linéaire, polynomiale, spline). Il est important de noter que chaque paramètre peut être initialisé ou redéfini
à différents endroits d’un programme SAS : dans les commandes spécifiques
symbol1 interpol=sm50s /* lissage spline */ (symbol, legend, axes, pattern, title, footnote), par la
value=diamond /* symbole */ commande goptions et dans chacune des procédures. Ceci impose de bien
height=3 /* taille du symb.*/ distinguer les paramètres globaux, applicables à tous les graphes, des para-
width=2; /* epaisseur */ mètres spécifiques à chaque graphe.
Réaliser les graphes suivants en exécutant les commande. Se reporter à l’an-
3.4 Options graphiques nexe pour expliciter et comprendre la syntaxe de chaque procédure.
Outre ceux décrits ci-dessus (hsize, vsize), cette commande redéfinit
les valeurs de plus de 80 paramètres affectant 4 Quantitatif
• les différents aspects du graphique :
– border cadre autour du graphique, Graphes usuels pour des variables quantitatives discrètes puis continues.
– gunit=cm|in|pct unité de mesure,
– rotate=landscape|portrait orientation du graphique, 4.1 Discret
• le texte : Création de la table
– ftext police du texte,
– ftitle police des titres, Lire les données de répartition en âge, à chaque âge est associé un effectif.
• texte, symboles, types de hachures, légendes.
• Les paramètres reprennent leurs valeurs par défaut à la suite de : data age;
– reset=all|global all concerne tous les paramètres tandis que que global input age effectif eff_cum;
n’affecte pas ceux définis dans la même commande. cards;
24 1 1
3.5 Titres et notes 26 2 3
29 3 6
Les commandes title et footnote définissent des lignes de texte au- 31 2 8
tour du graphique, elles suivent le même principe que celui décrit au para- 33 4 12
graphe I.1.4 et d’autres options sont disponibles : taille, couleur et police des 37 2 14
caractères, position, rotations de la ligne de texte et des caractéres, tracés de 38 4 18
lignes. 41 3 21
43 3 24
45 1 25
goptions reset=global gunit=pct border
46 6 31
ftext=swissb htext=3; 49 3 34
title1 height=5 ’Institut’; 50 1 35
title2 ’de’; 52 3 38
title3 height=5 ’Mathématiques’; 57 5 43
footnote1 font=script justify=left 59 2 45
’Universite de Toulouse’; 60 2 47
 4 Tutoriel 3 : Grahiques avec SAS/GRAPH

62 1 48 4.2 Continu
;
/* place impérative du ";" */ Répartition des exploitations agricoles par classes de surface agricole utile
run; (SAU). Repérer comment sont organisées les données, attention, les classes ne
proc print;run; sont pas d’amplitudes égales.
Fonction de répartition
Diagramme en bâton
data exploit;
Le graphe légitime est un diagramme bâton, pas un histogramme qui corres- input SAU dens_eff freq_cum;
pond à une variable continue. cards;
0 0 0
proc gplot data=age; 5 4.8 0.24
axis1 label=("Age" justify=right); 10 2.18 0.349
axis2 label=("Effectif") 20 1.78 0.527
order=(0 to 6 by 1) offset=(0,); 35 1.35 0.73
symbol1 interpol=needle value=dot; 50 0.68 0.832
plot effectif*age / haxis=axis1 200 0.112 1.
vaxis=axis2 hminor=4 vminor=0; ;
run; run;
quit; proc print;run;

proc gplot data=exploit;

De même, comme la variables est discrète, le diagramme cumulatif qui es- axis1 label=("SAU(ha)" justify=right)
time la fonciton de répartition n’est pas continu. order=(0 to 200 by 50) ;
axis2 label=("frequences" justify=right "Cumulees")
Diagramme cumulatif order=(0 to 1 by 0.5) offset=(0,);
symbol1 interpol=join value=dot;
proc gplot data=age; plot freq_cum*sau / haxis=axis1 vaxis=axis2
axis1 label=("Age" justify=right) hminor=0 vminor=4;
order=(20 to 70 by 10); run;
axis2 label=("Effectif" justify=right "Cumule") quit;
order=(0 to 50 by 10) offset=(0,);
symbol1 interpol=step value=none; Histogramme
plot eff_cum*age /
proc gplot data=exploit;
haxis=axis1 vaxis=axis2
axis1 label=("SAU(ha)" justify=right)
hminor=4 vminor=0;
order=(0 to 200 by 50);
run;
axis2 label=("densite de" justify=right "frequence")
quit;
length=6cm order=(0 to 5 by 1) offset=(0,);
pattern1 value=msolid ;
 5 Tutoriel 3 : Grahiques avec SAS/GRAPH

symbol1 interpol=steprj value=none; 5.3 Diagramme en colonne

plot dens_eff*sau / haxis=axis1 vaxis=axis2
hminor=4 vminor=0 areas=1; Les graphiques sont en couleur par défaut :
run; axis3 label=none ;
goptions reset=all; axis4 length=1cm;
quit; proc gchart data=csp;
hbar classe / noaxis subgroup=csp raxis=axis3
gaxis=axis4 sumvar=eff nostats;
5 Qualitatif run;
quit;
Graphiques destinées aux variables qualitatives. goptions reset=all;
5.1 Création de la table Graphique en noir pour impression avec hachurage automatique :
Répartition des catégories socio-professionnelles parmi les actifs. goptions colors=(black);
axis3 label=none ;
data csp;
axis4 length=1cm;
input csp $ 1-10 eff classe;
proc gchart data=csp;
cards;
hbar classe / noaxis subgroup=csp raxis=axis3
ar et com 1739 1
gaxis=axis4 sumvar=eff nostats;
agricult 1312 1
run;
cadres 2267 1
quit;
p inter 4327 1
goptions reset=all;
employes 5815 1
ouvriers 6049 1
; 5.4 Diagramme en secteur
run; proc print;
Le graphique est vide par défaut. on peut lui rajouter des hachures :
pattern1 v=p2n0 ;
5.2 Diagramme en barres pattern5 v=p2n90;
pattern3 v=P2x45;
axis2 label=("Effectif") pattern4 v=P4n0;
order=(0 to 8000 by 2000); pattern2 v=P4n90;
proc gchart data=csp; pattern6 v=P4x90;
vbar csp / sumvar=eff raxis=axis2 goptions colors=(black);
vminor=1; proc gchart data=csp;
run; pie csp /sumvar=eff noheading slice=outside
quit; percent=outside;
goptions reset=all; run;
quit;
Améliorer l’axe des ordonnées en changeant la commande order. goptions reset=all;
 6 Tutoriel 3 : Grahiques avec SAS/GRAPH

ou des couleurs : proc boxplot data=pento;

pattern1 v=solid; axis1 length=5cm ;

proc gchart data=csp; axis2 length=7cm;
pie csp /sumvar=eff noheading slice=outside plot rythme*facteur/notches vaxis=axis2 vminor=1
percent=outside; haxis=axis1 ;
run; run;
quit; goptions reset=all;
goptions reset=all; quit;

Un peu plus de travail permettrait de sélectionner des couleurs plus équili-

brées. Il est évidemment possible mais pas recommandé de remplacer pie par
6.2 Profils
pie3d. Voir aussi les graphiques de type “mosaïque” dans sas/insight. Les données
décrivent l’âge d’obtention du bac, la durée pour obtenir le DEUG (L2) et un
6 Bidimensionnel pourcentage calculé pour représenter des profils. Il est dans ce cas évidemment
plus simple d’utilise rle mosaÏc plot de SAS/Insght.
6.1 Boîtes parallèles data deug;
infile "[Link]";
Les données sont issues d’une expérience clinique. Le rythme cardiaque est input age $ duree $ prof;
mesuré en fonction de la concentration (facteur à 3 niveaux) d’une molécule run;
en cours de test. proc print;run;
axis1 length=5cm label=("Age au BAC");
data pento;
axis2 length=5cm label=("Pourcentage") /*
infile "[Link]";
order=(0 to 100 by 20)*/;
input code $ rythme facteur $;
goptions colors=(black);
run;
proc gchart data=deug;
proc print;run;
vbar age / subgroup=duree sumvar=prof
maxis=axis1 raxis=axis2 vminor=1
proc gplot data=pento;
midpoints="moyen" "inf18" "18ans"
symbol1 interpol=box;
"19ans" "sup19";
axis1 length=5cm offset=(1cm,1cm) ;
run;
axis2 length=5cm;
quit;
plot rythme*facteur=1/vaxis=axis2 vminor=1
goptions reset=all;
haxis=axis1 ;
run;
Refaire le même graphique en couleur.
goptions reset=all;
quit; 6.3 Nuage
La procédure boxplot permet de contrôler les paramètres graphiques des Il s’agit de représenter l’évolution du chiffre d’affaire des entreprises en
boîtes avec de très nombreuses options. fonction du nombre de salarié.
 7 Tutoriel 3 : Grahiques avec SAS/GRAPH

data entr; A quoi sert l’option group= ? Remplacez type=sum par type=mean. Es-
infile "[Link]"; sayer l’option patternid=group placée après hbar typeen. Tester de
input code $ nb ef ca; même les options gspace=5, noaxos, nostats, descending.
run;
proc print;run; Création d’un graphe représentant la production cumulée d’énergie pour
chaque année :
proc gplot data=entr;
goptions reset=global cback=white htitle=3
symbol1 interpol=r value=dot;
axis1 label=("Chiffre d affaire" justify=right);
htext=1 hpos=100 vpos=45;
axis2 label=("Effectif") offset=(0,); title1 "Production d""énergie 2001 à 2003";
plot ca*ef=1 / haxis=axis1 vaxis=axis2; pattern1 value=solid color=yellow;
run; pattern2 value=solid color=blue;
quit; pattern3 value=x3 color=green;
goptions reset=all; pattern4 value=solid color=red;
pattern5 value=x2 color=blue;
pattern6 value=solid color=green;
7 Énergies renouvelables axis1 label=("Energie produite (*)");
footnote1 justify=left "(*) exprimee en Kwh"
Le fichier de données [Link] contient la production d’énergie
justify=right;
d’origine renouvelable en France de 2001 à 2003, exprimée en Gwh. Les diffé-
axis2 label=none;
rentes sources sont l’énergie hydraulique, solaire, éolienne, l’énergie issue des
legend1 label=(position=(topo left)
déchets urbains solides, du bois et des déchets du bois, les biogaz (source ; Mi-
"Type d" justify=left "energie")
nistère de l’Économie, des Finances et de l’Industrie, 2004). On dispose donc
value=("biogaz" "bois" "dechets urbains solides"
en tout de trois variables annee, typeen et prod. Exécuter et commenter
"eolienne" "hydraulique" "solaire");
le code ci-dessous le plus précisément possible :
proc gchart data=enprod;
data enprod; vbar annee / sumvar=prod sum discrete raxis=axis1
infile "[Link]"; maxis=axis2 space=3 width=7
input annee typeen $ prod; subgroup=typeen legend=legend1;
proc print;run; run;
quit;
goptions reset=global vpos=45 goptions reset=all;
htitle=2 htext=1 hpos=100;
title1 "Production d energie 2001 a 2003"; A quoi sert l’option sum ?
proc gchart data=enprod; L’énergie hydraulique prend trop de place, celle-ci est retirée du graphique.
hbar3d typeen / sumvar=prod type=sum group=annee;
run; data enprod2;
quit; set enprod;
if typeen="hydrau" then delete;
 8 Tutoriel 3 : Grahiques avec SAS/GRAPH

run; <options statistiques> ;

proc format; star liste de variables
value $typeen "biogaz"="Biogaz" "eolien"="Eolienne" </<options d’apparence>
"sol"="Solaire" "dechsol"="Dechets solides" <options statistiques> > ;
"bois"="Bois";
run;
proc gchart date=enprod2;
pie typeen / sumvar=prod type=sum other=5 Options générales
coutline=black noheading group=annee across=2; • data=table sas indique le nom de la table ou, par défaut, la dernière crée,
format typeen $typeen.; • annotate= table contenant les compléments graphiques.
where annee in(2001 2003);
run; Options d’apparence
goptions reset=all; Elles spécifient les couleurs, les espacements et largeurs de colonnes ou
barres. Il est également possible d’adjoindre un cadre (frame), de suppri-
mer (nolegend) ou modifier la légende. Une option annotate peut être
Annexes : Syntaxe des commandes introduite au niveau de chaque commande.

procédure “gchart” Options statistiques

• sumvar= variable quantitative dont le cumul ou la moyenne est repré-
Cette procédure trace des diagrammes en barres (hbar), en colonnes et his-
senté,
togrammes (vbar), en secteurs (pie) et aréolaires (star). Elle peut traiter
• freq= variable de pondération des observations,
des variables quantitatives ou qualitatives ; les variables quantitatives sont co-
• midpoints= liste des bornes de classes,
dées explicitement ou automatiquement en classes ou, selon les besoins, som-
• levels= nombre de classes,
mées ou moyennées.
• type= spécifie ce que représente le graphique (par défaut une fréquence) :
Syntaxe cfreq (fréquence cumulée), cpt (pourcentage cumulé), pct (pourcen-
tage), sum ou mean (associées à sumvar=).
proc gchart <options générales> ; • group= représentation de plusieurs graphes côte à côte suivant les moda-
by <descending> variable ; lités de la variable spécifiée (hbar ou vbar),
vbar liste de variables • subgroup= découpage des barres ou colonnes selon la participation des
</<options d’apparence> modalités de la variable spécifiée (hbar ou vbar).
<options statistiques> <options d’axes> > > ;
hbar liste de variables Options d’axes
</<options d’apparence>
Deux options permettent de de définir les axes ou de leur assigner des décla-
<options statistiques>
rations antérieurs : gaxis=axisn pour l’axe des groupes et maxis=axisn
<options d’axes> > ;
pour celui des bornes où n caractérise la définition d’axe concernée (cf. para-
pie liste de variables
graphe V.2.1.).
</<options d’apparence>
 9 Tutoriel 3 : Grahiques avec SAS/GRAPH

procédure “gplot” plot liste des graphes sous la forme : y*x<=n|variable>, avec la même
syntaxe que précédemment pour désigner plusieurs graphes (a*(a b),. . . ).
Graphiques en haute résolution de nuages de points en deux dimensions. La variable y fournit les ordonnées et x les abscisses des points représen-
Syntaxe tés par des symboles définis dans la commande symboln ou par diffé-
rents symboles selon les valeurs de la variable spécifiée qui induit une
proc gplot <options générales> ; classification. Dans ce dernier cas, une légende est créée par défaut.
by <descending> variable ;
bubble liste des graphes sous la forme : y*x = size où size est une va-
plot liste de graphiques
riable indiquant la taille des bulles à tracer autour des centres de coordon-
</ < annotate=data-set >
nées x et y.
< options d’apparence>
< options d’axes>> ;
Annotate data set
bubble liste de graphiques
</ < annotate=data-set > Une table d’annotations, définie lors d’une étape data, est une table SAS
< options d’apparence> contenant les descriptifs d’un ensembles de graphiques qui viendront se
< options d’axes>> ; superposer aux résultats des procédures précédemment décrites (gchart,
gplot,...). Il est alors possible de positionner tout libellé ou toute figure
géométrique simple et ainsi de personnaliser ses graphiques.
Options générales Des applications immédiates sont, par exemple, la production de plans fac-
• data=table sas indique le nom de la table ou, par défaut, la dernière crée, toriels avec identifications des points (variables, individus, modalités) par des
• annotate=table sas table contenant les compléments graphiques. libellés explicites ou encore le tracé du cercle des corrélations en analyse en
• uniform impose les mêmes échelles aux axes des différents graphiques. composantes principales.
Par principe, chaque ligne ou “observation” d’une table d’annotations est
Options d’apparence
une commande de réalisation d’un graphique particulier. Les valeurs de cha-
Elles spécifient les couleurs, les polices de caractères, les tailles des bulles cune des “variables” spécifient comment réaliser ce graphique : type, empla-
(bubble), le hachurage d’aires, la définition de légendes, la superposition cement, couleur,. . . . Les variables de la table d’annotation ont des noms pré-
(plot). définis ; les plus usuelles sont :
Options d’axes function indique ce qu’il faut tracé : bar, draw, frame, pie, symbol, label,. . . ,
x positionnement en abscisses,
Deux options permettent de définir les axes ou de leur assigner des déclara-
tions antérieurs : vaxis=axisn pour l’axe vertical et haxis=axisn pour y positionnement en ordonnées,
l’axe horizontal où n caractérise la définition d’axe concernée (cf. paragraphe size hauteur des caractères,
5.2.1.). De plus, frame trace un cadre tandis que noaxis supprime les axes. xsys unité de mesure des abscisses,
Commandes ysys unité de mesure des ordonnées,
by suivi du nom d’une variable qualitative indique que les graphiques sont hsys unité de mesure des hauteurs,
tracés par groupe d’observations ; la table doit être triée. color couleur,
 10 Tutoriel 3 : Grahiques avec SAS/GRAPH

position d’un texte par rapport aux coordonnées (calé à gauche, centré,. . . ),
line type de ligne (par défaut, continue),
text texte du libellé
style police de caractères.
La mise en œuvre de ces fonctionnalités est un peu fastidieuse mais c’est
la seule façon de faire éditer par SAS certains types de graphiques dont les
fameux plans factoriels avec les libellés en clair de tous les points.
Création d’une table d’annotations :

data annocomp;
set outcomp;
x = prin1;
y = prin2;
xsys= ’2’;
ysys= ’3’;
text= lib_ind;
size= 0.8;
label x = ’axe1’;
label y = ’axe2’;
keep x y text xsys ysys size;
run;