Nuage direct et dual d'une ACP : MACRO SAS
Par Haache

Le , par Haache, Membre éprouvé
Encore une macro que je vous propose !
Celle ci fait le nuage des individus et celle des variables, elle trace le cercle de corrélation. Vous pouvez jouer sur les couleurs et tout! La macro est adaptée à la sortie ACP de la macro de l'INSEE. En fait, l'INSEE a créé une MACRO qui permet de réaliser l'ACP. Donc c'est la sortie des résultats que vous allez utiliser pour faire les nuages avec ma macro. Je vous suggère de télécharger la macro ACP de l'INSEE pour faire vos ACP. l'INSEE a aussi proposé une macro pour faire les nuages. Mais tester la mienne aussi pour voir.

La macro de l'INSEE qui permet de réaliser l'ACP est est téléchargeable ici : http://www.insee.fr/fr/methodes/defa...il_analyse.htm

Code : Sélectionner tout
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
%macro Plotacp(DATA= ,AXEH= ,AXEV= ,COND= ,CERCLE=N,BAREA=N)/ DES="Réprésentation des individus et variable en ACP";/*Je dupplique la table dans une autre _Graph en supprimant les données manquantes*/ 
DATA _Graph_; 
set &DATA; 
if not Missing(coord1); 
run; 
/*J'applique la condition. les observations et les variables que je veux afficher*/ 
DATA _Graph_; 
SET _Graph_; 
%if %length(&COND)^=0 %then %do; 
&COND; 
%end; 
run; 
/*Construction de la table Annotate*/ 
data ACPANNOT; 
length function $10.; 
set _Graph_; 
function="Label"; 
run; 
data ACPANNOT; 
length color $10.; 
length text $10.; 
set ACPANNOT; 
/*Les noms à écrire pour les observations et les variables dans le graphique*/ 
    text=_LABEL_; 
    position="3"; 
    x=coord&AXEH; 
    y=coord&AXEV; 
    xsys='2'; 
    ysys='2'; 
    if _TYPE_="OBSACT" then color="purple"; 
    else if _TYPE_="OBSSUP" then color="black"; 
    else if _TYPE_="VARACT" then color="blue"; 
    else color="red"; 
    if _TYPE_="OBSSUP" then size=1.7; 
    else size=0.9; 
    Output; 
/*Pour tracer les segments des variables. Les segments sont tracées du centre du cercle vers les points*/ 
    if _TYPE_ IN ("VARACT","VARSUP") then do;  
    x=0; 
    y=0; 
    function="Move"; 
    if _TYPE_="OBSACT" then color="blue"; 
    else if _TYPE_="OBSSUP" then color="red"; 
    else if _TYPE_="VARACT" then color="blue"; 
    else color="red"; 
    size=0.8; 
    when="B"; 
    end; 
    output; 
 
 
    if _TYPE_ IN ("VARACT","VARSUP") then do; 
    x=coord&AXEH; 
    y=coord&AXEV; 
    function="Draw"; 
    if _TYPE_="OBSACT" then color="blue"; 
    else if _TYPE_="OBSSUP" then color="red"; 
    else if _TYPE_="VARACT" then color="blue"; 
    else color="red"; 
    size=0.8; 
    when="A"; 
    end; 
output; 
/*Pour mettre les symboles*/ 
function="Symbol"; 
position="5"; 
when="B"; 
x=coord&AXEH; 
y=coord&AXEV; 
xsys='2'; 
ysys='2'; 
if _TYPE_="OBSACT" then text="triangle"; 
else if _TYPE_="OBSSUP" then text="dot"; 
else if _TYPE_="VARACT" then text="circle"; 
else text="circle"; 
if _TYPE_="OBSACT" then size=0.01*(co2_&AXEH+co2_&AXEV); 
else if _TYPE_="OBSSUP" then size=0.03*(co2_&AXEH+co2_&AXEV); 
else if _TYPE_="VARACT" then size=0.035*(co2_&AXEH+co2_&AXEV); 
else size=0.035*(co2_&AXEH+co2_&AXEV); 
if _TYPE_="OBSACT" then color="purple"; 
else if _TYPE_="OBSSUP" then color="yellow"; 
else if _TYPE_="VARACT" then color="blue"; 
else color="red"; 
output; 
run;  
 
 
%if &CERCLE=O %then %do; 
data ACPANNOT; 
set ACPANNOT end=fin; 
output ACPANNOT; 
if fin then do; 
    x=1; 
    y=0; 
    function="Move"; 
    color='black'; 
    size=0.2; 
    when="A"; 
