Discussion :

[Guyguy La Tambouille]

Bonjour,

Je rencontre actuellement un petit soucis et j'espère pouvoir trouver ici une réponse...
Je souhaite effectuer un test d'indépendance entre 2 variables : scoresSEUIL (5 modalités), et lib_FILIERE2 (6 modalités).
De plus, il me faudrait faire ce test d'un côté pour les femmes, puis en suivant pour les hommes.
Le problème que je rencontre est que ce test du khi2 devient un test exact de Fisher (car certains effectifs théoriques sont <5...et même <2!), et SAS me signale :
"WARNING: Computing exact p-values for this problem may require much time and memory. Press the system interrupt key to terminate exact computations."

Voici le code que j'utilise :

Code :

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

1
2
3
4
5
6
7
8
9
title2 "Test d'indépendance entre les scores qualitatifs (seuils) et les filières d'étude (par sexe)";
proc freq data=scores;
label lib_FILIERE2="Filière "
	scoresSEUIL="Seuils";
table scoresSEUIL*lib_FILIERE2 / fisher nopercent norow nocum relrisk expected;
*exact pchi or;
by sexe;
format SEXE f_sexe.;
run;

Du coup je ne parviens pas à obtenir de conclusion sur l'indépendance entre scoresSEUIL et FILIERE2 (quelques soit le sexe).
Auriez vous une petite astuce pour y parvenir ????

Vous remerciant par avance pour votre aide, je souhaite une agréable fin de journée aux lecteurs

Guyguy.

[Guyguy La Tambouille]

Personne n'a d'avis sur la question ??? Cela m'étonne...

De mon côté je suis finalement tombé sur ce lien : la méthode de Monte Carlo semble ainsi fournir une estimation (relativement précise et surtout largement plus rapide) de la p_value du test exact de Fisher.
Voici ci-joint le code associé :

Code :

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

1
2
3
4
5
6
7
8
9
title2 "Test d'indépendance entre les scores qualitatifs (seuils) et les filières d'étude (par sexe)";
proc freq data=scores;
label lib_FILIERE2="Filière "
	scoresSEUIL="Seuils";
table scoresSEUIL*lib_FILIERE2 / chisq nopercent norow nocum relrisk expected;
exact fisher / MC /*Option pour approximation de Monte Carlo*/;
by sexe;
format SEXE f_sexe.;
run;

Cependant, je me demande maintenant comment est ce que la p_value de Monte Carlo doit exactement être interprétée...
Voici ci-dessous l'une des sorties que j'obtiens pour le sexe=Homme :


Si certains savent interpréter correctement cette sortie, je suis preneur

Merci d'avance,

Guyguy.