Discussion :

[Noiram94]

Bonjour,
J'ai une question simple je pense (mais je ne trouve pas la réponse) :
Je voudrais savoir comment on obtient une p-value (chisq ou exact) non arrondie dans une proc freq. Je ne veux pas avoir des résultats comme : <.0001
Merci d'avance.
Marion

[Noiram94]

Bonjour, j'ai trouvé pour le chi² avec la macro suivante, mais je n'arrive pas à l'adapter pour un fisher exact test... Quelqu'un pourrait m'aider ?

macro FREQUENCE

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
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
%macro frequence (fichier,listeg,var,sortie);  
 
	/******************** FREQUENCE DE DATASET ************************ 
	* fichier 	: nom de la table à traiter			* 
	* listeg	: liste des variables				* 
	*  var		: variable à analyser 				*
	*  sortie 	: nom du fichier.rtf en sortie output		* 
 	 ******************************************************************/ 
 
	/* récupération du nombre de variables */ 
	%global nbvar; 
	data fichier; 
		set &fichier; 
		array v{*} &listeg; 
		call symput ("nbvar",dim(v)); 
	run; 
 
	/* on va créer une proc freq pour chaque variable */ 
	%do i=1 %to &nbvar; 
 
		/* récupération du nom de la variable à traiter */ 
		%global nomvar; 
		data fichier; 
			set fichier; 
			array lv{*} &listeg; 
			call symput("nomvar",vname(lv{&i})); 
		run; 
 
* frequence=tableau en liste des effectifs; 
* outpct=recuperation des frequences dans la table frequence; 
* test=recupere les tests du khi deux et le test de mantel haenszel; 
 
		proc freq data=fichier noprint ; 
			tables &nomvar*&var / chisq out=frequence(where=(&nomvar ne . and &var ne .)) outpct; 
		output out=test pchi ; 
 
		data tableau; 
			merge frequence test; 
		    eff=put(pct_col,numx6.1) || " ("|| left(put(count,5.) || ") "); 
 			rename &nomvar=Modalite;
			Variable="&nomvar";
			%let def=across "&var"; 
			format p_pchi numx6.4; 
 			 run;
 
/*Je colle tous les résultats pour obtenir un unique fichier &sortie*/
		data sortie;  
		%if &i=1 %then set tableau
		%else set sortie tableau;;	  
		run;  
 
/*J'efface le fichier test et j'en crée un nouveau contenant des données manquantes afin d'avoir des données manquantes pour les khi² incalculables*/
                proc delete data=test;
                run; 
 
                Data test;
                P_pchi=.;
                run;
%END; 
ods rtf body=&sortie;
 
proc report data=sortie nowd headline headskip; 
column variable modalite &var, eff p_pchi; 
define modalite/group; 
define variable/group;
define &var/&def; 
define eff/ '%col (Freq)'; 
* recuperation de la p_value du Khi Deux; 
define p_pchi/'P X2'; 
break after variable/skip; 
run; 
ods rtf close;
%mend frequence;
 
/**********************************************************************************************/

[hossward]

Bonjour,
J’ ai adapté ta macro pour quelle affiche la valeur « fisher exact test ».
Les variables à analyser doivent être numériques. Tu peux l’essayer :

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
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
 
%macro frequence (fichier,listeg,var,sortie);  
/******************** FREQUENCE DE DATASET ************************
* fichier : nom de la table à traiter *
* listeg : liste des variables *
* var : variable à analyser *
* sortie : nom du fichier.rtf en sortie output *
******************************************************************/
 
/* récupération du nombre de variables */
%global nbvar;
data fichier;
set &fichier;
array v{*} &listeg;
call symputx("nbvar", dim(v));	 
run;	 
 
/* on va créer une proc freq pour chaque variable */
%do i=1 %to &nbvar;
 
/* récupération du nom de la variable à traiter */
%global nomvar;
data fichier;
set fichier;
array lv{*} &listeg;
call symputx("nomvar", vname(lv{&i}));
run;	   
 
* frequence=tableau en liste des effectifs;
* outpct=recuperation des frequences dans la table frequence;
* test=recupere les tests du khi deux et le test de mantel haenszel;
 
proc freq data=fichier noprint ;
tables &nomvar*&var / FISHER out=frequence(where=(&nomvar ne . and &var ne .)) outpct;
output out=test FISHER ;
run;
 
data tableau;
merge frequence test;
eff=put(pct_col, numx6.1) || " ("|| left(put(count, 5.) || ") ") ;
rename &nomvar=Modalite;
Variable="&nomvar";
%let def=across "&var";
format XP2_FISH numx6.4;
run;
 
/*Je colle tous les résultats pour obtenir un unique fichier &sortie*/
data sortie;
%if &i=1 %then set tableau
%else set sortie tableau;;
run;
 
/*J'efface le fichier test et j'en crée un nouveau contenant des données manquantes afin d'avoir des données manquantes pour les khi² incalculables*/
proc delete data=test;
run;
 
Data test;
XP2_FISH=.;
run;
%END;
 
ods rtf body=&sortie;
proc report data=sortie nowd headline headskip;
column variable modalite &var, eff XP2_FISH;
define modalite/group;
define variable/group;
define &var/&def;
define eff/ '%col (Freq)';
* recuperation de la p_value du Khi Deux;
define XP2_FISH/'XP2_FISH';
break after variable/skip;
run;
ods rtf close;
%mend frequence;

Bon Courage
Ward

[Noiram94]

Ca marche! Merci beaucoup!