Discussion :

[étudiant11]

Bonjour ,
Est ce que c'est possible en présentant une ACM d'affecter à chaque variable une couleur différente.J'ai utilisé ce code pour avoir le graphique, simplement je l'ai eu en noir:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

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
%macro ACM(var,base,titre);
PROC corresp DATA=&base DIMENS=5  OUTC=ACM  MCA  ;
TABLES &var  ;
QUIT ;                       
DATA graphe;
   SET ACM;
   x = dim1;
   y = dim2;
   xsys = '2';
   ysys = '2';
   xmin=-1;
  xmax=1;
  ymax=1;
  ymin=-1;
   text = _name_;
   size = 0.8;
   LABEL x = 'dimension 1';
   LABEL y = 'dimension 2';
   KEEP x y text xsys ysys size;
RUN;
Proc print data=graphe ; RUN ; 
PROC GPLOT DATA = graphe;
TITLE &titre;
  PLOT y*x / ANNOTATE = graphe href=0 vref=0 haxis=-2 to 2 vaxis=-2 to 2;
  SYMBOL1 v=none;
RUN;
quit;
%mend;

Qu'est ce que je dois ajouter s'il vous plait à ce code pour différencier les variables en couleur.

[Haache]

C'est bien possible. Ton annotate est incomplète. Je ne vois pas la variable color ni la variable function . Avec color tu peux utiliser un if pour dire si telle variable alors color=yellow sinon color = red. J'espère que tu comprends

[étudiant11]

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
%macro PLOTCOR(DATA= ,ID= ,AXEH= ,AXEV= ,COND= )/ STORE DES="Réprésentation des individus et variable en AFC ou en ACM";
DATA _Graph_;
set &DATA;

Pour cette macro je n'ai pas compris qu'est ce que je dois mettre pour ID et COND??

[Haache]

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
%macro PLOTCOR(DATA= ,ID= ,AXEH= ,AXEV= ,COND= )/ STORE DES="Réprésentation des individus et variable en AFC ou en ACM";
DATA _Graph_;
set &DATA;

Pour cette macro je n'ai pas compris qu'est ce que je dois mettre pour ID et COND??

Attention! Ne pas confondre le code source de la macro et ton code.
Le premier code que j'ai mis qui va de

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

%macro PLOTCOR(DATA= ,ID= ,AXEH= ,AXEV= ,COND= )

jusqu'à

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

%mend;

est à compiler sans rien renseigner. Tu ne dois rien mettre dans DATA= ID= etc.... c'est pour créer la macro...
Si tu compiles cela, tu vas mettre ce code pour réaliser tes tâches

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

%PLOTCOR(DATA= ,ID= ,AXEH= ,AXEV= ,COND= )

C'est ici que tu renseignes tout.

Pour comprendre chaque macro variables ID et COND, j'ai l'impression que tu n'as pas lu tout le blog sinon tout est expliqué. En bas du graphique il y a "En Lire plus" il faut cliquer pour lire tout le texte. Comme c'est un long texte, c'est juste une partie qui est affiché et après ça il y a ....

Teste et confirme si ça marche

[étudiant11]

Merci de m'avoir répondu.
J'ai mis le code

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

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
%macro PLOTCOR(DATA= , AXEH= ,AXEV= ,COND= )/ STORE DES="Réprésentation des individus et variable en AFC ou en ACM";
DATA _Graph_;
set &DATA;
%if %length(&COND)^=0 %then;
&COND;
%end;
run;
 
 
data CORANNOT;
length color $10.;
length text $10.;
length function $10.;
set _Graph_;
function="Label";
position="1";
when="B";
x=Dim&AXEH;
y=Dim&AXEV;
xsys='2';
ysys='2';
if _TYPE_="OBS" then color="blue";
else if _TYPE_="VAR" then color="red";
else if _TYPE_="VARSUP" then color="green";
size=0.8;
Output;
function="Symbol";
position="5";
when="B";
x=Dim&AXEH;
y=Dim&AXEV;
xsys='2';
ysys='2';
text="triangle";
size=3*(SqCos&AXEH+SqCos&AXEV);
if _TYPE_="OBS" then color="blue";
else if _TYPE_="VAR" then color="red";
else if _TYPE_="SUPVAR" then color="green";
output;
run; 
 
 
 
 
goptions reset=all gunit=pct hsize=7 vsize= 7;
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=CORANNOT;
plot y*x=1/frame cframe=white href=0 vref=0 overlay annotate=CORANNOT grid haxis=axis1 vaxis=axis2;
run;
quit;
proc datasets lib=work memtype=DATA  nolist;
delete _Graph_ CORANNOT;
run; 
%mend Plotcor;

