Discussion :

[joyeux_lapin13]

Bonjour, j'essai de récupérer dans un output le warning de la proc FREQ quand il juge que la répartition n'est pas bonne pour le test du chi2.

Sur le site du support SAS ((http://support.sas.com/documentation...eq_sect010.htm )
ils proposent l'option WARN et donc la syntaxe suivante:

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
 
DATA TRY;
INPUT X Y;
CARDS;
1 1
1 1
0 0
1 0
;
RUN;
 
PROC FREQ DATA = TRY;
TABLES X * Y / CHISQ(WARN=OUTPUT) ;
OUTPUT OUT = ANALYSE FISHER CHISQ;
RUN;

Seulement ça ne marche pas... quelqu'un sait-il où cette option se place en fait?

En vous remerciant d'avance.

[MEGAMIND2]

Bonjour Julien,
Si tu veux tester la dépendance entre tes 2 variables et recupérer la proba (khi-2) à comparer au seuil de 5%:

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
DATA TRY;
INPUT X Y;
CARDS;
1 1
1 1
0 0
1 0
;
RUN;
ODS OUTPUT Chisq=CRAMER;
PROC FREQ DATA = TRY;
TABLES X * Y / CHISQ ;
RUN;
 
DATA CRAMER;
SET CRAMER;
WHERE Statistic='Khi-2';
rename Prob=Warning;
RUN;
 
DATA _NULL_;
SET CRAMER;
CALL SYMPUT("Warning",Warning);
RUN;
 
%PUT Warning=&Warning.;

[joyeux_lapin13]

Salut Brice, je pense m'être mal exprimé mais en fait ce que je cherche à faire c'est un test du chi2 ou de Fisher en fonction de ma répartition.

En outre j'arrive à récupérer les deux tests mais pour pouvoir décider si je vais choisir l'un ou l'autre il me faut pouvoir retourner le contenu du warning (vide si la répartition est ok, message

WARNING: 100% of the cells have expected counts less
than 5. Chi-Square may not be a valid test.

si elle ne l'est pas).

Ci-dessus j'ai indiqué la syntaxe que recommande SAS pour récupérer ce message d'erreur en OUTPUT, sauf que moi quand je m'en sers ça veut pas marcher:

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
 
14         GOPTIONS ACCESSIBLE;
15         PROC FREQ DATA = TRY;
16         TABLES X * Y / CHISQ(WARN=OUTPUT) ;
                               _
                               22
                               200
ERROR 22-322: Syntax error, expecting one of the following: ;, AGREE, ALL, ALPHA, BDT, BIN, BINC, BINOMIAL, BINOMIALC, CELLCHI2, 
              CHISQ, CL, CMH, CMH1, CMH2, CONTENTS, CONV, CONVERGE, CROSSLIST, CUMCOL, DEVIATION, EXACT, EXPECTED, FISHER, FORMAT, 
              JT, KAPPA, LIST, MAXITER, MEASURES, MISSING, MISSPRINT, NOCOL, NOCUM, NOFREQ, NOPERCENT, NOPRINT, NOROW, NOSPARSE, 
              NOWARN, OUT, OUTCUM, OUTEXPECT, OUTPCT, PDIFF, PLCORR, PLOTS, PRINTKWT, RDIFF, RELRISK, RISKDIFF, RISKDIFFC, SCORE, 
              SCORES, SCOROUT, SPARSE, TABLE, TESTF, TESTP, TOTPCT, TREND.  
ERROR 200-322: The symbol is not recognized and will be ignored.
17         OUTPUT OUT = ANALYSE FISHER CHISQ;
18         RUN;

Ce que je n'arrive pas à comprendre... à moins que cette commande soit propre à la toute dernière version de SAS (j'ai Entreprise Guide 4.2).

[MEGAMIND2]

C'est bizzare, le nowarn fonctionne mais pas warn
ça doit sûrement marcher en 9.2

[joyeux_lapin13]

C'est bien ce qui m'agace... surtout que j'aurais été curieux de voir le type d'information renvoyé par l'option, du genre s'il fournit également combien de cellule ont un nombre de valeurs inférieur à 5 pour vraiment opter pour un chi2 ou un fisher...

[MEGAMIND2]

Tu peux toutefois récupérer les effectifs du croisement des modalité de X avec Y et tu prends le max, tu stocks dans une macro variables et tu comparers cette valeur à 5 et si c'est > 5 -> fais un test de Khi-2 sinon Fischer ou Fisher je sais plus.

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
DATA TRY;
INPUT X Y;
CARDS;
1 1
1 1
0 0
1 0
;
RUN;
PROC FREQ DATA = TRY;
TABLES X * Y /chisq out=freq;
RUN;
 
proc transpose data=freq out=freq_t prefix=count;
var count;
run;
data freq_t;
set freq_t;
max=max (of count1-count3);
run;
DATA _NULL_;
SET freq_t;
call symput("max",max);
run;

[joyeux_lapin13]

LOL Brice, c'est ce que je viens de finir à l'instant tellement ça m'a saoulé:

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
 
/* Répartition des effectifs pour la variable en cours */
PROC FREQ DATA = ANALYSE;
TABLES Y * X/ LIST;
ODS OUTPUT LIST = FREQ;
RUN;
ODS OUTPUT CLOSE;
 
/* Etude des cellules pour voir si on va faire un test du chi2 ou de Fisher exact */
DATA FREQ;
SET FREQ;
IF (Frequency < 5) THEN Warning = 1;
ELSE Warning = 0;
CALL SYMPUT("Warning",Warning);
RUN;

Et avec un appel &Warning. dans un %IF je décide de si je garde le test du chi2 ou celui de Fisher.

L'avantage c'est qu'à long terme je vais pouvoir faire des scores avec mon warning pour mettre un seuil de cellule mal répartie et mieux décider le type de test.

Merci Brice pour ton intérêt à mon problème.

[MEGAMIND2]

JE me suis compliqué la vie avec une proc transpose et tout...
ta solution est plus sexy avec moins de lignes mais ça veut pas dire que la mienne est fausse

[olivier.decourt]

Je ne veux pas plomber l'ambiance, mais les Warnings sont déclenchés par la présence d'effectifs théoriques inférieurs à 5, pas d'effectifs observés. On devrait donc travailler avec l'option EXPECTED.

[joyeux_lapin13]

Et aprés il dit qu'il veut pas plomber l'ambiance.... lol merci Olivier je vais rectifier mon programme de suite!