output; 
do angle=0 to 44/7 by 0.01; 
        x=cos(angle); 
        y=sin(angle); 
        function="Draw"; 
        When="A"; 
        color='black';     
    output; 
end; 
output; 
    x=1; 
    y=0; 
    function="Draw"; 
    color='black'; 
    when="A"; 
output; 
end; 
run;  
%end;  
options nodate;     
goptions reset=all gunit=pct htext=2; 
%if &BAREA=O %then %do; 
goptions hsize=6 vsize= 6; 
%end; 
%else %do; 
goptions hsize=6 vsize= 6; 
%end; 
axis1 minor=none offset=(4,4) major=(h=0.1 w=1) label=("Axe&AXEH"); 
axis2  minor=none offset=(4,4) major=(h=0.1 w=1) label=("Axe&AXEV"); 
symbol v=none; 
proc gplot data=ACPANNOT; 
plot y*x=1/frame cframe=white href=0 vref=0 overlay annotate=ACPANNOT grid haxis=axis1 vaxis=axis2; 
run; 
quit; 
proc datasets lib=Work memtype=Data nolist; 
%if &CERCLE=O %then %do; 
delete Acpannot Acpannotm Acpannotn _graph_; 
%end; 
run; 
%Put Fin de macro; 
%put Ecris moi si tu rencontres de difficultés: ariel92and@gmail.com; 
%mend Plotacp;
Lorsque vous réalisé votre ACP avec la macro de l'INSEE voici mon code pour faire le graphique.

%Plotacp(DATA= ,AXEH= ,AXEV= ,COND= ,CERCLE=,BAREA=)

Dans DATA, vous spécifiez le nom de la table sortie que vous avez mise dans la macro de l'insee (ce que vous avez mis dans Out=).
Axev prend le numéro de l'axe que vous voulez mettre en vertical. soit 1, 2, 3 ....
Axeh prend le numéro de l'axe que vous voulez mettre en horizontal
Cond: Cette macro variable est très importante. Peut être que vous ne voulez pas afficher tous les points. Si vous voulez afficher les individus actifs seuls, vous allez mettre Cond= if _type_ in ("OBSACT"). Si c'est les variables supplémentaires seules, vous allez mettre Cond= if _type_ in ("VARSUP"). Si c'est les individus actifs et supplémentaires, vous allez mettre Cond= if _type_ in ("OBSACT" "OBSSUP").... A vous de voir les valeurs prises par la variable _type_ qui est dans la sortie pour pouvoir mettre ce qui est convenable
Dans la variable macro Cercle qui prend les valeurs O et N (O=OUI, N=NON), vous allez dire si vous voulez tracer le cercle de corrélation ou non
Dans la variable macro BAREA vous allez mettre O si vous voulez que SAS affiche un grand graphique. N si non.

Je vous laisse le code pour que vous pussiez modifier. Jouer sur les couleurs des variables actives, supplémentaires, individus et autres.
Voici des graphiques que j'ai sortis par la macro


Vous avez aimé cette actualité ? Alors partagez-la avec vos amis en cliquant sur les boutons ci-dessous :
Offres d'emploi IT
Ingénieur produit (FADEC militaire) H/F
Safran - Ile de France - 100 rue de Paris 91300 MASSY
Chef de projet technique H/F
Safran - Ile de France - Melun (77000)
Ingénieur développement logiciels temps réel embarqué H/F
Safran - Ile de France - Éragny (95610)

Voir plus d'offres Voir la carte des offres IT
Responsable bénévole de la rubrique SAS : fafabzh6 -