Le seul changement que j'ai effectué est de ne pas mettre "ID"
voila le message d'erreur que j'ai eu en faisant l'appel de la macro

[Haache]

Mais tu as touché au code source de la macro?????? C'est normal que ça ne marche pas. Le code de la macro est à exécuter sans retrancher ni augmenter une lettre

[étudiant11]

Je n'ai pas compris à quoi ça sert ce ID c'est pour cette raison que je l'ai pas mis.

[Haache]

Je n'ai pas compris à quoi ça sert se ID c'est pour cette raison que je l'ai pas mis.

Parfait!!! Allons-y pas-à-pas.
Lorsque tu as fait l'ACM avec proc corresp tu as mis une table de sortie dans OUT. Regarde dans le table et liste moi, s'il te plait, les noms des variables de cette table

[étudiant11]

voila le code que j'ai utilisé

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
PROC corresp DATA=Croix_q18_var_q19_var DIMENS=5 OUT=ACM1 MCA ;
tables col1 col2;
QUIT ;

et voila les deux premières colonnes de ma table sortie

[Haache]

OK Donc ce qui permet d'identifier les variables est _name_. Quand tu vas appeler la macro, tu mets ID=_name_.

[étudiant11]

ok merci
Mais j'ai eu le même message d'erreur

[Haache]

Ok il y a des erreurs. C'est justement ceux pourquoi je voulais qu'on teste. Chez moi, les macro ne sont pas Work. C'est pour quoi j'ai mis /Stored à la première ligne du code source. Mais toi, tu dois l'effacer.

Par ailleurs, à la quatrième ligne, il ya %do que j'ai oublié %if %length(&COND)^=0 %then %do;

Voici le code corrigé.

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

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
%macro PLOTCOR(DATA= , AXEH= ,AXEV= ,COND= )/ DES="Réprésentation des individus et variable en AFC ou en ACM";DATA _Graph_;
set &DATA;
%if %length(&COND)^=0 %then %do;
&COND;
%end;
run;
 
 
data CORANNOT;
length color $10.;
length text $10.;
length function $10.;
set _Graph_;
function="Label";
position="1";
when="B";
x=Dim&AXEH;
y=Dim&AXEV;
xsys='2';
ysys='2';
if _TYPE_="OBS" then color="blue";
else if _TYPE_="VAR" then color="red";
else if _TYPE_="VARSUP" then color="green";
size=0.8;
Output;
function="Symbol";
position="5";
when="B";
x=Dim&AXEH;
y=Dim&AXEV;
xsys='2';
ysys='2';
text="triangle";
size=3*(SqCos&AXEH+SqCos&AXEV);
if _TYPE_="OBS" then color="blue";
else if _TYPE_="VAR" then color="red";
else if _TYPE_="SUPVAR" then color="green";
output;
run; 
 
 
 
 
goptions reset=all gunit=pct hsize=7 vsize= 7;
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=CORANNOT;
plot y*x=1/frame cframe=white href=0 vref=0 overlay annotate=CORANNOT grid haxis=axis1 vaxis=axis2;
run;
quit;
proc datasets lib=work memtype=DATA  nolist;
delete _Graph_ CORANNOT;
run; 
%mend Plotcor;

Toutes mes excuses et je vais corriger dans le blog. C'est justement ces séances qui vont nous permettre d'améliorer la macro..
Teste et confirme si c'est bon.

Cordialement,

[étudiant11]

voila le résultat que j'ai eu

[Haache]

voila le résultat que j'ai eu

Tu n'as pas mis ID=_NAME_

[Haache]

On Avait mis le code source que tu as modifié.
C'est plutôt ça le code

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

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
%macro PLOTCOR(DATA= ,ID= ,AXEH= ,AXEV= ,COND= )/ DES="Réprésentation des individus et variable en AFC ou en ACM";DATA _Graph_;
set &DATA;
%if %length(&COND)^=0 %then %do;
&COND;
%end;
run;


