TP0

LIBNAME LIBRARY 'D:\Users\AdelinouxX\Desktop\Travail_SAS\FORMAT';
LIBNAME OUT 'D:\Users\AdelinouxX\Desktop\Travail_SAS\DATA';

ods pdf file = 'D:\Users\AdelinouxX\Desktop\Travail_SAS\RESULTS\reglog_polytomiale.pdf';
title1 'R?gression logistique polytomiale : X = EVALDEMO4, Y = polintr2, Z = reg3';
title2 'S?lection du mod?le';
proc logistic data = out.datalabel;
        class health3 (ref = 'A) Bonne santé') sblwcoa3 (ref = 'A) Oui')/param=ref;
    model happy3(event= 'B) Moyennement heureux') = health3 sblwcoa3 health3*sblwcoa3/link = glogit;
run;
/*Test de type 3 health3*sblwcoa3 p=0.8664  --> retrait , pas d'interaction*/

proc logistic data = out.datalabel;
        class health3 (ref = 'A) Bonne santé') sblwcoa3 (ref = 'A) Oui')/param=ref;
    model happy3(event= 'B) Moyennement heureux') = health3 sblwcoa3/link = glogit;
run;
/*Test de type 3 : health3  <.0001 sblwcoa3 0.0521 happy3 = health3 */

proc logistic data = out.datalabel;
        class hincfel2 (ref = 'A) Confortable') sblwcoa3 (ref = 'A) Oui')/param=ref;
    model happy3(event= 'B) Moyennement heureux') = hincfel2 sblwcoa3 hincfel2*sblwcoa3/link = glogit;
run;

/*p=0.0057  il y a interaction entre sblwcoa3 et hincfel2 sur happy3 */

proc logistic data = out.datalabel;
        class health3 (ref = 'A) Bonne santé') lknemny3 (ref = 'B) Plutôt pas')/param=ref;
    model happy3(event= 'B) Moyennement heureux') = health3 lknemny3 health3*lknemny3/link = glogit;
run;

/* health3*lknemny3 = 0.3759 pas d'interaction*/

proc logistic data = out.datalabel;
        class health3 (ref = 'A) Bonne santé') lknemny3 (ref = 'B) Plutôt pas')/param=ref;
    model happy3(event= 'B) Moyennement heureux') = health3 lknemny3/link = glogit;
run;

/* health3  <.0001 lknemny3 <.0001 happy3 = health3 + lknemny3 */

*/Modèles retenus :   happy3 = hincfel2*sblwcoa3 */

title2 "Examen des param?tres";
title3 "(4-1) r?gressions : choix 'i' Vs choix 'B) Moyennement heureux')";
proc logistic data = out.datalabel;
    class hincfel2 (ref = 'B) Confortable') sblwcoa3 (ref = 'A) Oui')/param=ref;
    model happy3(event= 'B) Moyennement heureux') = hincfel2 sblwcoa3/link = glogit;
run;

title;
ods pdf close;

ods pdf file = 'D:\Users\AdelinouxX\Desktop\Travail_SAS\FORMAT\RESULTS\cotes_prop_log.pdf';

title 'Modèle final';
proc logistic data = out.selection;
class hincfel2 (ref='A) Confortable') health3 (ref='A) Bonne santé') lknemny3(ref='B) Plutôt pas') sclmeet2 (ref='B) Beaucoup')
 sblwcoa3(ref='A) Oui') sclact2(ref='B) Beaucoup')
/ param=ref ref=first ;
model happy = hincfel2 health3 lknemny3 sclmeet2 sblwcoa3 sclact2

health3*sblwcoa3 hincfel2*sblwcoa3 health3*lknemny3
 / expb clparm=both ;
/*oddratio hincfel2;
/*oddratio health3;*/
/*oddratio lknemny3; */
/*oddratio sclmeet2; */
/*oddratio sblwcoa3;*/
/* oddratio sclact2; */

run;