proc tabulate : insérer une colonne [Résolu]

[polpo56]

Bonjour,

J'ai réalisé un tableau à l'aide de la proc tabulate avec les percentiles 25 50 et 75 de ma variable selon plusieurs catégorie. Cela fonctionne correctement mais je voudrais insérer une colonne avec les données d'un seul établissement pour pouvoir comparer aux percentiles.
N'hésitez pas à me poser d'autre question si je ne suis pas claire.

Merci beaucoup à celui qui pourra m'aider.

[MEGAMIND2]

Bonjour Pauline,
ça veut dire quoi: insérer une colonne avec les données d'un seul établissement?

[polpo56]

Alors en fait j'ai dans mon tableau trois colonnes avec les percentiles calculés sur tous les établissements et je veux avoir une quatrième colonne avec les chiffres associés à un seul établissement. En gros, je veux sélectionné une seule observation.

[MEGAMIND2]

Tu pourrais posté ton code stp! ça sera plus facile de comprendre

[polpo56]

Voila mon code :

Code :

1
2
3
4
5
6
7
8
proc tabulate DATA=pub.etab ;
    class discip;
    var ratio ;
    TABLE (discip=""),
    (ratio="référentiel CH")*
    (p25 p50 p75)
    /box="discipline";
run;

[MEGAMIND2]

Un échantillon de ton fichier si possible!

[polpo56]

FI DISCI ratio
010008407 21000 1.4704957307
010008407 24000 0.5282430699
010780054 21000 1.0618104477
010780054 24000 0.4615013663
010780062 21000 1.5356513169
010780062 24000 0.5067534247

Voila l'echantillon ca sera plus simple a expliquer c'est vrai. Donc pour l'instant avec mon code j'ai obtenu les percentiles pour l'ensemble des établissements (FI) et je veux une colonne avec les ratios de l'éatblissement 010008407 par exemple pour chacune des disciplines.
Merci d'essayer de m'aider!

[MEGAMIND2]

Si tu croises Discip avec Fi tu auras toutes tes combinaisons possibles en ligne et en colonne t'auras tes percentiles. C'est pas ce que tu veux? si oui, il y a forcément un moyen de filtrer par la suite.

[polpo56]

Le problème c'est que j'ai plus de 300 étblissement donc je ne veux pas avoir toutes les combinaisons possibles j'aimerais pouvoir faire le filtre directement

[MEGAMIND2]

Envoyé par polpo56

Le problème c'est que j'ai plus de 300 étblissement donc je ne veux pas avoir toutes les combinaisons possibles j'aimerais pouvoir faire le filtre directement

Et c'est quoi ton critère pour filtrer? comment tu as choisi 010008407?

[polpo56]

Ca je ne sais pas encore je vais peut être faire un regroupement par département et en choisir un je ne sais pas trop mais une fois que je sais comment avoir cette colonne dans mon tableau je pourrais choisir n'importe quel établissement.

[MEGAMIND2]

J'ai compris maintenant ton problème, en fait il faudrait réaliser deux proc tabulate et fusionner les deux tables résultantes. sachant que ça reste manuel, ie c'est à toi de choisir l'établissement sur lequel tu veux faire des comparaisons.

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
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
DATA TEST;
INPUT FI$1-10 DISCIP $11-16 ratio 17-29;
CARDS;
010008407 21000 1.4704957307
010008407 24000 0.5282430699
010780054 21000 1.0618104477
010780054 24000 0.4615013663
010780062 21000 1.5356513169
010780062 24000 0.5067534247
;
RUN;
 
proc tabulate DATA=test out=total;
	class discip FI;
	var ratio ;
	TABLE (discip=""),
	(ratio="référentiel CH")*
	(p25 p50 p75)
	/box="discipline";
run; 
 
proc tabulate DATA=test out= FI_single(DROP=_TYPE_ _PAGE_ _TABLE_);
	class discip FI;
	var ratio ;
	TABLE (discip=""),
	(ratio="référentiel CH")
	/box="discipline";
	WHERE FI='010008407';
run; 
 
PROC SORT DATA=total;
	BY DISCIP;
RUN;
PROC SORT DATA=FI_single;
	BY DISCIP;
RUN;
DATA COMPARE;
	MERGE total FI_single;
	BY DISCIP;
RUN;

[polpo56]

Merci beaucoup je vais essayer ça ça devrait marcher!!!

[olivier.decourt]

Bonjour Pauline.
Si je comprends bien, tu veux pouvoir tester l'appartenance de telle spécialité de tel établissement à un quartile précis.
Dans ce cas il n'est pas forcément nécessaire de calculer explicitement les valeurs des quartiles et on peut demander à la proc RANK d'affecter automatiquement chaque observation à son quartile.

Code :

1
2
3
4
5
6
7
8
PROC SORT DATA=pub.etab ;
BY discip ;
RUN ;
PROC RANK DATA=pub.etab OUT=work.etab_quartiles GROUPS=4 ;
  BY discip;
  VAR ratio ; 
  RANKS quartile ; /* contient 0, 1, 2 ou 3 selon le quartile auquel appartient l'observation */
RUN ;

Bon courage.
Olivier

[polpo56]

En fait ça m'intéresse d'avoir dans le même tableau à la fois les quartiles et l'établissement étudié.

[MEGAMIND2]

La première solution ne te satisfait pas?

[polpo56]

Si c'est bon j'ai réussi à faire ce que je voulais!!!