Proc Gpot: affichage de plusieurs graphs [Résolu]

[AstridG]

Bonjour,

Je voudrais afficher sur un seul graphique les courbes de valeurs prédites d'une variable mesurée sur 2 stations différentes entourées de leurs intervales de confiance à 95%. J'ai donc 3 instructions dans le statement plot (les variables x et y et l'intervalle de confiance inférieur et le supérieur).

Le problème étant que si j'utilise un by Station, j'ai deux graphs en sortie, et si j'utilise = station dans le statement plot, j'ai 3 graphs en sortie. Or je voudrais les courbes pour chaque station entourées de leurs intervalles de confiance respectifs sur le même graph.

Voici mon code:

Code :

1
2
3
4
5
6
7
 
proc gplot DATA= DB;
	BY Station;
	plot  (UpperIC LowerIC) * Var1 Var2 * Var1
               / overlay vaxis= axis1 haxis= axis2 Legend= Legend1;
run;
quit;

Merci!

Astrid

[olivier.decourt]

Bonjour Astrid.
Et avec

Code :

plot  (UpperIC LowerIC var2) * Var1 = station

ça ne fonctionne pas ?
Olivier

[AstridG]

Salut Olivier,

Non, c'est pareil, il sort 3 graphs. C'est énerant, il doit bien y avoir un moyen...

Merci pour ta réponse en tout cas!

Astrid

[Manoutz]

une esquisse de solution..

je plus gros de ton problème semble pris en compte, reste a jouer avec les options graphiques....

Code :

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
proc format;
value station 
	1="station 1"
	2="station 2"
;
run;
quit;
DATA station;
do station=1 TO 2;
		do i =1 TO 5;	
			varlow=-10*(rannorm(-3)+i );
			var=rannorm(-3)+i ;
			varhigh=10*(rannorm(-3)+i );
			output;
		end;
end;
format station station.;
run;
 
 
proc gplot DATA= station;
	plot  (varlow  var varhigh) * station
               / overlay vaxis= axis1 haxis= axis2 Legend= Legend1;
run;
quit;

reste que à voir si les données que j'ai simulé collent avec ta problématique... ce qui n'est pas forcément gagné

[AstridG]

Merci Manoutz,

Malheureusement ça ne répond pas tout à fait à mon problème parce que je veux faire les graphs par rapport à une variable Var1 définie. Or dans ta proposition les plots sont par rapport à la Station en X...

Astrid

[Manoutz]

ouais c'est bien ce que je me disais il manque une variable en considération...

j'ai également pensé "regrouper" deux variables avec plusieurs catégories (var, low, normal, high et station, 1 et 2) en 6 variables :

Code :

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
DATA station;
		do temps =1 TO 5;	
			varlowst1=-10*(rannorm(-3)+temps );
			varst1=rannorm(-3)+temps;
			varhighst1=10*(rannorm(-3)+temps);
			varlowst2=-10*(rannorm(-3)+temps );
			varst2=rannorm(-3)+temps;
			varhighst2=10*(rannorm(-3)+temps);
			output;
		end;
 
run;
 
 
proc gplot DATA= station;
	plot  (varlowst1  varst1 varhighst1 varlowst2  varst2 varhighst2)  *temps
               / overlay vaxis= axis1 haxis= axis2 Legend= Legend1;
run;
quit;

mais la encore je suis pas sur d'être dans le juste. Et puis il faut dire ce qui est, c'est pas super joli...

[AstridG]

Ça c'est l'idée qui m'a sauvée!!

C'est peut-être pas "super élégant", mais c'est malin et ça fait drôlement bien l'affaire! Puis sur le graph en sortie, on voit pas les bricolages qui sont derrière...

Merci millle fois!

Bonne soirée,

Astrid