data CORANNOT;
length color $10.;
length text $10.;
length function $10.;
set _Graph_;
function="Label";
position="1";
when="B";
x=Dim&AXEH;
y=Dim&AXEV;
xsys='2';
ysys='2';
text=&ID;
if _TYPE_="OBS" then color="blue";
else if _TYPE_="VAR" then color="red";
else if _TYPE_="VARSUP" then color="green";
size=0.8;
Output;
function="Symbol";
position="5";
when="B";
x=Dim&AXEH;
y=Dim&AXEV;
xsys='2';
ysys='2';
text="triangle";
size=3*(SqCos&AXEH+SqCos&AXEV);
if _TYPE_="OBS" then color="blue";
else if _TYPE_="VAR" then color="red";
else if _TYPE_="SUPVAR" then color="green";
output;
run; 




goptions reset=all gunit=pct hsize=7 vsize= 7;
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=CORANNOT;
plot y*x=1/frame cframe=white href=0 vref=0 overlay annotate=CORANNOT grid haxis=axis1 vaxis=axis2;
run;
quit;
proc datasets lib=work memtype=DATA  nolist;
delete _Graph_ CORANNOT;
run;  %mend Plotcor;

Pour l'appeler

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

%PLOTCOR(DATA= ,ID= ,AXEH= ,AXEV= ,COND= )

Cordialement,

[étudiant11]

Je l'ai ajouté et voila le résultat

[Haache]

C'est quand même bizarre, je viens d'essayer et ça a marché. Tu as vraiment des observations dans la table ou bien des variables? parce que tu as mis if type in ("OBS"). Si c'est des variables que tu as, tu doit mettre "VAR". Il faut voir dans la table les valeurs prises par la variable _type_.

"OBS" si c'est individu
"SUPOBS" Si c'est individu supplémentaire
"VAR" si c'est variable
"SUPVAR" si c'est variable supplémentaire

De plus, tu n'as pas fermé la parenthèse. tu as ouvert la parenthèse deux fois mais tu as fermé une fois. A la fin quand tu appelles la macro, tu dois fermer la parenthèse deux fois

[Haache]

Bonjour,
Où en es tu avec ton analyse?
Tu peux toujours poser des questions si tu n'est pas satisfaite.
Je voulais te dire aussi que axeh=2, axev=2 ne te donnera rien de bon mais tu auras des points sur la diagonale.
Normalement c'est axeh=1,axev=2 pour visualiser sur les deux premiers axes. EN tout cas, tu mets le numéro du premier axe que tu veux visualiser dans axeh et le second dans axev. Il ne faut pas mettre les mêmes numéros dans les deux.
Quand moi je fais une analyse et que je retiens 3 axes, je fais deux graphiques: le premier avec axeh=1,axev=2 et le second avec axeh=1,axev=3

[étudiant11]

Votre code m'a permis d'avoir le résultat mais le graphe que j'ai eu n'est pas trop clair et j n'ai pas de différents couleur pour chaque variable toutes les observations de toutes les variables sont en rouge

la macro que j'ai utilisée me donne aussi le résultat un peu plus clair mais tous en noir

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

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
 
%macro ACM(var,base,titre);
PROC corresp DATA=&base DIMENS=5  OUTC=ACM  MCA  ;
TABLES &var  ;
QUIT ;                       
DATA graphe;
   SET ACM;
   x = dim1;
   y = dim2;
   xsys = '2';
   ysys = '2';
   xmin=-1;
  xmax=1;
  ymax=1;
  ymin=-1;
   text = _name_;
   size = 1.2;
   LABEL x = 'dimension 1';
   LABEL y = 'dimension 2';
   KEEP x y text xsys ysys size;
RUN;
Proc print data=graphe ; RUN ; 
PROC GPLOT DATA = graphe;
TITLE &titre;
  PLOT y*x / ANNOTATE = graphe href=0 vref=0 haxis=-0.8 to 0.8 vaxis=-2 to 2;
  SYMBOL1 v=none;
RUN;
quit;
%mend;

[Haache]

Oui je vois.
C'est pour ça que j'ai laissé le code de la macro sinon j'aurais pu donner seulement les catalogues. C'est à toi de modifier les couleurs que j'ai mises au niveau de OBS VAR .... Moi j'ai mis red partout, mais toi en fonction de ce que tu veux, tu peux mettre green, yellow etc... Toujours dans mon blog tu verras que le nuage des individus que j'ai eu avec l'ACP a plusieurs couleurs.
Malheureusement, on ne peut pas créer une macro pour que SAS gère les couleurs automatiquement, c'est la seule possibilité. L'image ci-joint est issue du même code. Juste que j'ai changer les couleurs.

Bon courage